#include <BaseParticlePropagator.h>

Inheritance diagram for BaseParticlePropagator:

Public Member Functions
bool	backPropagate ()

	BaseParticlePropagator ()
	Default c'tor. More...

	BaseParticlePropagator (const RawParticle &myPart, double r, double z, double B)

	BaseParticlePropagator (const RawParticle &myPart, double r, double z, double B, double t)

double	getMagneticField () const
	Get the magnetic field. More...

int	getSuccess () const
	Has propagation been performed and was barrel or endcap reached ? More...

bool	hasDecayed () const
	Has the particle decayed while propagated ? More...

double	helixCentreDistToAxis () const
	The distance between the cylinder and the helix axes. More...

double	helixCentreDistToAxis (double xC, double yC) const

double	helixCentrePhi () const
	The azimuth if the vector joining the cylinder and the helix axes. More...

double	helixCentrePhi (double xC, double yC) const

double	helixCentreX () const
	The x coordinate of the helix axis. More...

double	helixCentreX (double radius, double phi) const

double	helixCentreY () const
	The y coordinate of the helix axis. More...

double	helixCentreY (double radius, double phi) const

double	helixRadius () const
	The helix Radius. More...

double	helixRadius (double pT) const

double	helixStartPhi () const
	The azimuth of the momentum at the vertex. More...

void	increaseRCyl (double delta)
	Just an internal trick. More...

void	init ()
	Initialize internal switches and quantities. More...

bool	inside () const
	Is the vertex inside the cylinder ? (stricly inside : true) More...

bool	inside (double rPos2) const

bool	onBarrel () const
	Is the vertex already on the cylinder barrel ? More...

bool	onBarrel (double rPos2) const

bool	onEndcap () const
	Is the vertex already on the cylinder endcap ? More...

bool	onEndcap (double rPos2) const

bool	onFiducial () const
	Is the vertex on some material ? More...

bool	onSurface () const
	Is the vertex already on the cylinder surface ? More...

bool	onSurface (double rPos2) const

RawParticle &	particle ()

RawParticle const &	particle () const
	The particle being propagated. More...

bool	propagate ()

BaseParticlePropagator	propagated () const

bool	propagateToBeamCylinder (const XYZTLorentzVector &v, double radius=0.)

bool	propagateToClosestApproach (double x0=0., double y0=0, bool first=true)

bool	propagateToEcal (bool first=true)

bool	propagateToEcalEntrance (bool first=true)

bool	propagateToHcalEntrance (bool first=true)

bool	propagateToHcalExit (bool first=true)

bool	propagateToHOLayer (bool first=true)

bool	propagateToNominalVertex (const XYZTLorentzVector &hit2=XYZTLorentzVector(0., 0., 0., 0.))

bool	propagateToPreshowerLayer1 (bool first=true)

bool	propagateToPreshowerLayer2 (bool first=true)

bool	propagateToVFcalEntrance (bool first=true)

void	resetDebug ()

void	setDebug ()
	Set the debug leve;. More...

void	setMagneticField (double b)
	Set the magnetic field. More...

void	setParticle (RawParticle const &iParticle)

void	setPropagationConditions (double r, double z, bool firstLoop=true)
	Set the propagation characteristics (rCyl, zCyl and first loop only) More...

void	setProperDecayTime (double t)
	Set the proper decay time. More...

double	xyImpactParameter (double x0=0., double y0=0.) const
	Transverse impact parameter. More...

double	zImpactParameter (double x0=0, double y0=0.) const
	Longitudinal impact parameter. More...

Protected Member Functions
double	c_light () const
	The speed of light in mm/ns (!) without clhep (yeaaahhh!) More...

Protected Attributes
bool	fiducial
	The particle traverses some real material. More...

int	success
	0:propagation still be done, 1:reached 'barrel', 2:reached 'endcaps' More...

Private Attributes
double	bField
	Magnetic field in the cylinder, oriented along the Z axis. More...

bool	debug
	The debug level. More...

bool	decayed
	The particle decayed while propagated ! More...

bool	firstLoop
	Do only the first half-loop. More...

RawParticle	particle_

int	propDir
	The propagation direction. More...

double	properDecayTime
	The proper decay time of the particle. More...

double	properTime
	The proper time of the particle. More...

double	rCyl
	Simulated particle that is to be resp has been propagated. More...

double	rCyl2

double	zCyl
	Half-height of the cylinder (centred at 0,0,0) to which propagation is done. More...

Detailed Description

This class is aimed at propagating charged and neutral particles (yet under the form of a RawParticle) from a given vertex to a cylinder, defined by a radius rCyl and a length 2*zCyl, centered in (0,0,0) and whose axis is parallel to the B field (B is oriented along z, by definition of the z axis).

The main method

bool propagate()

returns true if an intersection between the propagated RawParticle and the cylinder is found. The location of the intersection is given by the value of success:

success = 2 : an intersection is found forward in the cylinder endcaps. The RawParticle vertex and momentum are overwritten by the values at this intersection.
success = 1 : an intersection is found forward in the cylinder barrel. The RawParticle vertex and momentum are overwritten by the values at this intersection.
success =-2 : an intersection is found backward in the cylinder endcaps. The RawParticle vertex and momentum are NOT overwritten by the values at this intersection.
success =-1 : an intersection is found backward in the cylinder barrel. The RawParticle vertex and momentum are NOT overwritten by the values at this intersection.
success = 0 : No intersection is found after half a loop (only possible if firstLoop=true). The vertex and momentum are overwritten by the values after this half loop.

The method

propagated()

returns a new RawParticle with the propagated coordinates, if overwriting is not considered an advantage by the user.

Particles can be propagated backward with the method

backPropagate()

Member functions

o propagateToPreshowerLayer1(),
o propagateToPreshowerLayer2(),
o propagateToEcalEntrance(),
o propagateToHcalEntrance(),
o propagateToHcalExit(),
o propagateToClosestApproach(),

only need a momentum, a vertex and an electric charge to operate. Radii, half-lengths and default B field (4T) are set therein by default.

As of today, no average loss of energy (dE/dx, Brems) is considered in the propagation. No uncertainty (multiple scattering, dE/dx, Brems) is yet implemented.

Author: Patrick Janot $Date : 19-Aug-2002, with subsequent modification for FAMOS

Version: 15-Dec-2003

Definition at line 82 of file BaseParticlePropagator.h.

Constructor & Destructor Documentation

◆ BaseParticlePropagator() [1/3]

BaseParticlePropagator::BaseParticlePropagator ( )

Default c'tor.

Definition at line 5 of file BaseParticlePropagator.cc.

5 : rCyl(0.), rCyl2(0.), zCyl(0.), bField(0), properDecayTime(1E99) { ; }

◆ BaseParticlePropagator() [2/3]

BaseParticlePropagator::BaseParticlePropagator	(	const RawParticle &	myPart,
		double	r,
		double	z,
		double	B
	)

Constructors taking as arguments a RawParticle, as well as the radius, half-height and magnetic field defining the cylinder for which propagation is to be performed, and optionally, the proper decay time

Definition at line 12 of file BaseParticlePropagator.cc.

     : particle_(myPart), rCyl(R), rCyl2(R * R), zCyl(Z), bField(B), properDecayTime(1E99) {
   init();
 }

References init().

◆ BaseParticlePropagator() [3/3]

BaseParticlePropagator::BaseParticlePropagator	(	const RawParticle &	myPart,
		double	r,
		double	z,
		double	B,
		double	t
	)

Definition at line 7 of file BaseParticlePropagator.cc.

     : particle_(myPart), rCyl(R), rCyl2(R * R), zCyl(Z), bField(B), properDecayTime(t) {
   init();
 }

References init().

Member Function Documentation

◆ backPropagate()

bool BaseParticlePropagator::backPropagate ( )

Update the current instance, after the back-propagation of the particle to the surface of the cylinder

