|
|
|
#include <BaseCkfTrajectoryBuilder.h>
The component of track reconstruction that, strating from a seed, reconstructs all possible trajectories. The resulting trajectories may be mutually exclusive and require cleaning by a TrajectoryCleaner. The Trajectories are normally not smoothed.
Definition at line 46 of file BaseCkfTrajectoryBuilder.h.
typedef FreeTrajectoryState BaseCkfTrajectoryBuilder::FTS [protected] |
Definition at line 49 of file BaseCkfTrajectoryBuilder.h.
typedef std::pair<TSOS,std::vector<const DetLayer*> > BaseCkfTrajectoryBuilder::StateAndLayers [protected] |
Definition at line 52 of file BaseCkfTrajectoryBuilder.h.
| typedef std::vector<TempTrajectory> BaseCkfTrajectoryBuilder::TempTrajectoryContainer |
Reimplemented in CkfTrajectoryBuilder.
Definition at line 57 of file BaseCkfTrajectoryBuilder.h.
typedef TrajectoryMeasurement BaseCkfTrajectoryBuilder::TM [protected] |
Definition at line 51 of file BaseCkfTrajectoryBuilder.h.
| typedef std::vector<Trajectory> BaseCkfTrajectoryBuilder::TrajectoryContainer |
Reimplemented from TrajectoryBuilder.
Reimplemented in CkfTrajectoryBuilder.
Definition at line 56 of file BaseCkfTrajectoryBuilder.h.
| typedef TrajectoryContainer::iterator BaseCkfTrajectoryBuilder::TrajectoryIterator |
Reimplemented from TrajectoryBuilder.
Definition at line 58 of file BaseCkfTrajectoryBuilder.h.
typedef TrajectoryStateOnSurface BaseCkfTrajectoryBuilder::TSOS [protected] |
Definition at line 50 of file BaseCkfTrajectoryBuilder.h.
| BaseCkfTrajectoryBuilder::BaseCkfTrajectoryBuilder | ( | const edm::ParameterSet & | conf, |
| const TrajectoryStateUpdator * | updator, | ||
| const Propagator * | propagatorAlong, | ||
| const Propagator * | propagatorOpposite, | ||
| const Chi2MeasurementEstimatorBase * | estimator, | ||
| const TransientTrackingRecHitBuilder * | RecHitBuilder, | ||
| const MeasurementTracker * | measurementTracker, | ||
| const TrajectoryFilter * | filter, | ||
| const TrajectoryFilter * | inOutFilter = 0 |
||
| ) |
Definition at line 23 of file BaseCkfTrajectoryBuilder.cc.
References clustersToSkip_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and skipClusters_.
: theUpdator(updator), thePropagatorAlong(propagatorAlong),thePropagatorOpposite(propagatorOpposite), theEstimator(estimator),theTTRHBuilder(recHitBuilder), theMeasurementTracker(measurementTracker), theLayerMeasurements(new LayerMeasurements(theMeasurementTracker)), theForwardPropagator(0),theBackwardPropagator(0), theFilter(filter), theInOutFilter(inOutFilter) { if (conf.exists("clustersToSkip")){ skipClusters_=true; clustersToSkip_=conf.getParameter<edm::InputTag>("clustersToSkip"); } else skipClusters_=false; }
| BaseCkfTrajectoryBuilder::~BaseCkfTrajectoryBuilder | ( | ) | [virtual] |
Definition at line 49 of file BaseCkfTrajectoryBuilder.cc.
References theLayerMeasurements.
{
delete theLayerMeasurements;
}
| void BaseCkfTrajectoryBuilder::addToResult | ( | boost::shared_ptr< const TrajectorySeed > const & | seed, |
| TempTrajectory & | traj, | ||
| TrajectoryContainer & | result, | ||
| bool | inOut = false |
||
| ) | const [protected] |
Definition at line 160 of file BaseCkfTrajectoryBuilder.cc.
References Trajectory::empty(), Trajectory::foundHits(), Trajectory::lastMeasurement(), LogDebug, Trajectory::lostHits(), TempTrajectory::measurements(), Trajectory::measurements(), Trajectory::pop(), qualityFilter(), TrajectoryMeasurement::recHit(), Trajectory::setSharedSeed(), and TempTrajectory::toTrajectory().
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), and CkfTrajectoryBuilder::limitedCandidates().
{
// quality check
if ( !qualityFilter(tmptraj, inOut) ) return;
Trajectory traj = tmptraj.toTrajectory();
traj.setSharedSeed(seed);
// discard latest dummy measurements
while (!traj.empty() && !traj.lastMeasurement().recHit()->isValid()) traj.pop();
LogDebug("CkfPattern")<<inOut<<"=inOut option. pushing a Trajectory with: "<<traj.foundHits()<<" found hits. "<<traj.lostHits()
<<" lost hits. Popped :"<<(tmptraj.measurements().size())-(traj.measurements().size())<<" hits.";
result.push_back(std::move(traj));
}
| void BaseCkfTrajectoryBuilder::addToResult | ( | TempTrajectory const & | traj, |
| TempTrajectoryContainer & | result, | ||
| bool | inOut = false |
||
| ) | const [protected] |
Definition at line 177 of file BaseCkfTrajectoryBuilder.cc.
References TempTrajectory::empty(), TempTrajectory::foundHits(), TempTrajectory::lastMeasurement(), LogDebug, TempTrajectory::lostHits(), TempTrajectory::measurements(), TempTrajectory::pop(), qualityFilter(), and TrajectoryMeasurement::recHit().
{
// quality check
if ( !qualityFilter(tmptraj, inOut) ) return;
// discard latest dummy measurements
TempTrajectory traj = tmptraj;
while (!traj.empty() && !traj.lastMeasurement().recHit()->isValid()) traj.pop();
LogDebug("CkfPattern")<<inOut<<"=inOut option. pushing a TempTrajectory with: "<<traj.foundHits()<<" found hits. "<<traj.lostHits()
<<" lost hits. Popped :"<<(tmptraj.measurements().size())-(traj.measurements().size())<<" hits.";
result.push_back(std::move(traj));
}
| virtual bool BaseCkfTrajectoryBuilder::analyzeMeasurementsDebugger | ( | Trajectory & | traj, |
| std::vector< TrajectoryMeasurement > | meas, | ||
| const MeasurementTracker * | theMeasurementTracker, | ||
| const Propagator * | theForwardPropagator, | ||
| const Chi2MeasurementEstimatorBase * | theEstimator, | ||
| const TransientTrackingRecHitBuilder * | theTTRHBuilder | ||
| ) | const [inline, protected, virtual] |
Maximum number of lost hits per trajectory candidate. Maximum number of consecutive lost hits per trajectory candidate.
Reimplemented in CkfDebugTrajectoryBuilder.
Definition at line 95 of file BaseCkfTrajectoryBuilder.h.
Referenced by CkfTrajectoryBuilder::limitedCandidates().
{return true;}
| virtual bool BaseCkfTrajectoryBuilder::analyzeMeasurementsDebugger | ( | TempTrajectory & | traj, |
| std::vector< TrajectoryMeasurement > | meas, | ||
| const MeasurementTracker * | theMeasurementTracker, | ||
| const Propagator * | theForwardPropagator, | ||
| const Chi2MeasurementEstimatorBase * | theEstimator, | ||
| const TransientTrackingRecHitBuilder * | theTTRHBuilder | ||
| ) | const [inline, protected, virtual] |
Reimplemented in CkfDebugTrajectoryBuilder.
Definition at line 100 of file BaseCkfTrajectoryBuilder.h.
{return true;}
| virtual TempTrajectory BaseCkfTrajectoryBuilder::buildTrajectories | ( | const TrajectorySeed & | seed, |
| TrajectoryContainer & | ret, | ||
| const TrajectoryFilter * | |||
| ) | const [inline, virtual] |
Reimplemented in CkfTrajectoryBuilder, and GroupedCkfTrajectoryBuilder.
Definition at line 73 of file BaseCkfTrajectoryBuilder.h.
Referenced by cms::CkfTrackCandidateMakerBase::produceBase().
{ assert(0==1); return TempTrajectory();}
| TempTrajectory BaseCkfTrajectoryBuilder::createStartingTrajectory | ( | const TrajectorySeed & | seed | ) | const [protected] |
Definition at line 103 of file BaseCkfTrajectoryBuilder.cc.
References alongMomentum, TrajectorySeed::direction(), PrintoutHelper::dumpCandidate(), LogDebug, query::result, seedMeasurements(), theBackwardPropagator, and theForwardPropagator.
Referenced by CkfTrajectoryBuilder::buildTrajectories(), GroupedCkfTrajectoryBuilder::buildTrajectories(), and GroupedCkfTrajectoryBuilder::rebuildSeedingRegion().
{
TempTrajectory result(seed.direction());
if ( seed.direction() == alongMomentum) {
theForwardPropagator = &(*thePropagatorAlong);
theBackwardPropagator = &(*thePropagatorOpposite);
}
else {
theForwardPropagator = &(*thePropagatorOpposite);
theBackwardPropagator = &(*thePropagatorAlong);
}
seedMeasurements(seed, result);
LogDebug("CkfPattern")
<<" initial trajectory from the seed: "<<PrintoutHelper::dumpCandidate(result,true);
return result;
}
| virtual void BaseCkfTrajectoryBuilder::fillSeedHistoDebugger | ( | std::vector< TrajectoryMeasurement >::iterator | begin, |
| std::vector< TrajectoryMeasurement >::iterator | end | ||
| ) | const [inline, protected, virtual] |
Reimplemented in CkfDebugTrajectoryBuilder.
Definition at line 105 of file BaseCkfTrajectoryBuilder.h.
{;}
| BaseCkfTrajectoryBuilder::StateAndLayers BaseCkfTrajectoryBuilder::findStateAndLayers | ( | const TempTrajectory & | traj | ) | const [protected] |
Definition at line 231 of file BaseCkfTrajectoryBuilder.cc.
References TempTrajectory::direction(), TempTrajectory::empty(), TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), DetLayer::nextLayers(), and TrajectoryMeasurement::updatedState().
{
assert(!traj.empty());
TSOS const & currentState = traj.lastMeasurement().updatedState();
return StateAndLayers(currentState,traj.lastLayer()->nextLayers( *currentState.freeState(), traj.direction()) );
}
| BaseCkfTrajectoryBuilder::StateAndLayers BaseCkfTrajectoryBuilder::findStateAndLayers | ( | const TrajectorySeed & | seed, |
| const TempTrajectory & | traj | ||
| ) | const [protected] |
Definition at line 208 of file BaseCkfTrajectoryBuilder.cc.
References PTrajectoryStateOnDet::detId(), GeometricSearchTracker::detLayer(), TempTrajectory::direction(), TempTrajectory::empty(), TrajectoryStateOnSurface::freeState(), g, MeasurementTracker::geometricSearchTracker(), MeasurementTracker::geomTracker(), TrackingGeometry::idToDet(), TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), Propagator::magneticField(), DetLayer::nextLayers(), TrajectorySeed::startingState(), theForwardPropagator, theMeasurementTracker, trajectoryStateTransform::transientState(), and TrajectoryMeasurement::updatedState().
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), and CkfTrajectoryBuilder::findCompatibleMeasurements().
{
if (traj.empty())
{
//set the currentState to be the one from the trajectory seed starting point
PTrajectoryStateOnDet const & ptod = seed.startingState();
DetId id(ptod.detId());
const GeomDet * g = theMeasurementTracker->geomTracker()->idToDet(id);
const Surface * surface=&g->surface();
TSOS currentState(trajectoryStateTransform::transientState(ptod,surface,theForwardPropagator->magneticField()));
const DetLayer* lastLayer = theMeasurementTracker->geometricSearchTracker()->detLayer(id);
return StateAndLayers(currentState,lastLayer->nextLayers( *currentState.freeState(), traj.direction()) );
}
else
{
TSOS const & currentState = traj.lastMeasurement().updatedState();
return StateAndLayers(currentState,traj.lastLayer()->nextLayers( *currentState.freeState(), traj.direction()) );
}
}
| void BaseCkfTrajectoryBuilder::moveToResult | ( | TempTrajectory && | traj, |
| TempTrajectoryContainer & | result, | ||
| bool | inOut = false |
||
| ) | const [protected] |
Definition at line 192 of file BaseCkfTrajectoryBuilder.cc.
References LogDebug, and qualityFilter().
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer().
{
// quality check
if ( !qualityFilter(traj, inOut) ) return;
// discard latest dummy measurements
while (!traj.empty() && !traj.lastMeasurement().recHitR().isValid()) traj.pop();
LogDebug("CkfPattern")<<inOut<<"=inOut option. pushing a TempTrajectory with: "<<traj.foundHits()<<" found hits. "<<traj.lostHits();
// <<" lost hits. Popped :"<<(ttraj.measurements().size())-(traj.measurements().size())<<" hits.";
result.push_back(std::move(traj));
}
| bool BaseCkfTrajectoryBuilder::qualityFilter | ( | const TempTrajectory & | traj, |
| bool | inOut = false |
||
| ) | const [protected] |
Called at end of track building, to see if track should be kept
Definition at line 146 of file BaseCkfTrajectoryBuilder.cc.
References TrajectoryFilter::qualityFilter(), theFilter, and theInOutFilter.
Referenced by addToResult(), CkfTrajectoryBuilder::limitedCandidates(), and moveToResult().
{
// Called after building a trajectory is completed, to see if it is good enough
// to keep.
if (inOut) {
if (theInOutFilter == 0) edm::LogError("CkfPattern") << "CkfTrajectoryBuilder error: trying to use dedicated filter for in-out tracking phase, when none specified";
return theInOutFilter->qualityFilter(traj);
} else {
return theFilter->qualityFilter(traj);
}
}
| virtual void BaseCkfTrajectoryBuilder::rebuildTrajectories | ( | TempTrajectory const & | startingTraj, |
| const TrajectorySeed & | seed, | ||
| TrajectoryContainer & | result | ||
| ) | const [inline, virtual] |
Reimplemented in CkfTrajectoryBuilder, and GroupedCkfTrajectoryBuilder.
Definition at line 78 of file BaseCkfTrajectoryBuilder.h.
Referenced by cms::CkfTrackCandidateMakerBase::produceBase().
{ assert(0==1);}
| void BaseCkfTrajectoryBuilder::seedMeasurements | ( | const TrajectorySeed & | seed, |
| TempTrajectory & | result | ||
| ) | const [private] |
Definition at line 55 of file BaseCkfTrajectoryBuilder.cc.
References TransientTrackingRecHitBuilder::build(), GeometricSearchTracker::detLayer(), TempTrajectory::emplace(), MeasurementTracker::geometricSearchTracker(), MeasurementTracker::geomTracker(), TrackingGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), Propagator::magneticField(), Propagator::propagate(), TrajectorySeed::recHits(), TrajectorySeed::startingState(), GeomDet::surface(), theBackwardPropagator, theForwardPropagator, theMeasurementTracker, theTTRHBuilder, theUpdator, trajectoryStateTransform::transientState(), and TrajectoryStateUpdator::update().
Referenced by createStartingTrajectory().
{
TrajectorySeed::range hitRange = seed.recHits();
PTrajectoryStateOnDet pState( seed.startingState());
const GeomDet* gdet = theMeasurementTracker->geomTracker()->idToDet(pState.detId());
TSOS outerState = trajectoryStateTransform::transientState(pState, &(gdet->surface()),
theForwardPropagator->magneticField());
for (TrajectorySeed::const_iterator ihit = hitRange.first; ihit != hitRange.second; ihit++) {
TransientTrackingRecHit::RecHitPointer recHit = theTTRHBuilder->build(&(*ihit));
const GeomDet* hitGeomDet = recHit->det();
const DetLayer* hitLayer =
theMeasurementTracker->geometricSearchTracker()->detLayer(ihit->geographicalId());
TSOS invalidState( hitGeomDet->surface());
if (ihit == hitRange.second - 1) {
// the seed trajectory state should correspond to this hit
if (&gdet->surface() != &hitGeomDet->surface()) {
edm::LogError("CkfPattern") << "CkfTrajectoryBuilder error: the seed state is not on the surface of the detector of the last seed hit";
return; // FIXME: should throw exception
}
//TSOS updatedState = outerstate;
result.emplace(invalidState, outerState, recHit, 0, hitLayer);
}
else {
TSOS innerState = theBackwardPropagator->propagate(outerState,hitGeomDet->surface());
if(innerState.isValid()) {
TSOS innerUpdated = theUpdator->update(innerState,*recHit);
result.emplace(invalidState, innerUpdated, recHit, 0, hitLayer);
}
}
}
// method for debugging
// fix somehow
// fillSeedHistoDebugger(result.begin(),result.end());
}
| virtual void BaseCkfTrajectoryBuilder::setDebugger | ( | CkfDebugger * | dbg | ) | const [inline, virtual] |
Reimplemented in CkfDebugTrajectoryBuilder.
Definition at line 85 of file BaseCkfTrajectoryBuilder.h.
{;}
| void BaseCkfTrajectoryBuilder::setEvent | ( | const edm::Event & | event | ) | const [virtual] |
Implements TrajectoryBuilder.
Definition at line 240 of file BaseCkfTrajectoryBuilder.cc.
References clustersToSkip_, MeasurementTracker::setClusterToSkip(), skipClusters_, theMeasurementTracker, and MeasurementTracker::update().
Referenced by cms::CkfTrackCandidateMakerBase::produceBase().
{
theMeasurementTracker->update(event);
if (skipClusters_)
theMeasurementTracker->setClusterToSkip(clustersToSkip_,event);
}
| bool BaseCkfTrajectoryBuilder::toBeContinued | ( | TempTrajectory & | traj, |
| bool | inOut = false |
||
| ) | const [protected] |
Called after each new hit is added to the trajectory, to see if building this track should be continued
Definition at line 124 of file BaseCkfTrajectoryBuilder.cc.
References TempTrajectory::foundHits(), TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), LogTrace, TempTrajectory::lostHits(), TempTrajectory::measurements(), TrajectoryMeasurement::predictedState(), DetLayer::subDetector(), theFilter, theInOutFilter, and TrajectoryFilter::toBeContinued().
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), and CkfTrajectoryBuilder::limitedCandidates().
{
if (traj.measurements().size() > 400) {
edm::LogError("BaseCkfTrajectoryBuilder_InfiniteLoop");
LogTrace("BaseCkfTrajectoryBuilder_InfiniteLoop") <<
"Cropping Track After 400 Measurements:\n" <<
" Last predicted state: " << traj.lastMeasurement().predictedState() << "\n" <<
" Last layer subdetector: " << (traj.lastLayer() ? traj.lastLayer()->subDetector() : -1) << "\n" <<
" Found hits: " << traj.foundHits() << ", lost hits: " << traj.lostHits() << "\n\n";
return false;
}
// Called after each new hit is added to the trajectory, to see if it is
// worth continuing to build this track candidate.
if (inOut) {
if (theInOutFilter == 0) edm::LogError("CkfPattern") << "CkfTrajectoryBuilder error: trying to use dedicated filter for in-out tracking phase, when none specified";
return theInOutFilter->toBeContinued(traj);
} else {
return theFilter->toBeContinued(traj);
}
}
| void BaseCkfTrajectoryBuilder::unset | ( | ) | const [virtual] |
Reimplemented from TrajectoryBuilder.
Definition at line 247 of file BaseCkfTrajectoryBuilder.cc.
References skipClusters_, theMeasurementTracker, and MeasurementTracker::unsetClusterToSkip().
Referenced by cms::CkfTrackCandidateMakerBase::produceBase().
{
if (skipClusters_)
theMeasurementTracker->unsetClusterToSkip();
}
Definition at line 158 of file BaseCkfTrajectoryBuilder.h.
Referenced by BaseCkfTrajectoryBuilder(), and setEvent().
bool BaseCkfTrajectoryBuilder::skipClusters_ [private] |
Filter used at end of in-out tracking
Definition at line 157 of file BaseCkfTrajectoryBuilder.h.
Referenced by BaseCkfTrajectoryBuilder(), setEvent(), and unset().
const Propagator* BaseCkfTrajectoryBuilder::theBackwardPropagator [mutable, protected] |
Definition at line 143 of file BaseCkfTrajectoryBuilder.h.
Referenced by createStartingTrajectory(), GroupedCkfTrajectoryBuilder::rebuildSeedingRegion(), and seedMeasurements().
const Chi2MeasurementEstimatorBase* BaseCkfTrajectoryBuilder::theEstimator [protected] |
Definition at line 136 of file BaseCkfTrajectoryBuilder.h.
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), CkfDebugTrajectoryBuilder::analyzeMeasurementsDebugger(), MuonCkfTrajectoryBuilder::collectMeasurement(), GroupedCkfTrajectoryBuilder::estimator(), CkfTrajectoryBuilder::findCompatibleMeasurements(), CkfTrajectoryBuilder::limitedCandidates(), and MuonCkfTrajectoryBuilder::MuonCkfTrajectoryBuilder().
const TrajectoryFilter* BaseCkfTrajectoryBuilder::theFilter [private] |
Definition at line 154 of file BaseCkfTrajectoryBuilder.h.
Referenced by qualityFilter(), and toBeContinued().
const Propagator* BaseCkfTrajectoryBuilder::theForwardPropagator [mutable, protected] |
Definition at line 142 of file BaseCkfTrajectoryBuilder.h.
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), CkfDebugTrajectoryBuilder::analyzeMeasurementsDebugger(), GroupedCkfTrajectoryBuilder::buildTrajectories(), createStartingTrajectory(), MuonCkfTrajectoryBuilder::findCompatibleMeasurements(), CkfTrajectoryBuilder::findCompatibleMeasurements(), findStateAndLayers(), CkfTrajectoryBuilder::limitedCandidates(), and seedMeasurements().
const TrajectoryFilter* BaseCkfTrajectoryBuilder::theInOutFilter [private] |
Filter used at end of complete tracking
Definition at line 155 of file BaseCkfTrajectoryBuilder.h.
Referenced by qualityFilter(), and toBeContinued().
const LayerMeasurements* BaseCkfTrajectoryBuilder::theLayerMeasurements [protected] |
Definition at line 139 of file BaseCkfTrajectoryBuilder.h.
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), MuonCkfTrajectoryBuilder::collectMeasurement(), CkfTrajectoryBuilder::findCompatibleMeasurements(), and ~BaseCkfTrajectoryBuilder().
const MeasurementTracker* BaseCkfTrajectoryBuilder::theMeasurementTracker [protected] |
Definition at line 138 of file BaseCkfTrajectoryBuilder.h.
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), CkfDebugTrajectoryBuilder::analyzeMeasurementsDebugger(), MuonCkfTrajectoryBuilder::findCompatibleMeasurements(), findStateAndLayers(), CkfTrajectoryBuilder::limitedCandidates(), seedMeasurements(), setEvent(), and unset().
const Propagator* BaseCkfTrajectoryBuilder::thePropagatorAlong [protected] |
Definition at line 134 of file BaseCkfTrajectoryBuilder.h.
const Propagator* BaseCkfTrajectoryBuilder::thePropagatorOpposite [protected] |
Definition at line 135 of file BaseCkfTrajectoryBuilder.h.
const TransientTrackingRecHitBuilder* BaseCkfTrajectoryBuilder::theTTRHBuilder [protected] |
Definition at line 137 of file BaseCkfTrajectoryBuilder.h.
Referenced by CkfDebugTrajectoryBuilder::analyzeMeasurementsDebugger(), CkfTrajectoryBuilder::limitedCandidates(), and seedMeasurements().
const TrajectoryStateUpdator* BaseCkfTrajectoryBuilder::theUpdator [protected] |
Definition at line 133 of file BaseCkfTrajectoryBuilder.h.
Referenced by GroupedCkfTrajectoryBuilder::advanceOneLayer(), seedMeasurements(), CkfTrajectoryBuilder::updateTrajectory(), and GroupedCkfTrajectoryBuilder::updator().
1.7.3