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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 48 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(), magneticFieldToken_, muonCollectionTag_, muonCollectionToken_, skipMatchedMuons_, trackCollectionTag_, trackCollectionToken_, and trackerGeometryToken_.

                                                                                             {
  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>();
}

Member Function Documentation

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

override

Definition at line 76 of file TrajectorySeedFromMuonProducer.cc.

References RecoTauCleanerPlugins::charge, muonid::findOppositeTrack(), edm::Event::getByToken(), edm::EventSetup::getHandle(), 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(), edm::Event::put(), DetId::rawId(), dt_dqm_sourceclient_common_cff::reco, mps_fire::result, skipMatchedMuons_, HLT_FULL_cff::track, trackCollectionToken_, and trackerGeometryToken_.

                                                                                          {
  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));
}

Member Data Documentation

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

private

Definition at line 58 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

edm::InputTag TrajectorySeedFromMuonProducer::muonCollectionTag_

private

Definition at line 54 of file TrajectorySeedFromMuonProducer.cc.

Referenced by TrajectorySeedFromMuonProducer().

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

private

Definition at line 56 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

bool TrajectorySeedFromMuonProducer::skipMatchedMuons_

private

Definition at line 60 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

edm::InputTag TrajectorySeedFromMuonProducer::trackCollectionTag_

private

Definition at line 55 of file TrajectorySeedFromMuonProducer.cc.

Referenced by TrajectorySeedFromMuonProducer().

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

private

Definition at line 57 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().

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

private

Definition at line 59 of file TrajectorySeedFromMuonProducer.cc.

Referenced by produce(), and TrajectorySeedFromMuonProducer().