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 14 of file ThirdHitCorrection.h.

Member Typedef Documentation

typedef PixelRecoRange<float> ThirdHitCorrection::Range

Definition at line 17 of file ThirdHitCorrection.h.

Constructor & Destructor Documentation

ThirdHitCorrection::ThirdHitCorrection ( )

inline

Definition at line 19 of file ThirdHitCorrection.h.

{}

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 40 of file ThirdHitCorrection.h.

References init().

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

Member Function Documentation

void ThirdHitCorrection::correctRPhiRange ( Range &  range ) const

inline

Definition at line 58 of file ThirdHitCorrection.h.

References theMultScattCorrRPhi.

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

void ThirdHitCorrection::correctRZRange

(

Range &

range

)

const

Definition at line 106 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 40 of file ThirdHitCorrection.cc.

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

Referenced by ThirdHitCorrection().

                                 {

  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(), and EnergyCorrector::pt.

                                 {

  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 84 of file ThirdHitCorrection.cc.

References funct::abs(), HLT_FULL_cff::constraint, PixelRecoLineRZ::cotLine(), alignCSCRings::e, f, geometryCSVtoXML::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 77 of file ThirdHitCorrection.h.

bool ThirdHitCorrection::theBarrel

private

Definition at line 66 of file ThirdHitCorrection.h.

pixelrecoutilities::LongitudinalBendingCorrection ThirdHitCorrection::theBendingCorrection

private

Definition at line 76 of file ThirdHitCorrection.h.

PixelRecoLineRZ ThirdHitCorrection::theLine

private

Definition at line 71 of file ThirdHitCorrection.h.

float ThirdHitCorrection::theMScoeff =0

private

Definition at line 73 of file ThirdHitCorrection.h.

float ThirdHitCorrection::theMultScattCorrRPhi =0

private

Definition at line 72 of file ThirdHitCorrection.h.

Referenced by correctRPhiRange().

float ThirdHitCorrection::thePt

private

Definition at line 74 of file ThirdHitCorrection.h.

bool ThirdHitCorrection::theUseBendingCorrection

private

Definition at line 69 of file ThirdHitCorrection.h.

bool ThirdHitCorrection::theUseMultipleScattering

private

Definition at line 68 of file ThirdHitCorrection.h.