#include <OutInConversionTrackFinder.h>

Inheritance diagram for OutInConversionTrackFinder:

Public Member Functions
	OutInConversionTrackFinder (const edm::ParameterSet &config, const BaseCkfTrajectoryBuilder *trajectoryBuilder, edm::ConsumesCollector iC)

std::vector< Trajectory >	tracks (const TrajectorySeedCollection &seeds, TrackCandidateCollection &candidates) const override

	~OutInConversionTrackFinder () override

Public Member Functions inherited from ConversionTrackFinder
	ConversionTrackFinder (const edm::ParameterSet &config, const BaseCkfTrajectoryBuilder *trajectoryBuilder, edm::ConsumesCollector iC)

void	setEventSetup (const edm::EventSetup &es)
	Initialize EventSetup objects at each event. More...

virtual	~ConversionTrackFinder ()

Private Attributes
RedundantSeedCleaner *	theSeedCleaner_

TrajectoryCleanerBySharedHits *	theTrajectoryCleaner_

Additional Inherited Members
Protected Attributes inherited from ConversionTrackFinder
const BaseCkfTrajectoryBuilder *	theCkfTrajectoryBuilder_

std::unique_ptr< TransientInitialStateEstimator >	theInitialState_

const MeasurementTracker *	theMeasurementTracker_

std::string	theMeasurementTrackerName_

edm::ESGetToken< MeasurementTracker, CkfComponentsRecord >	theMeasurementTrackerToken_

const MagneticField *	theMF_

edm::ESHandle< Propagator >	thePropagator_

edm::ESGetToken< Propagator, TrackerDigiGeometryRecord >	thePropagatorToken_

const TrackerGeometry *	theTrackerGeom_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	theTrackerGeomToken_

KFUpdator *	theUpdator_

bool	useSplitHits_

Detailed Description

Author: Nancy Marinelli, U. of Notre Dame, US

Definition at line 30 of file OutInConversionTrackFinder.h.

Constructor & Destructor Documentation

◆ OutInConversionTrackFinder()

OutInConversionTrackFinder::OutInConversionTrackFinder	(	const edm::ParameterSet &	config,
		const BaseCkfTrajectoryBuilder *	trajectoryBuilder,
		edm::ConsumesCollector	iC
	)

Definition at line 20 of file OutInConversionTrackFinder.cc.

References Exception, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theSeedCleaner_, and theTrajectoryCleaner_.

     : ConversionTrackFinder(conf, trajectoryBuilder, iC) {
   theTrajectoryCleaner_ = new TrajectoryCleanerBySharedHits(conf);
 
   // get the seed cleaner
   std::string cleaner = conf.getParameter<std::string>("OutInRedundantSeedCleaner");
   if (cleaner == "CachingSeedCleanerBySharedInput") {
     theSeedCleaner_ = new CachingSeedCleanerBySharedInput();
   } else if (cleaner == "none") {
     theSeedCleaner_ = nullptr;
   } else {
     throw cms::Exception("OutInRedundantSeedCleaner not found, please use CachingSeedCleanerBySharedInput or none",
                          cleaner);
   }
 }

◆ ~OutInConversionTrackFinder()

OutInConversionTrackFinder::~OutInConversionTrackFinder ( )

override

Definition at line 38 of file OutInConversionTrackFinder.cc.

References theSeedCleaner_, and theTrajectoryCleaner_.

                                                         {
   delete theTrajectoryCleaner_;
   if (theSeedCleaner_)
     delete theSeedCleaner_;
 }

Member Function Documentation

◆ tracks()

std::vector< Trajectory > OutInConversionTrackFinder::tracks	(	const TrajectorySeedCollection &	seeds,
		TrackCandidateCollection &	candidates
	)		const

overridevirtual

Implements ConversionTrackFinder.

Definition at line 44 of file OutInConversionTrackFinder.cc.

References RedundantSeedCleaner::add(), alongMomentum, cms::cuda::assert(), TrajectoryCleanerBySharedHits::clean(), RedundantSeedCleaner::done(), PV3DBase< T, PVType, FrameType >::eta(), MeasurementTracker::geomTracker(), RedundantSeedCleaner::good(), BaseCkfTrajectoryBuilder::hitBuilder(), mps_fire::i, TrackerGeometry::idToDet(), RedundantSeedCleaner::init(), TrajectoryStateOnSurface::isValid(), LogDebug, PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), PV3DBase< T, PVType, FrameType >::phi(), precomputed_value_sort(), Propagator::propagate(), FastTrackerRecHitMaskProducer_cfi::recHits, mps_fire::result, ConversionTrackFinder::theCkfTrajectoryBuilder_, ConversionTrackFinder::theInitialState_, ConversionTrackFinder::theMeasurementTracker_, ConversionTrackFinder::thePropagator_, theSeedCleaner_, theTrajectoryCleaner_, TrajectoryBuilder::trajectories(), and ConversionTrackFinder::useSplitHits_.

