/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/RecoTracker/TkMSParametrization/interface/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  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);
  }


}
 
#endif