Definition at line 291 of file BaseParticlePropagator.cc.

                                            {
   // Backpropagate
   particle_.setMomentum(-particle_.Px(), -particle_.Py(), -particle_.Pz(), particle_.E());
   particle_.setCharge(-particle_.charge());
   propDir = -1;
   bool done = propagate();
   particle_.setMomentum(-particle_.Px(), -particle_.Py(), -particle_.Pz(), particle_.E());
   particle_.setCharge(-particle_.charge());
   propDir = +1;
  
   return done;
 }

References RawParticle::charge(), fileCollector::done, RawParticle::E(), particle_, propagate(), propDir, RawParticle::Px(), RawParticle::Py(), RawParticle::Pz(), RawParticle::setCharge(), and RawParticle::setMomentum().

Referenced by propagated(), and propagateToClosestApproach().

◆ c_light()

double BaseParticlePropagator::c_light ( ) const

inlineprotected

The speed of light in mm/ns (!) without clhep (yeaaahhh!)

Definition at line 150 of file BaseParticlePropagator.h.

150 { return 299.792458; }

Referenced by helixRadius(), and propagateToBeamCylinder().

◆ getMagneticField()

double BaseParticlePropagator::getMagneticField ( ) const

inline

Get the magnetic field.

Definition at line 302 of file BaseParticlePropagator.h.

302 { return bField; }

References bField.

Referenced by TrajectoryManager::createPSimHits(), and ConvBremSeedProducer::produce().

◆ getSuccess()

int BaseParticlePropagator::getSuccess ( ) const

inline

Has propagation been performed and was barrel or endcap reached ?

Definition at line 296 of file BaseParticlePropagator.h.

296 { return success; }

References success.

Referenced by PFTrackTransformer::addPoints(), PFTrackTransformer::addPointsAndBrems(), reco::tau::atECALEntrance(), FBaseSimEvent::fill(), ConvBremSeedProducer::GoodCluster(), ConvBremSeedProducer::produce(), TrajectoryManager::propagateToCalorimeters(), and TrajectoryManager::propagateToLayer().

◆ hasDecayed()

bool BaseParticlePropagator::hasDecayed ( ) const

inline

Has the particle decayed while propagated ?

Definition at line 293 of file BaseParticlePropagator.h.

293 { return decayed; }

References decayed.

Referenced by ParticlePropagator::propagateToBoundSurface(), and TrajectoryManager::propagateToCalorimeters().

◆ helixCentreDistToAxis() [1/2]

double BaseParticlePropagator::helixCentreDistToAxis ( ) const

inline

The distance between the cylinder and the helix axes.

Definition at line 235 of file BaseParticlePropagator.h.

                                               {
     // The distance between the cylinder and the helix axes
     double xC = helixCentreX();
     double yC = helixCentreY();
     return std::sqrt(xC * xC + yC * yC);
   }

References helixCentreX(), helixCentreY(), and mathSSE::sqrt().

Referenced by propagate(), propagateToBeamCylinder(), propagateToClosestApproach(), and xyImpactParameter().

◆ helixCentreDistToAxis() [2/2]

double BaseParticlePropagator::helixCentreDistToAxis	(	double	xC,
		double	yC
	)		const

inline

Definition at line 242 of file BaseParticlePropagator.h.

                                                                   {
     // Faster version of helixCentreDistToAxis
     return std::sqrt(xC * xC + yC * yC);
   }

References mathSSE::sqrt().

◆ helixCentrePhi() [1/2]

double BaseParticlePropagator::helixCentrePhi ( ) const

inline

The azimuth if the vector joining the cylinder and the helix axes.

Definition at line 248 of file BaseParticlePropagator.h.

                                        {
     // The azimuth if the vector joining the cylinder and the helix axes
     double xC = helixCentreX();
     double yC = helixCentreY();
     return xC == 0.0 && yC == 0.0 ? 0.0 : std::atan2(yC, xC);
   }

References helixCentreX(), and helixCentreY().

Referenced by propagate().

◆ helixCentrePhi() [2/2]

double BaseParticlePropagator::helixCentrePhi	(	double	xC,
		double	yC
	)		const

inline

Definition at line 255 of file BaseParticlePropagator.h.

                                                            {
     // Faster version of helixCentrePhi()
     return xC == 0.0 && yC == 0.0 ? 0.0 : std::atan2(yC, xC);
   }

◆ helixCentreX() [1/2]

double BaseParticlePropagator::helixCentreX ( ) const

inline

The x coordinate of the helix axis.

Definition at line 213 of file BaseParticlePropagator.h.

                                      {
     // The x coordinate of the helix axis
     return particle_.X() - helixRadius() * std::sin(helixStartPhi());
   }

References helixRadius(), helixStartPhi(), particle_, funct::sin(), and RawParticle::X().

Referenced by helixCentreDistToAxis(), helixCentrePhi(), propagate(), propagateToBeamCylinder(), propagateToClosestApproach(), and xyImpactParameter().

◆ helixCentreX() [2/2]

double BaseParticlePropagator::helixCentreX	(	double	radius,
		double	phi
	)		const

inline

Definition at line 218 of file BaseParticlePropagator.h.

                                                               {
     // Fast version of helixCentreX()
     return particle_.X() - radius * std::sin(phi);
   }

References particle_, phi, CosmicsPD_Skims::radius, funct::sin(), and RawParticle::X().

◆ helixCentreY() [1/2]

double BaseParticlePropagator::helixCentreY ( ) const

inline

The y coordinate of the helix axis.

Definition at line 224 of file BaseParticlePropagator.h.

                                      {
     // The y coordinate of the helix axis
     return particle_.Y() + helixRadius() * std::cos(helixStartPhi());
   }

References funct::cos(), helixRadius(), helixStartPhi(), particle_, and RawParticle::Y().

Referenced by helixCentreDistToAxis(), helixCentrePhi(), propagate(), propagateToBeamCylinder(), propagateToClosestApproach(), and xyImpactParameter().

◆ helixCentreY() [2/2]

double BaseParticlePropagator::helixCentreY	(	double	radius,
		double	phi
	)		const

inline

Definition at line 229 of file BaseParticlePropagator.h.

                                                               {
     // Fast version of helixCentreX()
     return particle_.Y() + radius * std::cos(phi);
   }

References funct::cos(), particle_, phi, CosmicsPD_Skims::radius, and RawParticle::Y().

◆ helixRadius() [1/2]

double BaseParticlePropagator::helixRadius ( ) const

inline

The helix Radius.

Definition at line 189 of file BaseParticlePropagator.h.

                                     {
     // The helix Radius
     //
     // The helix' Radius sign accounts for the orientation of the magnetic field
     // (+ = along z axis) and the sign of the electric charge of the particle.
     // It signs the rotation of the (charged) particle around the z axis:
     // Positive means anti-clockwise, negative means clockwise rotation.
     //
     // The radius is returned in cm to match the units in RawParticle.
     return particle_.charge() == 0 ? 0.0 : -particle_.Pt() / (c_light() * 1e-5 * bField * particle_.charge());
   }

References bField, c_light(), RawParticle::charge(), MillePedeFileConverter_cfg::e, particle_, and RawParticle::Pt().

Referenced by helixCentreX(), helixCentreY(), propagate(), propagateToClosestApproach(), propagateToNominalVertex(), and xyImpactParameter().

◆ helixRadius() [2/2]

double BaseParticlePropagator::helixRadius ( double pT ) const

inline

Definition at line 201 of file BaseParticlePropagator.h.

                                              {
     // a faster version of helixRadius, once Perp() has been computed
     return particle_.charge() == 0 ? 0.0 : -pT / (c_light() * 1e-5 * bField * particle_.charge());
   }

References bField, c_light(), RawParticle::charge(), MillePedeFileConverter_cfg::e, particle_, and PVValHelper::pT.

◆ helixStartPhi()

double BaseParticlePropagator::helixStartPhi ( ) const

inline

The azimuth of the momentum at the vertex.

Definition at line 207 of file BaseParticlePropagator.h.

                                       {
     // The azimuth of the momentum at the vertex
     return particle_.Px() == 0.0 && particle_.Py() == 0.0 ? 0.0 : std::atan2(particle_.Py(), particle_.Px());
   }

References particle_, RawParticle::Px(), and RawParticle::Py().

