Inheritance diagram for TrajectorySeedFromMuonProducer:

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

	TrajectorySeedFromMuonProducer (const edm::ParameterSet &)

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
edm::InputTag	muonCollectionTag_

edm::EDGetTokenT< edm::View< reco::Muon > >	muonCollectionToken_

bool	skipMatchedMuons_

edm::InputTag	trackCollectionTag_

edm::EDGetTokenT< reco::TrackCollection >	trackCollectionToken_

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

Detailed Description

Definition at line 47 of file TrajectorySeedFromMuonProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 63 of file TrajectorySeedFromMuonProducer.cc.

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

 {
   muonCollectionTag_  = iConfig.getParameter<edm::InputTag>("muonCollectionTag");
   trackCollectionTag_ = iConfig.getParameter<edm::InputTag>("trackCollectionTag");
   skipMatchedMuons_   = iConfig.getParameter<bool>("skipMatchedMuons");
   
   muonCollectionToken_ = consumes<edm::View<reco::Muon> >(muonCollectionTag_);
   trackCollectionToken_ = consumes<reco::TrackCollection> (trackCollectionTag_);
 
   produces<TrajectorySeedCollection>();
 }

Member Function Documentation

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

override

Definition at line 76 of file TrajectorySeedFromMuonProducer.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, DEFINE_FWK_MODULE, muonid::findOppositeTrack(), edm::EventSetup::get(), edm::Event::getByToken(), TrackerGeometry::idToDet(), edm::Ref< C, T, F >::isNonnull(), LogTrace, seedCreatorFromRegionConsecutiveHitsEDProducer_cff::magneticField, eostools::move(), metsig::muon, muonCollectionToken_, oppositeToMomentum, trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), edm::Event::put(), DetId::rawId(), mps_fire::result, skipMatchedMuons_, GeomDet::surface(), HiIsolationCommonParameters_cff::track, and trackCollectionToken_.

 {
   using namespace edm;
   using namespace reco;
   using namespace std;
 
   // Product
   auto result = std::make_unique<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.getByToken(muonCollectionToken_, muonCollectionHandle);
   
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   iEvent.getByToken(trackCollectionToken_, 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(std::move(result));
 
 }