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"
#include "FWCore/Utilities/interface/GCC11Compatibility.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