MultipleScatteringParametrisation Class Reference

#include <MultipleScatteringParametrisation.h>

Public Types
enum	Consecutive { notAssumeConsecutive, useConsecutive }
enum	X0Source { useDetLayer =0, useX0AtEta =1, useX0DataAveraged =2 }

Public Member Functions
void	init (const DetLayer *layer, const edm::EventSetup &iSetup, X0Source x0source=useX0AtEta)
	MultipleScatteringParametrisation ()
	MultipleScatteringParametrisation (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
float	operator() (float pT, const PixelRecoPointRZ &pointV, const PixelRecoPointRZ &pointO, int ol) const

Static Public Member Functions
static void	initKeepers (const edm::EventSetup &iSetup)

Private Attributes
MSLayer	theLayer
MSLayersKeeper const *	theLayerKeeper = 0

Static Private Attributes
static const float	x0ToSigma = 0.0136f

Detailed Description

Parametrisation of multiple scattering sigma in tracker.

Definition at line 20 of file MultipleScatteringParametrisation.h.

Member Enumeration Documentation

enum MultipleScatteringParametrisation::Consecutive

Enumerator
notAssumeConsecutive
useConsecutive

Definition at line 27 of file MultipleScatteringParametrisation.h.

{ notAssumeConsecutive, useConsecutive };

enum MultipleScatteringParametrisation::X0Source

Enumerator
useDetLayer
useX0AtEta
useX0DataAveraged

Definition at line 26 of file MultipleScatteringParametrisation.h.

{ useDetLayer=0, useX0AtEta=1, useX0DataAveraged=2 };

Constructor & Destructor Documentation

MultipleScatteringParametrisation::MultipleScatteringParametrisation ( )

inline

Definition at line 29 of file MultipleScatteringParametrisation.h.

{}

MultipleScatteringParametrisation::MultipleScatteringParametrisation ( const DetLayer *  layer, const edm::EventSetup &  iSetup, X0Source  x0source = useX0AtEta  )

inline

Definition at line 31 of file MultipleScatteringParametrisation.h.

References init().

                                                                     {
       init(layer,iSetup,x0source);
  }

Member Function Documentation

void MultipleScatteringParametrisation::init	(	const DetLayer *	layer,
		const edm::EventSetup &	iSetup,
		X0Source	x0source = useX0AtEta
	)

Definition at line 56 of file MultipleScatteringParametrisation.cc.

Referenced by MultipleScatteringParametrisation().

                                                                            {

  theLayerKeeper = keepers(x0Source);

  // FIXME not thread safe: move elsewhere...
  initKeepers(iSetup);

  if (!layer) return;
  theLayer = theLayerKeeper->layer(layer);
} 

void MultipleScatteringParametrisation::initKeepers ( const edm::EventSetup &  iSetup )

static

Definition at line 49 of file MultipleScatteringParametrisation.cc.

                                                                              {
  keepers.init(iSetup);
}

float MultipleScatteringParametrisation::operator()	(	float	pt,
		float	cotTheta,
		float	transverseIP = 0.
	)		const

MS sigma at the layer for which parametrisation is initialised; particle assumed to come from nominal vertex, "fast" methods called

Definition at line 68 of file MultipleScatteringParametrisation.cc.

{
  float sumX0D = theLayer.sumX0D(cotTheta); 
  return x0ToSigma * sumX0D /pT;
}

float MultipleScatteringParametrisation::operator()	(	float	pt,
		float	cotTheta,
		const PixelRecoPointRZ &	point,
		float	transverseIP = 0.
	)		const

MS sigma at the layer for which parametrisation is initialised; particle assumed to come from constraint point (inner to layer). layer by layer contribution is calculated

Definition at line 76 of file MultipleScatteringParametrisation.cc.

References MSLayersAtAngle::sumX0D().

{

  PixelRecoLineRZ lineIO(pointI, cotTheta, tip);
  PixelRecoPointRZ pointO = theLayer.crossing(lineIO).first;

  const MSLayersAtAngle & layersAtEta = theLayerKeeper->layers(cotTheta);
  
  float sumX0D = layersAtEta.sumX0D(pointI, pointO);
  return x0ToSigma * sumX0D /pT;
}

float MultipleScatteringParametrisation::operator()	(	float	pt,
		float	cotTheta,
		const PixelRecoPointRZ &	point,
		int	ol
	)		const

Definition at line 91 of file MultipleScatteringParametrisation.cc.

References MSLayersAtAngle::sumX0D().

                                                                                                                      {

  PixelRecoLineRZ lineIO(pointI, cotTheta);
  PixelRecoPointRZ pointO = theLayer.crossing(lineIO).first;

  const MSLayersAtAngle & layersAtEta = theLayerKeeper->layers(cotTheta);
  
  float sumX0D = layersAtEta.sumX0D(il, theLayer.seqNum(), pointI, pointO);
  return x0ToSigma * sumX0D /pT;
}

float MultipleScatteringParametrisation::operator()	(	float	pt,
		const PixelRecoPointRZ &	point1,
		const PixelRecoPointRZ &	point2,
		Consecutive	consecutive = notAssumeConsecutive,
		float	transverseIP = 0.
	)		const

MS sigma at the layer for which parametrisation is initialised; particle assumed to be measured at point1 and point2, it is assumed that layer is between point1 and point2. layer by layer contribution is calculated

Definition at line 104 of file MultipleScatteringParametrisation.cc.

References PixelRecoLineRZ::cotLine(), PixelRecoPointRZ::r(), mathSSE::sqrt(), and MSLayersAtAngle::sumX0D().

{   


  PixelRecoLineRZ lineIO(pointI, pointO, tip);
  PixelRecoPointRZ pointM = theLayer.crossing(lineIO).first;
  float cotTheta = lineIO.cotLine();

  if (consecutive==useConsecutive) {
    float dist = fabs(  (pointO.r()-pointM.r())
                      * (pointM.r()-pointI.r())
                      / (pointO.r()-pointI.r()) );
    return  x0ToSigma * sqrt(theLayer.x0(cotTheta)) * dist /pT;
  } else {
    const MSLayersAtAngle & layersAtEta = theLayerKeeper->layers(cotTheta);
    float sumX0D = layersAtEta.sumX0D(pointI, pointM, pointO);
    return x0ToSigma * sumX0D /pT;
  }
}

float MultipleScatteringParametrisation::operator()	(	float	pT,
		const PixelRecoPointRZ &	pointV,
		const PixelRecoPointRZ &	pointO,
		int	ol
	)		const

Definition at line 129 of file MultipleScatteringParametrisation.cc.

References PixelRecoLineRZ::cotLine(), MSLayersAtAngle::sumX0D(), and PixelRecoPointRZ::z().

{   

  PixelRecoLineRZ lineIO(pointV, pointO);
  PixelRecoPointRZ pointI = theLayer.crossing(lineIO).first;
  float cotTheta = lineIO.cotLine();

  const MSLayersAtAngle & layersAtEta = theLayerKeeper->layers(cotTheta);
  float sumX0D = layersAtEta.sumX0D(pointV.z(), theLayer.seqNum(), ol, pointI, pointO);
  return x0ToSigma * sumX0D /pT;
}

Member Data Documentation

MSLayer MultipleScatteringParametrisation::theLayer

private

Definition at line 78 of file MultipleScatteringParametrisation.h.

MSLayersKeeper const* MultipleScatteringParametrisation::theLayerKeeper = 0

private

Definition at line 79 of file MultipleScatteringParametrisation.h.

const float MultipleScatteringParametrisation::x0ToSigma = 0.0136f

staticprivate

Definition at line 80 of file MultipleScatteringParametrisation.h.