TrajectorySeedFromMuonProducer Class Reference

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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magneticFieldToken_
edm::InputTag	muonCollectionTag_
edm::EDGetTokenT< edm::View< reco::Muon > >	muonCollectionToken_
bool	skipMatchedMuons_
edm::InputTag	trackCollectionTag_
edm::EDGetTokenT< reco::TrackCollection >	trackCollectionToken_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_

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

Constructor & Destructor Documentation

TrajectorySeedFromMuonProducer()

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

explicit

Definition at line 63 of file TrajectorySeedFromMuonProducer.cc.

                                                                                             {
  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_);
  magneticFieldToken_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
  trackerGeometryToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
 
  produces<TrajectorySeedCollection>();
}

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

Member Function Documentation

produce()

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

override

Definition at line 76 of file TrajectorySeedFromMuonProducer.cc.

                                                                                          {
  using namespace edm;
  using namespace reco;
  using namespace std;
 
  // Product
  auto result = std::make_unique<TrajectorySeedCollection>();
 
  edm::ESHandle<MagneticField> magneticField = iSetup.getHandle(magneticFieldToken_);
 
  edm::ESHandle<TrackerGeometry> trackerGeometry = iSetup.getHandle(trackerGeometryToken_);
 
  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));
}

References ALCARECOTkAlJpsiMuMu_cff::charge, muonid::findOppositeTrack(), edm::EventSetup::getHandle(), TrackerGeometry::idToDet(), iEvent, edm::Ref< C, T, F >::isNonnull(), LogTrace, HLT_FULL_cff::magneticField, magneticFieldToken_, eostools::move(), HLT_FULL_cff::muon, muonCollectionToken_, oppositeToMomentum, trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), DetId::rawId(), mps_fire::result, skipMatchedMuons_, GeomDet::surface(), HLT_FULL_cff::track, trackCollectionToken_, and trackerGeometryToken_.

Member Data Documentation

magneticFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TrajectorySeedFromMuonProducer::magneticFieldToken_

private

Definition at line 58 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

muonCollectionTag_

edm::InputTag TrajectorySeedFromMuonProducer::muonCollectionTag_

private

Definition at line 54 of file TrajectorySeedFromMuonProducer.cc.

Referenced by TrajectorySeedFromMuonProducer().

muonCollectionToken_

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

private

Definition at line 56 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

skipMatchedMuons_

bool TrajectorySeedFromMuonProducer::skipMatchedMuons_

private

Definition at line 60 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

trackCollectionTag_

edm::InputTag TrajectorySeedFromMuonProducer::trackCollectionTag_

private

Definition at line 55 of file TrajectorySeedFromMuonProducer.cc.

Referenced by TrajectorySeedFromMuonProducer().

trackCollectionToken_

edm::EDGetTokenT<reco::TrackCollection> TrajectorySeedFromMuonProducer::trackCollectionToken_

private

Definition at line 57 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

trackerGeometryToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TrajectorySeedFromMuonProducer::trackerGeometryToken_

private

Definition at line 59 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().