PixelRecoUtilities.h

Go to the documentation of this file.

#ifndef PixelRecoUtilities_H
#define PixelRecoUtilities_H
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "MagneticField/Engine/interface/MagneticField.h"

namespace PixelRecoUtilities {
  template <typename T>
  T bendingRadius(T pt, const MagneticField& field) {
    return pt * field.inverseBzAtOriginInGeV();
  }

  template <typename T>
  T curvature(T InversePt, const MagneticField& field) {
    return InversePt / field.inverseBzAtOriginInGeV();
  }

  template <typename T>
  T inversePt(T curvature, const MagneticField& field) {
    return curvature * field.inverseBzAtOriginInGeV();
  }

  inline double longitudinalBendingCorrection(double radius, double pt, const MagneticField& field) {
    double invCurv = bendingRadius(pt, field);
    if (invCurv == 0.)
      return 0.;
    return radius / 6. * radius * radius / (2. * invCurv * 2. * invCurv);
  }

}  // namespace PixelRecoUtilities

#endif