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

Inheritance diagram for TriggerMatcherToHLTDebug:

Public Member Functions
void	beginRun (edm::Run &run, const edm::EventSetup &eventSetup)
void	produce (edm::Event &event, const edm::EventSetup &eventSetup)
	TriggerMatcherToHLTDebug (const edm::ParameterSet &pset)
virtual	~TriggerMatcherToHLTDebug ()
	Destructor.
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.
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

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.

{}

Member Function Documentation

void TriggerMatcherToHLTDebug::beginRun	(	edm::Run &	run,
		const edm::EventSetup &	eventSetup
	)		`[virtual]`

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
	)		`[virtual]`

Implements edm::EDProducer.

Definition at line 200 of file TriggerMatcherToHLTDebug.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, beamspotTag_, deltaR(), deltaR_, 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(), TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), edm::Ref< C, T, F >::key(), edm::Ptr< T >::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);
}