#include <MuonAnalysis/MuonAssociators/plugins/TriggerMatcherToHLTDebug.cc>

Inheritance diagram for TriggerMatcherToHLTDebug:

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

void	produce (edm::Event &event, const edm::EventSetup &eventSetup) override

	TriggerMatcherToHLTDebug (const edm::ParameterSet &pset)

virtual	~TriggerMatcherToHLTDebug ()
	Destructor. More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef edm::AssociationMap < edm::OneToMany< std::vector < L2MuonTrajectorySeed > , std::vector < L2MuonTrajectorySeed > > >	SeedMap

Private Member Functions
template<typename T >
void	storeValueMap (edm::Event &iEvent, const edm::Handle< edm::View< reco::Muon > > &handle, const std::vector< T > &values, const std::string &label) const
	Store extra information in a ValueMap. More...

Private Attributes
edm::InputTag	beamspotTag_

double	deltaR_

PropagateToMuon	l1matcher_

edm::InputTag	l1Tag_

double	max_Dr_L2

double	max_Dr_L3

double	max_Dz_L2

double	max_Dz_L3

double	max_Eta_L2

double	max_Eta_L3

std::string	metname

int	min_N_L2

int	min_N_L3

int	min_Nhits_L2

int	min_Nhits_L3

double	min_Pt_L2

double	min_Pt_L3

int	minL1Quality_

double	nsigma_Pt_L2

double	nsigma_Pt_L3

edm::InputTag	seedMapTag_

edm::InputTag	tagTag_

edm::InputTag	theL2MuonsLabel

edm::InputTag	theL2SeedsLabel

edm::InputTag	theL3MuonsLabel

edm::InputTag	theL3SeedsLabel

edm::InputTag	theL3TkTracksLabel

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: Matches RECO muons to Trigger ones using HLTDEBUG information. Muon is first matched to L1 using the PropagateToMuon tool from this same package, then all compatible L1s are examined and the corresponding L2 and L3 objects are searched using the references inside those objects.

Definition at line 74 of file TriggerMatcherToHLTDebug.cc.

Member Typedef Documentation

typedef edm::AssociationMap<edm::OneToMany<std::vector<L2MuonTrajectorySeed>, std::vector<L2MuonTrajectorySeed> > > TriggerMatcherToHLTDebug::SeedMap

private

Definition at line 126 of file TriggerMatcherToHLTDebug.cc.

Constructor & Destructor Documentation

TriggerMatcherToHLTDebug::TriggerMatcherToHLTDebug ( const edm::ParameterSet & pset )

explicit

Definition at line 144 of file TriggerMatcherToHLTDebug.cc.

