#include <CosmicTrackFinder.h>

Inheritance diagram for cms::CosmicTrackFinder:

Public Member Functions
	CosmicTrackFinder (const edm::ParameterSet &conf)

virtual void	produce (edm::Event &e, const edm::EventSetup &c)

virtual	~CosmicTrackFinder ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef TrajectoryStateOnSurface	TSOS

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 67 of file CosmicTrackFinder.h.

Member Typedef Documentation

typedef TrajectoryStateOnSurface cms::CosmicTrackFinder::TSOS

private

Definition at line 70 of file CosmicTrackFinder.h.

Constructor & Destructor Documentation

cms::CosmicTrackFinder::CosmicTrackFinder ( const edm::ParameterSet & conf )

explicit

Definition at line 30 of file CosmicTrackFinder.cc.

References geometry, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), matchedrecHitsToken_, pixelRecHitsToken_, rphirecHitsToken_, seedToken_, stereorecHitsToken_, AlCaHLTBitMon_QueryRunRegistry::string, and useHitsSplitting_.

                                                                   : 
     cosmicTrajectoryBuilder_(conf) ,
     crackTrajectoryBuilder_(conf)
   {
     geometry=conf.getUntrackedParameter<std::string>("GeometricStructure","STANDARD");
     useHitsSplitting_=conf.getParameter<bool>("useHitsSplitting");
     matchedrecHitsToken_ = consumes<SiStripMatchedRecHit2DCollection>(
         conf.getParameter<edm::InputTag>("matchedRecHits"));
     rphirecHitsToken_ = consumes<SiStripRecHit2DCollection>(
         conf.getParameter<edm::InputTag>("rphirecHits"));
     stereorecHitsToken_ = consumes<SiStripRecHit2DCollection>(
         conf.getParameter<edm::InputTag>("stereorecHits"));
     pixelRecHitsToken_ = consumes<SiPixelRecHitCollection>(
         conf.getParameter<edm::InputTag>("pixelRecHits"));
     seedToken_ = consumes<TrajectorySeedCollection>(
         conf.getParameter<edm::InputTag>("cosmicSeeds"));
 
     produces<TrackCandidateCollection>();
   }

cms::CosmicTrackFinder::~CosmicTrackFinder ( )

virtual

Definition at line 52 of file CosmicTrackFinder.cc.

52 { }

Member Function Documentation

void cms::CosmicTrackFinder::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

.....

Implements edm::stream::EDProducerBase.

Definition at line 55 of file CosmicTrackFinder.cc.

