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 (const edm::EventSetup &es, float pt, const DetLayer &layer1, const DetLayer &layer2, const DetLayer &layer3, bool useMultipleScattering, bool useBendingCorrection)
void	init (const edm::EventSetup &es, float pt, const DetLayer &layer3, bool useMultipleScattering, bool useBendingCorrection)
void	init (const PixelRecoLineRZ &line, const PixelRecoPointRZ &constraint, int ol)
	ThirdHitCorrection ()
	ThirdHitCorrection (const edm::EventSetup &es, float pt, const DetLayer *layer, const PixelRecoLineRZ &line, const PixelRecoPointRZ &constraint, int ol, bool useMultipleScattering, bool useBendingCorrection)

Private Attributes
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 16 of file ThirdHitCorrection.h.

Member Typedef Documentation

typedef PixelRecoRange<float> ThirdHitCorrection::Range

Definition at line 18 of file ThirdHitCorrection.h.

Constructor & Destructor Documentation

ThirdHitCorrection::ThirdHitCorrection ( )

inline

Definition at line 20 of file ThirdHitCorrection.h.

References init, DiDispStaMuonMonitor_cfi::pt, HLT_2018_cff::useBendingCorrection, and HLT_2018_cff::useMultipleScattering.

{}

ThirdHitCorrection::ThirdHitCorrection ( const edm::EventSetup &  es, float  pt, const DetLayer *  layer, const PixelRecoLineRZ &  line, const PixelRecoPointRZ &  constraint, int  ol, bool  useMultipleScattering, bool  useBendingCorrection  )

inline

Definition at line 36 of file ThirdHitCorrection.h.

References HLT_2018_cff::constraint, init, and mps_splice::line.

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

Member Function Documentation

void ThirdHitCorrection::correctRPhiRange ( Range &  range ) const

inline

Definition at line 50 of file ThirdHitCorrection.h.

References FastTimerService_cff::range.

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

void ThirdHitCorrection::correctRZRange

(

Range &

range

)

const

Definition at line 97 of file ThirdHitCorrection.cc.

References funct::abs(), corr, PixelRecoRange< T >::max(), PixelRecoRange< T >::min(), CosmicsPD_Skims::radius, 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;
    }
  }
}

void ThirdHitCorrection::init	(	const edm::EventSetup &	es,
		float	pt,
		const DetLayer &	layer1,
		const DetLayer &	layer2,
		const DetLayer &	layer3,
		bool	useMultipleScattering,
		bool	useBendingCorrection
	)

Definition at line 37 of file ThirdHitCorrection.cc.

References Surface::bounds(), SimpleLineRZ::cotLine(), init, DetLayer::isBarrel(), Bounds::length(), mps_splice::line, GeometricSearchDet::position(), SimpleLineRZ::rAtZ(), DetLayer::seqNum(), ForwardDetLayer::specificSurface(), BarrelDetLayer::specificSurface(), GeometricSearchDet::surface(), PV3DBase< T, PVType, FrameType >::z(), SimpleLineRZ::zAtR(), and SiStripMonitorCluster_cfi::zmax.

                                                         {
  init(es, pt, layer3, useMultipleScattering, useBendingCorrection);

  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());
}

void ThirdHitCorrection::init	(	const edm::EventSetup &	es,
		float	pt,
		const DetLayer &	layer3,
		bool	useMultipleScattering,
		bool	useBendingCorrection
	)

Definition at line 17 of file ThirdHitCorrection.cc.

References DetLayer::isBarrel(), DiDispStaMuonMonitor_cfi::pt, HLT_2018_cff::useBendingCorrection, and HLT_2018_cff::useMultipleScattering.

                                                         {
  theUseMultipleScattering = useMultipleScattering;
  theUseBendingCorrection = useBendingCorrection;
  if (useBendingCorrection)
    theBendingCorrection.init(pt, es);

  theMultScattCorrRPhi = 0;
  theMScoeff = 0;

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

  if (theUseMultipleScattering)
    sigmaRPhi.init(&layer3, es);
}

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

Definition at line 80 of file ThirdHitCorrection.cc.

References funct::abs(), HLT_2018_cff::constraint, PixelRecoLineRZ::cotLine(), 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

MultipleScatteringParametrisation ThirdHitCorrection::sigmaRPhi

private

Definition at line 69 of file ThirdHitCorrection.h.

bool ThirdHitCorrection::theBarrel

private

Definition at line 58 of file ThirdHitCorrection.h.

pixelrecoutilities::LongitudinalBendingCorrection ThirdHitCorrection::theBendingCorrection

private

Definition at line 68 of file ThirdHitCorrection.h.

PixelRecoLineRZ ThirdHitCorrection::theLine

private

Definition at line 63 of file ThirdHitCorrection.h.

float ThirdHitCorrection::theMScoeff = 0

private

Definition at line 65 of file ThirdHitCorrection.h.

float ThirdHitCorrection::theMultScattCorrRPhi = 0

private

Definition at line 64 of file ThirdHitCorrection.h.

float ThirdHitCorrection::thePt

private

Definition at line 66 of file ThirdHitCorrection.h.

bool ThirdHitCorrection::theUseBendingCorrection

private

Definition at line 61 of file ThirdHitCorrection.h.

bool ThirdHitCorrection::theUseMultipleScattering

private

Definition at line 60 of file ThirdHitCorrection.h.