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	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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 60 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 71 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, HLT_2018_cff::magneticField, eostools::move(), HLT_2018_cff::muon, muonCollectionToken_, oppositeToMomentum, trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), edm::Event::put(), DetId::rawId(), mps_fire::result, skipMatchedMuons_, GeomDet::surface(), HLT_2018_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));
 }