#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...

virtual	~FastTSGFromIOHit ()
	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	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &	iC
	)

constructor

Definition at line 18 of file FastTSGFromIOHit.cc.

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

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

FastTSGFromIOHit::~FastTSGFromIOHit ( )

virtual

destructor

Definition at line 34 of file FastTSGFromIOHit.cc.

References LogTrace, and theCategory.

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

Member Function Documentation

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

private

Definition at line 101 of file FastTSGFromIOHit.cc.

References hiPixelPairStep_cff::deltaPhi, TrackingRegion::direction(), MuonErrorMatrixValues_cff::etaRange, RectangularEtaPhiTrackingRegion::etaRange(), PixelRecoRange< T >::max(), PixelRecoRange< T >::mean(), CoreSimTrack::momentum(), PV3DBase< T, PVType, FrameType >::phi(), RectangularEtaPhiTrackingRegion::phiMargin(), ALCARECOTkAlBeamHalo_cff::ptMin, TrackingRegion::ptMin(), TkTrackingRegionsMargin< T >::right(), and mathSSE::sqrt().

Referenced by trackerSeeds().

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

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 41 of file FastTSGFromIOHit.cc.

References clean(), edm::Event::getByToken(), TrackerSeedGenerator::getEvent(), SurveyInfoScenario_cff::seed, HiRegitMuonMixedTripletStep_cff::seedCollections, seedsTks, simTrackMatching_cfi::simTrack, FastTrackerRecHit::simTrackId(), tkConvValidator_cfi::simTracks, simTracksTk, and thePtCut.

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