Referenced by helixCentreX(), helixCentreY(), propagate(), propagateToClosestApproach(), and xyImpactParameter().

◆ increaseRCyl()

void BaseParticlePropagator::increaseRCyl ( double delta )

inline

Just an internal trick.

Definition at line 172 of file BaseParticlePropagator.h.

                                          {
     rCyl = rCyl + delta;
     rCyl2 = rCyl * rCyl;
   }

References dumpMFGeometry_cfg::delta, rCyl, and rCyl2.

◆ init()

void BaseParticlePropagator::init ( void )

Initialize internal switches and quantities.

Definition at line 17 of file BaseParticlePropagator.cc.

                                   {
   // The status of the propagation
   success = 0;
   // Propagate only for the first half-loop
   firstLoop = true;
   //
   fiducial = true;
   // The particle has not yet decayed
   decayed = false;
   // The proper time is set to zero
   properTime = 0.;
   // The propagation direction is along the momentum
   propDir = 1;
   //
   debug = false;
 }

References debug, decayed, fiducial, firstLoop, propDir, properTime, and success.

Referenced by BaseParticlePropagator().

◆ inside() [1/2]

bool BaseParticlePropagator::inside ( ) const

inline

Is the vertex inside the cylinder ? (stricly inside : true)

Definition at line 261 of file BaseParticlePropagator.h.

                              {
     return (particle_.R2() < rCyl2 - 0.00001 * rCyl && fabs(particle_.Z()) < zCyl - 0.00001);
   }

References particle_, RawParticle::R2(), rCyl, rCyl2, RawParticle::Z(), and zCyl.

Referenced by propagate().

◆ inside() [2/2]

bool BaseParticlePropagator::inside ( double rPos2 ) const

inline

Definition at line 265 of file BaseParticlePropagator.h.

                                          {
     return (rPos2 < rCyl2 - 0.00001 * rCyl && fabs(particle_.Z()) < zCyl - 0.00001);
   }

References particle_, rCyl, rCyl2, RawParticle::Z(), and zCyl.

◆ onBarrel() [1/2]

bool BaseParticlePropagator::onBarrel ( ) const

inline

Is the vertex already on the cylinder barrel ?

Definition at line 275 of file BaseParticlePropagator.h.

                                {
     double rPos2 = particle_.R2();
     return (fabs(rPos2 - rCyl2) < 0.00001 * rCyl && fabs(particle_.Z()) <= zCyl);
   }

References particle_, RawParticle::R2(), rCyl, rCyl2, RawParticle::Z(), and zCyl.

Referenced by onSurface(), and propagate().

◆ onBarrel() [2/2]

bool BaseParticlePropagator::onBarrel ( double rPos2 ) const

inline

Definition at line 280 of file BaseParticlePropagator.h.

                                            {
     return (fabs(rPos2 - rCyl2) < 0.00001 * rCyl && fabs(particle_.Z()) <= zCyl);
   }

References particle_, rCyl, rCyl2, RawParticle::Z(), and zCyl.

◆ onEndcap() [1/2]

bool BaseParticlePropagator::onEndcap ( ) const

inline

Is the vertex already on the cylinder endcap ?

Definition at line 285 of file BaseParticlePropagator.h.

285 { return (fabs(fabs(particle_.Z()) - zCyl) < 0.00001 && particle_.R2() <= rCyl2); }

References particle_, RawParticle::R2(), rCyl2, RawParticle::Z(), and zCyl.

Referenced by onSurface(), and propagate().

◆ onEndcap() [2/2]

bool BaseParticlePropagator::onEndcap ( double rPos2 ) const

inline

Definition at line 287 of file BaseParticlePropagator.h.

287 { return (fabs(fabs(particle_.Z()) - zCyl) < 0.00001 && rPos2 <= rCyl2); }

References particle_, rCyl2, RawParticle::Z(), and zCyl.

◆ onFiducial()

bool BaseParticlePropagator::onFiducial ( ) const

inline

Is the vertex on some material ?

Definition at line 290 of file BaseParticlePropagator.h.

290 { return fiducial; }

References fiducial.

Referenced by TrajectoryManager::propagateToLayer().

◆ onSurface() [1/2]

bool BaseParticlePropagator::onSurface ( ) const

inline

Is the vertex already on the cylinder surface ?

Definition at line 270 of file BaseParticlePropagator.h.

270 { return (onBarrel() || onEndcap()); }

References onBarrel(), and onEndcap().

Referenced by propagate().

◆ onSurface() [2/2]

bool BaseParticlePropagator::onSurface ( double rPos2 ) const

inline

Definition at line 272 of file BaseParticlePropagator.h.

272 { return (onBarrel(rPos2) || onEndcap(rPos2)); }

References onBarrel(), and onEndcap().

◆ particle() [1/2]

RawParticle& BaseParticlePropagator::particle ( )

inline

Definition at line 165 of file BaseParticlePropagator.h.

165 { return particle_; }

References particle_.

◆ particle() [2/2]

RawParticle const& BaseParticlePropagator::particle ( ) const

inline

The particle being propagated.

Definition at line 164 of file BaseParticlePropagator.h.

164 { return particle_; }

References particle_.

Referenced by PFTrackTransformer::addPoints(), PFTrackTransformer::addPointsAndBrems(), MuonSimHitProducer::applyMaterialEffects(), reco::tau::atECALEntrance(), TrajectoryManager::createPSimHits(), ParticlePropagator::fieldMap(), FBaseSimEvent::fill(), ConvBremSeedProducer::GoodCluster(), ParticlePropagator::initProperDecayTime(), MaterialEffects::interact(), MaterialEffects::normalVector(), PythiaDecays::particleDaughters(), ParticlePropagator::ParticlePropagator(), ConvBremSeedProducer::produce(), ParticlePropagator::propagateToBoundSurface(), l1tpf::propagateToCalo(), TrajectoryManager::propagateToCalorimeters(), MaterialEffects::radLengths(), and TrajectoryManager::updateWithDaughters().

◆ propagate()

bool BaseParticlePropagator::propagate ( )

Update the current instance, after the propagation of the particle to the surface of the cylinder

