#include <TempTrajectory.h>

Public Types
typedef TrackingRecHit::ConstRecHitContainer	ConstRecHitContainer

typedef cmsutils::bqueue < TrajectoryMeasurement >	DataContainer

typedef ConstRecHitContainer	RecHitContainer

Public Member Functions
int	cccBadHits () const

float	chiSquared () const
	Value of the raw Chi2 of the trajectory, not normalised to the N.D.F. More...

PropagationDirection	direction () const

float	dPhiCacheForLoopersReconstruction () const

template<typename... Args>
void	emplace (Args &&...args)

template<typename... Args>
void	emplace (double chi2Increment, Args &&...args)

bool	empty () const
	True if trajectory has no measurements. More...

const TrajectoryMeasurement &	firstMeasurement () const

int	foundHits () const

void	incrementLoops ()

void	invalidate ()
	Method to invalidate a trajectory. Useful during ambiguity resolution. More...

bool	isLooper () const

bool	isValid () const

void	join (TempTrajectory &segment)

const DetLayer *	lastLayer () const
	Redundant method, returns the layer of lastMeasurement() . More...

const TrajectoryMeasurement &	lastMeasurement () const

int	lostHits () const

const DataContainer &	measurements () const

signed char	nLoops () const

int	numberOfCCCBadHits (float ccc_threshold)

TempTrajectory &	operator= (TempTrajectory const &rh)=default

TempTrajectory &	operator= (TempTrajectory &&rh)

void	pop ()

void	popInvalidTail ()
	Pops out all the invalid hits on the tail. More...

void	push (const TrajectoryMeasurement &tm)

void	push (TrajectoryMeasurement &&tm)

void	push (TempTrajectory const &segment)

void	push (const TrajectoryMeasurement &tm, double chi2Increment)

void	push (TrajectoryMeasurement &&tm, double chi2Increment)

unsigned int	seedNHits () const

void	setDPhiCacheForLoopersReconstruction (float dphi)

void	setNLoops (signed char value)

void	setStopReason (StopReason s)

StopReason	stopReason () const

	TempTrajectory ()

	TempTrajectory (PropagationDirection dir, unsigned char nhseed)

	TempTrajectory (TempTrajectory const &rh)=default

	TempTrajectory (TempTrajectory &&rh)

	TempTrajectory (Trajectory &&traj)
	construct TempTrajectory from standard Trajectory More...

Trajectory	toTrajectory () const
	Convert to a standard Trajectory. More...

	~TempTrajectory ()
	destruct a TempTrajectory More...

Static Public Member Functions
static bool	inactive ()

Private Member Functions
bool	badForCCC (const TrajectoryMeasurement &tm)

void	check () const

void	pushAux (double chi2Increment)

void	updateBadForCCC (float ccc_threshold)

Static Private Member Functions
static bool	lost (const TrackingRecHit &hit)

Private Attributes
StopReason	stopReason_

float	theCCCThreshold_

float	theChiSquared

DataContainer	theData

signed char	theDirection

float	theDPhiCache

unsigned char	theNHseed

signed char	theNLoops

signed short	theNumberOfCCCBadHits_

signed short	theNumberOfFoundHits

signed short	theNumberOfLostHits

bool	theValid

Detailed Description

A class for detailed particle trajectory representation. It is used during trajectory building to "grow" a trajectory. The trajectory is represented as an ordered sequence of TrajectoryMeasurement objects with a stack-like interface. The measurements are added to the Trajectory in the order of increasing precision: each new TrajectoryMeasurement is assumed to improve the precision of the last one, normally by adding a constraint from a new RecHit. However the Trajectory class does not have the means to verify that measurements are added in the correct order, and thus cannot guarantee the order, which is the responsibility of the TrajectoryBuilder. The Trajectory provides some security by allowing to add or remove measurements only on one of it's ends, with push(TM) and pop() methods. The last measurement in a Trajectory can thus be either the innermost (closest to the interaction point) or the outermost, depending on the way the Trajectory was built. The direction of building is represented as a PropagationDirection, which has two possible values: alongMomentum (outwards) and oppositeToMomentum (inwards), and is accessed with the direction() method.

