#include <TSGFromL1Muon.h>

Inheritance diagram for TSGFromL1Muon:

Public Member Functions
void	beginRun (const edm::Run &run, const edm::EventSetup &es) override

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

	TSGFromL1Muon (const edm::ParameterSet &cfg)

	~TSGFromL1Muon () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
edm::ParameterSet	theConfig

edm::EDGetTokenT< PixelTrackFilter >	theFilterToken

std::unique_ptr< L1MuonPixelTrackFitter >	theFitter

OrderedHitsGenerator *	theHitGenerator

L1MuonSeedsMerger *	theMerger

std::unique_ptr< L1MuonRegionProducer >	theRegionProducer

edm::InputTag	theSourceTag

edm::EDGetTokenT< l1extra::L1MuonParticleCollection >	theSourceToken

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

Description: EDPRoducer to generate L3MuonTracjectorySeed from L1MuonParticles

Author: Marcin Konecki

Definition at line 24 of file TSGFromL1Muon.h.

Constructor & Destructor Documentation

TSGFromL1Muon::TSGFromL1Muon ( const edm::ParameterSet & cfg )

Definition at line 38 of file TSGFromL1Muon.cc.

References reco::get(), edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theConfig, theFilterToken, theFitter, theHitGenerator, theRegionProducer, theSourceTag, and theSourceToken.

   : theConfig(cfg),theHitGenerator(nullptr),theMerger(nullptr)
 {
   produces<L3MuonTrajectorySeedCollection>();
   theSourceTag = cfg.getParameter<edm::InputTag>("L1MuonLabel");
 
   edm::ConsumesCollector iC = consumesCollector();
   theFilterToken = consumes<PixelTrackFilter>(cfg.getParameter<edm::InputTag>("Filter"));
 
   edm::ParameterSet hitsfactoryPSet =
       theConfig.getParameter<edm::ParameterSet>("OrderedHitsFactoryPSet");
   std::string hitsfactoryName = hitsfactoryPSet.getParameter<std::string>("ComponentName");
   theHitGenerator = OrderedHitsGeneratorFactory::get()->create( hitsfactoryName, hitsfactoryPSet, iC);
 
   theSourceToken=iC.consumes<L1MuonParticleCollection>(theSourceTag);
 
   theRegionProducer = std::make_unique<L1MuonRegionProducer>(theConfig.getParameter<edm::ParameterSet>("RegionFactoryPSet"));
   theFitter = std::make_unique<L1MuonPixelTrackFitter>(cfg.getParameter<edm::ParameterSet>("FitterPSet"));
 }

TSGFromL1Muon::~TSGFromL1Muon ( )

override

Definition at line 58 of file TSGFromL1Muon.cc.

References theHitGenerator, and theMerger.

 {
   delete theMerger;
   delete theHitGenerator;
 }

Member Function Documentation

void TSGFromL1Muon::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	es
	)

override

Definition at line 64 of file TSGFromL1Muon.cc.

References edm::ParameterSet::getParameter(), theConfig, and theMerger.

 {
   edm::ParameterSet cleanerPSet = theConfig.getParameter<edm::ParameterSet>("CleanerPSet");
   theMerger = new L1MuonSeedsMerger(cleanerPSet);
 }

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

override

Definition at line 71 of file TSGFromL1Muon.cc.

References objects.IsoTrackAnalyzer::candidates, ALCARECOTkAlBeamHalo_cff::filter, edm::Event::getByToken(), hfClusterShapes_cfi::hits, mps_fire::i, SeedFromProtoTrack::isValid(), LogDebug, eostools::move(), edm::Event::put(), L1MuonSeedsMerger::resolve(), mps_fire::result, OrderedHitsGenerator::run(), SurveyInfoScenario_cff::seed, OrderedSeedingHits::size(), SeedingHitSet::size(), theFilterToken, theFitter, theHitGenerator, theMerger, theRegionProducer, theSourceToken, HiIsolationCommonParameters_cff::track, l1t::tracks, and SeedFromProtoTrack::trajectorySeed().

 {
   auto result = std::make_unique<L3MuonTrajectorySeedCollection>();
 
   edm::Handle<L1MuonParticleCollection> l1muon;
   ev.getByToken(theSourceToken, l1muon);
 
   edm::Handle<PixelTrackFilter> hfilter;
   ev.getByToken(theFilterToken, hfilter);
   const PixelTrackFilter& filter = *hfilter;
 
   LogDebug("TSGFromL1Muon")<<l1muon->size()<<" l1 muons to seed from.";
 
   L1MuonParticleCollection::const_iterator muItr = l1muon->begin(); 
   L1MuonParticleCollection::const_iterator muEnd = l1muon->end(); 
   for  ( size_t iL1 = 0;  muItr < muEnd; ++muItr, ++iL1) {
        
     if (muItr->gmtMuonCand().empty()) continue;
 
     const L1MuGMTCand & muon = muItr->gmtMuonCand();
     l1extra::L1MuonParticleRef l1Ref(l1muon, iL1);
 
     theRegionProducer->setL1Constraint(muon);
     theFitter->setL1Constraint(muon);
 
     typedef std::vector<std::unique_ptr<TrackingRegion> > Regions;
     Regions regions = theRegionProducer->regions();
     for (Regions::const_iterator ir=regions.begin(); ir != regions.end(); ++ir) {
 
       L1MuonSeedsMerger::TracksAndHits tracks;
       const TrackingRegion & region = **ir;
       const OrderedSeedingHits & candidates = theHitGenerator->run(region,ev,es);
 
       unsigned int nSets = candidates.size();
       for (unsigned int ic= 0; ic <nSets; ic++) {
 
         const SeedingHitSet & hits =  candidates[ic]; 
         std::vector<const TrackingRecHit *> trh;
         for (unsigned int i= 0, nHits = hits.size(); i< nHits; ++i) trh.push_back( hits[i]->hit() );
 
         theFitter->setPxConstraint(hits);
         reco::Track* track = theFitter->run(es, trh, region);
         if (!track) continue;
 
         if (!filter(track, trh) ) { delete track; continue; }
         tracks.push_back(L1MuonSeedsMerger::TrackAndHits(track, hits));
       }
   
       if(theMerger) theMerger->resolve(tracks);
       for (L1MuonSeedsMerger::TracksAndHits::const_iterator it = tracks.begin();
         it != tracks.end(); ++it) {
 
         SeedFromProtoTrack seed( *(it->first), it->second, es);
         if (seed.isValid()) (*result).push_back(L3MuonTrajectorySeed(seed.trajectorySeed(),l1Ref));
 
 //      GlobalError vtxerr( sqr(region->originRBound()), 0, sqr(region->originRBound()),
 //                                               0, 0, sqr(region->originZBound()));
 //      SeedFromConsecutiveHits seed( candidates[ic],region->origin(), vtxerr, es);
 //      if (seed.isValid()) (*result).push_back( seed.TrajSeed() );
         delete it->first;
       }
     }
   }
 
   LogDebug("TSGFromL1Muon")<<result->size()<<" seeds to the event.";
   ev.put(std::move(result));
 }