Definition at line 34 of file BaseParticlePropagator.cc.

                                        {
   //
   // Check that the particle is not already on the cylinder surface
   //
   double rPos2 = particle_.R2();
  
   if (onBarrel(rPos2)) {
     success = 1;
     return true;
   }
   //
   if (onEndcap(rPos2)) {
     success = 2;
     return true;
   }
   //
   // Treat neutral particles (or no magnetic field) first
   //
   if (fabs(particle_.charge()) < 1E-12 || bField < 1E-12) {
     //
     // First check that the particle crosses the cylinder
     //
     double pT2 = particle_.Perp2();
     //    double pT2 = pT*pT;
     //    double b2 = rCyl * pT;
     double b2 = rCyl2 * pT2;
     double ac = fabs(particle_.X() * particle_.Py() - particle_.Y() * particle_.Px());
     double ac2 = ac * ac;
     //
     // The particle never crosses (or never crossed) the cylinder
     //
     //    if ( ac > b2 ) return false;
     if (ac2 > b2)
       return false;
  
     //    double delta  = std::sqrt(b2*b2 - ac*ac);
     double delta = std::sqrt(b2 - ac2);
     double tplus = -(particle_.X() * particle_.Px() + particle_.Y() * particle_.Py()) + delta;
     double tminus = -(particle_.X() * particle_.Px() + particle_.Y() * particle_.Py()) - delta;
     //
     // Find the first (time-wise) intersection point with the cylinder
     //
  
     double solution = tminus < 0 ? tplus / pT2 : tminus / pT2;
     if (solution < 0.)
       return false;
     //
     // Check that the particle does (did) not exit from the endcaps first.
     //
     double zProp = particle_.Z() + particle_.Pz() * solution;
     if (fabs(zProp) > zCyl) {
       tplus = (zCyl - particle_.Z()) / particle_.Pz();
       tminus = (-zCyl - particle_.Z()) / particle_.Pz();
       solution = tminus < 0 ? tplus : tminus;
       if (solution < 0.)
         return false;
       success = 2;
     } else {
       success = 1;
     }
  
     //
     // Check the decay time
     //
     double delTime = propDir * particle_.mass() * solution;
     double factor = 1.;
     properTime += delTime;
     if (properTime > properDecayTime) {
       factor = 1. - (properTime - properDecayTime) / delTime;
       properTime = properDecayTime;
       decayed = true;
     }
  
     //
     // Compute the coordinates of the RawParticle after propagation
     //
     double xProp = particle_.X() + particle_.Px() * solution * factor;
     double yProp = particle_.Y() + particle_.Py() * solution * factor;
     zProp = particle_.Z() + particle_.Pz() * solution * factor;
     double tProp = particle_.T() + particle_.E() * solution * factor;
  
     //
     // Last check : Although propagated to the endcaps, R could still be
     // larger than rCyl because the cylinder was too short...
     //
     if (success == 2 && xProp * xProp + yProp * yProp > rCyl2) {
       success = -1;
       return true;
     }
  
     //
     // ... and update the particle with its propagated coordinates
     //
     particle_.setVertex(XYZTLorentzVector(xProp, yProp, zProp, tProp));
  
     return true;
   }
   //
   // Treat charged particles with magnetic field next.
   //
   else {
     //
     // First check that the particle can cross the cylinder
     //
     double pT = particle_.Pt();
     double radius = helixRadius(pT);
     double phi0 = helixStartPhi();
     double xC = helixCentreX(radius, phi0);
     double yC = helixCentreY(radius, phi0);
     double dist = helixCentreDistToAxis(xC, yC);
     //
     // The helix either is outside or includes the cylinder -> no intersection
     //
     if (fabs(fabs(radius) - dist) > rCyl)
       return false;
     //
     // The particle is already away from the endcaps, and will never come back
     // Could be back-propagated, so warn the user that it is so.
     //
     if (particle_.Z() * particle_.Pz() > zCyl * fabs(particle_.Pz())) {
       success = -2;
       return true;
     }
  
     double pZ = particle_.Pz();
     double phiProp, zProp;
     //
     // Does the helix cross the cylinder barrel ? If not, do the first half-loop
     //
     double rProp = std::min(fabs(radius) + dist - 0.000001, rCyl);
  
     // Check for rounding errors in the ArcSin.
     double sinPhiProp = (rProp * rProp - radius * radius - dist * dist) / (2. * dist * radius);
  
     //
     double deltaPhi = 1E99;
  
     // Taylor development up to third order for large momenta
     if (1. - fabs(sinPhiProp) < 1E-9) {
       double cphi0 = std::cos(phi0);
       double sphi0 = std::sin(phi0);
       double r0 = (particle_.X() * cphi0 + particle_.Y() * sphi0) / radius;
       double q0 = (particle_.X() * sphi0 - particle_.Y() * cphi0) / radius;
       double rcyl2 = (rCyl2 - particle_.X() * particle_.X() - particle_.Y() * particle_.Y()) / (radius * radius);
       double delta = r0 * r0 + rcyl2 * (1. - q0);
  
       // This is a solution of a second order equation, assuming phi = phi0 + epsilon
       // and epsilon is small.
       deltaPhi = radius > 0 ? (-r0 + std::sqrt(delta)) / (1. - q0) : (-r0 - std::sqrt(delta)) / (1. - q0);
     }
  
     // Use complete calculation otherwise, or if the delta phi is large anyway.
     if (fabs(deltaPhi) > 1E-3) {
       //    phiProp =  fabs(sinPhiProp) < 0.99999999  ?
       //      std::asin(sinPhiProp) : M_PI/2. * sinPhiProp;
       phiProp = std::asin(sinPhiProp);
  
       //
       // Compute phi after propagation (two solutions, but the asin definition
       // (between -pi/2 and +pi/2) takes care of the wrong one.
       //
       phiProp = helixCentrePhi(xC, yC) + (inside(rPos2) || onSurface(rPos2) ? phiProp : M_PI - phiProp);
  
       //
       // Solve the obvious two-pi ambiguities: more than one turn!
       //
       if (fabs(phiProp - phi0) > 2. * M_PI)
         phiProp += (phiProp > phi0 ? -2. * M_PI : +2. * M_PI);
  
       //
       // Check that the phi difference has the right sign, according to Q*B
       // (Another two pi ambuiguity, slightly more subtle)
       //
       if ((phiProp - phi0) * radius < 0.0)
         radius > 0.0 ? phiProp += 2 * M_PI : phiProp -= 2 * M_PI;
  
     } else {
       // Use Taylor
       phiProp = phi0 + deltaPhi;
     }
  
     //
     // Check that the particle does not exit from the endcaps first.
     //
     zProp = particle_.Z() + (phiProp - phi0) * pZ * radius / pT;
     if (fabs(zProp) > zCyl) {
       zProp = zCyl * fabs(pZ) / pZ;
       phiProp = phi0 + (zProp - particle_.Z()) / radius * pT / pZ;
       success = 2;
  
     } else {
       //
       // If the particle does not cross the barrel, either process the
       // first-half loop only or propagate all the way down to the endcaps
       //
       if (rProp < rCyl) {
         if (firstLoop || fabs(pZ) / pT < 1E-10) {
           success = 0;
  
         } else {
           zProp = zCyl * fabs(pZ) / pZ;
           phiProp = phi0 + (zProp - particle_.Z()) / radius * pT / pZ;
           success = 2;
         }
         //
         // The particle crossed the barrel
         //
       } else {
         success = 1;
       }
     }
     //
     // Compute the coordinates of the RawParticle after propagation
     //
     //
     // Check the decay time
     //
     double delTime = propDir * (phiProp - phi0) * radius * particle_.mass() / pT;
     double factor = 1.;
     properTime += delTime;
     if (properTime > properDecayTime) {
       factor = 1. - (properTime - properDecayTime) / delTime;
       properTime = properDecayTime;
       decayed = true;
     }
  
     zProp = particle_.Z() + (zProp - particle_.Z()) * factor;
  
     phiProp = phi0 + (phiProp - phi0) * factor;
  
     double sProp = std::sin(phiProp);
     double cProp = std::cos(phiProp);
     double xProp = xC + radius * sProp;
     double yProp = yC - radius * cProp;
     double tProp = particle_.T() + (phiProp - phi0) * radius * particle_.E() / pT;
     double pxProp = pT * cProp;
     double pyProp = pT * sProp;
  
     //
     // Last check : Although propagated to the endcaps, R could still be
     // larger than rCyl because the cylinder was too short...
     //
     if (success == 2 && xProp * xProp + yProp * yProp > rCyl2) {
       success = -1;
       return true;
     }
     //
     // ... and update the particle vertex and momentum
     //
     particle_.setVertex(XYZTLorentzVector(xProp, yProp, zProp, tProp));
     particle_.setMomentum(pxProp, pyProp, pZ, particle_.E());
     //    particle_.SetPx(pxProp);
     //    particle_.SetPy(pyProp);
     return true;
   }
 }

References b2, bField, RawParticle::charge(), funct::cos(), decayed, dumpMFGeometry_cfg::delta, SiPixelRawToDigiRegional_cfi::deltaPhi, RawParticle::E(), DQMScaleToClient_cfi::factor, firstLoop, helixCentreDistToAxis(), helixCentrePhi(), helixCentreX(), helixCentreY(), helixRadius(), helixStartPhi(), inside(), M_PI, RawParticle::mass(), min(), onBarrel(), onEndcap(), onSurface(), particle_, RawParticle::Perp2(), propDir, properDecayTime, properTime, PVValHelper::pT, RawParticle::Pt(), HLT_2018_cff::pT2, RawParticle::Px(), RawParticle::Py(), RawParticle::Pz(), SiStripMonitorCluster_cfi::q0, RawParticle::R2(), CosmicsPD_Skims::radius, rCyl, rCyl2, RawParticle::setMomentum(), RawParticle::setVertex(), funct::sin(), mathSSE::sqrt(), success, RawParticle::T(), RawParticle::X(), RawParticle::Y(), RawParticle::Z(), and zCyl.

Referenced by PFTrackTransformer::addPoints(), PFTrackTransformer::addPointsAndBrems(), backPropagate(), ConvBremSeedProducer::produce(), propagated(), ParticlePropagator::propagateToBoundSurface(), propagateToEcal(), propagateToEcalEntrance(), propagateToHcalEntrance(), propagateToHcalExit(), propagateToHOLayer(), propagateToPreshowerLayer1(), propagateToPreshowerLayer2(), and propagateToVFcalEntrance().

◆ propagated()

BaseParticlePropagator BaseParticlePropagator::propagated ( ) const

Return a new instance, corresponding to the particle propagated to the surface of the cylinder

Definition at line 304 of file BaseParticlePropagator.cc.

                                                                 {
   // Copy the input particle in a new instance
   BaseParticlePropagator myPropPart(*this);
   // Allow for many loops
   myPropPart.firstLoop = false;
   // Propagate the particle !
   myPropPart.propagate();
   // Back propagate if the forward propagation was not a success
   if (myPropPart.success < 0)
     myPropPart.backPropagate();
   // Return the new instance
   return myPropPart;
 }

References backPropagate(), firstLoop, propagate(), and success.

Referenced by ParticlePropagator::propagated().

◆ propagateToBeamCylinder()

bool BaseParticlePropagator::propagateToBeamCylinder	(	const XYZTLorentzVector &	v,
		double	radius = `0.`
	)

Definition at line 586 of file BaseParticlePropagator.cc.

                                                                                               {
   // For neutral (or BField = 0, simply check that the track passes through the cylinder
   if (particle_.charge() == 0. || bField == 0.)
     return fabs(particle_.Px() * particle_.Y() - particle_.Py() * particle_.X()) / particle_.Pt() < radius;
  
   // Now go to the charged particles
  
   // A few needed intermediate variables (to make the code understandable)
   // The inner cylinder
   double r = radius;
   double r2 = r * r;
   double r4 = r2 * r2;
   // The two hits
   double dx = particle_.X() - v.X();
   double dy = particle_.Y() - v.Y();
   double dz = particle_.Z() - v.Z();
   double Sxy = particle_.X() * v.X() + particle_.Y() * v.Y();
   double Dxy = particle_.Y() * v.X() - particle_.X() * v.Y();
   double Dxy2 = Dxy * Dxy;
   double Sxy2 = Sxy * Sxy;
   double SDxy = dx * dx + dy * dy;
   double SSDxy = std::sqrt(SDxy);
   double ATxy = std::atan2(dy, dx);
   // The coefficients of the second order equation for the trajectory radius
   double a = r2 - Dxy2 / SDxy;
   double b = r * (r2 - Sxy);
   double c = r4 - 2. * Sxy * r2 + Sxy2 + Dxy2;
  
   // Here are the four possible solutions for
   // 1) The trajectory radius
   std::array<double, 4> helixRs = {{0.}};
   helixRs[0] = (b - std::sqrt(b * b - a * c)) / (2. * a);
   helixRs[1] = (b + std::sqrt(b * b - a * c)) / (2. * a);
   helixRs[2] = -helixRs[0];
   helixRs[3] = -helixRs[1];
   // 2) The azimuthal direction at the second point
   std::array<double, 4> helixPhis = {{0.}};
   helixPhis[0] = std::asin(SSDxy / (2. * helixRs[0]));
   helixPhis[1] = std::asin(SSDxy / (2. * helixRs[1]));
   helixPhis[2] = -helixPhis[0];
   helixPhis[3] = -helixPhis[1];
   // Only two solutions are valid though
   double solution1 = 0.;
   double solution2 = 0.;
   double phi1 = 0.;
   double phi2 = 0.;
   // Loop over the four possibilities
   for (unsigned int i = 0; i < 4; ++i) {
     helixPhis[i] = ATxy + helixPhis[i];
     // The centre of the helix
     double xC = helixCentreX(helixRs[i], helixPhis[i]);
     double yC = helixCentreY(helixRs[i], helixPhis[i]);
     // The distance of that centre to the origin
     double dist = helixCentreDistToAxis(xC, yC);
     /*
     std::cout << "Solution "<< i << " = " << helixRs[i] << " accepted ? " 
           << fabs(fabs(helixRs[i]) - dist ) << " - " << radius << " = " 
           << fabs(fabs(helixRs[i]) - dist ) - fabs(radius) << " " 
           << ( fabs( fabs(fabs(helixRs[i]) - dist) -fabs(radius) ) < 1e-6)  << std::endl;
     */
     // Check that the propagation will lead to a trajectroy tangent to the inner cylinder
     if (fabs(fabs(fabs(helixRs[i]) - dist) - fabs(radius)) < 1e-6) {
       // Fill first solution
       if (solution1 == 0.) {
         solution1 = helixRs[i];
         phi1 = helixPhis[i];
         // Fill second solution (ordered)
       } else if (solution2 == 0.) {
         if (helixRs[i] < solution1) {
           solution2 = solution1;
           solution1 = helixRs[i];
           phi2 = phi1;
           phi1 = helixPhis[i];
         } else {
           solution2 = helixRs[i];
           phi2 = helixPhis[i];
         }
         // Must not happen!
       } else {
         std::cout << "warning !!! More than two solutions for this track !!! " << std::endl;
       }
     }
   }
  
   // Find the largest possible pT compatible with coming from within the inne cylinder
   double pT = 0.;
   double helixPhi = 0.;
   double helixR = 0.;
   // This should not happen
   if (solution1 * solution2 == 0.) {
     std::cout << "warning !!! Less than two solution for this track! " << solution1 << " " << solution2 << std::endl;
     return false;
     // If the two solutions have opposite sign, it means that pT = infinity is a possibility.
   } else if (solution1 * solution2 < 0.) {
     pT = 1000.;
     double norm = pT / SSDxy;
     particle_.setCharge(+1.);
     particle_.setMomentum(dx * norm, dy * norm, dz * norm, 0.);
     particle_.SetE(std::sqrt(particle_.momentum().Vect().Mag2()));
     // Otherwise take the solution that gives the largest transverse momentum
   } else {
     if (solution1 < 0.) {
       helixR = solution1;
       helixPhi = phi1;
       particle_.setCharge(+1.);
     } else {
       helixR = solution2;
       helixPhi = phi2;
       particle_.setCharge(-1.);
     }
     pT = fabs(helixR) * 1e-5 * c_light() * bField;
     double norm = pT / SSDxy;
     particle_.setMomentum(pT * std::cos(helixPhi), pT * std::sin(helixPhi), dz * norm, 0.);
     particle_.SetE(std::sqrt(particle_.momentum().Vect().Mag2()));
   }
  
   // Propagate to closest approach to get the Z value (a bit of an overkill)
   return propagateToClosestApproach();
 }

References a, b, bField, HltBtagPostValidation_cff::c, c_light(), RawParticle::charge(), funct::cos(), gather_cfg::cout, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, MillePedeFileConverter_cfg::e, helixCentreDistToAxis(), helixCentreX(), helixCentreY(), mps_fire::i, RawParticle::momentum(), particle_, propagateToClosestApproach(), PVValHelper::pT, RawParticle::Pt(), RawParticle::Px(), RawParticle::Py(), alignCSCRings::r, diffTwoXMLs::r2, CosmicsPD_Skims::radius, RawParticle::setCharge(), RawParticle::SetE(), RawParticle::setMomentum(), funct::sin(), mathSSE::sqrt(), findQualityFiles::v, RawParticle::X(), RawParticle::Y(), and RawParticle::Z().

◆ propagateToClosestApproach()

bool BaseParticlePropagator::propagateToClosestApproach	(	double	x0 = `0.`,
		double	y0 = `0`,
		bool	first = `true`
	)

Update the particle after propagation to the closest approach from Z axis, to the preshower layer 1 & 2, to the ECAL entrance, to the HCAL entrance, the HCAL 2nd and 3rd layer (not coded yet), the VFCAL entrance, or any BoundSurface(disk or cylinder)