References alongMomentum, assert(), funct::C, cosmicTrajectoryBuilder_, gather_cfg::cout, crackTrajectoryBuilder_, edm::OwnVector< T, P >::front(), TrackingRecHit::geographicalId(), geometry, edm::EventSetup::get(), edm::Event::getByToken(), CosmicTrajectoryBuilder::hitBuilder(), CRackTrajectoryBuilder::hitBuilder(), TrajectoryStateOnSurface::isValid(), cmsLHEtoEOSManager::l, edm::EDConsumerBase::labelsForToken(), TrajectoryStateOnSurface::localParameters(), LogDebug, TrajectoryStateOnSurface::magneticField(), matchedrecHitsToken_, edm::ProductLabels::module, convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), beampixel_dqm_sourceclient-live_cfg::pixelHitCollection, pixelRecHitsToken_, edm::Event::put(), DetId::rawId(), TrackingRecHit::rawId(), HLT_FULL_cff::recHits, edm::OwnVector< T, P >::reverse(), rphirecHitsToken_, CosmicTrajectoryBuilder::run(), CRackTrajectoryBuilder::run(), fileCollector::seed, seedToken_, stereorecHitsToken_, TrajectoryStateOnSurface::surface(), patCandidatesForDimuonsSequences_cff::tracker, and useHitsSplitting_.

   {
     using namespace std;
 
     // retrieve seeds
     edm::Handle<TrajectorySeedCollection> seed;
     e.getByToken(seedToken_, seed);
 
     //retrieve PixelRecHits
     static const SiPixelRecHitCollection s_empty;
     const SiPixelRecHitCollection *pixelHitCollection = &s_empty;
     edm::Handle<SiPixelRecHitCollection> pixelHits;
     if (geometry!="MTCC" && (geometry!="CRACK" )) {
       if( e.getByToken(pixelRecHitsToken_, pixelHits)) {
     pixelHitCollection = pixelHits.product();
       } else {
         Labels l;
         labelsForToken(pixelRecHitsToken_, l);
     edm::LogWarning("CosmicTrackFinder")
             << "Collection SiPixelRecHitCollection with InputTag "
             << l.module
             << " cannot be found, using empty collection of same type.";
       }
     }
     
    
 
 
  //retrieve StripRecHits
     edm::Handle<SiStripMatchedRecHit2DCollection> matchedrecHits;
     e.getByToken(matchedrecHitsToken_, matchedrecHits);
     edm::Handle<SiStripRecHit2DCollection> rphirecHits;
     e.getByToken(rphirecHitsToken_, rphirecHits);
     edm::Handle<SiStripRecHit2DCollection> stereorecHits;
     e.getByToken(stereorecHitsToken_, stereorecHits);
 
     // Step B: create empty output collection
     std::auto_ptr<TrackCandidateCollection> output(new TrackCandidateCollection);
 
     edm::ESHandle<TrackerGeometry> tracker;
     es.get<TrackerDigiGeometryRecord>().get(tracker);
     edm::LogVerbatim("CosmicTrackFinder") << "========== Cosmic Track Finder Info ==========";
     edm::LogVerbatim("CosmicTrackFinder") << " Numbers of Seeds " << (*seed).size();
     if((*seed).size()>0){
 
       std::vector<Trajectory> trajoutput;
       
 
       const TransientTrackingRecHitBuilder * hitBuilder = nullptr;
       if(geometry!="CRACK" ) {
         cosmicTrajectoryBuilder_.run(*seed,
                                    *stereorecHits,
                                    *rphirecHits,
                                    *matchedrecHits,
                    *pixelHitCollection,
                                    es,
                                    e,
                                    trajoutput);
         hitBuilder = cosmicTrajectoryBuilder_.hitBuilder();
       } else {
         crackTrajectoryBuilder_.run(*seed,
                                    *stereorecHits,
                                    *rphirecHits,
                                    *matchedrecHits,
                    *pixelHitCollection,
                                    es,
                                    e,
                                    trajoutput);
     hitBuilder = crackTrajectoryBuilder_.hitBuilder();
       }
       assert(hitBuilder);
       Traj2TrackHits t2t(hitBuilder,true);
 
       edm::LogVerbatim("CosmicTrackFinder") << " Numbers of Temp Trajectories " << trajoutput.size();
       edm::LogVerbatim("CosmicTrackFinder") << "========== END Info ==========";
       if(trajoutput.size()>0){
 
         // crazyness...
     std::vector<Trajectory*> tmpTraj;
     std::vector<Trajectory>::iterator itr;
     for (itr=trajoutput.begin();itr!=trajoutput.end();itr++)tmpTraj.push_back(&(*itr));
 
     //The best track is selected
     //FOR MTCC the criteria are:
     //1)# of layers,2) # of Hits,3)Chi2
     if (geometry=="MTCC")  stable_sort(tmpTraj.begin(),tmpTraj.end(),CompareTrajLay());
     else  stable_sort(tmpTraj.begin(),tmpTraj.end(),CompareTrajChi());
 
 
     const Trajectory  theTraj = *(*tmpTraj.begin());
     bool seedplus=(theTraj.seed().direction()==alongMomentum);
 
         // std::cout << "CosmicTrackFinder " <<"Reconstruction " <<  (seedplus ? "along" : "opposite to") << " momentum" << std::endl;
         LogDebug("CosmicTrackFinder")<<"Reconstruction " <<  (seedplus ? "along" : "opposite to") << " momentum";
 
     /*
     // === the convention is to save always final tracks with hits sorted *along* momentum
         // --- this is NOT what was necessaraly happening before and not consistent with what done in standard CKF (this is a candidate not a track)    
     */
          edm::OwnVector<TrackingRecHit> recHits;
          if(theTraj.direction() == alongMomentum) std::cout << "cosmic: along momentum... " << std::endl;
          t2t(theTraj,recHits,useHitsSplitting_);
          recHits.reverse();  // according to original code
 
         /*
     Trajectory::RecHitContainer thits;
     //it->recHitsV(thits);
         theTraj.recHitsV(thits,useHitsSplitting_);
     edm::OwnVector<TrackingRecHit> recHits;
     recHits.reserve(thits.size());
     // reverse hit order
     for (Trajectory::RecHitContainer::const_iterator hitIt = thits.end()-1;
          hitIt >= thits.begin(); hitIt--) {
       recHits.push_back( (**hitIt).hit()->clone());
     }
        */
 
     TSOS firstState;
     unsigned int firstId;
 
         // assume not along momentum....
     firstState=theTraj.lastMeasurement().updatedState();
         firstId = theTraj.lastMeasurement().recHitR().rawId();
     //firstId = recHits.front().rawId();
 
     /*
     cout << "firstState y, z: " << firstState.globalPosition().y() 
          << " , " << firstState.globalPosition().z() <<  endl;
 
     */
 
     AlgebraicSymMatrix55 C = AlgebraicMatrixID();
     TSOS startingState( firstState.localParameters(), LocalTrajectoryError(C),
                 firstState.surface(),
                 firstState.magneticField() );
 
     // protection againt invalid initial states
     if (! firstState.isValid()) {
       edm::LogWarning("CosmicTrackFinder") << "invalid innerState, will not make TrackCandidate";
       edm::OrphanHandle<TrackCandidateCollection> rTrackCand = e.put(output);  
       return;
     }
 
         // FIXME in case of slitting this can happen see CkfTrackCandidateMakerBase
     if(firstId != recHits.front().rawId()){
       edm::LogWarning("CosmicTrackFinder") <<"Mismatch in DetID of first hit: firstID= " <<firstId
                            << "   DetId= " << recHits.front().geographicalId().rawId();
       edm::OrphanHandle<TrackCandidateCollection> rTrackCand = e.put(output);  
       return;
     }
 
     PTrajectoryStateOnDet const &  state = trajectoryStateTransform::persistentState( startingState, firstId);
     
     
     output->push_back(TrackCandidate(recHits,theTraj.seed(),state,theTraj.seedRef() ) );
     
       }
 
     }
     edm::OrphanHandle<TrackCandidateCollection> rTrackCand = e.put(output);  
   }

