ThirdHitCorrection Class Reference

#include <ThirdHitCorrection.h>

Public Types
typedef PixelRecoRange< float >	Range

Public Member Functions
void	correctRPhiRange (Range &range) const
void	correctRZRange (Range &range) const
void	init (float pt, const DetLayer &layer1, const DetLayer &layer2, const DetLayer &layer3, bool useMultipleScattering, const MultipleScatteringParametrisationMaker *msmaker, bool useBendingCorrection, const MagneticField *bfield)
void	init (float pt, const DetLayer &layer3, bool useMultipleScattering, const MultipleScatteringParametrisationMaker *msmaker, bool useBendingCorrection, const MagneticField *bfield)
void	init (const PixelRecoLineRZ &line, const PixelRecoPointRZ &constraint, int ol)
	ThirdHitCorrection ()
	ThirdHitCorrection (float pt, const DetLayer *layer, const PixelRecoLineRZ &line, const PixelRecoPointRZ &constraint, int ol, bool useMultipleScattering, const MultipleScatteringParametrisationMaker *msmaker, bool useBendingCorrection, const MagneticField *bfield)

Private Attributes
std::optional< MultipleScatteringParametrisation >	sigmaRPhi
bool	theBarrel
pixelrecoutilities::LongitudinalBendingCorrection	theBendingCorrection
PixelRecoLineRZ	theLine
float	theMScoeff = 0
float	theMultScattCorrRPhi = 0
float	thePt
bool	theUseBendingCorrection
bool	theUseMultipleScattering

Detailed Description

Definition at line 17 of file ThirdHitCorrection.h.

Member Typedef Documentation

Range

typedef PixelRecoRange<float> ThirdHitCorrection::Range

Definition at line 19 of file ThirdHitCorrection.h.

Constructor & Destructor Documentation

ThirdHitCorrection() [1/2]

ThirdHitCorrection::ThirdHitCorrection ( )

inline

Definition at line 21 of file ThirdHitCorrection.h.

{}

ThirdHitCorrection() [2/2]

ThirdHitCorrection::ThirdHitCorrection ( float  pt, const DetLayer *  layer, const PixelRecoLineRZ &  line, const PixelRecoPointRZ &  constraint, int  ol, bool  useMultipleScattering, const MultipleScatteringParametrisationMaker *  msmaker, bool  useBendingCorrection, const MagneticField *  bfield  )

inline

Definition at line 39 of file ThirdHitCorrection.h.

References HLT_2023v12_cff::constraint, init(), nano_mu_digi_cff::layer, mps_splice::line, DiDispStaMuonMonitor_cfi::pt, HLT_2023v12_cff::useBendingCorrection, and highPtTripletStepTrackingRegions_cfi::useMultipleScattering.

                                                  {
    init(pt, *layer, useMultipleScattering, msmaker, useBendingCorrection, bfield);
    init(line, constraint, ol);
  }

Member Function Documentation

correctRPhiRange()

void ThirdHitCorrection::correctRPhiRange ( Range &  range ) const

inline

Definition at line 54 of file ThirdHitCorrection.h.

References FastTimerService_cff::range, and theMultScattCorrRPhi.

                                            {
    range.first -= theMultScattCorrRPhi;
    range.second += theMultScattCorrRPhi;
  }

correctRZRange()

void ThirdHitCorrection::correctRZRange

(

Range &

range

)

const

Definition at line 101 of file ThirdHitCorrection.cc.

References funct::abs(), alignCSCRings::corr, CosmicsPD_Skims::radius, FastTimerService_cff::range, and crabWrapper::theLine.

                                                          {
  range.first -= theMScoeff;
  range.second += theMScoeff;

  if (theUseBendingCorrection) {
    if (theBarrel) {
      float cotTheta = theLine.cotLine();
      if (cotTheta > 0) {
        float radius = theLine.rAtZ(range.max());
        float corr = theBendingCorrection(radius) * cotTheta;
        range.second += corr;
      } else {
        float radius = theLine.rAtZ(range.min());
        float corr = theBendingCorrection(radius) * std::abs(cotTheta);
        range.first -= corr;
      }
    } else {
      float radius = range.max();
      float corr = theBendingCorrection(radius);
      range.first -= corr;
    }
  }
}

init() [1/3]

void ThirdHitCorrection::init	(	float	pt,
		const DetLayer &	layer1,
		const DetLayer &	layer2,
		const DetLayer &	layer3,
		bool	useMultipleScattering,
		const MultipleScatteringParametrisationMaker *	msmaker,
		bool	useBendingCorrection,
		const MagneticField *	bfield
	)

Definition at line 40 of file ThirdHitCorrection.cc.

References Surface::bounds(), init, DetLayer::isBarrel(), Bounds::length(), mps_splice::line, GeometricSearchDet::position(), DiDispStaMuonMonitor_cfi::pt, DetLayer::seqNum(), BarrelDetLayer::specificSurface(), ForwardDetLayer::specificSurface(), GeometricSearchDet::surface(), HLT_2023v12_cff::useBendingCorrection, highPtTripletStepTrackingRegions_cfi::useMultipleScattering, PV3DBase< T, PVType, FrameType >::z(), SiPixelPI::zero, and SiStripMonitorCluster_cfi::zmax.