Definition at line 325 of file BaseParticlePropagator.cc.

                                                                                         {
   // Pre-computed quantities
   double pT = particle_.Pt();
   double radius = helixRadius(pT);
   double phi0 = helixStartPhi();
  
   // The Helix centre coordinates are computed wrt to the z axis
   // Actually the z axis might be centred on 0,0: it's the beam spot position (x0,y0)!
   double xC = helixCentreX(radius, phi0);
   double yC = helixCentreY(radius, phi0);
   double distz = helixCentreDistToAxis(xC - x0, yC - y0);
   double dist0 = helixCentreDistToAxis(xC, yC);
  
   //
   // Propagation to point of clostest approach to z axis
   //
   double rz, r0, z;
   if (particle_.charge() != 0.0 && bField != 0.0) {
     rz = fabs(fabs(radius) - distz) + std::sqrt(x0 * x0 + y0 * y0) + 0.0000001;
     r0 = fabs(fabs(radius) - dist0) + 0.0000001;
   } else {
     rz = fabs(particle_.Px() * (particle_.Y() - y0) - particle_.Py() * (particle_.X() - x0)) / particle_.Pt();
     r0 = fabs(particle_.Px() * particle_.Y() - particle_.Py() * particle_.X()) / particle_.Pt();
   }
  
   z = 999999.;
  
   // Propagate to the first interesection point
   // with cylinder of radius sqrt(x0**2+y0**2)
   setPropagationConditions(rz, z, first);
   bool done = backPropagate();
  
   // Unsuccessful propagation - should not happen!
   if (!done)
     return done;
  
   // The z axis is (0,0) - no need to go further
   if (fabs(rz - r0) < 1E-10)
     return done;
   double dist1 = (particle_.X() - x0) * (particle_.X() - x0) + (particle_.Y() - y0) * (particle_.Y() - y0);
  
   // We are already at closest approach - no need to go further
   if (dist1 < 1E-10)
     return done;
  
   // Keep for later if it happens to be the right solution
   XYZTLorentzVector vertex1 = particle_.vertex();
   XYZTLorentzVector momentum1 = particle_.momentum();
  
   // Propagate to the distance of closest approach to (0,0)
   setPropagationConditions(r0, z, first);
   done = backPropagate();
   if (!done)
     return done;
  
   // Propagate to the first interesection point
   // with cylinder of radius sqrt(x0**2+y0**2)
   setPropagationConditions(rz, z, first);
   done = backPropagate();
   if (!done)
     return done;
   double dist2 = (particle_.X() - x0) * (particle_.X() - x0) + (particle_.Y() - y0) * (particle_.Y() - y0);
  
   // Keep the good solution.
   if (dist2 > dist1) {
     particle_.setVertex(vertex1);
     particle_.setMomentum(momentum1.X(), momentum1.Y(), momentum1.Z(), momentum1.E());
   }
  
   // Done
   return done;
 }

References backPropagate(), bField, RawParticle::charge(), fileCollector::done, dqmdumpme::first, helixCentreDistToAxis(), helixCentreX(), helixCentreY(), helixRadius(), helixStartPhi(), RawParticle::momentum(), particle_, PVValHelper::pT, RawParticle::Pt(), RawParticle::Px(), RawParticle::Py(), CosmicsPD_Skims::radius, RawParticle::setMomentum(), setPropagationConditions(), RawParticle::setVertex(), mathSSE::sqrt(), RawParticle::vertex(), RawParticle::X(), RawParticle::Y(), and z.

Referenced by propagateToBeamCylinder(), ParticlePropagator::propagateToClosestApproach(), and propagateToNominalVertex().

◆ propagateToEcal()

bool BaseParticlePropagator::propagateToEcal ( bool first = true )

Definition at line 398 of file BaseParticlePropagator.cc.

                                                        {
   //
   // Propagation to Ecal (including preshower) after the
   // last Tracker layer
   // TODO: include proper geometry
   // Geometry taken from CMS ECAL TDR
   //
   // First propagate to global barrel / endcap cylinder
   //  setPropagationConditions(1290. , 3045 , first);
   //  New position of the preshower as of CMSSW_1_7_X
   setPropagationConditions(129.0, 303.353, first);
   return propagate();
 }

References dqmdumpme::first, propagate(), and setPropagationConditions().

◆ propagateToEcalEntrance()

bool BaseParticlePropagator::propagateToEcalEntrance ( bool first = true )

Definition at line 450 of file BaseParticlePropagator.cc.

                                                                {
   //
   // Propagation to ECAL entrance
   // TODO: include proper geometry
   // Geometry taken from CMS ECAL TDR
   //
   // First propagate to global barrel / endcap cylinder
   setPropagationConditions(129.0, 320.9, first);
   bool done = propagate();
  
   // Go to endcap cylinder in the "barrel cut corner"
   // eta = 1.479 -> cos^2(theta) = 0.81230
   //  if ( done && eta > 1.479 && success == 1 ) {
   if (done && particle_.cos2ThetaV() > 0.81230 && success == 1) {
     setPropagationConditions(152.6, 320.9, first);
     done = propagate();
   }
  
   // We are not in the ECAL acceptance
   // eta = 3.0 -> cos^2(theta) = 0.99013
   if (particle_.cos2ThetaV() > 0.99014)
     success = 0;
  
   return done;
 }

References RawParticle::cos2ThetaV(), fileCollector::done, dqmdumpme::first, particle_, propagate(), setPropagationConditions(), and success.

Referenced by PFTrackTransformer::addPointsAndBrems(), reco::tau::atECALEntrance(), FBaseSimEvent::fill(), ConvBremSeedProducer::GoodCluster(), reco::tau::PFRecoTauChargedHadronFromGenericTrackPlugin< TrackClass >::operator()(), GoodSeedProducer::produce(), l1tpf::propagateToCalo(), TrajectoryManager::propagateToCalorimeters(), and LowPtGsfElectronSeedProducer::propagateTrackToCalo().

◆ propagateToHcalEntrance()

bool BaseParticlePropagator::propagateToHcalEntrance ( bool first = true )

Definition at line 476 of file BaseParticlePropagator.cc.

                                                                {
   //
   // Propagation to HCAL entrance
   // TODO: include proper geometry
   // Geometry taken from CMSSW_3_1_X xml (Sunanda)
   //
  
   // First propagate to global barrel / endcap cylinder
   setPropagationConditions(177.5, 335.0, first);
   propDir = 0;
   bool done = propagate();
   propDir = 1;
  
   // If went through the bottom of HB cylinder -> re-propagate to HE surface
   if (done && success == 2) {
     setPropagationConditions(300.0, 400.458, first);
     propDir = 0;
     done = propagate();
     propDir = 1;
   }
  
   // out of the HB/HE acceptance
   // eta = 3.0 -> cos^2(theta) = 0.99014
   if (done && particle_.cos2ThetaV() > 0.99014)
     success = 0;
  
   return done;
 }

References RawParticle::cos2ThetaV(), fileCollector::done, dqmdumpme::first, particle_, propagate(), propDir, setPropagationConditions(), and success.

Referenced by PFTrackTransformer::addPointsAndBrems(), FBaseSimEvent::fill(), and TrajectoryManager::propagateToCalorimeters().

◆ propagateToHcalExit()

bool BaseParticlePropagator::propagateToHcalExit ( bool first = true )

Definition at line 528 of file BaseParticlePropagator.cc.

                                                            {
   //
   // Propagation to HCAL exit
   // TODO: include proper geometry
   // Geometry taken from CMSSW_3_1_X xml (Sunanda)
   //
  
   // Approximate it to a single cylinder as it is not that crucial.
   setPropagationConditions(263.9, 554.1, first);
   //  this->rawPart().particle_.setCharge(0.0); ?? Shower Propagation ??
   propDir = 0;
   bool done = propagate();
   propDir = 1;
  
   return done;
 }

References fileCollector::done, dqmdumpme::first, propagate(), propDir, and setPropagationConditions().

Referenced by PFTrackTransformer::addPointsAndBrems(), and FBaseSimEvent::fill().

◆ propagateToHOLayer()

bool BaseParticlePropagator::propagateToHOLayer ( bool first = true )

Definition at line 545 of file BaseParticlePropagator.cc.

                                                           {
   //
   // Propagation to Layer0 of HO (Ring-0)
   // TODO: include proper geometry
   // Geometry taken from CMSSW_3_1_X xml (Sunanda)
   //
  
   // Approximate it to a single cylinder as it is not that crucial.
   setPropagationConditions(387.6, 1000.0, first);  //Dia is the middle of HO \sqrt{384.8**2+46.2**2}
   //  setPropagationConditions(410.2, 1000.0, first); //sqrt{406.6**2+54.2**2} //for outer layer
   //  this->rawPart().setCharge(0.0); ?? Shower Propagation ??
   propDir = 0;
   bool done = propagate();
   propDir = 1;
  
   if (done && std::abs(particle_.Z()) > 700.25)
     success = 0;
  
   return done;
 }

References funct::abs(), fileCollector::done, dqmdumpme::first, particle_, propagate(), propDir, setPropagationConditions(), success, and RawParticle::Z().

Referenced by PFTrackTransformer::addPointsAndBrems(), and FBaseSimEvent::fill().

◆ propagateToNominalVertex()

bool BaseParticlePropagator::propagateToNominalVertex ( const XYZTLorentzVector & hit2 = XYZTLorentzVector(0., 0., 0., 0.) )

Definition at line 566 of file BaseParticlePropagator.cc.

                                                                                 {
   // Not implemented for neutrals (used for electrons only)
   if (particle_.charge() == 0. || bField == 0.)
     return false;
  
   // Define the proper pT direction to meet the point (vx,vy) in (x,y)
   double dx = particle_.X() - v.X();
   double dy = particle_.Y() - v.Y();
   double phi = std::atan2(dy, dx) + std::asin(std::sqrt(dx * dx + dy * dy) / (2. * helixRadius()));
  
   // The absolute pT (kept to the original value)
   double pT = particle_.pt();
  
   // Set the new pT
   particle_.SetPx(pT * std::cos(phi));
   particle_.SetPy(pT * std::sin(phi));
  
   return propagateToClosestApproach(v.X(), v.Y());
 }

References bField, RawParticle::charge(), funct::cos(), PVValHelper::dx, PVValHelper::dy, helixRadius(), particle_, phi, propagateToClosestApproach(), PVValHelper::pT, RawParticle::pt(), RawParticle::SetPx(), RawParticle::SetPy(), funct::sin(), mathSSE::sqrt(), findQualityFiles::v, RawParticle::X(), and RawParticle::Y().

Referenced by ParticlePropagator::propagateToNominalVertex().

◆ propagateToPreshowerLayer1()

bool BaseParticlePropagator::propagateToPreshowerLayer1 ( bool first = true )

Definition at line 412 of file BaseParticlePropagator.cc.

                                                                   {
   //
   // Propagation to Preshower Layer 1
   // TODO: include proper geometry
   // Geometry taken from CMS ECAL TDR
   //
   // First propagate to global barrel / endcap cylinder
   //  setPropagationConditions(1290., 3045 , first);
   //  New position of the preshower as of CMSSW_1_7_X
   setPropagationConditions(129.0, 303.353, first);
   bool done = propagate();
  
   // Check that were are on the Layer 1
   if (done && (particle_.R2() > 125.0 * 125.0 || particle_.R2() < 45.0 * 45.0))
     success = 0;
  
   return done;
 }

References fileCollector::done, dqmdumpme::first, particle_, propagate(), RawParticle::R2(), setPropagationConditions(), and success.

Referenced by PFTrackTransformer::addPointsAndBrems(), FBaseSimEvent::fill(), and TrajectoryManager::propagateToCalorimeters().

◆ propagateToPreshowerLayer2()

bool BaseParticlePropagator::propagateToPreshowerLayer2 ( bool first = true )

Definition at line 431 of file BaseParticlePropagator.cc.

                                                                   {
   //
   // Propagation to Preshower Layer 2
   // TODO: include proper geometry
   // Geometry taken from CMS ECAL TDR
   //
   // First propagate to global barrel / endcap cylinder
   //  setPropagationConditions(1290. , 3090 , first);
   //  New position of the preshower as of CMSSW_1_7_X
   setPropagationConditions(129.0, 307.838, first);
   bool done = propagate();
  
   // Check that we are on Layer 2
   if (done && (particle_.R2() > 125.0 * 125.0 || particle_.R2() < 45.0 * 45.0))
     success = 0;
  
   return done;
 }

References fileCollector::done, dqmdumpme::first, particle_, propagate(), RawParticle::R2(), setPropagationConditions(), and success.

Referenced by PFTrackTransformer::addPointsAndBrems(), FBaseSimEvent::fill(), and TrajectoryManager::propagateToCalorimeters().

◆ propagateToVFcalEntrance()

bool BaseParticlePropagator::propagateToVFcalEntrance ( bool first = true )

Definition at line 505 of file BaseParticlePropagator.cc.

                                                                 {
   //
   // Propagation to VFCAL (HF) entrance
   // Geometry taken from xml: Geometry/CMSCommonData/data/cms.xml
  
   setPropagationConditions(400.0, 1114.95, first);
   propDir = 0;
   bool done = propagate();
   propDir = 1;
  
   if (!done)
     success = 0;
  
   // We are not in the VFCAL acceptance
   // eta = 3.0  -> cos^2(theta) = 0.99014
   // eta = 5.2  -> cos^2(theta) = 0.9998755
   double c2teta = particle_.cos2ThetaV();
   if (done && (c2teta < 0.99014 || c2teta > 0.9998755))
     success = 0;
  
   return done;
 }

References RawParticle::cos2ThetaV(), fileCollector::done, dqmdumpme::first, particle_, propagate(), propDir, setPropagationConditions(), and success.

Referenced by FBaseSimEvent::fill(), and TrajectoryManager::propagateToCalorimeters().

◆ resetDebug()

void BaseParticlePropagator::resetDebug ( )

inline

Definition at line 306 of file BaseParticlePropagator.h.

306 { debug = false; }

References debug.

◆ setDebug()

void BaseParticlePropagator::setDebug ( )

inline

Set the debug leve;.

Definition at line 305 of file BaseParticlePropagator.h.

305 { debug = true; }

References debug.

◆ setMagneticField()

void BaseParticlePropagator::setMagneticField ( double b )

inline

Set the magnetic field.

Definition at line 299 of file BaseParticlePropagator.h.

299 { bField = b; }

References b, and bField.

Referenced by PFTrackTransformer::addPointsAndBrems(), FBaseSimEvent::fill(), ParticlePropagator::ParticlePropagator(), ParticlePropagator::propagateToBoundSurface(), ParticlePropagator::propagateToClosestApproach(), and ParticlePropagator::propagateToNominalVertex().

◆ setParticle()

void BaseParticlePropagator::setParticle ( RawParticle const & iParticle )

inline

Definition at line 166 of file BaseParticlePropagator.h.

166 { particle_ = iParticle; }

References particle_.

◆ setPropagationConditions()

void BaseParticlePropagator::setPropagationConditions	(	double	r,
		double	z,
		bool	firstLoop = `true`
	)

Set the propagation characteristics (rCyl, zCyl and first loop only)

Definition at line 318 of file BaseParticlePropagator.cc.

                                                                                 {
   rCyl = R;
   rCyl2 = R * R;
   zCyl = Z;
   firstLoop = f;
 }

References f, firstLoop, dttmaxenums::R, rCyl, rCyl2, DOFs::Z, and zCyl.

Referenced by PFTrackTransformer::addPoints(), PFTrackTransformer::addPointsAndBrems(), propagateToClosestApproach(), propagateToEcal(), propagateToEcalEntrance(), propagateToHcalEntrance(), propagateToHcalExit(), propagateToHOLayer(), propagateToPreshowerLayer1(), propagateToPreshowerLayer2(), propagateToVFcalEntrance(), and ParticlePropagator::setPropagationConditions().

◆ setProperDecayTime()

void BaseParticlePropagator::setProperDecayTime ( double t )

inline