References edm::ParameterSet::getParameter(), metname, theL2MuonsLabel, theL2SeedsLabel, theL3MuonsLabel, theL3SeedsLabel, and theL3TkTracksLabel.

                                                                              :
   tagTag_(pset.getParameter<edm::InputTag>("tags")),
   l1Tag_(pset.getParameter<edm::InputTag>("l1s")),
   l1matcher_(pset.getParameter<edm::ParameterSet>("l1matcherConfig")),
   deltaR_(pset.getParameter<double>("deltaR")),
   minL1Quality_(pset.getParameter<int32_t>("MinL1Quality")),
   beamspotTag_(pset.getParameter<edm::InputTag>("BeamSpotTag")),
   min_N_L2(pset.getParameter<int> ("MinN_L2")),
   max_Eta_L2(pset.getParameter<double> ("MaxEta_L2")),
   min_Nhits_L2(pset.getParameter<int> ("MinNhits_L2")),
   max_Dr_L2(pset.getParameter<double> ("MaxDr_L2")),
   max_Dz_L2(pset.getParameter<double> ("MaxDz_L2")),
   min_Pt_L2(pset.getParameter<double> ("MinPt_L2")),
   nsigma_Pt_L2(pset.getParameter<double> ("NSigmaPt_L2")),
   min_N_L3(pset.getParameter<int> ("MinN_L3")),
   max_Eta_L3(pset.getParameter<double> ("MaxEta_L3")),
   min_Nhits_L3(pset.getParameter<int> ("MinNhits_L3")),
   max_Dr_L3(pset.getParameter<double> ("MaxDr_L3")),
   max_Dz_L3(pset.getParameter<double> ("MaxDz_L3")),
   min_Pt_L3(pset.getParameter<double> ("MinPt_L3")),
   nsigma_Pt_L3(pset.getParameter<double> ("NSigmaPt_L3")),
   seedMapTag_(pset.getParameter<edm::InputTag >("SeedMapTag"))
 {
 
 
   theL2SeedsLabel = pset.getParameter<InputTag>("L2Seeds_Collection");
   theL2MuonsLabel = pset.getParameter<InputTag>("L2Muons_Collection");
   theL3SeedsLabel = pset.getParameter<InputTag>("L3Seeds_Collection");
   theL3TkTracksLabel = pset.getParameter<InputTag>("L3TkTracks_Collection");
   theL3MuonsLabel = pset.getParameter<InputTag>("L3Muons_Collection");
 
   
   metname = "TriggerMatcherToHLTDebug";
 
   produces<edm::ValueMap<int> > ("propagatesToM2"); 
   produces<edm::ValueMap<int> > ("hasL1Particle"); 
   produces<edm::ValueMap<int> > ("hasL1Filtered"); 
   produces<edm::ValueMap<int> > ("hasL2Seed"); 
   produces<edm::ValueMap<int> > ("hasL2Muon"); 
   produces<edm::ValueMap<int> > ("hasL2MuonFiltered"); 
   produces<edm::ValueMap<int> > ("hasL3Seed"); 
   produces<edm::ValueMap<int> > ("hasL3Track"); 
   produces<edm::ValueMap<int> > ("hasL3Muon"); 
   produces<edm::ValueMap<int> > ("hasL3MuonFiltered"); 
   
   produces<edm::ValueMap<reco::CandidatePtr> > ("l1Candidate"); 
   produces<edm::ValueMap<reco::CandidatePtr> > ("l2Candidate"); 
   produces<edm::ValueMap<reco::CandidatePtr> > ("l3Candidate"); 
 } 

TriggerMatcherToHLTDebug::~TriggerMatcherToHLTDebug ( )

virtual

Destructor.

Definition at line 197 of file TriggerMatcherToHLTDebug.cc.

197 {}

Member Function Documentation

void TriggerMatcherToHLTDebug::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 408 of file TriggerMatcherToHLTDebug.cc.

References PropagateToMuon::init(), and l1matcher_.

                                                                                   {
   l1matcher_.init(iSetup);
 }

