#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::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) 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

Private Attributes
edm::ParameterSet	conf_

CosmicTrajectoryBuilder	cosmicTrajectoryBuilder_

CRackTrajectoryBuilder	crackTrajectoryBuilder_

std::string	geometry

bool	trinevents

bool	useHitsSplitting_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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 conf_, geometry, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and useHitsSplitting_.

                                                                   : 
     cosmicTrajectoryBuilder_(conf) ,
     crackTrajectoryBuilder_(conf) ,
     conf_(conf)
   {
     geometry=conf_.getUntrackedParameter<std::string>("GeometricStructure","STANDARD");
     useHitsSplitting_=conf.getParameter<bool>("useHitsSplitting");
     produces<TrackCandidateCollection>();
   }

cms::CosmicTrackFinder::~CosmicTrackFinder ( )

virtual

Definition at line 42 of file CosmicTrackFinder.cc.

42 { }

Member Function Documentation

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

virtual

.....

Implements edm::EDProducer.

Definition at line 45 of file CosmicTrackFinder.cc.

References alongMomentum, funct::C, conf_, cosmicTrajectoryBuilder_, gather_cfg::cout, crackTrajectoryBuilder_, edm::OwnVector< T, P >::front(), TrackingRecHit::geographicalId(), geometry, edm::EventSetup::get(), edm::Event::getByLabel(), edm::ParameterSet::getParameter(), CosmicTrajectoryBuilder::hitBuilder(), CRackTrajectoryBuilder::hitBuilder(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localParameters(), LogDebug, TrajectoryStateOnSurface::magneticField(), convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), edm::Handle< T >::product(), edm::Event::put(), DetId::rawId(), TrackingRecHit::rawId(), edm::OwnVector< T, P >::reverse(), CosmicTrajectoryBuilder::run(), CRackTrajectoryBuilder::run(), fileCollector::seed, TrajectoryStateOnSurface::surface(), patCandidatesForDimuonsSequences_cff::tracker, and useHitsSplitting_.

   {
     using namespace std  ;
     edm::InputTag matchedrecHitsTag = conf_.getParameter<edm::InputTag>("matchedRecHits");
     edm::InputTag rphirecHitsTag = conf_.getParameter<edm::InputTag>("rphirecHits");
     edm::InputTag stereorecHitsTag = conf_.getParameter<edm::InputTag>("stereorecHits");
     edm::InputTag pixelRecHitsTag = conf_.getParameter<edm::InputTag>("pixelRecHits");  
 
 
     edm::InputTag seedTag = conf_.getParameter<edm::InputTag>("cosmicSeeds");
     // retrieve seeds
     edm::Handle<TrajectorySeedCollection> seed;
     e.getByLabel(seedTag,seed);  
 
   //retrieve PixelRecHits
     static const SiPixelRecHitCollection s_empty;
     const SiPixelRecHitCollection *pixelHitCollection = &s_empty;
     edm::Handle<SiPixelRecHitCollection> pixelHits;
     if (geometry!="MTCC" && (geometry!="CRACK" )) {
       if( e.getByLabel(pixelRecHitsTag, pixelHits)) {
     pixelHitCollection = pixelHits.product();
       } else {
     edm::LogWarning("CosmicTrackFinder") << "Collection SiPixelRecHitCollection with InputTag " << pixelRecHitsTag << " cannot be found, using empty collection of same type.";
       }
     }
     
    
 
 
  //retrieve StripRecHits
     edm::Handle<SiStripMatchedRecHit2DCollection> matchedrecHits;
     e.getByLabel( matchedrecHitsTag ,matchedrecHits);
     edm::Handle<SiStripRecHit2DCollection> rphirecHits;
     e.getByLabel( rphirecHitsTag ,rphirecHits);
     edm::Handle<SiStripRecHit2DCollection> stereorecHits;
     e.getByLabel( stereorecHitsTag, 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.

Referenced by CosmicTrackFinder(), and produce().

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().

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().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation