MuTriggerAnalyzer Class Reference

Inheritance diagram for MuTriggerAnalyzer:

Public Member Functions
	MuTriggerAnalyzer (const edm::ParameterSet &pset)
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup)
virtual void	endJob ()
bool	IsMuMatchedToHLTMu (const reco::Muon &, std::vector< reco::Particle >, double, double)

Private Attributes
TH1D *	deltaPtOverPt_
TH1D *	deltaR_
double	etaMuCut_
edm::Handle < trigger::TriggerEvent >	handleTriggerEvent_
std::string	hltPath_
TH1D *	hTrigMuonPtDenS_
TH1D *	hTrigMuonPtNumS_
std::string	L3FilterName_
double	maxDeltaR_
double	maxDPtRel_
edm::InputTag	muonTag_
const int	nbins_
const double	ptMax_
double	ptMuCut_
edm::InputTag	trigEv_
edm::Handle< edm::TriggerResults >	triggerResults_
edm::TriggerNames const *	trigNames_
edm::InputTag	trigTag_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 49 of file MuTriggerAnalyzer.cc.

Constructor & Destructor Documentation

MuTriggerAnalyzer::MuTriggerAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 112 of file MuTriggerAnalyzer.cc.

References deltaPtOverPt_, deltaR_, hTrigMuonPtDenS_, hTrigMuonPtNumS_, TFileDirectory::make(), maxDeltaR_, maxDPtRel_, nbins_, and ptMax_.

                                                                :
  trigTag_(cfg.getParameter<edm::InputTag> ("TrigTag")),
  trigEv_(cfg.getParameter<edm::InputTag> ("triggerEvent")),
  muonTag_(cfg.getUntrackedParameter<edm::InputTag>("muons")),
  ptMuCut_(cfg.getUntrackedParameter<double>("ptMuCut")),
  etaMuCut_(cfg.getUntrackedParameter<double>("etaMuCut")),
  hltPath_(cfg.getParameter<std::string >("hltPath")),
  L3FilterName_(cfg.getParameter<std::string >("L3FilterName")),
  maxDPtRel_(cfg.getParameter<double>("maxDPtRel")),
  maxDeltaR_(cfg.getParameter<double>("maxDeltaR")),
  nbins_(cfg.getParameter<double>("ptMax_")), 
  ptMax_(cfg.getParameter<double>("ptMax_")){
   Service<TFileService> fs;
 
 hTrigMuonPtNumS_ = fs->make<TH1D>("hTrigMuonPtNumS", "hTrigMuonPtNumS", nbins_ + 1, 0, ptMax_); 
  hTrigMuonPtDenS_ = fs->make<TH1D>("hTrigMuonPtDenS", "hTrigMuonPtDenS", nbins_ +1 , 0, ptMax_); 
  deltaR_ = fs->make<TH1D>("deltaR", "deltaR", nbins_+1, 0, maxDeltaR_); 
  deltaPtOverPt_ = fs->make<TH1D>("deltaPtOverPt", "deltaPtOverPt", nbins_ + 1, 0, maxDPtRel_); 

}

Member Function Documentation