void TriggerMatcherToHLTDebug::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 200 of file TriggerMatcherToHLTDebug.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, beamspotTag_, deltaR(), deltaR_, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::TrackBase::dxy(), reco::TrackBase::dz(), reco::TrackBase::error(), PV3DBase< T, PVType, FrameType >::eta(), reco::TrackBase::eta(), event(), PropagateToMuon::extrapolate(), TrajectoryStateOnSurface::globalPosition(), i, edm::HandleBase::id(), edm::Ref< C, T, F >::id(), edm::HandleBase::isValid(), TrajectoryStateOnSurface::isValid(), edm::Ptr< T >::key(), edm::Ref< C, T, F >::key(), l1matcher_, l1Tag_, L2Muons_cfi::L2Muons, L3Muons_cfi::L3Muons, LogTrace, max_Dr_L2, max_Dr_L3, max_Dz_L2, max_Dz_L3, max_Eta_L2, max_Eta_L3, metname, min_Nhits_L2, min_Nhits_L3, min_Pt_L2, min_Pt_L3, RPCpg::mu, patZpeak::muons, nsigma_Pt_L2, nsigma_Pt_L3, reco::TrackBase::numberOfValidHits(), reco::TrackBase::parameter(), PV3DBase< T, PVType, FrameType >::phi(), reco::BeamSpot::position(), reco::TrackBase::pt(), L1MuGMTCand::quality(), seedMapTag_, edm::RefVector< C, T, F >::size(), tagTag_, theL2MuonsLabel, theL2SeedsLabel, theL3MuonsLabel, theL3SeedsLabel, and theL3TkTracksLabel.

                                                                                  {
   
   Handle<View<reco::Muon> > muons;
   event.getByLabel(tagTag_,muons);
   
   Handle<l1extra::L1MuonParticleCollection> L1Muons;
   event.getByLabel(l1Tag_,L1Muons);
   
   Handle<L2MuonTrajectorySeedCollection> L2Seeds;
   event.getByLabel(theL2SeedsLabel,L2Seeds);
   
   Handle<RecoChargedCandidateCollection> L2Muons;
   event.getByLabel(theL2MuonsLabel,L2Muons);
   
   Handle<L3MuonTrajectorySeedCollection> L3Seeds;
   event.getByLabel(theL3SeedsLabel,L3Seeds);
   
   Handle<reco::TrackCollection> L3TkTracks;
   event.getByLabel(theL3TkTracksLabel,L3TkTracks);
   
   Handle<RecoChargedCandidateCollection> L3Muons;
   event.getByLabel(theL3MuonsLabel,L3Muons);
   
   //beam spot
   BeamSpot beamSpot;
   Handle<BeamSpot> recoBeamSpotHandle;
   event.getByLabel(beamspotTag_,recoBeamSpotHandle);
   beamSpot = *recoBeamSpotHandle;
   
   //new for the MAP!!!!
   edm::Handle<SeedMap> seedMapHandle;
   event.getByLabel(seedMapTag_, seedMapHandle);
   
   
   size_t nmu = muons->size();
   std::vector<int> propagatesToM2(nmu), hasL1Particle(nmu), hasL1Filtered(nmu);
   std::vector<int> hasL2Seed(nmu), hasL2Muon(nmu), hasL2MuonFiltered(nmu);
   std::vector<int> hasL3Seed(nmu), hasL3Track(nmu), hasL3TrackFiltered(nmu), hasL3Muon(nmu), hasL3MuonFiltered(nmu);
   std::vector<reco::CandidatePtr> l1ptr(nmu), l2ptr(nmu), l3ptr(nmu);
   
   for (size_t i = 0; i < nmu; ++i) {
     const reco::Muon &mu = (*muons)[i];
     
     // Propagate to muon station (using the L1 tool)
     TrajectoryStateOnSurface stateAtMB2 = l1matcher_.extrapolate(mu);
     if (!stateAtMB2.isValid())  continue;
     propagatesToM2[i] = 1;
     
     double etaTk = stateAtMB2.globalPosition().eta();
     double phiTk = stateAtMB2.globalPosition().phi();
     l1extra::L1MuonParticleCollection::const_iterator it;
     vector<l1extra::L1MuonParticleRef>::const_iterator itMu3;
     L2MuonTrajectorySeedCollection::const_iterator iSeed;
     L3MuonTrajectorySeedCollection::const_iterator iSeedL3;
     RecoChargedCandidateCollection::const_iterator iL2Muon;
     reco::TrackCollection::const_iterator tktrackL3;
     RecoChargedCandidateCollection::const_iterator iL3Muon;
     
     reco::CandidatePtr thisL1, thisL2, thisL3;
     for(it = L1Muons->begin(); it != L1Muons->end(); ++it) {
       
       const L1MuGMTExtendedCand muonCand = (*it).gmtMuonCand();
       unsigned int quality =  muonCand.quality();
       
       double L1phi =(*it).phi();
       double L1eta =(*it).eta();
       double L1pt =(*it).pt();
       double dR=deltaR(etaTk,phiTk,L1eta,L1phi);
       
       //CONDIZIONE-> CE NE E' UNA ASSOCIATA?
       if (dR >= deltaR_) continue;
       thisL1 = reco::CandidatePtr(L1Muons, it - L1Muons->begin());
       if (!hasL1Particle[i]) l1ptr[i] = thisL1; // if nobody reached L1 before, then we're the best L1 found up to now.
       hasL1Particle[i]++;
       
       if ((quality <= 3) || (L1pt<7)) continue;
       if (!hasL1Filtered[i]) l1ptr[i] = thisL1; // if nobody reached L1 before, then we're the best L1 found up to now.
       hasL1Filtered[i]++;
       
       if(!L2Seeds.isValid()) continue;
       //LOOP SULLA COLLEZIONE DEI SEED
       for( iSeed = L2Seeds->begin(); iSeed != L2Seeds->end(); ++iSeed) {
     
     l1extra::L1MuonParticleRef l1FromSeed = iSeed->l1Particle();
     if (l1FromSeed.id() != L1Muons.id()) throw cms::Exception("CorruptData") << "You're using a different L1 collection than the one used by L2 seeds.\n";
     if (l1FromSeed.key() != thisL1.key()) continue;
     if (!hasL2Seed[i]) l1ptr[i] = thisL1; // if nobody reached here before, we're the best L1
     hasL2Seed[i]++;
     
     if(!L2Muons.isValid()) continue;
     //LOOP SULLA COLLEZIONE L2MUON
     for( iL2Muon = L2Muons->begin(); iL2Muon != L2Muons->end(); ++iL2Muon) {
       
       
       //MI FACCIO DARE REF E GUARDO SE E' UGUALE AL L2SEED ASSOCIATO
       //BEFORE THE ASSOCIATION MAP!!!!! 
       //edm::Ref<L2MuonTrajectorySeedCollection> l2seedRef = iL2Muon->track()->seedRef().castTo<edm::Ref<L2MuonTrajectorySeedCollection> >();
       //l1extra::L1MuonParticleRef l1FromL2 = l2seedRef->l1Particle();
       
       //if (l1FromL2.id() != l1FromSeed.id()) throw cms::Exception("CorruptData") << "You're using L2s with a different L1 collection than the one used by L2 seeds.\n";
       //if (l1FromL2 != l1FromSeed) continue;
       
       //AFTER THE ASSOCIATION MAP
       const edm::RefVector<L2MuonTrajectorySeedCollection>& seeds = (*seedMapHandle)[iL2Muon->track()->seedRef().castTo<edm::Ref<L2MuonTrajectorySeedCollection> >()];
       //      bool isTriggered = false;
       for(size_t jjj=0; jjj<seeds.size(); jjj++){
         
         if(seeds[jjj]->l1Particle()!= l1FromSeed) continue;
         
       }
           
       
       thisL2 = reco::CandidatePtr(L2Muons, iL2Muon - L2Muons->begin()) ;
       if (!hasL2Muon[i]) { l1ptr[i] = thisL1; l2ptr[i] = thisL2; } // if nobody reached here before, we're the best L1 and L2)
       hasL2Muon[i]++;
       
       LogTrace(metname) <<"L2MUON TROVATO!"<<endl;
       const reco::Track & L2Track = *iL2Muon->track();
       double Eta_L2= L2Track.eta();
       double Pt_L2= L2Track.pt(); 
       int nValidHits_L2= L2Track.numberOfValidHits();
       double BSPos_L2 = L2Track.dxy(beamSpot.position());
       double dz_L2 =L2Track.dz();
       double err0_L2 = L2Track.error(0);
       double abspar0_L2 = fabs(L2Track.parameter(0));
       double ptLx_L2 = Pt_L2;
       if (abspar0_L2>0) ptLx_L2 += nsigma_Pt_L2*err0_L2/abspar0_L2*Pt_L2;
       
       //GUARDO SE L2MUON ASSOCIATO AVREBBE PASSATO IL FILTRO
       bool passFilter = (((fabs(Eta_L2))<=max_Eta_L2)&&(nValidHits_L2>=min_Nhits_L2)&&((fabs(BSPos_L2))<=max_Dr_L2)&&((fabs(dz_L2))<=max_Dz_L2)&&(ptLx_L2>=min_Pt_L2));
       if (!passFilter) continue;
       if (!hasL2MuonFiltered[i]) { l1ptr[i] = thisL1; l2ptr[i] = thisL2; } // if nobody reached here before, we're the best L1 and L2)
       hasL2MuonFiltered[i]++;
       
       const reco::TrackRef L2FilteredRef = iL2Muon->track();
       
       //########L3 PART##############
       if (!L3Seeds.isValid()) continue;
       for (iSeedL3 = L3Seeds->begin(); iSeedL3!= L3Seeds->end(); ++iSeedL3){
         
         TrackRef staTrack = iSeedL3->l2Track();
         if (staTrack!=L2FilteredRef) continue;
         if (!hasL3Seed[i]) { l1ptr[i] = thisL1; l2ptr[i] = thisL2; } // if nobody reached here before, we're the best L1 and L2)
         hasL3Seed[i]++;
         
         if (!L3TkTracks.isValid()) continue;
         for (tktrackL3 = L3TkTracks->begin(); tktrackL3!= L3TkTracks->end(); ++tktrackL3){
           
           edm::Ref<L3MuonTrajectorySeedCollection> l3seedRef = tktrackL3->seedRef().castTo<edm::Ref<L3MuonTrajectorySeedCollection> >();
           TrackRef staTrack2 = l3seedRef->l2Track();
           
           if (staTrack2!=L2FilteredRef) continue;
           if (!hasL3Track[i]) { l1ptr[i] = thisL1; l2ptr[i] = thisL2; } // if nobody reached here before, we're the best L1 and L2)
           hasL3Track[i]++;
           
           if (!L3Muons.isValid()) continue;
           for (iL3Muon = L3Muons->begin(); iL3Muon != L3Muons->end(); ++iL3Muon) {
         
         edm::Ref<L3MuonTrajectorySeedCollection> l3seedRef2 = iL3Muon->track()->seedRef().castTo<edm::Ref<L3MuonTrajectorySeedCollection> >();
         TrackRef staTrack3 = l3seedRef2->l2Track();
         
         if (staTrack3!=L2FilteredRef) continue;
         thisL3 = reco::CandidatePtr(L3Muons, iL3Muon - L3Muons->begin());
         
         if (!hasL3Muon[i]) { l1ptr[i] = thisL1; l2ptr[i] = thisL2; l3ptr[i] = thisL3; } // if nobody reached here before, we're the best L1, L2, L3
         hasL3Muon[i]++;
         
         const reco::Track &L3Track = *iL3Muon->track();
         double Eta_L3= L3Track.eta();
         double Pt_L3= L3Track.pt(); 
         int nValidHits_L3= L3Track.numberOfValidHits();
         double BSPos_L3 = L3Track.dxy(beamSpot.position());
         double dz_L3 =L3Track.dz();
         double err0_L3 = L3Track.error(0);
         double abspar0_L3 = fabs(L3Track.parameter(0));
         double ptLx_L3 = Pt_L3;
         
         if (abspar0_L3>0) ptLx_L3 += nsigma_Pt_L3*err0_L3/abspar0_L3*Pt_L3;
 
         if(((fabs(Eta_L3))<=max_Eta_L3)&&(nValidHits_L3>=min_Nhits_L3)&&((fabs(BSPos_L3))<=max_Dr_L3)&&((fabs(dz_L3))<=max_Dz_L3)&&(ptLx_L3>=min_Pt_L3)){
           
           if (!hasL3MuonFiltered[i]) { l1ptr[i] = thisL1; l2ptr[i] = thisL2; l3ptr[i] = thisL3; } // if nobody reached here before, we're the best L1, L2, L3
           hasL3MuonFiltered[i]++;
           
         }//L3MUON FILTERED ASSOCIATO TROVATO
           }//L3MUON LOOP
         }// L3 TRACKS
       }// L3 SEEDS
     }//T L2 MUONS
       }// L2 SEEDS
     }//L1 MUONS
   } // RECO MUONS
   storeValueMap<int>(event, muons, propagatesToM2,    "propagatesToM2");
   storeValueMap<int>(event, muons, hasL1Particle,     "hasL1Particle");
   storeValueMap<int>(event, muons, hasL1Filtered,     "hasL1Filtered");
   storeValueMap<int>(event, muons, hasL2Seed,         "hasL2Seed");
   storeValueMap<int>(event, muons, hasL2Muon,         "hasL2Muon");
   storeValueMap<int>(event, muons, hasL2MuonFiltered, "hasL2MuonFiltered");
   storeValueMap<int>(event, muons, hasL3Seed,         "hasL3Seed");
   storeValueMap<int>(event, muons, hasL3Track,        "hasL3Track");
   storeValueMap<int>(event, muons, hasL3Muon,         "hasL3Muon");
   storeValueMap<int>(event, muons, hasL3MuonFiltered, "hasL3MuonFiltered");
   storeValueMap<reco::CandidatePtr>(event, muons, l1ptr, "l1Candidate");
   storeValueMap<reco::CandidatePtr>(event, muons, l2ptr, "l2Candidate");
   storeValueMap<reco::CandidatePtr>(event, muons, l3ptr, "l3Candidate");
 } // METHOD

template<typename T >

void TriggerMatcherToHLTDebug::storeValueMap	(	edm::Event &	iEvent,
		const edm::Handle< edm::View< reco::Muon > > &	handle,
		const std::vector< T > &	values,
		const std::string &	label
	)		const

private

Store extra information in a ValueMap.

Definition at line 414 of file TriggerMatcherToHLTDebug.cc.

References edm::helper::Filler< Map >::fill(), patZpeak::handle, edm::helper::Filler< Map >::insert(), and edm::Event::put().

                                                       {
   using namespace edm; using namespace std;
   auto_ptr<ValueMap<T> > valMap(new ValueMap<T>());
   typename edm::ValueMap<T>::Filler filler(*valMap);
   filler.insert(handle, values.begin(), values.end());
   filler.fill();
   iEvent.put(valMap, label);
 }