#include <FastTSGFromIOHit.h>

Inheritance diagram for FastTSGFromIOHit:

Public Member Functions
	FastTSGFromIOHit (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
	constructor More...

void	trackerSeeds (const TrackCand &, const TrackingRegion &, const TrackerTopology *tTopo, std::vector< TrajectorySeed > &) override
	generate seed(s) for a track More...

	~FastTSGFromIOHit () override
	destructor More...

Public Member Functions inherited from TrackerSeedGenerator
const edm::Event *	getEvent () const

virtual void	init (const MuonServiceProxy *service)

virtual void	setEvent (const edm::Event &)

	TrackerSeedGenerator ()

virtual void	trackerSeeds (const TrackCand &, const TrackingRegion &, const TrackerTopology *, BTSeedCollection &)

virtual	~TrackerSeedGenerator ()
	destructor More...

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

Private Attributes
std::vector< edm::EDGetTokenT< TrajectorySeedCollection > >	seedsTks

edm::EDGetTokenT< edm::SimTrackContainer >	simTracksTk

std::string	theCategory

double	thePtCut

Additional Inherited Members
Public Types inherited from TrackerSeedGenerator
typedef std::vector< TrajectorySeed >	BTSeedCollection

typedef std::pair< const Trajectory *, reco::TrackRef >	TrackCand

Protected Attributes inherited from TrackerSeedGenerator
const edm::Event *	theEvent

const MuonServiceProxy *	theProxyService

Detailed Description

Emulate TSGFromIOHit in RecoMuon

Author: Adam Everett - Purdue University

Tracker Seed Generator by propagating and updating a standAlone muon to the first pixel hits in tracker system

Emulate TSGFromIOHit in RecoMuon

Author: Adam Everett - Purdue University

Definition at line 25 of file FastTSGFromIOHit.h.

Constructor & Destructor Documentation

◆ FastTSGFromIOHit()

FastTSGFromIOHit::FastTSGFromIOHit	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &	iC
	)

constructor

Definition at line 18 of file FastTSGFromIOHit.cc.

                                                                                            {
   theCategory = "FastSimulation|Muons||FastTSGFromIOHit";
  
   thePtCut = iConfig.getParameter<double>("PtCut");
  
   simTracksTk = iC.consumes<edm::SimTrackContainer>(iConfig.getParameter<edm::InputTag>("SimTrackCollectionLabel"));
   const auto& seedLabels = iConfig.getParameter<std::vector<edm::InputTag> >("SeedCollectionLabels");
   for (const auto& seedLabel : seedLabels) {
     seedsTks.push_back(iC.consumes<TrajectorySeedCollection>(seedLabel));
   }
 }

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), seedsTks, simTracksTk, theCategory, and thePtCut.

◆ ~FastTSGFromIOHit()

FastTSGFromIOHit::~FastTSGFromIOHit ( )

override

destructor

Definition at line 30 of file FastTSGFromIOHit.cc.

30 { LogTrace(theCategory) << " FastTSGFromIOHit dtor called "; }

References LogTrace, and theCategory.

Member Function Documentation

◆ clean()

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

private

Definition at line 85 of file FastTSGFromIOHit.cc.

                                                           {
   // 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());
   if (!inPtRange) {
     return false;
   }
   return true;
 }

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

Referenced by trackerSeeds().

◆ trackerSeeds()

void FastTSGFromIOHit::trackerSeeds	(	const TrackCand &	staMuon,
		const TrackingRegion &	region,
		const TrackerTopology *	tTopo,
		std::vector< TrajectorySeed > &	result
	)

override

generate seed(s) for a track

Definition at line 32 of file FastTSGFromIOHit.cc.

                                                                        {
   // Make a ref to l2 muon
   reco::TrackRef muRef(staMuon.second);
  
   // Cut on muons with low momenta
   if (muRef->pt() < thePtCut || muRef->innerMomentum().Rho() < thePtCut || muRef->innerMomentum().R() < 2.5) {
     return;
   }
  
   // Retrieve the Monte Carlo truth (SimTracks)
   edm::Handle<edm::SimTrackContainer> simTracks;
   getEvent()->getByToken(simTracksTk, simTracks);
  
   // Retrieve Seed collection
   std::vector<edm::Handle<TrajectorySeedCollection> > seedCollections;
   seedCollections.resize(seedsTks.size());
   for (unsigned iSeed = 0; iSeed < seedsTks.size(); iSeed++) {
     getEvent()->getByToken(seedsTks[iSeed], seedCollections[iSeed]);
   }
  
   // cast the tracking region
   const RectangularEtaPhiTrackingRegion& regionRef = dynamic_cast<const RectangularEtaPhiTrackingRegion&>(region);
  
   // select and store seeds
   std::set<unsigned> simTrackIds;
   for (const auto& seeds : seedCollections) {
     for (const auto& seed : *seeds) {
       // Find the simtrack corresponding to the seed
       TrajectorySeed::range recHitRange = seed.recHits();
       const FastTrackerRecHit* firstRecHit = (const FastTrackerRecHit*)(&(*(recHitRange.first)));
       int simTrackId = firstRecHit->simTrackId(0);
       const SimTrack& simTrack = (*simTracks)[simTrackId];
  
       // skip if simTrack already associated to a seed
       if (simTrackIds.find(simTrackId) != simTrackIds.end()) {
         continue;
       }
       simTrackIds.insert(simTrackId);
  
       // filter seed
       if (!clean(muRef, regionRef, &seed, simTrack)) {
         continue;
       }
  
       // store results
       result.push_back(L3MuonTrajectorySeed(seed, muRef));
     }
   }
 }