SegmentationUtil.h

Go to the documentation of this file.

/*
 * SegmentationUtil.h
 *
 *  Created on: Oct 31, 2013
 *      Author: Christian Grefe, CERN
 */

#ifndef DDSegmentation_SEGMENTATIONUTIL_H_
#define DDSegmentation_SEGMENTATIONUTIL_H_

#include <cmath>
#include <vector>

namespace DD4hep {
namespace DDSegmentation {
namespace Util {





double magFromXYZ(const Vector3D& position) {
        return std::sqrt(position.X * position.X + position.Y * position.Y + position.Z * position.Z);
}

double radiusFromXYZ(const Vector3D& position) {
        return std::sqrt(position.X * position.X + position.Y * position.Y);
}

double thetaFromXYZ(const Vector3D& position) {
        return std::acos(position.Z / radiusFromXYZ(position));
}

double phiFromXYZ(const Vector3D& position) {
        return std::atan2(position.Y, position.X);
}


Vector3D positionFromRPhiZ(double r, double phi, double z) {
        return Vector3D(r * std::cos(phi), r * std::sin(phi), z);
}

double magFromRPhiZ(double r, double /* phi */, double z) {
        return std::sqrt(r * r + z * z);
}

double xFromRPhiZ(double r, double phi, double /* z */) {
        return r * std::cos(phi);
}

double yFromRPhiZ(double r, double phi, double /* z */) {
        return r * std::sin(phi);
}

double thetaFromRPhiZ(double r, double /* phi */, double z) {
        return r * std::atan(z / r);
}


Vector3D positionFromRThetaPhi(double r, double theta, double phi) {
        return Vector3D(r * std::cos(phi), r * std::sin(phi), r * std::tan(theta));
}

Vector3D positionFromMagThetaPhi(double mag, double theta, double phi) {
        double r = mag * sin(theta);
        return Vector3D(r * std::cos(phi), r * std::sin(phi), mag * std::cos(theta));
}

} /* namespace Util */
} /* namespace DDSegmentation */
} /* namespace DD4hep */

#endif /* DDSegmentation_SEGMENTATIONUTIL_H_ */