Referenced by ThirdHitCorrection().

                                                           {
  init(pt, layer3, useMultipleScattering, msmaker, useBendingCorrection, bfield);

  if (!theUseMultipleScattering)
    return;

  auto point3 = [&]() -> PixelRecoPointRZ {
    if (theBarrel) {
      const BarrelDetLayer &bl = static_cast<const BarrelDetLayer &>(layer3);
      float rLayer = bl.specificSurface().radius();
      auto zmax = 0.5f * layer3.surface().bounds().length();
      return PixelRecoPointRZ(rLayer, zmax);
    } else {
      const ForwardDetLayer &fl = static_cast<const ForwardDetLayer &>(layer3);
      auto maxR = fl.specificSurface().outerRadius();
      auto layerZ = layer3.position().z();
      return PixelRecoPointRZ(maxR, layerZ);
    }
  };

  PixelRecoPointRZ zero(0., 0.);
  SimpleLineRZ line(zero, point3());

  auto point2 = [&]() -> PixelRecoPointRZ {
    if (layer2.isBarrel()) {
      const BarrelDetLayer &bl = static_cast<const BarrelDetLayer &>(layer2);
      float rLayer = bl.specificSurface().radius();
      return PixelRecoPointRZ(rLayer, line.zAtR(rLayer));
    } else {
      auto layerZ = layer2.position().z();
      return PixelRecoPointRZ(line.rAtZ(layerZ), layerZ);
    }
  };

  theMultScattCorrRPhi = 3.f * (*sigmaRPhi)(pt, line.cotLine(), point2(), layer2.seqNum());
}

init() [2/3]

void ThirdHitCorrection::init	(	float	pt,
		const DetLayer &	layer3,
		bool	useMultipleScattering,
		const MultipleScatteringParametrisationMaker *	msmaker,
		bool	useBendingCorrection,
		const MagneticField *	bfield
	)

Definition at line 18 of file ThirdHitCorrection.cc.

References DetLayer::isBarrel(), MultipleScatteringParametrisationMaker::parametrisation(), DiDispStaMuonMonitor_cfi::pt, HLT_2023v12_cff::useBendingCorrection, and highPtTripletStepTrackingRegions_cfi::useMultipleScattering.

                                                           {
  theUseMultipleScattering = useMultipleScattering;
  theUseBendingCorrection = useBendingCorrection;
  if (useBendingCorrection) {
    theBendingCorrection.init(pt, *bfield);
  }

  theMultScattCorrRPhi = 0;
  theMScoeff = 0;

  theBarrel = layer3.isBarrel();
  thePt = pt;

  if (theUseMultipleScattering)
    sigmaRPhi = msmaker->parametrisation(&layer3);
}

init() [3/3]

void ThirdHitCorrection::init	(	const PixelRecoLineRZ &	line,
		const PixelRecoPointRZ &	constraint,
		int	ol
	)

Definition at line 84 of file ThirdHitCorrection.cc.

References funct::abs(), HLT_2023v12_cff::constraint, MillePedeFileConverter_cfg::e, f, mps_splice::line, funct::sqr(), mathSSE::sqrt(), and crabWrapper::theLine.

                                                                                                     {
  theLine = line;
  if (!theUseMultipleScattering)
    return;

  // auto newCorr = theMultScattCorrRPhi;
  theMultScattCorrRPhi = 3.f * (*sigmaRPhi)(thePt, line.cotLine(), constraint, il);
  // std::cout << "ThirdHitCorr " << (theBarrel ? "B " : "F " )<< theMultScattCorrRPhi << ' ' << newCorr << ' ' << newCorr/theMultScattCorrRPhi << std::endl;
  float overSinTheta = std::sqrt(1.f + sqr(line.cotLine()));
  if (theBarrel) {
    theMScoeff = theMultScattCorrRPhi * overSinTheta;
  } else {
    float overCosTheta = std::abs(line.cotLine()) < 1.e-4f ? 1.e4f : overSinTheta / std::abs(line.cotLine());
    theMScoeff = theMultScattCorrRPhi * overCosTheta;
  }
}

Member Data Documentation

sigmaRPhi

std::optional<MultipleScatteringParametrisation> ThirdHitCorrection::sigmaRPhi

private

Definition at line 73 of file ThirdHitCorrection.h.

theBarrel

bool ThirdHitCorrection::theBarrel

private

Definition at line 62 of file ThirdHitCorrection.h.

theBendingCorrection

pixelrecoutilities::LongitudinalBendingCorrection ThirdHitCorrection::theBendingCorrection

private

Definition at line 72 of file ThirdHitCorrection.h.

theLine

PixelRecoLineRZ ThirdHitCorrection::theLine

private

Definition at line 67 of file ThirdHitCorrection.h.

theMScoeff

float ThirdHitCorrection::theMScoeff = 0

private

Definition at line 69 of file ThirdHitCorrection.h.

theMultScattCorrRPhi

float ThirdHitCorrection::theMultScattCorrRPhi = 0

private

Definition at line 68 of file ThirdHitCorrection.h.

Referenced by correctRPhiRange().

thePt

float ThirdHitCorrection::thePt

private

Definition at line 70 of file ThirdHitCorrection.h.

theUseBendingCorrection

bool ThirdHitCorrection::theUseBendingCorrection

private

Definition at line 65 of file ThirdHitCorrection.h.

theUseMultipleScattering

bool ThirdHitCorrection::theUseMultipleScattering

private

Definition at line 64 of file ThirdHitCorrection.h.