CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
TempTrajectory Class Reference

#include <TempTrajectory.h>

Public Types

typedef
TransientTrackingRecHit::ConstRecHitContainer 
ConstRecHitContainer
 
typedef cmsutils::bqueue
< TrajectoryMeasurement
DataContainer
 
typedef ConstRecHitContainer RecHitContainer
 

Public Member Functions

double chiSquared () const
 Value of the raw Chi2 of the trajectory, not normalised to the N.D.F. More...
 
PropagationDirection direction () const
 
bool empty () const
 True if trajectory has no measurements. More...
 
const TrajectoryMeasurementfirstMeasurement () const
 
int foundHits () const
 obsolete name, use measurements() instead. More...
 
void invalidate ()
 Method to invalidate a trajectory. Useful during ambiguity resolution. More...
 
bool isValid () const
 
void join (TempTrajectory &segment)
 
const DetLayerlastLayer () const
 Redundant method, returns the layer of lastMeasurement() . More...
 
const TrajectoryMeasurementlastMeasurement () const
 
int lostHits () const
 
const DataContainermeasurements () const
 
void pop ()
 
void popInvalidTail ()
 Pops out all the invalid hits on the tail. More...
 
void push (const TrajectoryMeasurement &tm)
 
void push (const TempTrajectory &segment)
 
void push (const TrajectoryMeasurement &tm, double chi2Increment)
 