Member Data Documentation

edm::ParameterSet cms::CosmicTrackFinder::conf_

private

Definition at line 82 of file CosmicTrackFinder.h.

CosmicTrajectoryBuilder cms::CosmicTrackFinder::cosmicTrajectoryBuilder_

private

Definition at line 80 of file CosmicTrackFinder.h.

Referenced by produce().

CRackTrajectoryBuilder cms::CosmicTrackFinder::crackTrajectoryBuilder_

private

Definition at line 81 of file CosmicTrackFinder.h.

Referenced by produce().

std::string cms::CosmicTrackFinder::geometry

private

Definition at line 83 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), Vispa.Gui.ConnectableWidget.ConnectableWidget::leaveEvent(), and produce().

edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> cms::CosmicTrackFinder::matchedrecHitsToken_

private

Definition at line 86 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiPixelRecHitCollection> cms::CosmicTrackFinder::pixelRecHitsToken_

private

Definition at line 89 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiStripRecHit2DCollection> cms::CosmicTrackFinder::rphirecHitsToken_

private

Definition at line 87 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<TrajectorySeedCollection> cms::CosmicTrackFinder::seedToken_

private

Definition at line 90 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiStripRecHit2DCollection> cms::CosmicTrackFinder::stereorecHitsToken_

private

Definition at line 88 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

bool cms::CosmicTrackFinder::trinevents

private

Definition at line 84 of file CosmicTrackFinder.h.

bool cms::CosmicTrackFinder::useHitsSplitting_

private

Definition at line 85 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

Private Attributes
edm::ParameterSet	conf_

CosmicTrajectoryBuilder	cosmicTrajectoryBuilder_

CRackTrajectoryBuilder	crackTrajectoryBuilder_

std::string	geometry

edm::EDGetTokenT < SiStripMatchedRecHit2DCollection >	matchedrecHitsToken_

edm::EDGetTokenT < SiPixelRecHitCollection >	pixelRecHitsToken_

edm::EDGetTokenT < SiStripRecHit2DCollection >	rphirecHitsToken_

edm::EDGetTokenT < TrajectorySeedCollection >	seedToken_

edm::EDGetTokenT < SiStripRecHit2DCollection >	stereorecHitsToken_

bool	trinevents

bool	useHitsSplitting_

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation