PixelRecoUtilities.h

Go to the documentation of this file.

#ifndef PixelRecoUtilities_H
#define PixelRecoUtilities_H
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include <string>

namespace PixelRecoUtilities {

  struct FieldAt0 {
    FieldAt0(const edm::EventSetup& es);
    float fieldInInvGev;
  };

  inline float fieldInInvGev(const edm::EventSetup& iSetup) {
    static const FieldAt0 fieldAt0(iSetup);
    return fieldAt0.fieldInInvGev;
  }
  //  void MaginTesla(const edm::EventSetup& c);
  template <typename T>
  T bendingRadius(T pt, const edm::EventSetup& iSetup) {
    return pt * fieldInInvGev(iSetup);
  }

  template <typename T>
  T curvature(T InversePt, const edm::EventSetup& iSetup) {
    return InversePt / fieldInInvGev(iSetup);
  }

  template <typename T>
  T inversePt(T curvature, const edm::EventSetup& iSetup) {
    return curvature * fieldInInvGev(iSetup);
  }

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

}  // namespace PixelRecoUtilities

#endif