const TrajectorySeedseed () const
 Access to the seed used to reconstruct the Trajectory. More...
 
 TempTrajectory ()
 
 TempTrajectory (const TrajectorySeed &seed)
 
 TempTrajectory (const TrajectorySeed &seed, PropagationDirection dir)
 
 TempTrajectory (const boost::shared_ptr< const TrajectorySeed > &seed, PropagationDirection dir)
 
 TempTrajectory (PropagationDirection dir)
 
 TempTrajectory (const 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 ()
 
static bool lost (const TransientTrackingRecHit &hit)
 

Private Member Functions

void check () const
 
void pushAux (const TrajectoryMeasurement &tm, double chi2Increment)
 

Private Attributes

double theChiSquared
 
DataContainer theData
 
PropagationDirection theDirection
 
bool theDirectionValidity
 
int theNumberOfFoundHits
 
int theNumberOfLostHits
 
boost::shared_ptr< const
TrajectorySeed
theSeed
 
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 40 of file TempTrajectory.h.

Member Typedef Documentation

Definition at line 45 of file TempTrajectory.h.

Definition at line 44 of file TempTrajectory.h.

Definition at line 46 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 54 of file TempTrajectory.h.

54  : theChiSquared(0), theValid(true),
57  theSeed()
58  {}
PropagationDirection theDirection
boost::shared_ptr< const TrajectorySeed > theSeed
TempTrajectory::TempTrajectory ( const TrajectorySeed seed)
inline

Constructor of an empty trajectory with undefined direction. The direction will be defined at the moment of the push of a second measurement, from the relative radii of the first and second measurements.

Definition at line 67 of file TempTrajectory.h.

67  :
68  theChiSquared(0), theValid(true),
71  theSeed( new TrajectorySeed(seed) )
72  {}
PropagationDirection theDirection
boost::shared_ptr< const TrajectorySeed > theSeed
TempTrajectory::TempTrajectory ( const TrajectorySeed seed,
PropagationDirection  dir 
)
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 78 of file TempTrajectory.h.

78  :
79  theChiSquared(0), theValid(true),
82  theSeed( new TrajectorySeed(seed) )
83  {}
PropagationDirection theDirection
boost::shared_ptr< const TrajectorySeed > theSeed
dbl *** dir
Definition: mlp_gen.cc:35
TempTrajectory::TempTrajectory ( const boost::shared_ptr< const TrajectorySeed > &  seed,
PropagationDirection  dir 
)
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 89 of file TempTrajectory.h.

89  :
90  theChiSquared(0), theValid(true),
93  theSeed( seed )
94  {}
const TrajectorySeed & seed() const
Access to the seed used to reconstruct the Trajectory.
PropagationDirection theDirection
boost::shared_ptr< const TrajectorySeed > theSeed
dbl *** dir
Definition: mlp_gen.cc:35
TempTrajectory::TempTrajectory ( PropagationDirection  dir)
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 101 of file TempTrajectory.h.

101  :
102  theChiSquared(0), theValid(true),
105  theSeed()
106  {}
PropagationDirection theDirection
boost::shared_ptr< const TrajectorySeed > theSeed
dbl *** dir
Definition: mlp_gen.cc:35
TempTrajectory::TempTrajectory ( const Trajectory traj)

construct TempTrajectory from standard Trajectory

Definition at line 7 of file TempTrajectory.cc.

References begin, end, Trajectory::measurements(), and push().

7  :
8  theChiSquared(0), theValid(traj.isValid()),
11  theSeed( traj.sharedSeed() ){
12 
13  Trajectory::DataContainer::const_iterator begin=traj.measurements().begin();
14  Trajectory::DataContainer::const_iterator end=traj.measurements().end();
15 
16  for(Trajectory::DataContainer::const_iterator it=begin; it!=end; ++it){
17  push(*it);
18  }
19 
20 
21 }
PropagationDirection theDirection
const boost::shared_ptr< const TrajectorySeed > & sharedSeed() const
Definition: Trajectory.h:274
PropagationDirection const & direction() const
Definition: Trajectory.cc:195
DataContainer const & measurements() const
Definition: Trajectory.h:169
#define end
Definition: vmac.h:38
boost::shared_ptr< const TrajectorySeed > theSeed
bool isValid() const
Definition: Trajectory.h:225
#define begin
Definition: vmac.h:31
void push(const TrajectoryMeasurement &tm)
TempTrajectory::~TempTrajectory ( )

destruct a TempTrajectory

Definition at line 23 of file TempTrajectory.cc.

23  {
24 }

Member Function Documentation

void TempTrajectory::check ( void  ) const
private

Definition at line 133 of file TempTrajectory.cc.

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

Referenced by firstMeasurement(), lastLayer(), and lastMeasurement().

133  {
134  if ( theData.size() == 0)
135  throw cms::Exception("TrackingTools/PatternTools","Trajectory::check() - information requested from empty Trajectory");
136 }
DataContainer theData
size_type size() const
Definition: bqueue.h:143
double TempTrajectory::chiSquared ( ) const
inline

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

Definition at line 226 of file TempTrajectory.h.

References theChiSquared.

Referenced by HICTrajectoryBuilder::qualityFilter(), and GroupedTrajCandLess::score().

226 { return theChiSquared;}
PropagationDirection TempTrajectory::direction ( ) const

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

Definition at line 128 of file TempTrajectory.cc.

References edm::hlt::Exception, theDirection, and theDirectionValidity.

Referenced by GroupedCkfTrajectoryBuilder::backwardFit(), MuonCkfTrajectoryBuilder::findCompatibleMeasurements(), HICTrajectoryBuilder::findCompatibleMeasurements(), BaseCkfTrajectoryBuilder::findStateAndLayers(), join(), and push().

128  {
130  else throw cms::Exception("TrackingTools/PatternTools","Trajectory::direction() requested but not set");
131 }
PropagationDirection theDirection
bool TempTrajectory::empty ( void  ) const
inline
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 193 of file TempTrajectory.h.

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

193  {
194  check(); return theData.front();
195  }
DataContainer theData
void check() const
int TempTrajectory::foundHits ( ) const
inline

obsolete name, use measurements() instead.

Return all RecHits in a container.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 213 of file TempTrajectory.h.

References theNumberOfFoundHits.

Referenced by BaseCkfTrajectoryBuilder::addToResult(), GroupedCkfTrajectoryBuilder::advanceOneLayer(), GroupedCkfTrajectoryBuilder::backwardFit(), TrajectoryLessByFoundHits::operator()(), ThresholdPtTrajectoryFilter::qualityFilter(), MinPtTrajectoryFilter::qualityFilter(), HICTrajectoryBuilder::qualityFilter(), GroupedTrajCandLess::score(), MinPtTrajectoryFilter::toBeContinued(), ThresholdPtTrajectoryFilter::toBeContinued(), and BaseCkfTrajectoryBuilder::toBeContinued().

213 { return theNumberOfFoundHits;}
static bool TempTrajectory::inactive ( )
inlinestatic

Definition of inactive Det from the Trajectory point of view.

Definition at line 248 of file TempTrajectory.h.

249  {return false;}//FIXME
void TempTrajectory::invalidate ( )
inline

Method to invalidate a trajectory. Useful during ambiguity resolution.

Definition at line 240 of file TempTrajectory.h.

References theValid.

240 { theValid = false;}
bool TempTrajectory::isValid ( void  ) const
inline
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 90 of file TempTrajectory.cc.

References alongMomentum, cmsutils::bqueue< T >::back(), cmsutils::bqueue< T >::clear(), direction(), cmsutils::bqueue< T >::front(), TrajectoryStateOnSurface::globalPosition(), cmsutils::bqueue< T >::join(), lost(), measurements(), oppositeToMomentum, PV3DBase< T, PVType, FrameType >::perp(), push(), cmsutils::bqueue< T >::rbegin(), cmsutils::bqueue< T >::rend(), cmsutils::bqueue< T >::shared(), cmsutils::bqueue< T >::size(), theChiSquared, theData, theDirection, theDirectionValidity, theNumberOfFoundHits, theNumberOfLostHits, and TrajectoryMeasurement::updatedState().

90  {
91  assert (segment.direction() == theDirection) ;
92  if (segment.theData.shared()) {
93  push(segment);
94  segment.theData.clear(); // obey the contract, and increase the chances it will be not shared one day
95  } else {
96  for (DataContainer::const_iterator it = segment.measurements().rbegin(), ed = segment.measurements().rend(); it != ed; --it) {
97  if ( it->recHit()->isValid()) theNumberOfFoundHits++;
98  else if (lost( *(it->recHit()) ) ) theNumberOfLostHits++;
99  theChiSquared += it->estimate();
100  }
101  theData.join(segment.theData);
102 
103  if ( !theDirectionValidity && theData.size() >= 2) {
108  theDirectionValidity = true;
109  }
110  }
111 }
const_iterator rend() const
Definition: bqueue.h:142
T perp() const
Definition: PV3DBase.h:66
const DataContainer & measurements() const
GlobalPoint globalPosition() const
PropagationDirection theDirection
static bool lost(const TransientTrackingRecHit &hit)
PropagationDirection direction() const
DataContainer theData
bool shared()
Definition: bqueue.h:151
void clear()
Definition: bqueue.h:172
TrajectoryStateOnSurface updatedState() const
const _bqueue_itr< value_type > const_iterator
Definition: bqueue.h:88
iterator rbegin()
Definition: bqueue.h:140
void join(bqueue< T > &other)
Definition: bqueue.h:158
size_type size() const
Definition: bqueue.h:143
void push(const TrajectoryMeasurement &tm)
const DetLayer* TempTrajectory::lastLayer ( ) const
inline
const TrajectoryMeasurement& TempTrajectory::lastMeasurement ( ) const
inline
bool TempTrajectory::lost ( const TransientTrackingRecHit hit)
static

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

Definition at line 138 of file TempTrajectory.cc.

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

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

139 {
140  if ( hit.isValid()) return false;
141  else {
142  // // A DetLayer is always inactive in this logic.
143  // // The DetLayer is the Det of an invalid RecHit only if no DetUnit
144  // // is compatible with the predicted state, so we don't really expect
145  // // a hit in this case.
146 
147  if(hit.geographicalId().rawId() == 0) {return false;}
148  else{
149  return hit.getType() == TrackingRecHit::missing;
150  }
151  }
152 }
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
Type getType() const
bool isValid() const
DetId geographicalId() const
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 220 of file TempTrajectory.h.

References theNumberOfLostHits.

Referenced by BaseCkfTrajectoryBuilder::addToResult(), GroupedTrajCandLess::score(), BaseCkfTrajectoryBuilder::toBeContinued(), and HICTrajectoryBuilder::toBeContinued().

220 { return theNumberOfLostHits;}
const DataContainer& TempTrajectory::measurements ( ) const
inline
void TempTrajectory::pop ( )

Remove the last measurement from the trajectory.

Definition at line 26 of file TempTrajectory.cc.

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

Referenced by BaseCkfTrajectoryBuilder::addToResult().

26  {
27  if (!empty()) {
28  if (theData.back().recHit()->isValid()) theNumberOfFoundHits--;
29  else if(lost(* (theData.back().recHit()) )) theNumberOfLostHits--;
30  theData.pop_back();
31  }
32 }
bool empty() const
True if trajectory has no measurements.
ConstRecHitPointer recHit() const
static bool lost(const TransientTrackingRecHit &hit)
DataContainer theData
void pop_back()
Definition: bqueue.h:130
void TempTrajectory::popInvalidTail ( )

Pops out all the invalid hits on the tail.

void TempTrajectory::push ( const TrajectoryMeasurement tm)
void TempTrajectory::push ( const TempTrajectory 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)

Definition at line 79 of file TempTrajectory.cc.

References direction(), list(), measurements(), push(), cmsutils::bqueue< T >::rbegin(), cmsutils::bqueue< T >::rend(), and theDirection.

79  {
80  assert (segment.direction() == theDirection) ;
81  __gnu_cxx::slist<const TrajectoryMeasurement*> list;
82  for (DataContainer::const_iterator it = segment.measurements().rbegin(), ed = segment.measurements().rend(); it != ed; --it) {
83  list.push_front(&(*it));
84  }
85  for(__gnu_cxx::slist<const TrajectoryMeasurement*>::const_iterator it = list.begin(), ed = list.end(); it != ed; ++it) {
86  push(**it);
87  }
88 }
const_iterator rend() const
Definition: bqueue.h:142
const DataContainer & measurements() const
PropagationDirection theDirection
PropagationDirection direction() const
const _bqueue_itr< value_type > const_iterator
Definition: bqueue.h:88
iterator rbegin()
Definition: bqueue.h:140
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger list("!*","!HLTx*"if it matches 2 triggers or more) will accept the event if all the matching triggers are FAIL.It will reject the event if any of the triggers are PASS or EXCEPTION(this matches the behavior of"!*"before the partial wildcard feature was incorporated).Triggers which are in the READY state are completely ignored.(READY should never be returned since the trigger paths have been run
void push(const TrajectoryMeasurement &tm)
void TempTrajectory::push ( const TrajectoryMeasurement tm,
double  chi2Increment 
)

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

Definition at line 44 of file TempTrajectory.cc.

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

44  {
45  pushAux(tm,chi2Increment);
46  theData.push_back(tm);
47 }
DataContainer theData
void pushAux(const TrajectoryMeasurement &tm, double chi2Increment)
void push_back(const T &val)
Definition: bqueue.h:119
void TempTrajectory::pushAux ( const TrajectoryMeasurement tm,
double  chi2Increment 
)
private

Definition at line 56 of file TempTrajectory.cc.

References alongMomentum, cmsutils::bqueue< T >::back(), cmsutils::bqueue< T >::front(), TrajectoryStateOnSurface::globalPosition(), lost(), oppositeToMomentum, PV3DBase< T, PVType, FrameType >::perp(), TrajectoryMeasurement::recHit(), cmsutils::bqueue< T >::size(), theChiSquared, theData, theDirection, theDirectionValidity, theNumberOfFoundHits, theNumberOfLostHits, and TrajectoryMeasurement::updatedState().

Referenced by push().

57 {
58  if ( tm.recHit()->isValid()) {
60  }
61  //else if (lost( tm.recHit()) && !inactive(tm.recHit().det())) theNumberOfLostHits++;
62  else if (lost( *(tm.recHit()) ) ) theNumberOfLostHits++;
63 
64 
65  theChiSquared += chi2Increment;
66 
67  // in case of a Trajectory constructed without direction,
68  // determine direction from the radii of the first two measurements
69 
70  if ( !theDirectionValidity && theData.size() >= 2) {
75  theDirectionValidity = true;
76  }
77 }
T perp() const
Definition: PV3DBase.h:66
GlobalPoint globalPosition() const
ConstRecHitPointer recHit() const
PropagationDirection theDirection
static bool lost(const TransientTrackingRecHit &hit)
DataContainer theData
TrajectoryStateOnSurface updatedState() const
size_type size() const
Definition: bqueue.h:143
const TrajectorySeed& TempTrajectory::seed ( ) const
inline

Access to the seed used to reconstruct the Trajectory.

Definition at line 243 of file TempTrajectory.h.

References theSeed.

Referenced by GroupedCkfTrajectoryBuilder::backwardFit(), MuonCkfTrajectoryBuilder::findCompatibleMeasurements(), BaseCkfTrajectoryBuilder::findStateAndLayers(), and GroupedCkfTrajectoryBuilder::rebuildSeedingRegion().

243 { return *theSeed;}
boost::shared_ptr< const TrajectorySeed > theSeed
Trajectory TempTrajectory::toTrajectory ( ) const

Convert to a standard Trajectory.

Definition at line 154 of file TempTrajectory.cc.

References Trajectory::push(), cmsutils::bqueue< T >::rbegin(), cmsutils::bqueue< T >::rend(), Trajectory::reserve(), cmsutils::bqueue< T >::size(), theData, theDirection, and theSeed.

Referenced by BaseCkfTrajectoryBuilder::addToResult(), HICTrajectoryBuilder::addToResult(), CkfDebugTrajectoryBuilder::analyzeMeasurementsDebugger(), HICTrajectoryBuilder::findCompatibleMeasurements(), HICTrajectoryBuilder::qualityFilter(), and HICTrajectoryBuilder::updateTrajectory().

154  {
156 
157  traj.reserve(theData.size());
158  static std::vector<const TrajectoryMeasurement*> work;
159  work.resize(theData.size(), 0);
160  std::vector<const TrajectoryMeasurement*>::iterator workend = work.end(), itwork = workend;
161  for (TempTrajectory::DataContainer::const_iterator it = theData.rbegin(), ed = theData.rend(); it != ed; --it) {
162  --itwork; *itwork = (&(*it));
163  }
164  for (; itwork != workend; ++itwork) {
165  traj.push(**itwork);
166  }
167  return traj;
168 }
const_iterator rend() const
Definition: bqueue.h:142
PropagationDirection theDirection
DataContainer theData
boost::shared_ptr< const TrajectorySeed > theSeed
iterator rbegin()
Definition: bqueue.h:140
size_type size() const
Definition: bqueue.h:143

Member Data Documentation

double TempTrajectory::theChiSquared
private

Definition at line 274 of file TempTrajectory.h.

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

DataContainer TempTrajectory::theData
private
PropagationDirection TempTrajectory::theDirection
private

Definition at line 280 of file TempTrajectory.h.

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

bool TempTrajectory::theDirectionValidity
private

Definition at line 281 of file TempTrajectory.h.

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

int TempTrajectory::theNumberOfFoundHits
private

Definition at line 277 of file TempTrajectory.h.

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

int TempTrajectory::theNumberOfLostHits
private

Definition at line 278 of file TempTrajectory.h.

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

boost::shared_ptr<const TrajectorySeed> TempTrajectory::theSeed
private

Definition at line 283 of file TempTrajectory.h.

Referenced by seed(), and toTrajectory().

bool TempTrajectory::theValid
private

Definition at line 275 of file TempTrajectory.h.

Referenced by invalidate(), and isValid().