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