Definition at line 41 of file TempTrajectory.h.

Member Typedef Documentation

typedef TrackingRecHit::ConstRecHitContainer TempTrajectory::ConstRecHitContainer

Definition at line 46 of file TempTrajectory.h.

typedef cmsutils::bqueue<TrajectoryMeasurement> TempTrajectory::DataContainer

Definition at line 45 of file TempTrajectory.h.

typedef ConstRecHitContainer TempTrajectory::RecHitContainer

Definition at line 47 of file TempTrajectory.h.

Constructor & Destructor Documentation

TempTrajectory::TempTrajectory ( )

inline

Default constructor of an empty trajectory with undefined seed and undefined direction. This constructor is necessary in order to transiently copy vector<Trajectory> in the edm::Event

Definition at line 55 of file TempTrajectory.h.

TempTrajectory::TempTrajectory	(	PropagationDirection	dir,
		unsigned char	nhseed
	)

inline

Constructor of an empty trajectory with defined direction. No check is made in the push method that measurements are added in the correct direction.

Definition at line 75 of file TempTrajectory.h.

TempTrajectory::TempTrajectory ( TempTrajectory const & rh )

default

TempTrajectory::TempTrajectory ( TempTrajectory && rh )

inline

Definition at line 95 of file TempTrajectory.h.

TempTrajectory::TempTrajectory ( Trajectory && traj )

explicit

construct TempTrajectory from standard Trajectory

Definition at line 11 of file TempTrajectory.cc.

References begin, end, eostools::move(), and push().

                                                 :
   theChiSquared(0),
   theNumberOfFoundHits(0),
   theNumberOfLostHits(0),
   theNumberOfCCCBadHits_(0),
   theDirection(traj.direction()),
   theValid(traj.isValid()),
   theNHseed(traj.seedNHits()),
   theNLoops(traj.nLoops()),
   theDPhiCache(traj.dPhiCacheForLoopersReconstruction()),
   theCCCThreshold_(traj.cccThreshold()),
   stopReason_(traj.stopReason()) {
 
   Trajectory::DataContainer::const_iterator begin=traj.measurements().begin();
   Trajectory::DataContainer::const_iterator end=traj.measurements().end();
 
   for(Trajectory::DataContainer::const_iterator it=begin; it!=end; ++it){
     push(std::move(*it));
   }
 
 }

TempTrajectory::~TempTrajectory ( )

inline

destruct a TempTrajectory

Definition at line 133 of file TempTrajectory.h.

Member Function Documentation

bool TempTrajectory::badForCCC ( const TrajectoryMeasurement & tm )

private

Definition at line 115 of file TempTrajectory.cc.

