OutInConversionTrackFinder Class Reference

#include <OutInConversionTrackFinder.h>

Inheritance diagram for OutInConversionTrackFinder:

Public Member Functions
	OutInConversionTrackFinder (const edm::EventSetup &es, const edm::ParameterSet &config)
virtual std::vector< Trajectory >	tracks (const TrajectorySeedCollection &seeds, TrackCandidateCollection &candidates) const
virtual	~OutInConversionTrackFinder ()
Public Member Functions inherited from ConversionTrackFinder
	ConversionTrackFinder (const edm::EventSetup &es, const edm::ParameterSet &config)
void	setEvent (const edm::Event &e)
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
edm::ParameterSet	conf_
const TrajectoryBuilder *	theCkfTrajectoryBuilder_
TransientInitialStateEstimator *	theInitialState_
const MeasurementTracker *	theMeasurementTracker_
std::string	theMeasurementTrackerName_
const MagneticField *	theMF_
edm::ESHandle< Propagator >	thePropagator_
const TrackerGeometry *	theTrackerGeom_
KFUpdator *	theUpdator_
bool	useSplitHits_

Detailed Description

Id:

OutInConversionTrackFinder.h,v 1.10 2013/05/23 17:57:02 gartung Exp

Date:

2013/05/23 17:57:02

Revision:

1.10

Author

Nancy Marinelli, U. of Notre Dame, US

Definition at line 33 of file OutInConversionTrackFinder.h.

Constructor & Destructor Documentation

OutInConversionTrackFinder::OutInConversionTrackFinder	(	const edm::EventSetup &	es,
		const edm::ParameterSet &	config
	)

Definition at line 22 of file OutInConversionTrackFinder.cc.

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

                                                             :  ConversionTrackFinder( es, conf )
{


  theTrajectoryCleaner_ = new TrajectoryCleanerBySharedHits(conf);

  // get the seed cleaner
  std::string cleaner = conf_.getParameter<std::string>("OutInRedundantSeedCleaner");
  if (cleaner == "SeedCleanerByHitPosition") {
    theSeedCleaner_ = new SeedCleanerByHitPosition();
  } else if (cleaner == "CachingSeedCleanerByHitPosition") {
    theSeedCleaner_ = new CachingSeedCleanerByHitPosition();
  } else if (cleaner == "CachingSeedCleanerBySharedInput") {
    theSeedCleaner_ = new CachingSeedCleanerBySharedInput();
  } else if (cleaner == "none") {
    theSeedCleaner_ = 0;
  } else {
    throw cms::Exception("OutInRedundantSeedCleaner not found", cleaner);
  }
  

 
 
  
}

OutInConversionTrackFinder::~OutInConversionTrackFinder ( )

virtual

Definition at line 52 of file OutInConversionTrackFinder.cc.

References theSeedCleaner_, and theTrajectoryCleaner_.

                                                        {

  delete theTrajectoryCleaner_;
  if (theSeedCleaner_) delete theSeedCleaner_;

}

Member Function Documentation

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

virtual

Implements ConversionTrackFinder.

Definition at line 60 of file OutInConversionTrackFinder.cc.

References RedundantSeedCleaner::add(), TrajectoryCleanerBySharedHits::clean(), RedundantSeedCleaner::done(), PV3DBase< T, PVType, FrameType >::eta(), MeasurementTracker::geomTracker(), RedundantSeedCleaner::good(), i, TrackingGeometry::idToDet(), RedundantSeedCleaner::init(), TransientInitialStateEstimator::innerState(), TrajectoryStateOnSurface::isValid(), LogDebug, PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), PV3DBase< T, PVType, FrameType >::phi(), precomputed_value_sort(), edm::OwnVector< T, P >::push_back(), edm::OwnVector< T, P >::reserve(), query::result, edm::OwnVector< T, P >::size(), evf::utils::state, 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 == 0) {
      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.size() ) {
    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"; 
  //}
  
  
  // Converted to track candidates  
  for (std::vector<Trajectory>::const_iterator it =  result.begin(); it != result.end(); it++) {
    //    if( !it->isValid() ) continue;

    edm::OwnVector<TrackingRecHit> recHits;
    Trajectory::RecHitContainer thits;
    it->recHitsV(thits,useSplitHits_);
    recHits.reserve(thits.size());
    for (Trajectory::RecHitContainer::const_iterator hitIt = thits.begin(); hitIt != thits.end(); hitIt++) {
      recHits.push_back( (**hitIt).hit()->clone());
    }
    
    
    std::pair<TrajectoryStateOnSurface, const GeomDet*> initState =  theInitialState_->innerState( *it);
    // temporary protection againt invalid initial states
    if (! initState.first.isValid() || initState.second == 0) {
      //cout << "invalid innerState, will not make TrackCandidate" << endl;
      continue;
    }

    PTrajectoryStateOnDet state;
    if(useSplitHits_ && (initState.second != thits.front()->det()) && thits.front()->det() ){ 
      TrajectoryStateOnSurface propagated = thePropagator_->propagate(initState.first,thits.front()->det()->surface());
      if (!propagated.isValid()) continue;
      state = trajectoryStateTransform::persistentState(propagated,
                             thits.front()->det()->geographicalId().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;
  


}

Member Data Documentation

RedundantSeedCleaner* OutInConversionTrackFinder::theSeedCleaner_

private

Definition at line 52 of file OutInConversionTrackFinder.h.

Referenced by OutInConversionTrackFinder(), tracks(), and ~OutInConversionTrackFinder().

TrajectoryCleanerBySharedHits* OutInConversionTrackFinder::theTrajectoryCleaner_

private

Definition at line 51 of file OutInConversionTrackFinder.h.

Referenced by OutInConversionTrackFinder(), tracks(), and ~OutInConversionTrackFinder().