FastTSGFromL2Muon Class Reference

#include <FastTSGFromL2Muon.h>

Inheritance diagram for FastTSGFromL2Muon:

Public Member Functions
	FastTSGFromL2Muon (const edm::ParameterSet &cfg)
void	produce (edm::Event &ev, const edm::EventSetup &es) override
	~FastTSGFromL2Muon () override=default
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 Member Functions
bool	clean (reco::TrackRef muRef, RectangularEtaPhiTrackingRegion *region, const BasicTrajectorySeed *aSeed, const SimTrack &theSimTrack)

Private Attributes
const edm::EDGetTokenT< reco::TrackCollection >	l2TrackToken_
std::vector< edm::EDGetTokenT< edm::View< TrajectorySeed > > >	seedToken_
const edm::EDGetTokenT< edm::SimTrackContainer >	simTrackToken_
const edm::InputTag	theL2CollectionLabel
const double	thePtCut
std::unique_ptr< MuonTrackingRegionBuilder >	theRegionBuilder
const std::vector< edm::InputTag >	theSeedCollectionLabels
const edm::InputTag	theSimTrackCollectionLabel

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 26 of file FastTSGFromL2Muon.h.

Constructor & Destructor Documentation

FastTSGFromL2Muon()

FastTSGFromL2Muon::FastTSGFromL2Muon

(

const edm::ParameterSet &

cfg

)

Definition at line 16 of file FastTSGFromL2Muon.cc.

References looper::cfg, fileCollector::seed, seedToken_, theRegionBuilder, and theSeedCollectionLabels.

    : thePtCut(cfg.getParameter<double>("PtCut")),
      theL2CollectionLabel(cfg.getParameter<edm::InputTag>("MuonCollectionLabel")),
      theSeedCollectionLabels(cfg.getParameter<std::vector<edm::InputTag> >("SeedCollectionLabels")),
      theSimTrackCollectionLabel(cfg.getParameter<edm::InputTag>("SimTrackCollectionLabel")),
      simTrackToken_(consumes<edm::SimTrackContainer>(theSimTrackCollectionLabel)),
      l2TrackToken_(consumes<reco::TrackCollection>(theL2CollectionLabel)) {
  produces<L3MuonTrajectorySeedCollection>();

  for (auto& seed : theSeedCollectionLabels)
    seedToken_.emplace_back(consumes<edm::View<TrajectorySeed> >(seed));

  edm::ParameterSet regionBuilderPSet = cfg.getParameter<edm::ParameterSet>("MuonTrackingRegionBuilder");
  theRegionBuilder = std::make_unique<MuonTrackingRegionBuilder>(regionBuilderPSet, consumesCollector());
}

~FastTSGFromL2Muon()

FastTSGFromL2Muon::~FastTSGFromL2Muon ( )

overridedefault

Member Function Documentation

clean()

bool FastTSGFromL2Muon::clean ( reco::TrackRef  muRef, RectangularEtaPhiTrackingRegion *  region, const BasicTrajectorySeed *  aSeed, const SimTrack &  theSimTrack  )

private

Definition at line 114 of file FastTSGFromL2Muon.cc.

References srCondWrite_cfg::deltaPhi, generateTowerEtThresholdLUT::etaRange, CoreSimTrack::momentum(), ptMin, nano_mu_digi_cff::region, TkTrackingRegionsMargin< T >::right(), and mathSSE::sqrt().

Referenced by produce().

                                                           {
  // Eta cleaner
  const PixelRecoRange<float>& etaRange = region->etaRange();
  double etaSeed = theSimTrack.momentum().Eta();
  double etaLimit = (fabs(fabs(etaRange.max()) - fabs(etaRange.mean())) < 0.05)
                        ? 0.05
                        : fabs(fabs(etaRange.max()) - fabs(etaRange.mean()));
  bool inEtaRange = etaSeed >= (etaRange.mean() - etaLimit) && etaSeed <= (etaRange.mean() + etaLimit);
  if (!inEtaRange)
    return false;

  // Phi cleaner
  const TkTrackingRegionsMargin<float>& phiMargin = region->phiMargin();
  double phiSeed = theSimTrack.momentum().Phi();
  double phiLimit = (phiMargin.right() < 0.05) ? 0.05 : phiMargin.right();
  bool inPhiRange = (fabs(deltaPhi(phiSeed, double(region->direction().phi()))) < phiLimit);
  if (!inPhiRange)
    return false;

  // pt cleaner
  double ptSeed = std::sqrt(theSimTrack.momentum().Perp2());
  double ptMin = (region->ptMin() > 3.5) ? 3.5 : region->ptMin();
  bool inPtRange = ptSeed >= ptMin && ptSeed <= 2 * (muRef->pt());
  return inPtRange;
}

produce()

void FastTSGFromL2Muon::produce ( edm::Event &  ev, const edm::EventSetup &  es  )

override

Definition at line 32 of file FastTSGFromL2Muon.cc.