Referenced by ConversionTrackCandidateProducer::buildCollections().

                                                                                                      {
   //  std::cout  << "OutInConversionTrackFinder::tracks getting " <<  outInSeeds.size() << " Out-In seeds " << "\n";;
 
   std::vector<Trajectory> tmpO;
   tmpO.erase(tmpO.begin(), tmpO.end());
 
   std::vector<Trajectory> result;
   result.erase(result.begin(), result.end());
 
   std::vector<Trajectory> rawResult;
   if (theSeedCleaner_)
     theSeedCleaner_->init(&rawResult);
 
   /*
   for(TrajectorySeedCollection::const_iterator iSeed=outInSeeds.begin(); iSeed!=outInSeeds.end();iSeed++){
     DetId tmpId = DetId( iSeed->startingState().detId());
     const GeomDet* tmpDet  = theMeasurementTracker_->geomTracker()->idToDet( tmpId );
     GlobalVector gv = tmpDet->surface().toGlobal( iSeed->startingState().parameters().momentum() );
     
     LogDebug("OutInConversionTrackFinder") << " OutInConversionTrackFinder::tracks hits in the seed " << iSeed->nHits() << "\n";
     LogDebug("OutInConversionTrackFinder")<< " OutInConversionTrackFinder::tracks seed starting state position  " << iSeed->startingState().parameters().position() << " momentum " <<  iSeed->startingState().parameters().momentum() << " charge " << iSeed->startingState().parameters().charge() << " R " << gv.perp() << " eta " << gv.eta() << " phi " << gv.phi() << "\n";
     
     TrajectorySeed::range hitRange = iSeed->recHits();
     for (TrajectorySeed::const_iterator ihit = hitRange.first; ihit != hitRange.second; ihit++) {
       
       if ( ihit->isValid() ) {
     
     LogDebug("OutInConversionTrackFinder")  << " Valid hit global position " << theMeasurementTracker_->geomTracker()->idToDet((ihit)->geographicalId())->surface().toGlobal((ihit)->localPosition()) << " R " << theMeasurementTracker_->geomTracker()->idToDet((ihit)->geographicalId())->surface().toGlobal((ihit)->localPosition()).perp() << " phi " << theMeasurementTracker_->geomTracker()->idToDet((ihit)->geographicalId())->surface().toGlobal((ihit)->localPosition()).phi() << " eta " << theMeasurementTracker_->geomTracker()->idToDet((ihit)->geographicalId())->surface().toGlobal((ihit)->localPosition()).eta() <<    "\n" ;
     
       }
     }
   } 
   
   */
 
   int goodSeed = 0;
   std::vector<Trajectory> theTmpTrajectories;
   for (TrajectorySeedCollection::const_iterator iSeed = outInSeeds.begin(); iSeed != outInSeeds.end(); iSeed++) {
     theTmpTrajectories.clear();
 
     if (!theSeedCleaner_ || theSeedCleaner_->good(&(*iSeed))) {
       goodSeed++;
 
       DetId tmpId = DetId(iSeed->startingState().detId());
       const GeomDet* tmpDet = theMeasurementTracker_->geomTracker()->idToDet(tmpId);
       GlobalVector gv = tmpDet->surface().toGlobal(iSeed->startingState().parameters().momentum());
 
       //      std::cout << " OutInConversionTrackFinder::tracks hits in the seed " << iSeed->nHits() << "\n";
       LogDebug("OutInConversionTrackFinder") << " OutInConversionTrackFinder::tracks seed starting state position  "
                                              << iSeed->startingState().parameters().position() << " momentum "
                                              << iSeed->startingState().parameters().momentum() << " charge "
                                              << iSeed->startingState().parameters().charge() << " R " << gv.perp()
                                              << " eta " << gv.eta() << " phi " << gv.phi() << "\n";
 
       theCkfTrajectoryBuilder_->trajectories(*iSeed, theTmpTrajectories);
 
       LogDebug("OutInConversionTrackFinder")
           << "OutInConversionTrackFinder::track returned " << theTmpTrajectories.size() << " trajectories"
           << "\n";
 
       theTrajectoryCleaner_->clean(theTmpTrajectories);
 
       for (std::vector<Trajectory>::const_iterator it = theTmpTrajectories.begin(); it != theTmpTrajectories.end();
            it++) {
         if (it->isValid()) {
           rawResult.push_back(*it);
           if (theSeedCleaner_)
             theSeedCleaner_->add(&(*it));
         }
       }
     }
   }  // end loop over the seeds
   LogDebug("OutInConversionTrackFinder") << " OutInConversionTrackFinder::track Good seeds " << goodSeed << "\n";
   LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder::track rawResult size after cleaning "
                                          << rawResult.size() << "\n";
 
   if (theSeedCleaner_)
     theSeedCleaner_->done();
 
   std::vector<Trajectory> unsmoothedResult;
   theTrajectoryCleaner_->clean(rawResult);
 
   for (std::vector<Trajectory>::const_iterator itraw = rawResult.begin(); itraw != rawResult.end(); itraw++) {
     if ((*itraw).isValid()) {
       tmpO.push_back(*itraw);
       LogDebug("OutInConversionTrackFinder") << " rawResult num hits " << (*itraw).foundHits() << "\n";
     }
   }
 
   LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  tmpO size " << tmpO.size()
                                          << " before sorting "
                                          << "\n";
   //  for (std::vector<Trajectory>::const_iterator it =tmpO.begin(); it != tmpO.end(); it++) {
   // LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  tmpO num of hits " << (*it).foundHits() << " before ordering " << "\n";
   // }
 
   precomputed_value_sort(tmpO.begin(), tmpO.end(), ExtractNumOfHits());
 
   LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  tmpO after sorting "
                                          << "\n";
   //  for (std::vector<Trajectory>::const_iterator it =tmpO.begin(); it != tmpO.end(); it++) {
   // LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  tmpO  num of hits " << (*it).foundHits() << "\n";
   // }
 
   for (int i = tmpO.size() - 1; i >= 0; i--) {
     unsmoothedResult.push_back(tmpO[i]);
   }
   LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  unsmoothedResult size  "
                                          << unsmoothedResult.size() << "\n";
 
   // for (std::vector<Trajectory>::const_iterator it =  unsmoothedResult.begin(); it !=  unsmoothedResult.end(); it++) {
   // LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  unsmoothedResult  after reordering " <<(*it).foundHits() <<  "\n";
   //  }
 
   // Check if the inner state is valid
   tmpO.clear();
   LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  tmpO size " << tmpO.size()
                                          << " after clearing "
                                          << "\n";
   for (std::vector<Trajectory>::const_iterator it = unsmoothedResult.begin(); it != unsmoothedResult.end(); it++) {
     if (!it->isValid())
       continue;
 
     std::pair<TrajectoryStateOnSurface, const GeomDet*> initState = theInitialState_->innerState(*it);
     //  LogDebug("OutInConversionTrackFinder") << " Initial state parameters " << initState.first << "\n";
 
     // temporary protection againt invalid initial states
     if (!initState.first.isValid() || initState.second == nullptr) {
       LogDebug("OutInConversionTrackFinder") << "invalid innerState, will not make TrackCandidate"
                                              << "\n";
       ;
       continue;
     }
     tmpO.push_back(*it);
   }
 
   LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  tmpO size " << tmpO.size() << " after filling "
                                          << "\n";
   if (!tmpO.empty()) {
     std::vector<Trajectory>::iterator it = tmpO.begin();
 
     // only send out the two best tracks
     result.push_back(*it);
     if (tmpO.size() > 1)
       result.push_back(*(++it));
   }
 
   //  for (std::vector<Trajectory>::const_iterator it =result.begin(); it != result.end(); it++) {
   // LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  Result  num of hits " << (*it).foundHits() << "\n";
   //}
 
   // Convert to TrackCandidates and fill in the output_p
   Traj2TrackHits t2t(theCkfTrajectoryBuilder_->hitBuilder(), true);
   for (std::vector<Trajectory>::const_iterator it = result.begin(); it != result.end(); it++) {
     edm::OwnVector<TrackingRecHit> recHits;
     if (it->direction() == alongMomentum)
       LogDebug("OutInConversionTrackFinder") << "OutInConv along momentum... " << std::endl;
     t2t(*it, recHits, useSplitHits_);
 
     std::pair<TrajectoryStateOnSurface, const GeomDet*> initState = theInitialState_->innerState(*it);
 
     assert(initState.second == recHits.front().det());
 
     // temporary protection againt invalid initial states
     if (!initState.first.isValid() || initState.second == nullptr) {
       LogDebug("OutInConversionTrackFinder") << "invalid innerState, will not make TrackCandidate" << std::endl;
       continue;
     }
 
     PTrajectoryStateOnDet state;
     if (useSplitHits_ && (initState.second != recHits.front().det()) && recHits.front().det()) {
       TrajectoryStateOnSurface propagated =
           thePropagator_->propagate(initState.first, recHits.front().det()->surface());
       if (!propagated.isValid())
         continue;
       state = trajectoryStateTransform::persistentState(propagated, recHits.front().rawId());
     } else
       state = trajectoryStateTransform::persistentState(initState.first, initState.second->geographicalId().rawId());
 
     LogDebug("OutInConversionTrackFinder") << "OutInConversionTrackFinder  Number of hits for the track candidate "
                                            << recHits.size() << " TSOS charge " << initState.first.charge() << "\n";
 
     output_p.push_back(TrackCandidate(recHits, it->seed(), state));
   }
 
   //  std::cout << "  Returning " << result.size() << "Out In Trajectories  " << "\n";
 
   return result;
 }