MultipleScatteringParametrisation.h

Go to the documentation of this file.

#ifndef MultipleScatteringParametrisation_H
#define MultipleScatteringParametrisation_H
 
#include "TrackingTools/DetLayers/interface/DetLayer.h"
#include "RecoTracker/TkMSParametrization/interface/PixelRecoPointRZ.h"
#include "RecoTracker/TkMSParametrization/interface/MSLayer.h"
#include "FWCore/Framework/interface/EventSetup.h"
 
class MSLayersKeeper;
class PixelRecoPointRZ;
class DetLayer;
 
class MultipleScatteringParametrisation {
public:
  static void initKeepers(const edm::EventSetup &iSetup);
 
  enum X0Source { useDetLayer = 0, useX0AtEta = 1, useX0DataAveraged = 2 };
  enum Consecutive { notAssumeConsecutive, useConsecutive };
 
  MultipleScatteringParametrisation() {}
 
  MultipleScatteringParametrisation(const DetLayer *layer,
                                    const edm::EventSetup &iSetup,
                                    X0Source x0source = useX0AtEta) {
    init(layer, iSetup, x0source);
  }
 
  void init(const DetLayer *layer, const edm::EventSetup &iSetup, X0Source x0source = useX0AtEta);
 
  float operator()(float pt, float cotTheta, float transverseIP = 0.) const;
 
  float operator()(float pt, float cotTheta, const PixelRecoPointRZ &point, float transverseIP = 0.) const;
  float operator()(float pt, float cotTheta, const PixelRecoPointRZ &point, int ol) const;
 
  float operator()(float pt,
                   const PixelRecoPointRZ &point1,
                   const PixelRecoPointRZ &point2,
                   Consecutive consecutive = notAssumeConsecutive,
                   float transverseIP = 0.) const;
 
  // as above, pointV is at vertex and pointO is on layer ol
  float operator()(float pT, const PixelRecoPointRZ &pointV, const PixelRecoPointRZ &pointO, int ol) const;
 
private:
  MSLayer theLayer;
  MSLayersKeeper const *theLayerKeeper = nullptr;
  static const float x0ToSigma;
};
#endif