References edm::Range< T >::begin(), clean(), Options::const, makeMEIFBenchmarkPlots::ev, iseed, l2TrackToken_, eostools::move(), TrajectorySeed::recHits(), nano_mu_digi_cff::region, mps_fire::result, seedmultiplicitymonitor_cfi::seedCollections, seedToken_, simTrackToken_, thePtCut, theRegionBuilder, and theSeedCollectionLabels.

                                                                     {
  // Initialize the output product
  std::unique_ptr<L3MuonTrajectorySeedCollection> result(new L3MuonTrajectorySeedCollection());

  // Region builder
  theRegionBuilder->setEvent(ev, es);

  // Retrieve the Monte Carlo truth (SimTracks)
  const edm::Handle<edm::SimTrackContainer>& theSimTracks = ev.getHandle(simTrackToken_);

  // Retrieve L2 muon collection
  const edm::Handle<reco::TrackCollection>& l2muonH = ev.getHandle(l2TrackToken_);

  // Retrieve Seed collection
  unsigned seedCollections = theSeedCollectionLabels.size();
  std::vector<edm::Handle<edm::View<TrajectorySeed> > > theSeeds;
  theSeeds.resize(seedCollections);
  unsigned seed_size = 0;
  for (unsigned iseed = 0; iseed < seedCollections; ++iseed) {
    ev.getByToken(seedToken_[iseed], theSeeds[iseed]);
    seed_size += theSeeds[iseed]->size();
  }

  // Loop on L2 muons
  unsigned int imu = 0;
  unsigned int imuMax = l2muonH->size();
  edm::LogVerbatim("FastTSGFromL2Muon") << "Found " << imuMax << " L2 muons";
  for (; imu != imuMax; ++imu) {
    // Make a ref to l2 muon
    reco::TrackRef muRef(l2muonH, imu);

    // Cut on muons with low momenta
    if (muRef->pt() < thePtCut || muRef->innerMomentum().Rho() < thePtCut || muRef->innerMomentum().R() < 2.5)
      continue;

    // Define the region of interest
    std::unique_ptr<RectangularEtaPhiTrackingRegion> region = theRegionBuilder->region(muRef);

    // Copy the collection of seeds (ahem, this is time consuming!)
    std::vector<TrajectorySeed> tkSeeds;
    std::set<unsigned> tkIds;
    tkSeeds.reserve(seed_size);
    for (unsigned iseed = 0; iseed < seedCollections; ++iseed) {
      edm::Handle<edm::View<TrajectorySeed> > aSeedCollection = theSeeds[iseed];
      unsigned nSeeds = aSeedCollection->size();
      for (unsigned seednr = 0; seednr < nSeeds; ++seednr) {
        // The seed
        const BasicTrajectorySeed* aSeed = &((*aSeedCollection)[seednr]);

        // The SimTrack id associated to the first hit of the Seed
        int simTrackId = static_cast<FastTrackerRecHit const&>(*aSeed->recHits().begin()).simTrackId(0);

        // Track already associated to a seed
        std::set<unsigned>::iterator tkId = tkIds.find(simTrackId);
        if (tkId != tkIds.end())
          continue;

        const SimTrack& theSimTrack = (*theSimTracks)[simTrackId];

        if (clean(muRef, region.get(), aSeed, theSimTrack))
          tkSeeds.push_back(*aSeed);
        tkIds.insert(simTrackId);

      }  // End loop on seeds

    }  // End loop on seed collections

    // Now create the Muon Trajectory Seed
    unsigned int is = 0;
    unsigned int isMax = tkSeeds.size();
    for (; is != isMax; ++is) {
      result->push_back(L3MuonTrajectorySeed(tkSeeds[is], muRef));
    }  // End of tk seed loop

  }  // End of l2 muon loop

  edm::LogVerbatim("FastTSGFromL2Muon") << "Found " << result->size() << " seeds for muons";

  //put in the event
  ev.put(std::move(result));
}

Member Data Documentation

l2TrackToken_

const edm::EDGetTokenT<reco::TrackCollection> FastTSGFromL2Muon::l2TrackToken_

private

Definition at line 45 of file FastTSGFromL2Muon.h.

Referenced by produce().

seedToken_

std::vector<edm::EDGetTokenT<edm::View<TrajectorySeed> > > FastTSGFromL2Muon::seedToken_

private

Definition at line 46 of file FastTSGFromL2Muon.h.

Referenced by FastTSGFromL2Muon(), and produce().

simTrackToken_

const edm::EDGetTokenT<edm::SimTrackContainer> FastTSGFromL2Muon::simTrackToken_

private

Definition at line 44 of file FastTSGFromL2Muon.h.

Referenced by produce().

theL2CollectionLabel

const edm::InputTag FastTSGFromL2Muon::theL2CollectionLabel

private

Definition at line 40 of file FastTSGFromL2Muon.h.

thePtCut

const double FastTSGFromL2Muon::thePtCut

private

Definition at line 39 of file FastTSGFromL2Muon.h.

Referenced by produce().

theRegionBuilder

std::unique_ptr<MuonTrackingRegionBuilder> FastTSGFromL2Muon::theRegionBuilder

private

Definition at line 48 of file FastTSGFromL2Muon.h.

Referenced by FastTSGFromL2Muon(), and produce().

theSeedCollectionLabels

const std::vector<edm::InputTag> FastTSGFromL2Muon::theSeedCollectionLabels

private

Definition at line 41 of file FastTSGFromL2Muon.h.

Referenced by FastTSGFromL2Muon(), and produce().

theSimTrackCollectionLabel

const edm::InputTag FastTSGFromL2Muon::theSimTrackCollectionLabel

private

Definition at line 42 of file FastTSGFromL2Muon.h.