void MuTriggerAnalyzer::analyze ( const edm::Event &  event, const edm::EventSetup &  setup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 142 of file MuTriggerAnalyzer.cc.

References gather_cfg::cout, edm::InputTag::encode(), reco::LeafCandidate::eta(), eta(), etaMuCut_, edm::Event::getByLabel(), hltPath_, i, reco::Muon::isGlobalMuon(), IsMuMatchedToHLTMu(), gen::k, L3FilterName_, maxDeltaR_, maxDPtRel_, patZpeak::muons, muonTag_, mergeVDriftHistosByStation::name, L1TEmulatorMonitor_cff::p, reco::LeafCandidate::p4(), reco::LeafCandidate::pt(), ptMuCut_, trigEv_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), trigNames, trigNames_, and trigTag_.

                                                                     {
      bool singleTrigFlag1 = false;
          //      bool overlap = false;
  
      // Trigger
      Handle<TriggerResults> triggerResults;
      TriggerNames trigNames;
      if (!ev.getByLabel(trigTag_, triggerResults)) {
    LogWarning("") << ">>> TRIGGER collection does not exist !!!";
    return;
      }
      ev.getByLabel(trigTag_, triggerResults); 
      // trigNames.init(*triggerResults);
      trigNames_ = &ev.triggerNames(*triggerResults);
      bool trigger_fired = false;
      for (unsigned int i=0; i<triggerResults->size(); i++) {
    

        std::string trigName = trigNames_->triggerName(i);
    if ( trigName == hltPath_ && triggerResults->accept(i)) {
      trigger_fired = true;
    }
      }   
      edm::Handle< trigger::TriggerEvent > handleTriggerEvent;
      //  LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << hltPath_ << ")";
          if ( ! ev.getByLabel( trigEv_, handleTriggerEvent ))  {
    LogWarning( "errorTriggerEventValid" ) << "trigger::TriggerEvent product with InputTag " << trigEv_.encode() << " not in event";
    return;
      }
      ev.getByLabel( trigEv_, handleTriggerEvent );
      const trigger::TriggerObjectCollection & toc(handleTriggerEvent->getObjects());
      size_t nMuHLT =0;
      std::vector<reco::Particle>  HLTMuMatched; 
      for ( size_t ia = 0; ia < handleTriggerEvent->sizeFilters(); ++ ia) {
    std::string fullname = handleTriggerEvent->filterTag(ia).encode();
    std::string name;
    size_t p = fullname.find_first_of(':');
    if ( p != std::string::npos) {
      name = fullname.substr(0, p);
    }
    else {
      name = fullname;
    }
    if ( &toc !=0 ) {
      const trigger::Keys & k = handleTriggerEvent->filterKeys(ia);
      for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
        if (name == L3FilterName_  ) { 
          HLTMuMatched.push_back(toc[*ki].particle());
          nMuHLT++;     
        }
      }    
    }
      }
      
    //  looping on muon....
      Handle<View<Muon> >   muons;     
      if (!ev.getByLabel(muonTag_, muons)) {           
    LogError("") << ">>> muon collection does not exist !!!";     
    return;     
      }
      
      ev.getByLabel(muonTag_, muons);  
      //saving only muons with pt> ptMuCut and eta<etaMuCut  
      std::vector<reco::Muon>  highPtGlbMuons; 
    
      for (unsigned int i=0; i<muons->size(); i++ ){    
        const reco::Muon & mu = muons->at(i);
    double pt = mu.pt();
    double eta = mu.eta();
    if (pt> ptMuCut_ && fabs(eta)< etaMuCut_) {
      if (mu.isGlobalMuon()) highPtGlbMuons.push_back(mu);    
    }
      }
      unsigned int nHighPtGlbMu = highPtGlbMuons.size();
      std::cout << "I've got " << nHighPtGlbMu << " nHighPtGlbMu" << std::endl;
      // unsigned int nHighPtStaMu = highPtStaMuons.size();

    // stop the loop after 10 cicles....  
        (nHighPtGlbMu> 10)?   nHighPtGlbMu=10 : 1; 

    if (nHighPtGlbMu>0 ){
      
           for(unsigned int i =0 ; i < nHighPtGlbMu ; i++) {
        reco::Muon muon1 = highPtGlbMuons[i];
        math::XYZTLorentzVector mu1(muon1.p4());
        //      double pt1= muon1.pt();
        
        singleTrigFlag1 = IsMuMatchedToHLTMu ( muon1,  HLTMuMatched ,maxDeltaR_, maxDPtRel_ );
        
      }
      
    }      
      
}

void MuTriggerAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 133 of file MuTriggerAnalyzer.cc.

References gather_cfg::cout, hTrigMuonPtDenS_, hTrigMuonPtNumS_, i, and nbins_.

                               {
  for (int i = 0; i < nbins_+1; ++i){
    std::cout << "number of reco muon in bin " << i << " = " << hTrigMuonPtDenS_->GetBinContent(i) << std::endl;
    std::cout << "number of hlt muon in bin " << i << " = " << hTrigMuonPtNumS_->GetBinContent(i) << std::endl;

 
  }
}

