Inheritance diagram for TrajectorySeedFromMuonProducer:

Public Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)

	TrajectorySeedFromMuonProducer (const edm::ParameterSet &)

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::InputTag	muonCollectionTag_

bool	skipMatchedMuons_

edm::InputTag	trackCollectionTag_

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

typedef WorkerT< EDProducer >	WorkerType

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::EDProducer
CurrentProcessingContext const *	currentContext () const

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 48 of file TrajectorySeedFromMuonProducer.cc.

Constructor & Destructor Documentation

TrajectorySeedFromMuonProducer::TrajectorySeedFromMuonProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 61 of file TrajectorySeedFromMuonProducer.cc.

References edm::ParameterSet::getParameter(), muonCollectionTag_, skipMatchedMuons_, and trackCollectionTag_.

 {
   muonCollectionTag_  = iConfig.getParameter<edm::InputTag>("muonCollectionTag");
   trackCollectionTag_ = iConfig.getParameter<edm::InputTag>("trackCollectionTag");
   skipMatchedMuons_   = iConfig.getParameter<bool>("skipMatchedMuons");
   produces<TrajectorySeedCollection>();
 }

Member Function Documentation

void TrajectorySeedFromMuonProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDProducer.

Definition at line 70 of file TrajectorySeedFromMuonProducer.cc.

References DeDxDiscriminatorTools::charge(), muonid::findOppositeTrack(), edm::EventSetup::get(), edm::Event::getByLabel(), edm::Ref< C, T, F >::isNonnull(), LogTrace, metsig::muon, muonCollectionTag_, oppositeToMomentum, trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), edm::Event::put(), DetId::rawId(), L1Trigger_dataformats::reco, query::result, skipMatchedMuons_, and trackCollectionTag_.

 {
   using namespace edm;
   using namespace reco;
   using namespace std;
 
   // Product
   std::auto_ptr<TrajectorySeedCollection> result(new TrajectorySeedCollection());
   
   edm::ESHandle<MagneticField> magneticField;
   iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
   
   edm::ESHandle<TrackerGeometry> trackerGeometry;
   iSetup.get<TrackerDigiGeometryRecord>().get(trackerGeometry);
   
   
   
   edm::Handle<edm::View<Muon> > muonCollectionHandle;
   iEvent.getByLabel(muonCollectionTag_, muonCollectionHandle);
   
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   iEvent.getByLabel(trackCollectionTag_, trackCollectionHandle);
 
   // Loop over the muon track
   for (edm::View<Muon>::const_iterator muon = muonCollectionHandle->begin();  muon != muonCollectionHandle->end(); ++muon) {
     // muon must have a tracker track
     if( muon->innerTrack().isNull() ) continue;
     edm::RefToBase<reco::Track> track(muon->innerTrack());
     // check if there is a back-to-back track
     if( skipMatchedMuons_ &&
     muonid::findOppositeTrack(trackCollectionHandle, *track).isNonnull() ) continue;
     if( (!track->innerOk()) || (!track->recHit(0)->isValid())) continue;
     GlobalPoint innerPosition(track->innerPosition().x(), track->innerPosition().y(), track->innerPosition().z());
     GlobalVector innerMomentum(track->innerMomentum().x(), track->innerMomentum().y(), track->innerMomentum().z());
     int charge = track->charge();
     const reco::Track::CovarianceMatrix innerStateCovariance = track->innerStateCovariance();
     DetId innerDetId = DetId(track->innerDetId());
     // Construct the information necessary to make a TrajectoryStateOnSurface
     GlobalTrajectoryParameters globalTrajParams(innerPosition, innerMomentum, charge, &(*magneticField));
     CurvilinearTrajectoryError curviError(innerStateCovariance);
     FreeTrajectoryState tracker_state(globalTrajParams, curviError);
     LogTrace("MuonIdentification") << "Track Inner FTS: "<<tracker_state;
     
     TrajectoryStateOnSurface tracker_tsos  = TrajectoryStateOnSurface(globalTrajParams, curviError, trackerGeometry->idToDet(innerDetId)->surface());
     
     // Make Hits, push back the innermost Hit
     edm::OwnVector<TrackingRecHit> trackHits;
     trackHits.push_back(track->recHit(0)->clone() );
     
     // Make TrajectorySeed
     PTrajectoryStateOnDet const & PTraj = trajectoryStateTransform::persistentState(tracker_tsos, innerDetId.rawId());
     TrajectorySeed trajectorySeed(PTraj, trackHits, oppositeToMomentum);
     LogTrace("MuonIdentification")<<"Trajectory Seed Direction: "<< trajectorySeed.direction()<<endl;
     result->push_back(trajectorySeed);
   }
   
   iEvent.put(result);
 
 }