Base class for particle flow input reconstructed tracks and simulated particles. More...

#include <PFTrack.h>

Inheritance diagram for reco::PFTrack:

Public Member Functions
void	addPoint (const reco::PFTrajectoryPoint &trajPt)

void	calculatePositionREP ()

double	charge () const

int	color () const

const reco::PFTrajectoryPoint &	extrapolatedPoint (unsigned layerid) const

std::vector< reco::PFTrajectoryPoint >::const_iterator	innermostMeasurement () const
	iterator on innermost tracker measurement More...

unsigned int	nTrajectoryMeasurements () const

unsigned int	nTrajectoryPoints () const

std::vector< reco::PFTrajectoryPoint >::const_iterator	outermostMeasurement () const
	iterator on outermost tracker measurement More...

	PFTrack ()

	PFTrack (double charge)

	PFTrack (const PFTrack &other)

void	setColor (int color)

void	setPoint (unsigned int index, const reco::PFTrajectoryPoint &measurement)
	set a trajectory point More...

const reco::PFTrajectoryPoint &	trajectoryPoint (unsigned index) const

const std::vector< reco::PFTrajectoryPoint > &	trajectoryPoints () const

Protected Attributes
double	charge_
	charge More...

int	color_
	color (transient) More...

unsigned int	indexInnermost_
	index innermost tracker measurement More...

unsigned int	indexOutermost_
	index outermost tracker measurement More...

std::vector< reco::PFTrajectoryPoint >	trajectoryPoints_
	vector of trajectory points More...

Static Protected Attributes
static const unsigned int	nMaxTrackingLayers_ = 17
	maximal number of tracking layers More...

Detailed Description

Base class for particle flow input reconstructed tracks and simulated particles.

A PFTrack contains a vector of PFTrajectoryPoint objects. These points are stored in a vector to benefit from the random access. One must take care of storing the points in the right order, and it might even be necessary to insert dummy points.

For a PFRecTrack, the ordering of the points is the following:

point 1: Closest approach
point 2: Beam Pipe
point 3 to n: Trajectory measurements (from the tracker layers)
point n+1: Preshower Layer 1, or dummy point if not in the preshower zone
point n+2: Preshower Layer 2, or dummy point if not in the preshower zone
point n+3: ECAL Entrance
point n+4: ECAL expected position of the shower maximum, assuming the track is an electron track.
point n+5: HCAL Entrance

For a PFSimParticle, the ordering of the points is the following.

If the particle decays before ECAL:
point 1: start point
point 2: end point
If the particle does not decay before ECAL:
point 1: start point
point 2: PS1 or dummy
point 3: PS2 or dummy
point 4: ECAL entrance
point 5: HCAL entrance

PFRecTracks and PFSimParticles are created in the PFTrackProducer module.

Author: Renaud Bruneliere

Date: July 2006

Definition at line 63 of file PFTrack.h.

Constructor & Destructor Documentation

PFTrack::PFTrack ( )

Definition at line 11 of file PFTrack.cc.

References reco::PFTrajectoryPoint::NLayers, nMaxTrackingLayers_, and trajectoryPoints_.

                  :
   charge_(0.),
   indexInnermost_(0),
   indexOutermost_(0),
   color_(1) {
 
   // prepare vector of trajectory points for propagated positions
   trajectoryPoints_.reserve(PFTrajectoryPoint::NLayers + nMaxTrackingLayers_);
 }

PFTrack::PFTrack ( double charge )

Definition at line 22 of file PFTrack.cc.

References reco::PFTrajectoryPoint::NLayers, nMaxTrackingLayers_, and trajectoryPoints_.

                               : 
   charge_(charge), 
   indexInnermost_(0),
   indexOutermost_(0),
   color_(1) {
 
   // prepare vector of trajectory points for propagated positions
   trajectoryPoints_.reserve(PFTrajectoryPoint::NLayers + nMaxTrackingLayers_);
 }

PFTrack::PFTrack ( const PFTrack & other )

Definition at line 33 of file PFTrack.cc.

                                      :
   charge_(other.charge_), 
   trajectoryPoints_(other.trajectoryPoints_),
   indexInnermost_(other.indexInnermost_),
   indexOutermost_(other.indexOutermost_),
   color_(other.color_)
 {}

Member Function Documentation

void PFTrack::addPoint ( const reco::PFTrajectoryPoint & trajPt )

add a trajectory measurement

Definition at line 42 of file PFTrack.cc.

References reco::PFTrajectoryPoint::BeamPipeOrEndVertex, indexInnermost_, indexOutermost_, reco::PFTrajectoryPoint::isTrackerLayer(), and trajectoryPoints_.

Referenced by PFTrackTransformer::addPoints(), and PFTrackTransformer::addPointsAndBrems().

                                                       {
   
   //   cout<<"adding "<<trajPt<<endl;
 
   if (trajPt.isTrackerLayer()) {
     if (!indexOutermost_) { // first time a measurement is added
       if (trajectoryPoints_.size() < PFTrajectoryPoint::BeamPipeOrEndVertex + 1) {
         PFTrajectoryPoint dummyPt;
         for (unsigned iPt = trajectoryPoints_.size(); iPt < PFTrajectoryPoint::BeamPipeOrEndVertex + 1; iPt++)
           trajectoryPoints_.push_back(dummyPt);
       } else if (trajectoryPoints_.size() > PFTrajectoryPoint::BeamPipeOrEndVertex + 1) {
         // throw an exception here
         //      edm::LogError("PFTrack")<<"trajectoryPoints_.size() is too large = " 
         //                              <<trajectoryPoints_.size()<<"\n";
       }
       indexOutermost_ = indexInnermost_ = PFTrajectoryPoint::BeamPipeOrEndVertex + 1;
     } else 
       indexOutermost_++;
   }
   // Use push_back instead of insert in order to gain time
   trajectoryPoints_.push_back(trajPt);
 
   //   cout<<"adding point "<<*this<<endl;
 }

void PFTrack::calculatePositionREP ( )

calculate posrep_ once and for all for each point

Definition at line 68 of file PFTrack.cc.

Referenced by PFTrackTransformer::addPointsAndBrems(), reco::GsfPFRecTrack::calculateBremPositionREP(), ConvBremPFTrackFinder::runConvBremFinder(), and setPoint().

                                    {
   
   //for(unsigned i=0; i<trajectoryPoints_.size(); i++) {
   //  trajectoryPoints_[i].calculatePositionREP();
   //}
 }

double reco::PFTrack::charge ( ) const

inline

Returns: electric charge

Definition at line 87 of file PFTrack.h.

References charge_.

Referenced by extrapolatedPoint(), reco::GsfPFRecTrack::GsfPFRecTrack(), reco::operator<<(), reco::PFSimParticle::PFSimParticle(), and reco::PFRecTrack::~PFRecTrack().

87 { return charge_; }

reco::PFTrack::charge_

double charge_

charge

Definition: PFTrack.h:130

int reco::PFTrack::color ( ) const

inline

Definition at line 122 of file PFTrack.h.

References color_.

Referenced by setColor().

122 { return color_; }

reco::PFTrack::color_

int color_

color (transient)

Definition: PFTrack.h:142

const reco::PFTrajectoryPoint & PFTrack::extrapolatedPoint ( unsigned layerid ) const

Returns: an extrapolated point

Definition at line 76 of file PFTrack.cc.

References charge(), reco::PFTrajectoryPoint::ClosestApproach, Exception, i, indexInnermost_, reco::PFTrajectoryPoint::momentum(), reco::PFTrajectoryPoint::NLayers, nTrajectoryMeasurements(), reco::operator<<(), MillePedeFileConverter_cfg::out, reco::PFTrajectoryPoint::position(), HiIsolationCommonParameters_cff::track, trajectoryPoint(), trajectoryPoints(), and trajectoryPoints_.

Referenced by PFElecTkProducer::isSharingEcalEnergyWithEgSC(), ConvBremPFTrackFinder::runConvBremFinder(), HGCalAndBREMLinker::testLink(), GSFAndHCALLinker::testLink(), HCALAndBREMLinker::testLink(), GSFAndECALLinker::testLink(), ECALAndBREMLinker::testLink(), GSFAndHGCalLinker::testLink(), LinkByRecHit::testTrackAndClusterByRecHit(), and trajectoryPoint().

                                                                               {
   const unsigned offset_layerid = nTrajectoryMeasurements() + layerid;
   if( layerid >= reco::PFTrajectoryPoint::NLayers ||
       offset_layerid >= trajectoryPoints_.size() ) {
 
     // cout<<(*this)<<endl;
     // cout<<"lid "<<layerid<<" "<<nTrajectoryMeasurements()<<" "<<trajectoryPoints_.size()<<endl;
     
     throw cms::Exception("SizeError")<<"PFRecTrack::extrapolatedPoint: cannot access "
                      <<layerid
                      <<" #traj meas = "<<nTrajectoryMeasurements()
                      <<" #traj points = "<<trajectoryPoints_.size()
                      <<endl
                      <<(*this);
     // assert(0);
   }
   if (layerid < indexInnermost_)
     return trajectoryPoints_[ layerid ];
   else
     return trajectoryPoints_[ offset_layerid ];  
 }