bool MuTriggerAnalyzer::IsMuMatchedToHLTMu ( const reco::Muon &  mu, std::vector< reco::Particle >  HLTMu, double  DR, double  DPtRel  )

private

Definition at line 79 of file MuTriggerAnalyzer.cc.

References gather_cfg::cout, deltaR(), deltaR_, reco::LeafCandidate::eta(), HcalObjRepresent::Fill(), gen::k, and reco::LeafCandidate::pt().

Referenced by analyze().

                                                                                                                           {
  size_t dim =  HLTMu.size();
  size_t nPass=0;
  
  // filling the denumerator;
  double muRecoPt= mu.pt();
  hTrigMuonPtDenS_-> Fill(muRecoPt);         
 
  if (dim==0) return false;
  for (size_t k =0; k< dim; k++ ) {
    if (  (deltaR(HLTMu[k], mu) < DR)   && (fabs(HLTMu[k].pt() - mu.pt())/ HLTMu[k].pt()<DPtRel)){     
      nPass++ ;
      std::cout<< " matching a muon " << std::endl;  
      std::cout  << "muon reco  pt : " << muRecoPt<< std::endl;     
      std::cout  << "muon reco  eta " <<  mu.eta() << std::endl;   
      std::cout << "muon trigger  pt "<< HLTMu[k].pt() << std::endl; 
      // filling the numerator, at the same bin as the denum..... 
      hTrigMuonPtNumS_-> Fill(muRecoPt);
      deltaR_-> Fill(deltaR(HLTMu[k], mu));
      deltaPtOverPt_-> Fill(fabs(HLTMu[k].pt() - mu.pt())/ HLTMu[k].pt());

         
      std::cout << "muon trigger  eta : "<< HLTMu[k].eta() << std::endl;   
      std::cout  <<"deltaR((HLTMu[k], mu)): "<< deltaR(HLTMu[k], mu) << std::endl;
      std::cout  <<"deltaPtOverPt: "<< fabs(HLTMu[k].pt() - mu.pt())/ HLTMu[k].pt() << std::endl;
    }
  }
  
  return (nPass>0);
}

Member Data Documentation

TH1D* MuTriggerAnalyzer::deltaPtOverPt_

private

Definition at line 74 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

TH1D* MuTriggerAnalyzer::deltaR_

private

Definition at line 73 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

double MuTriggerAnalyzer::etaMuCut_

private

Definition at line 61 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::Handle< trigger::TriggerEvent > MuTriggerAnalyzer::handleTriggerEvent_

private

Definition at line 66 of file MuTriggerAnalyzer.cc.

std::string MuTriggerAnalyzer::hltPath_

private

Definition at line 62 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

TH1D* MuTriggerAnalyzer::hTrigMuonPtDenS_

private

Definition at line 72 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

TH1D* MuTriggerAnalyzer::hTrigMuonPtNumS_

private

Definition at line 71 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

std::string MuTriggerAnalyzer::L3FilterName_

private

Definition at line 63 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

double MuTriggerAnalyzer::maxDeltaR_

private

Definition at line 68 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

double MuTriggerAnalyzer::maxDPtRel_

private

Definition at line 68 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

edm::InputTag MuTriggerAnalyzer::muonTag_

private

Definition at line 60 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

const int MuTriggerAnalyzer::nbins_

private

Definition at line 69 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

const double MuTriggerAnalyzer::ptMax_

private

Definition at line 70 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

double MuTriggerAnalyzer::ptMuCut_

private

Definition at line 61 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::InputTag MuTriggerAnalyzer::trigEv_

private

Definition at line 59 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::Handle<edm::TriggerResults> MuTriggerAnalyzer::triggerResults_

private

Definition at line 64 of file MuTriggerAnalyzer.cc.

edm::TriggerNames const* MuTriggerAnalyzer::trigNames_

private

Definition at line 65 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::InputTag MuTriggerAnalyzer::trigTag_

private

Definition at line 58 of file MuTriggerAnalyzer.cc.

Referenced by analyze().