Set the proper decay time.

Definition at line 169 of file BaseParticlePropagator.h.

169 { properDecayTime = t; }

References properDecayTime, and OrderedSet::t.

Referenced by ParticlePropagator::initProperDecayTime(), and ParticlePropagator::ParticlePropagator().

◆ xyImpactParameter()

double BaseParticlePropagator::xyImpactParameter	(	double	x0 = `0.`,
		double	y0 = `0.`
	)		const

Transverse impact parameter.

Definition at line 706 of file BaseParticlePropagator.cc.

                                                                            {
   double ip = 0.;
   double pT = particle_.Pt();
  
   if (particle_.charge() != 0.0 && bField != 0.0) {
     double radius = helixRadius(pT);
     double phi0 = helixStartPhi();
  
     // The Helix centre coordinates are computed wrt to the z axis
     double xC = helixCentreX(radius, phi0);
     double yC = helixCentreY(radius, phi0);
     double distz = helixCentreDistToAxis(xC - x0, yC - y0);
     ip = distz - fabs(radius);
   } else {
     ip = fabs(particle_.Px() * (particle_.Y() - y0) - particle_.Py() * (particle_.X() - x0)) / pT;
   }
  
   return ip;
 }

References bField, RawParticle::charge(), helixCentreDistToAxis(), helixCentreX(), helixCentreY(), helixRadius(), helixStartPhi(), particle_, PVValHelper::pT, RawParticle::Pt(), RawParticle::Px(), RawParticle::Py(), CosmicsPD_Skims::radius, RawParticle::X(), and RawParticle::Y().

◆ zImpactParameter()

double BaseParticlePropagator::zImpactParameter	(	double	x0 = `0`,
		double	y0 = `0.`
	)		const

inline

Longitudinal impact parameter.

Definition at line 181 of file BaseParticlePropagator.h.

                                                                       {
     // Longitudinal impact parameter
     return particle_.Z() - particle_.Pz() * std::sqrt(((particle_.X() - x0) * (particle_.X() - x0) +
                                                        (particle_.Y() - y0) * (particle_.Y() - y0)) /
                                                       particle_.Perp2());
   }

References particle_, RawParticle::Perp2(), RawParticle::Pz(), mathSSE::sqrt(), RawParticle::X(), RawParticle::Y(), and RawParticle::Z().

Member Data Documentation

◆ bField

double BaseParticlePropagator::bField

private

Magnetic field in the cylinder, oriented along the Z axis.

Definition at line 137 of file BaseParticlePropagator.h.

Referenced by getMagneticField(), helixRadius(), propagate(), propagateToBeamCylinder(), propagateToClosestApproach(), propagateToNominalVertex(), setMagneticField(), and xyImpactParameter().

◆ debug

bool BaseParticlePropagator::debug

private

The debug level.

Definition at line 141 of file BaseParticlePropagator.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), init(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), resetDebug(), runTauIdMVA.TauIDEmbedder::runTauID(), setDebug(), and pkg.AbstractPkg::write().

◆ decayed

bool BaseParticlePropagator::decayed

private

The particle decayed while propagated !

Definition at line 156 of file BaseParticlePropagator.h.

Referenced by hasDecayed(), init(), and propagate().

◆ fiducial

bool BaseParticlePropagator::fiducial

protected

The particle traverses some real material.

Definition at line 147 of file BaseParticlePropagator.h.

Referenced by init(), onFiducial(), and ParticlePropagator::propagateToBoundSurface().

◆ firstLoop

bool BaseParticlePropagator::firstLoop

private

Do only the first half-loop.

Definition at line 154 of file BaseParticlePropagator.h.

Referenced by init(), propagate(), propagated(), ParticlePropagator::setPropagationConditions(), and setPropagationConditions().

◆ particle_

RawParticle BaseParticlePropagator::particle_

private

◆ propDir

int BaseParticlePropagator::propDir

private

The propagation direction.

Definition at line 160 of file BaseParticlePropagator.h.

Referenced by backPropagate(), init(), propagate(), propagateToHcalEntrance(), propagateToHcalExit(), propagateToHOLayer(), and propagateToVFcalEntrance().

◆ properDecayTime

double BaseParticlePropagator::properDecayTime

private

The proper decay time of the particle.

Definition at line 139 of file BaseParticlePropagator.h.

Referenced by ParticlePropagator::initProperDecayTime(), propagate(), and setProperDecayTime().

◆ properTime

double BaseParticlePropagator::properTime

private

The proper time of the particle.

Definition at line 158 of file BaseParticlePropagator.h.

Referenced by init(), and propagate().

◆ rCyl

double BaseParticlePropagator::rCyl

private

Simulated particle that is to be resp has been propagated.

Radius of the cylinder (centred at 0,0,0) to which propagation is done

Definition at line 132 of file BaseParticlePropagator.h.

Referenced by increaseRCyl(), inside(), onBarrel(), propagate(), and setPropagationConditions().

◆ rCyl2

double BaseParticlePropagator::rCyl2

private

Definition at line 133 of file BaseParticlePropagator.h.

Referenced by increaseRCyl(), inside(), onBarrel(), onEndcap(), propagate(), and setPropagationConditions().

◆ success

int BaseParticlePropagator::success

protected

0:propagation still be done, 1:reached 'barrel', 2:reached 'endcaps'

Definition at line 145 of file BaseParticlePropagator.h.

Referenced by ParticlePropagator::fieldMap(), getSuccess(), init(), propagate(), propagated(), ParticlePropagator::propagateToBoundSurface(), propagateToEcalEntrance(), propagateToHcalEntrance(), propagateToHOLayer(), propagateToPreshowerLayer1(), propagateToPreshowerLayer2(), and propagateToVFcalEntrance().

◆ zCyl

double BaseParticlePropagator::zCyl

private

Half-height of the cylinder (centred at 0,0,0) to which propagation is done.

Definition at line 135 of file BaseParticlePropagator.h.

Referenced by inside(), onBarrel(), onEndcap(), propagate(), and setPropagationConditions().

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BaseParticlePropagator() [1/3]

◆ BaseParticlePropagator() [2/3]

◆ BaseParticlePropagator() [3/3]

Member Function Documentation

◆ backPropagate()

◆ c_light()

◆ getMagneticField()

◆ getSuccess()

◆ hasDecayed()

◆ helixCentreDistToAxis() [1/2]

◆ helixCentreDistToAxis() [2/2]

◆ helixCentrePhi() [1/2]

◆ helixCentrePhi() [2/2]

◆ helixCentreX() [1/2]

◆ helixCentreX() [2/2]

◆ helixCentreY() [1/2]

◆ helixCentreY() [2/2]

◆ helixRadius() [1/2]

◆ helixRadius() [2/2]

◆ helixStartPhi()

◆ increaseRCyl()

◆ init()

◆ inside() [1/2]

◆ inside() [2/2]

◆ onBarrel() [1/2]

◆ onBarrel() [2/2]

◆ onEndcap() [1/2]

◆ onEndcap() [2/2]

◆ onFiducial()

◆ onSurface() [1/2]

◆ onSurface() [2/2]

◆ particle() [1/2]

◆ particle() [2/2]

◆ propagate()

◆ propagated()

◆ propagateToBeamCylinder()

◆ propagateToClosestApproach()

◆ propagateToEcal()

◆ propagateToEcalEntrance()

◆ propagateToHcalEntrance()

◆ propagateToHcalExit()

◆ propagateToHOLayer()

◆ propagateToNominalVertex()

◆ propagateToPreshowerLayer1()

◆ propagateToPreshowerLayer2()

◆ propagateToVFcalEntrance()

◆ resetDebug()

◆ setDebug()

◆ setMagneticField()

◆ setParticle()

◆ setPropagationConditions()

◆ setProperDecayTime()

◆ xyImpactParameter()

◆ zImpactParameter()

Member Data Documentation

◆ bField

◆ debug

◆ decayed

◆ fiducial

◆ firstLoop

◆ particle_

◆ propDir

◆ properDecayTime

◆ properTime

◆ rCyl

◆ rCyl2

◆ success

◆ zCyl