References siStripClusterTools::chargePerCM(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localParameters(), TrajectoryMeasurement::recHit(), theCCCThreshold_, and TrajectoryMeasurement::updatedState().

Referenced by pop(), pushAux(), and updateBadForCCC().

                                                               {
   auto const * thit = dynamic_cast<const BaseTrackerRecHit*>( tm.recHit()->hit() );
   if (!thit)
     return false;
   if (thit->isPixel())
     return false;
   if (!tm.updatedState().isValid())
     return false;
   return siStripClusterTools::chargePerCM(thit->rawId(),
                                           thit->firstClusterRef().stripCluster(),
                                           tm.updatedState().localParameters()) < theCCCThreshold_;
 }

int TempTrajectory::cccBadHits ( ) const

inline

Number of hits that are not compatible with the CCC used during patter recognition. Used mainly as a criteria for abandoning a trajectory candidate during trajectory building.

Definition at line 236 of file TempTrajectory.h.

void TempTrajectory::check ( void ) const

private

Definition at line 97 of file TempTrajectory.cc.

References cmsutils::bqueue< T >::size(), and theData.

                                  {
   if ( theData.size() == 0) 
     throw cms::Exception("TrackingTools/PatternTools","Trajectory::check() - information requested from empty Trajectory");
 }

float TempTrajectory::chiSquared ( ) const

inline

Value of the raw Chi2 of the trajectory, not normalised to the N.D.F.

Definition at line 245 of file TempTrajectory.h.

PropagationDirection TempTrajectory::direction ( ) const

Direction of "growing" of the trajectory. Possible values are alongMomentum (outwards) and oppositeToMomentum (inwards).

Definition at line 93 of file TempTrajectory.cc.

References theDirection.

                                                      {
   return PropagationDirection(theDirection);
 }

float TempTrajectory::dPhiCacheForLoopersReconstruction ( ) const

inline

accessor to the delta phi angle betweem the directions of the two measurements on the last two layers crossed by the trajectory

Definition at line 282 of file TempTrajectory.h.

template<typename... Args>

void TempTrajectory::emplace ( Args &&... args )

inline

Definition at line 148 of file TempTrajectory.h.

template<typename... Args>

void TempTrajectory::emplace	(	double	chi2Increment,
		Args &&...	args
	)

inline

Definition at line 185 of file TempTrajectory.h.

bool TempTrajectory::empty ( void ) const

inline

True if trajectory has no measurements.

Definition at line 242 of file TempTrajectory.h.

Referenced by pop().

const TrajectoryMeasurement& TempTrajectory::firstMeasurement ( ) const

inline

Access to the first measurement. It is the least precise one in a trajectory before smoothing. It is precise in a smoothed trajectory. It's the innermost measurement if direction() == alongMomentum, the outermost one if direction() == oppositeToMomentum.

Definition at line 210 of file TempTrajectory.h.

int TempTrajectory::foundHits ( ) const

inline

Number of valid RecHits used to determine the trajectory. Can be less than the number of measurements in data() since detector layers crossed without using RecHits from them are also stored as measurements.

Definition at line 224 of file TempTrajectory.h.

Referenced by lessByFoundHits(), TrajectoryLessByFoundHits::operator()(), ThresholdPtTrajectoryFilter::qualityFilter(), MinPtTrajectoryFilter::qualityFilter(), ThresholdPtTrajectoryFilter::toBeContinued(), and MinPtTrajectoryFilter::toBeContinued().

static bool TempTrajectory::inactive ( )

inlinestatic

Definition of inactive Det from the Trajectory point of view.

Definition at line 264 of file TempTrajectory.h.

void TempTrajectory::incrementLoops ( )

inline

Definition at line 292 of file TempTrajectory.h.

void TempTrajectory::invalidate ( )

inline

Method to invalidate a trajectory. Useful during ambiguity resolution.

Definition at line 259 of file TempTrajectory.h.

bool TempTrajectory::isLooper ( ) const

inline

Definition at line 288 of file TempTrajectory.h.

bool TempTrajectory::isValid ( void ) const

inline

Returns true if the Trajectory is valid. Trajectories are invalidated e.g. during ambiguity resolution.

Definition at line 256 of file TempTrajectory.h.

Referenced by IntermediateTrajectoryCleaner::clean(), ChargeSignificanceTrajectoryFilter::qualityFilter(), and core.AutoHandle.AutoHandle::ReallyLoad().

void TempTrajectory::join ( TempTrajectory & segment )

Add a new sets of measurements to a Trajectory Exactly like push(TempTrajectory), but it doesn't copy the data (the input segment will be reset to an empty one)

Definition at line 75 of file TempTrajectory.cc.

References assert(), cmsutils::bqueue< T >::clear(), cmsutils::bqueue< T >::join(), push(), cmsutils::bqueue< T >::shared(), theCCCThreshold_, theChiSquared, theData, theDirection, theNumberOfCCCBadHits_, theNumberOfFoundHits, theNumberOfLostHits, and updateBadForCCC().

                                                   {
   assert (segment.theDirection == theDirection) ;
 
   if (theCCCThreshold_ != segment.theCCCThreshold_)
     segment.updateBadForCCC(theCCCThreshold_);
   if (segment.theData.shared()) {
     push(segment);
     segment.theData.clear(); // obey the contract, and increase the chances it will be not shared one day
   } else {
     theData.join(segment.theData);
     theNumberOfFoundHits+= segment.theNumberOfFoundHits;
     theNumberOfLostHits += segment.theNumberOfLostHits;
     theNumberOfCCCBadHits_ += segment.theNumberOfCCCBadHits_;
     theChiSquared += segment.theChiSquared;
   }
 }

const DetLayer* TempTrajectory::lastLayer ( ) const

inline

Redundant method, returns the layer of lastMeasurement() .

Definition at line 269 of file TempTrajectory.h.

const TrajectoryMeasurement& TempTrajectory::lastMeasurement ( ) const

inline

Access to the last measurement. It's the most precise one in a trajectory before smoothing. It's the outermost measurement if direction() == alongMomentum, the innermost one if direction() == oppositeToMomentum.

Definition at line 200 of file TempTrajectory.h.

Referenced by ThresholdPtTrajectoryFilter::qualityFilter(), MinPtTrajectoryFilter::qualityFilter(), ThresholdPtTrajectoryFilter::toBeContinued(), and MinPtTrajectoryFilter::toBeContinued().

bool TempTrajectory::lost ( const TrackingRecHit & hit )

staticprivate

Definition of what it means for a hit to be "lost". This definition is also used by the TrajectoryBuilder.

Definition at line 102 of file TempTrajectory.cc.

References TrackingRecHit::geographicalId(), TrackingRecHit::getType(), TrackingRecHit::isValid(), likely, TrackingRecHit::missing, and DetId::rawId().

Referenced by pop(), and pushAux().

 {
   if  likely(hit.isValid()) return false;
 
   //     // A DetLayer is always inactive in this logic.
   //     // The DetLayer is the Det of an invalid RecHit only if no DetUnit 
   //     // is compatible with the predicted state, so we don't really expect
   //     // a hit in this case.
 
   if(hit.geographicalId().rawId() == 0) {return false;}
   return hit.getType() == TrackingRecHit::missing;
 }

int TempTrajectory::lostHits ( ) const

inline

Number of detector layers crossed without valid RecHits. Used mainly as a criteria for abandoning a trajectory candidate during trajectory building.

Definition at line 230 of file TempTrajectory.h.

const DataContainer& TempTrajectory::measurements ( ) const

inline

Return all measurements in a container.

Definition at line 216 of file TempTrajectory.h.

Referenced by push().

signed char TempTrajectory::nLoops ( ) const

inline

Definition at line 289 of file TempTrajectory.h.

int TempTrajectory::numberOfCCCBadHits ( float ccc_threshold )

Definition at line 143 of file TempTrajectory.cc.

References theNumberOfCCCBadHits_, and updateBadForCCC().

                                                           {
   updateBadForCCC(ccc_threshold);
   return theNumberOfCCCBadHits_;
 }

TempTrajectory& TempTrajectory::operator= ( TempTrajectory const & rh )

default

TempTrajectory& TempTrajectory::operator= ( TempTrajectory && rh )

inline

Definition at line 109 of file TempTrajectory.h.

void TempTrajectory::pop ( )

Remove the last measurement from the trajectory.

Definition at line 34 of file TempTrajectory.cc.

References cmsutils::bqueue< T >::back(), badForCCC(), empty(), lost(), cmsutils::bqueue< T >::pop_back(), TrajectoryMeasurement::recHit(), theData, theNumberOfCCCBadHits_, theNumberOfFoundHits, and theNumberOfLostHits.

                          { 
   if (!empty()) {
     if (theData.back().recHit()->isValid()) theNumberOfFoundHits--;
     else if(lost(* (theData.back().recHit()) )) theNumberOfLostHits--;
     else if(badForCCC(theData.back())) theNumberOfCCCBadHits_--;
     theData.pop_back();
   }
 }

void TempTrajectory::popInvalidTail ( )

Pops out all the invalid hits on the tail.

void TempTrajectory::push ( const TrajectoryMeasurement & tm )

inline

Add a new measurement to a Trajectory. The Chi2 of the trajectory is incremented by the value of tm.estimate() .

Definition at line 139 of file TempTrajectory.h.

Referenced by join(), and TempTrajectory().

void TempTrajectory::push ( TrajectoryMeasurement && tm )

inline

Definition at line 143 of file TempTrajectory.h.

void TempTrajectory::push ( TempTrajectory const & segment )

Add a new sets of measurements to a Trajectory The sorting of hits in the other trajectory must match the one inside this trajectory (that is, both along or both opposite to momentum) (the input segment will be reset to an empty one)

Definition at line 58 of file TempTrajectory.cc.

References assert(), i, measurements(), N, cmsutils::bqueue< T >::push_back(), cmsutils::bqueue< T >::size(), theCCCThreshold_, theChiSquared, theData, theDirection, theNumberOfCCCBadHits_, theNumberOfFoundHits, theNumberOfLostHits, and tmp.

                                                        {
   assert (segment.theDirection == theDirection) ;
   assert (segment.theCCCThreshold_ == theCCCThreshold_);
 
   const int N = segment.measurements().size();
   TrajectoryMeasurement const * tmp[N];
   int i=0;
   //for (DataContainer::const_iterator it = segment.measurements().rbegin(), ed = segment.measurements().rend(); it != ed; --it)
   for ( auto const & tm : segment.measurements())
     tmp[i++] =&tm;
   while(i!=0) theData.push_back(*tmp[--i]);
   theNumberOfFoundHits+= segment.theNumberOfFoundHits;
   theNumberOfLostHits += segment.theNumberOfLostHits;
   theNumberOfCCCBadHits_ += segment.theNumberOfCCCBadHits_;
   theChiSquared += segment.theChiSquared;
 }

void TempTrajectory::push	(	const TrajectoryMeasurement &	tm,
		double	chi2Increment
	)

inline

same as the one-argument push, but the trajectory Chi2 is incremented by chi2Increment. Useful e.g. in trajectory smoothing.

Definition at line 173 of file TempTrajectory.h.

void TempTrajectory::push	(	TrajectoryMeasurement &&	tm,
		double	chi2Increment
	)

inline

Definition at line 178 of file TempTrajectory.h.

void TempTrajectory::pushAux ( double chi2Increment )

private

Definition at line 45 of file TempTrajectory.cc.

References cmsutils::bqueue< T >::back(), badForCCC(), lost(), TrajectoryMeasurement::recHit(), theChiSquared, theData, theNumberOfCCCBadHits_, theNumberOfFoundHits, and theNumberOfLostHits.

                                                  {
   const TrajectoryMeasurement& tm = theData.back();
   if ( tm.recHit()->isValid()) {
     theNumberOfFoundHits++;
    }
   //else if (lost( tm.recHit()) && !inactive(tm.recHit().det())) theNumberOfLostHits++;
   else if (lost( *(tm.recHit()) ) )   theNumberOfLostHits++;
   else if (badForCCC(tm)) theNumberOfCCCBadHits_++;
 
   theChiSquared += chi2Increment;
 
 }

unsigned int TempTrajectory::seedNHits ( ) const

inline

Definition at line 239 of file TempTrajectory.h.

void TempTrajectory::setDPhiCacheForLoopersReconstruction ( float dphi )

inline

method to set the delta phi angle betweem the directions of the two measurements on the last two layers crossed by the trajectory

Definition at line 286 of file TempTrajectory.h.

void TempTrajectory::setNLoops ( signed char value )

inline

Definition at line 291 of file TempTrajectory.h.

void TempTrajectory::setStopReason ( StopReason s )

inline

Definition at line 295 of file TempTrajectory.h.

Referenced by MinPtTrajectoryFilter::toBeContinued().

StopReason TempTrajectory::stopReason ( ) const

inline

Definition at line 294 of file TempTrajectory.h.

Trajectory TempTrajectory::toTrajectory ( ) const

Convert to a standard Trajectory.

Definition at line 148 of file TempTrajectory.cc.

References i, AlCaHLTBitMon_ParallelJobs::p, Trajectory::push(), cmsutils::bqueue< T >::rbegin(), cmsutils::bqueue< T >::rend(), Trajectory::reserve(), Trajectory::setNLoops(), Trajectory::setStopReason(), cmsutils::bqueue< T >::size(), stopReason_, theData, theDirection, theNLoops, and tmp.

Referenced by CkfDebugTrajectoryBuilder::analyzeMeasurementsDebugger().

                                               {
   PropagationDirection p=PropagationDirection(theDirection);
   Trajectory traj(p);
   traj.setNLoops(theNLoops);
   traj.setStopReason(stopReason_);
 
   traj.reserve(theData.size());
   const TrajectoryMeasurement* tmp[theData.size()];
   int i=0;
   for (DataContainer::const_iterator it = theData.rbegin(), ed = theData.rend(); it != ed; --it)
     tmp[i++] = &(*it);
   while(i!=0) traj.push(*tmp[--i]);
   return traj;
 }

void TempTrajectory::updateBadForCCC ( float ccc_threshold )

private

Definition at line 128 of file TempTrajectory.cc.

References badForCCC(), h, theCCCThreshold_, theData, and theNumberOfCCCBadHits_.

Referenced by join(), and numberOfCCCBadHits().

                                                         {
   // If the supplied threshold is the same as the currently cached
   // one, then return the current number of bad hits for CCC,
   // otherwise do a new full rescan.
   if (ccc_threshold == theCCCThreshold_)
     return;
 
   theCCCThreshold_ = ccc_threshold;
   theNumberOfCCCBadHits_ = 0;
   for (auto const & h : theData) {
     if (badForCCC(h))
       theNumberOfCCCBadHits_++;
   }
 }

Member Data Documentation

StopReason TempTrajectory::stopReason_

private

Definition at line 331 of file TempTrajectory.h.

Referenced by toTrajectory().

float TempTrajectory::theCCCThreshold_

private

Definition at line 330 of file TempTrajectory.h.

Referenced by badForCCC(), join(), push(), and updateBadForCCC().

float TempTrajectory::theChiSquared

private

Definition at line 316 of file TempTrajectory.h.

Referenced by join(), push(), and pushAux().

DataContainer TempTrajectory::theData

private

Definition at line 314 of file TempTrajectory.h.

Referenced by check(), join(), pop(), push(), pushAux(), toTrajectory(), and updateBadForCCC().

signed char TempTrajectory::theDirection

private

Definition at line 323 of file TempTrajectory.h.

Referenced by direction(), join(), push(), and toTrajectory().

float TempTrajectory::theDPhiCache

private

Definition at line 329 of file TempTrajectory.h.

unsigned char TempTrajectory::theNHseed

private

Definition at line 326 of file TempTrajectory.h.

signed char TempTrajectory::theNLoops

private

Definition at line 328 of file TempTrajectory.h.

Referenced by toTrajectory().

signed short TempTrajectory::theNumberOfCCCBadHits_

private

Definition at line 320 of file TempTrajectory.h.

Referenced by join(), numberOfCCCBadHits(), pop(), push(), pushAux(), and updateBadForCCC().

signed short TempTrajectory::theNumberOfFoundHits

private

Definition at line 318 of file TempTrajectory.h.

Referenced by join(), pop(), push(), and pushAux().

signed short TempTrajectory::theNumberOfLostHits

private

Definition at line 319 of file TempTrajectory.h.

Referenced by join(), pop(), push(), and pushAux().

bool TempTrajectory::theValid

private

Definition at line 324 of file TempTrajectory.h.

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation