#include <HLTJetMETValidation.h>

Inheritance diagram for HLTJetMETValidation:

Public Member Functions
	HLTJetMETValidation (const edm::ParameterSet &)

	~HLTJetMETValidation ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginRun (const edm::Run &, const edm::EventSetup &)

virtual void	endJob ()

void	getHLTResults (const edm::TriggerResults &, const edm::TriggerNames &triggerNames)

Private Attributes
std::vector< MonitorElement * >	_meGenJetEta

std::vector< MonitorElement * >	_meGenJetEtaTrg

std::vector< MonitorElement * >	_meGenJetEtaTrgLow

std::vector< MonitorElement * >	_meGenJetEtaTrgMC

std::vector< MonitorElement * >	_meGenJetPhi

std::vector< MonitorElement * >	_meGenJetPhiTrg

std::vector< MonitorElement * >	_meGenJetPhiTrgLow

std::vector< MonitorElement * >	_meGenJetPhiTrgMC

std::vector< MonitorElement * >	_meGenJetPt

std::vector< MonitorElement * >	_meGenJetPtTrg

std::vector< MonitorElement * >	_meGenJetPtTrgLow

std::vector< MonitorElement * >	_meGenJetPtTrgMC

std::vector< MonitorElement * >	_meGenMET

std::vector< MonitorElement * >	_meGenMETTrg

std::vector< MonitorElement * >	_meGenMETTrgLow

std::vector< MonitorElement * >	_meGenMETTrgMC

std::vector< MonitorElement * >	_meRecoJetEta

std::vector< MonitorElement * >	_meRecoJetEtaTrg

std::vector< MonitorElement * >	_meRecoJetEtaTrgLow

std::vector< MonitorElement * >	_meRecoJetEtaTrgMC

std::vector< MonitorElement * >	_meRecoJetPhi

std::vector< MonitorElement * >	_meRecoJetPhiTrg

std::vector< MonitorElement * >	_meRecoJetPhiTrgLow

std::vector< MonitorElement * >	_meRecoJetPhiTrgMC

std::vector< MonitorElement * >	_meRecoJetPt

std::vector< MonitorElement * >	_meRecoJetPtTrg

std::vector< MonitorElement * >	_meRecoJetPtTrgLow

std::vector< MonitorElement * >	_meRecoJetPtTrgMC

std::vector< MonitorElement * >	_meRecoMET

std::vector< MonitorElement * >	_meRecoMETTrg

std::vector< MonitorElement * >	_meRecoMETTrgLow

std::vector< MonitorElement * >	_meRecoMETTrgMC

MonitorElement *	_triggerResults

edm::EDGetTokenT< PFJetCollection >	CaloJetAlgorithm

edm::EDGetTokenT < CaloMETCollection >	CaloMETColl

int	evtCnt

edm::EDGetTokenT < GenJetCollection >	GenJetAlgorithm

edm::EDGetTokenT < GenMETCollection >	GenMETColl

HLTConfigProvider	hltConfig_

bool	HLTinit_

edm::EDGetTokenT < edm::TriggerResults >	HLTriggerResults

std::vector< std::string >	hltTrgJet

std::vector< std::string >	hltTrgJetLow

std::vector< std::string >	hltTrgMet

std::vector< std::string >	hltTrgMetLow

std::map< std::string, bool >	hltTriggerMap

std::vector< bool >	hlttrigs

std::string	MyTrigger

std::string	outFile_

std::string	patternJetTrg_

std::string	patternMetTrg_

std::string	patternMuTrg_

DQMStore *	store

std::map< std::string, bool > ::iterator	trig_iter

edm::EDGetTokenT < TriggerEventWithRefs >	triggerEventObject_
	InputTag of TriggerEventWithRefs to analyze. More...

std::string	triggerTag_

bool	writeFile_

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

Definition at line 67 of file HLTJetMETValidation.h.

Constructor & Destructor Documentation

HLTJetMETValidation::HLTJetMETValidation ( const edm::ParameterSet & ps )

explicit

Definition at line 5 of file HLTJetMETValidation.cc.

References evtCnt, and store.

                                                                   : 
   triggerEventObject_(consumes<TriggerEventWithRefs>(ps.getUntrackedParameter<edm::InputTag>("triggerEventObject"))),
   CaloJetAlgorithm( consumes<PFJetCollection>(ps.getUntrackedParameter<edm::InputTag>( "CaloJetAlgorithm" ) )),
   GenJetAlgorithm( consumes<GenJetCollection>(ps.getUntrackedParameter<edm::InputTag>( "GenJetAlgorithm" ) )),
   CaloMETColl( consumes<CaloMETCollection>(ps.getUntrackedParameter<edm::InputTag>( "CaloMETCollection" ) )),
   GenMETColl( consumes<GenMETCollection>(ps.getUntrackedParameter<edm::InputTag>( "GenMETCollection" ) )),
   HLTriggerResults(consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>( "HLTriggerResults" ) )),
   triggerTag_(ps.getUntrackedParameter<std::string>("DQMFolder","SingleJet")),
   patternJetTrg_(ps.getUntrackedParameter<std::string>("PatternJetTrg","")),
   patternMetTrg_(ps.getUntrackedParameter<std::string>("PatternMetTrg","")),
   patternMuTrg_(ps.getUntrackedParameter<std::string>("PatternMuTrg","")),
   outFile_(ps.getUntrackedParameter<std::string>("OutputFileName","")),
   HLTinit_(false),
   //JL
   writeFile_(ps.getUntrackedParameter<bool>("WriteFile",false))
 {
   evtCnt=0;
 
   store = &*edm::Service<DQMStore>();
 
 }

HLTJetMETValidation::~HLTJetMETValidation ( )

Definition at line 27 of file HLTJetMETValidation.cc.

28 {

29 }

Member Function Documentation

void HLTJetMETValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 193 of file HLTJetMETValidation.cc.

References _meGenJetEta, _meGenJetEtaTrg, _meGenJetEtaTrgLow, _meGenJetEtaTrgMC, _meGenJetPhi, _meGenJetPhiTrg, _meGenJetPhiTrgLow, _meGenJetPhiTrgMC, _meGenJetPt, _meGenJetPtTrg, _meGenJetPtTrgLow, _meGenJetPtTrgMC, _meGenMET, _meGenMETTrg, _meGenMETTrgLow, _meGenMETTrgMC, _meRecoJetEta, _meRecoJetEtaTrg, _meRecoJetEtaTrgLow, _meRecoJetEtaTrgMC, _meRecoJetPhi, _meRecoJetPhiTrg, _meRecoJetPhiTrgLow, _meRecoJetPhiTrgMC, _meRecoJetPt, _meRecoJetPtTrg, _meRecoJetPtTrgLow, _meRecoJetPtTrgMC, _meRecoMET, _meRecoMETTrg, _meRecoMETTrgLow, _meRecoMETTrgMC, CaloJetAlgorithm, FamosSequences_cff::caloJets, CaloMETColl, evtCnt, HcalObjRepresent::Fill(), relval_steps::gen(), GenJetAlgorithm, HepMCValidationHelper::genMet(), GenMETColl, edm::Event::getByToken(), getHLTResults(), HLTriggerResults, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, hltTriggerMap, i, dt_dqm_sourceclient_common_cff::reco, trig_iter, triggerEventObject_, and edm::Event::triggerNames().

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace l1extra;
   using namespace trigger;
   
   evtCnt++;
   //get The triggerEvent
 
   Handle<TriggerEventWithRefs> trigEv;
   iEvent.getByToken(triggerEventObject_,trigEv);
 
 // get TriggerResults object
 
   bool gotHLT=true;
   //bool myTrig=false;
   //bool myTrigLow=false;
   std::vector<bool> myTrigJ;
   for (size_t it=0;it<hltTrgJet.size();it++) myTrigJ.push_back(false);
   std::vector<bool> myTrigJLow;
   for (size_t it=0;it<hltTrgJetLow.size();it++) myTrigJLow.push_back(false);
   std::vector<bool> myTrigM;
   for (size_t it=0;it<hltTrgMet.size();it++) myTrigM.push_back(false);
   std::vector<bool> myTrigMLow;
   for (size_t it=0;it<hltTrgMetLow.size();it++) myTrigMLow.push_back(false);
 
 
   Handle<TriggerResults> hltresults,hltresultsDummy;
   iEvent.getByToken(HLTriggerResults,hltresults);
   if (! hltresults.isValid() ) { 
     //std::cout << "  -- No HLTRESULTS"; 
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No HLTRESULTS";    
     gotHLT=false;
   }
 
   if (gotHLT) {
     const edm::TriggerNames & triggerNames = iEvent.triggerNames(*hltresults);
     getHLTResults(*hltresults, triggerNames);
     //    trig_iter=hltTriggerMap.find(MyTrigger);
     //trig_iter=hltTriggerMap.find(_HLTPath.label());
 
     //---pick-up the jet trigger decisions
     for (size_t it=0;it<hltTrgJet.size();it++) {
       trig_iter=hltTriggerMap.find(hltTrgJet[it]);
       if (trig_iter==hltTriggerMap.end()){
     //std::cout << "Could not find trigger path with name: " << _probefilter.label() << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "Could not find trigger path with name: " << _probefilter.label(); 
       }else{
     myTrigJ[it]=trig_iter->second;
       }
       //std::cout<<hltTrgJet[it].c_str()<<" "<<myTrigJ[it]<<std::endl;
     }
     for (size_t it=0;it<hltTrgJetLow.size();it++) {
       trig_iter=hltTriggerMap.find(hltTrgJetLow[it]);
       if (trig_iter==hltTriggerMap.end()){
     //std::cout << "Could not find trigger path with name: " << _probefilter.label() << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "Could not find trigger path with name: " << _probefilter.label(); 
       }else{
     myTrigJLow[it]=trig_iter->second;
       }
       //std::cout<<hltTrgJetLow[it].c_str()<<" "<<myTrigJLow[it]<<std::endl;
     }
     //---pick-up the met trigger decisions
     for (size_t it=0;it<hltTrgMet.size();it++) {
       trig_iter=hltTriggerMap.find(hltTrgMet[it]);
       if (trig_iter==hltTriggerMap.end()){
     //std::cout << "Could not find trigger path with name: " << _probefilter.label() << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "Could not find trigger path with name: " << _probefilter.label(); 
       }else{
     myTrigM[it]=trig_iter->second;
       }
       //std::cout<<hltTrgMet[it].c_str()<<" "<<myTrigM[it]<<std::endl;
     }
     for (size_t it=0;it<hltTrgMetLow.size();it++) {
       trig_iter=hltTriggerMap.find(hltTrgMetLow[it]);
       if (trig_iter==hltTriggerMap.end()){
     //std::cout << "Could not find trigger path with name: " << _probefilter.label() << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "Could not find trigger path with name: " << _probefilter.label(); 
       }else{
     myTrigMLow[it]=trig_iter->second;
       }
       //std::cout<<hltTrgMetLow[it].c_str()<<" "<<myTrigMLow[it]<<std::endl;
     }
   }
 
   Handle<PFJetCollection> caloJets,caloJetsDummy;
   iEvent.getByToken( CaloJetAlgorithm, caloJets );
   double calJetPt=-1.;
   double calJetEta=-999.;
   double calJetPhi=-999.;
   //double calHT=0;
   if (caloJets.isValid()) { 
     //Loop over the CaloJets and fill some histograms
     int jetInd = 0;
     for( PFJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
       //std::cout << "CALO JET #" << jetInd << std::endl << cal->print() << std::endl;
       if (jetInd == 0){
     calJetPt=cal->pt();
     calJetEta=cal->eta();
     calJetPhi=cal->phi();
     for (size_t it=0;it<hltTrgJet.size();it++) {
       _meRecoJetPt[it]->Fill( calJetPt );
       _meRecoJetEta[it]->Fill( calJetEta );
       _meRecoJetPhi[it]->Fill( calJetPhi );
       if (myTrigJ[it]) _meRecoJetPtTrgMC[it]->Fill( calJetPt );
       if (myTrigJ[it]) _meRecoJetEtaTrgMC[it]->Fill( calJetEta );
       if (myTrigJ[it]) _meRecoJetPhiTrgMC[it]->Fill( calJetPhi );
       if (myTrigJ[it] && myTrigJLow[it]) _meRecoJetPtTrg[it]->Fill( calJetPt );
       if (myTrigJ[it] && myTrigJLow[it]) _meRecoJetEtaTrg[it]->Fill( calJetEta );
       if (myTrigJ[it] && myTrigJLow[it]) _meRecoJetPhiTrg[it]->Fill( calJetPhi );
       if (myTrigJLow[it]) _meRecoJetPtTrgLow[it]->Fill( calJetPt );
       if (myTrigJLow[it]) _meRecoJetEtaTrgLow[it]->Fill( calJetEta );
       if (myTrigJLow[it]) _meRecoJetPhiTrgLow[it]->Fill( calJetPhi );
     }
     jetInd++;
       }
       /*
       if (cal->pt()>30) {
     calHT+=cal->pt();
       }
       */
     }
     /*
     _meRecoHT->Fill( calHT );
     if (myTrig) _meRecoHTTrg->Fill( calHT );
     if (myTrigLow) _meRecoHTTrgLow->Fill( calHT );
     */
   }else{
     //std::cout << "  -- No CaloJets" << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No CaloJets"; 
   }
 
   Handle<GenJetCollection> genJets,genJetsDummy;
   iEvent.getByToken( GenJetAlgorithm, genJets );
   double genJetPt=-1.;
   double genJetEta=-999.;
   double genJetPhi=-999.;
   //double genHT=0;
   if (genJets.isValid()) { 
     //Loop over the GenJets and fill some histograms
     int jetInd = 0;
     for( GenJetCollection::const_iterator gen = genJets->begin(); gen != genJets->end(); ++ gen ) {
       // std::cout << "CALO JET #" << jetInd << std::endl << cal->print() << std::endl;
       if (jetInd == 0){
     genJetPt=gen->pt();
     genJetEta=gen->eta();
     genJetPhi=gen->phi();
     for (size_t it=0;it<hltTrgJet.size();it++) {
       _meGenJetPt[it]->Fill( genJetPt );
       _meGenJetEta[it]->Fill( genJetEta );
       _meGenJetPhi[it]->Fill( genJetPhi );
       if (myTrigJ[it]) _meGenJetPtTrgMC[it]->Fill( genJetPt );
       if (myTrigJ[it]) _meGenJetEtaTrgMC[it]->Fill( genJetEta );
       if (myTrigJ[it]) _meGenJetPhiTrgMC[it]->Fill( genJetPhi );
       if (myTrigJ[it] && myTrigJLow[it]) _meGenJetPtTrg[it]->Fill( genJetPt );
       if (myTrigJ[it] && myTrigJLow[it]) _meGenJetEtaTrg[it]->Fill( genJetEta );
       if (myTrigJ[it] && myTrigJLow[it]) _meGenJetPhiTrg[it]->Fill( genJetPhi );
       if (myTrigJLow[it]) _meGenJetPtTrgLow[it]->Fill( genJetPt );
       if (myTrigJLow[it]) _meGenJetEtaTrgLow[it]->Fill( genJetEta );
       if (myTrigJLow[it]) _meGenJetPhiTrgLow[it]->Fill( genJetPhi );
     }
     jetInd++;
       }
       /*
       if (gen->pt()>30) {
     genHT+=gen->pt();
       }
       */
     }
     /*
     _meGenHT->Fill( genHT );
     if (myTrig) _meGenHTTrg->Fill( genHT );
     if (myTrigLow) _meGenHTTrgLow->Fill( genHT );
     */
   }else{
     //std::cout << "  -- No GenJets" << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No GenJets"; 
   }
   
 
   edm::Handle<CaloMETCollection> recmet, recmetDummy;
   iEvent.getByToken(CaloMETColl,recmet);
 
   double calMet=-1;
   if (recmet.isValid()) { 
     typedef CaloMETCollection::const_iterator cmiter;
     //std::cout << "Size of MET collection" <<  recmet.size() << std::endl;
     for ( cmiter i=recmet->begin(); i!=recmet->end(); i++) {
       calMet = i->pt();
       for (size_t it=0;it<hltTrgMet.size();it++) {
     _meRecoMET[it] -> Fill(calMet);
     if (myTrigM[it]) _meRecoMETTrgMC[it] -> Fill(calMet);
     if (myTrigM[it] && myTrigMLow[it]) _meRecoMETTrg[it] -> Fill(calMet);
     if (myTrigMLow[it]) _meRecoMETTrgLow[it] -> Fill(calMet);
       }
     }
   }else{
     //std::cout << "  -- No MET Collection with name: " << CaloMETColl << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No MET Collection with name: "<< CaloMETColl; 
   }
   
   edm::Handle<GenMETCollection> genmet, genmetDummy;
   iEvent.getByToken(GenMETColl,genmet);
 
   double genMet=-1;
   if (genmet.isValid()) { 
     typedef GenMETCollection::const_iterator cmiter;
     //std::cout << "Size of GenMET collection" <<  recmet.size() << std::endl;
     for ( cmiter i=genmet->begin(); i!=genmet->end(); i++) {
       genMet = i->pt();
       for (size_t it=0;it<hltTrgMet.size();it++) {
     _meGenMET[it] -> Fill(genMet);
     if (myTrigM[it]) _meGenMETTrgMC[it] -> Fill(genMet);
     if (myTrigM[it] && myTrigMLow[it]) _meGenMETTrg[it] -> Fill(genMet);
     if (myTrigMLow[it]) _meGenMETTrgLow[it] -> Fill(genMet);
       }
     }
   }else{
     //std::cout << "  -- No GenMET Collection with name: " << GenMETColl << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No GenMET Collection with name: "<< GenMETColl; 
   }
   
 }

void HLTJetMETValidation::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 37 of file HLTJetMETValidation.cc.

References _meGenJetEta, _meGenJetEtaTrg, _meGenJetEtaTrgLow, _meGenJetEtaTrgMC, _meGenJetPhi, _meGenJetPhiTrg, _meGenJetPhiTrgLow, _meGenJetPhiTrgMC, _meGenJetPt, _meGenJetPtTrg, _meGenJetPtTrgLow, _meGenJetPtTrgMC, _meGenMET, _meGenMETTrg, _meGenMETTrgLow, _meGenMETTrgMC, _meRecoJetEta, _meRecoJetEtaTrg, _meRecoJetEtaTrgLow, _meRecoJetEtaTrgMC, _meRecoJetPhi, _meRecoJetPhiTrg, _meRecoJetPhiTrgLow, _meRecoJetPhiTrgMC, _meRecoJetPt, _meRecoJetPtTrg, _meRecoJetPtTrgLow, _meRecoJetPtTrgMC, _meRecoMET, _meRecoMETTrg, _meRecoMETTrgLow, _meRecoMETTrgMC, DQMStore::book1D(), hltConfig_, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, HLTConfigProvider::init(), j, patternJetTrg_, patternMetTrg_, patternMuTrg_, HLTConfigProvider::removeVersion(), DQMStore::setCurrentFolder(), store, AlCaHLTBitMon_QueryRunRegistry::string, HLTConfigProvider::triggerNames(), and triggerTag_.

                                                                                {
 
   
   bool foundMuTrg = false;
   std::string trgMuNm;
   bool changedConfig;
   //--define search patterns
   /*
   TPRegexp patternJet("HLT_Jet([0-9]*)?(_v[0-9]*)?$");
   TPRegexp patternMet("HLT_(PF*)?M([E,H]*)?T([0-9]*)?(_v[0-9]*)?$");
   TPRegexp patternMu("HLT_Mu([0-9]*)?(_v[0-9]*)?$");
   */
   TPRegexp patternJet(patternJetTrg_);
   TPRegexp patternMet(patternMetTrg_);
   TPRegexp patternMu(patternMuTrg_);
 
   if (!hltConfig_.init(iRun, iSetup, "HLT", changedConfig)) {
     edm::LogError("HLTJetMET") << "Initialization of HLTConfigProvider failed!!"; 
     return;
   }
   std::vector<std::string> validTriggerNames = hltConfig_.triggerNames();
   for (size_t j = 0; j < validTriggerNames.size(); j++) {
     //---find the muon path
     if (TString(validTriggerNames[j]).Contains(patternMu)) {
       //std::cout <<validTriggerNames[j].c_str()<<std::endl;
       if (!foundMuTrg) trgMuNm = validTriggerNames[j].c_str();
       foundMuTrg = true;
     }
     //---find the jet paths
     if (TString(validTriggerNames[j]).Contains(patternJet)) {
       hltTrgJet.push_back(validTriggerNames[j]);
     }
     //---find the met paths
     if (TString(validTriggerNames[j]).Contains(patternMet)) {
       hltTrgMet.push_back(validTriggerNames[j]);
     }
   }
   
   //----set the denominator paths
   for (size_t it=0;it<hltTrgJet.size();it++) {
     if (it==0 && foundMuTrg) hltTrgJetLow.push_back(trgMuNm);//--lowest threshold uses muon path
     if (it==0 && !foundMuTrg) hltTrgJetLow.push_back(hltTrgJet[it]);//---if no muon then itself
     if (it!=0) hltTrgJetLow.push_back(hltTrgJet[it-1]);
     //std::cout<<hltTrgJet[it].c_str()<<" "<<hltTrgJetLow[it].c_str()<<std::endl;
   }
   int itm(0), itpm(0), itmh(0), itpmh(0);
   for (size_t it=0;it<hltTrgMet.size();it++) {
     if (TString(hltTrgMet[it]).Contains("PF")) {
       if (TString(hltTrgMet[it]).Contains("MHT")) {
     if( 0 == itpmh ) {
       if( foundMuTrg ) hltTrgMetLow.push_back(trgMuNm);
       else hltTrgMetLow.push_back(hltTrgMet[it]);
     }
     else hltTrgMetLow.push_back(hltTrgMet[it-1]);
     itpmh++;
       }
       if (TString(hltTrgMet[it]).Contains("MET")) {
     if( 0 == itpm ) {
       if( foundMuTrg ) hltTrgMetLow.push_back(trgMuNm);
       else hltTrgMetLow.push_back(hltTrgMet[it]);
     }
     else hltTrgMetLow.push_back(hltTrgMet[it-1]);
     itpm++;
       }
     }
     else {
       if (TString(hltTrgMet[it]).Contains("MHT")) {
     if( 0 == itmh ) {
       if( foundMuTrg ) hltTrgMetLow.push_back(trgMuNm);
       else hltTrgMetLow.push_back(hltTrgMet[it]);
     }
     else hltTrgMetLow.push_back(hltTrgMet[it-1]);
     itmh++;
       }
       if (TString(hltTrgMet[it]).Contains("MET")) {
     if( 0 == itm ) {
       if( foundMuTrg ) hltTrgMetLow.push_back(trgMuNm);
       else hltTrgMetLow.push_back(hltTrgMet[it]);
     }
     else hltTrgMetLow.push_back(hltTrgMet[it-1]);
     itm++;
       }
     }
     //std::cout<<hltTrgMet[it].c_str()<<" "<<hltTrgMetLow[it].c_str()<<std::endl;
   }
 
   //----define dqm folders and elements
   for (size_t it=0;it<hltTrgJet.size();it++) {
     //std::cout<<hltTrgJet[it].c_str()<<" "<<hltTrgJetLow[it].c_str()<<std::endl;
     //store->setCurrentFolder(triggerTag_+hltTrgJet[it].c_str());
     std::string trgPathName = HLTConfigProvider::removeVersion(triggerTag_+hltTrgJet[it].c_str());
     //std::cout << "str = " << triggerTag_+hltTrgJet[it].c_str() << std::endl;
     //std::cout << "trgPathName = " << trgPathName << std::endl;
     store->setCurrentFolder(trgPathName);
     //_meRecoJetPt= store->book1D("_meRecoJetPt","Single Reconstructed Jet Pt",100,0,500);
     _meRecoJetPt.push_back(store->book1D("_meRecoJetPt","Single Reconstructed Jet Pt",100,0,500));
     _meRecoJetPtTrgMC.push_back(store->book1D("_meRecoJetPtTrgMC","Single Reconstructed Jet Pt -- HLT Triggered",100,0,500));
     _meRecoJetPtTrg.push_back(store->book1D("_meRecoJetPtTrg","Single Reconstructed Jet Pt -- HLT Triggered",100,0,500));
     _meRecoJetPtTrgLow.push_back(store->book1D("_meRecoJetPtTrgLow","Single Reconstructed Jet Pt -- HLT Triggered Low",100,0,500));
     
     _meRecoJetEta.push_back(store->book1D("_meRecoJetEta","Single Reconstructed Jet Eta",100,-10,10));
     _meRecoJetEtaTrgMC.push_back(store->book1D("_meRecoJetEtaTrgMC","Single Reconstructed Jet Eta -- HLT Triggered",100,-10,10));
     _meRecoJetEtaTrg.push_back(store->book1D("_meRecoJetEtaTrg","Single Reconstructed Jet Eta -- HLT Triggered",100,-10,10));
     _meRecoJetEtaTrgLow.push_back(store->book1D("_meRecoJetEtaTrgLow","Single Reconstructed Jet Eta -- HLT Triggered Low",100,-10,10));
     
     _meRecoJetPhi.push_back(store->book1D("_meRecoJetPhi","Single Reconstructed Jet Phi",100,-4.,4.));
     _meRecoJetPhiTrgMC.push_back(store->book1D("_meRecoJetPhiTrgMC","Single Reconstructed Jet Phi -- HLT Triggered",100,-4.,4.));
     _meRecoJetPhiTrg.push_back(store->book1D("_meRecoJetPhiTrg","Single Reconstructed Jet Phi -- HLT Triggered",100,-4.,4.));
     _meRecoJetPhiTrgLow.push_back(store->book1D("_meRecoJetPhiTrgLow","Single Reconstructed Jet Phi -- HLT Triggered Low",100,-4.,4.));
     
     _meGenJetPt.push_back(store->book1D("_meGenJetPt","Single Generated Jet Pt",100,0,500));
     _meGenJetPtTrgMC.push_back(store->book1D("_meGenJetPtTrgMC","Single Generated Jet Pt -- HLT Triggered",100,0,500));
     _meGenJetPtTrg.push_back(store->book1D("_meGenJetPtTrg","Single Generated Jet Pt -- HLT Triggered",100,0,500));
     _meGenJetPtTrgLow.push_back(store->book1D("_meGenJetPtTrgLow","Single Generated Jet Pt -- HLT Triggered Low",100,0,500));
     
     _meGenJetEta.push_back(store->book1D("_meGenJetEta","Single Generated Jet Eta",100,-10,10));
     _meGenJetEtaTrgMC.push_back(store->book1D("_meGenJetEtaTrgMC","Single Generated Jet Eta -- HLT Triggered",100,-10,10));
     _meGenJetEtaTrg.push_back(store->book1D("_meGenJetEtaTrg","Single Generated Jet Eta -- HLT Triggered",100,-10,10));
     _meGenJetEtaTrgLow.push_back(store->book1D("_meGenJetEtaTrgLow","Single Generated Jet Eta -- HLT Triggered Low",100,-10,10));
     
     _meGenJetPhi.push_back(store->book1D("_meGenJetPhi","Single Generated Jet Phi",100,-4.,4.));
     _meGenJetPhiTrgMC.push_back(store->book1D("_meGenJetPhiTrgMC","Single Generated Jet Phi -- HLT Triggered",100,-4.,4.));
     _meGenJetPhiTrg.push_back(store->book1D("_meGenJetPhiTrg","Single Generated Jet Phi -- HLT Triggered",100,-4.,4.));
     _meGenJetPhiTrgLow.push_back(store->book1D("_meGenJetPhiTrgLow","Single Generated Jet Phi -- HLT Triggered Low",100,-4.,4.));
     
   }
   for (size_t it=0;it<hltTrgMet.size();it++) {
     //std::cout<<hltTrgMet[it].c_str()<<" "<<hltTrgMetLow[it].c_str()<<std::endl;
     //store->setCurrentFolder(triggerTag_+hltTrgMet[it].c_str());
     std::string trgPathName = HLTConfigProvider::removeVersion(triggerTag_+hltTrgMet[it].c_str());
     store->setCurrentFolder(trgPathName);
     _meRecoMET.push_back(store->book1D("_meRecoMET","Reconstructed Missing ET",100,0,500));
     _meRecoMETTrgMC.push_back(store->book1D("_meRecoMETTrgMC","Reconstructed Missing ET -- HLT Triggered",100,0,500));
     _meRecoMETTrg.push_back(store->book1D("_meRecoMETTrg","Reconstructed Missing ET -- HLT Triggered",100,0,500));
     _meRecoMETTrgLow.push_back(store->book1D("_meRecoMETTrgLow","Reconstructed Missing ET -- HLT Triggered Low",100,0,500));
     
     _meGenMET.push_back(store->book1D("_meGenMET","Generated Missing ET",100,0,500));
     _meGenMETTrgMC.push_back(store->book1D("_meGenMETTrgMC","Generated Missing ET -- HLT Triggered",100,0,500));
     _meGenMETTrg.push_back(store->book1D("_meGenMETTrg","Generated Missing ET -- HLT Triggered",100,0,500));
     _meGenMETTrgLow.push_back(store->book1D("_meGenMETTrgLow","Generated Missing ET -- HLT Triggered Low",100,0,500));
   }
 }

void HLTJetMETValidation::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 181 of file HLTJetMETValidation.cc.

References outFile_, cuy::save, and writeFile_.

 {
 
   //Write DQM thing..
   if(outFile_.size()>0)
     if (&*edm::Service<DQMStore>() && writeFile_) {
       edm::Service<DQMStore>()->save (outFile_);
     }
   
 }

void HLTJetMETValidation::getHLTResults	(	const edm::TriggerResults &	hltresults,
		const edm::TriggerNames &	triggerNames
	)

private

Definition at line 419 of file HLTJetMETValidation.cc.

References accept(), edm::HLTGlobalStatus::accept(), HLTinit_, hltTriggerMap, edm::HLTGlobalStatus::size(), AlCaHLTBitMon_QueryRunRegistry::string, trig_iter, and edm::TriggerNames::triggerName().

Referenced by analyze().

                                                                               {
 
   int ntrigs=hltresults.size();
   if (! HLTinit_){
     HLTinit_=true;
     
     //if (writeFile_) std::cout << "Number of HLT Paths: " << ntrigs << std::endl;
     for (int itrig = 0; itrig != ntrigs; ++itrig){
       std::string trigName = triggerNames.triggerName(itrig);
       // std::cout << "trigger " << itrig << ": " << trigName << std::endl; 
     }
   }
   
   for (int itrig = 0; itrig != ntrigs; ++itrig){
     std::string trigName = triggerNames.triggerName(itrig);
      bool accept=hltresults.accept(itrig);
 
      //if (accept) _triggerResults->Fill(float(itrig));
 
      // fill the trigger map
      typedef std::map<std::string,bool>::value_type valType;
      trig_iter=hltTriggerMap.find(trigName);
      if (trig_iter==hltTriggerMap.end())
        hltTriggerMap.insert(valType(trigName,accept));
      else
        trig_iter->second=accept;
   }
 }

Member Data Documentation

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetEta

private

Definition at line 115 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetEtaTrg

private

Definition at line 117 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetEtaTrgLow

private

Definition at line 118 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetEtaTrgMC

private

Definition at line 116 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPhi

private

Definition at line 119 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPhiTrg

private

Definition at line 121 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPhiTrgLow

private

Definition at line 122 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPhiTrgMC

private

Definition at line 120 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPt

private

Definition at line 111 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPtTrg

private

Definition at line 113 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPtTrgLow

private

Definition at line 114 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenJetPtTrgMC

private

Definition at line 112 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenMET

private

Definition at line 128 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenMETTrg

private

Definition at line 130 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenMETTrgLow

private

Definition at line 131 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meGenMETTrgMC

private

Definition at line 129 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetEta

private

Definition at line 102 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetEtaTrg

private

Definition at line 104 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetEtaTrgLow

private

Definition at line 105 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetEtaTrgMC

private

Definition at line 103 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPhi

private

Definition at line 106 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPhiTrg

private

Definition at line 108 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPhiTrgLow

private

Definition at line 109 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPhiTrgMC

private

Definition at line 107 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPt

private

Definition at line 98 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPtTrg

private

Definition at line 100 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPtTrgLow

private

Definition at line 101 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoJetPtTrgMC

private

Definition at line 99 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoMET

private

Definition at line 124 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoMETTrg

private

Definition at line 126 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoMETTrgLow

private

Definition at line 127 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<MonitorElement*> HLTJetMETValidation::_meRecoMETTrgMC

private

Definition at line 125 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

MonitorElement* HLTJetMETValidation::_triggerResults

private

Definition at line 135 of file HLTJetMETValidation.h.

edm::EDGetTokenT<PFJetCollection> HLTJetMETValidation::CaloJetAlgorithm

private

Definition at line 82 of file HLTJetMETValidation.h.

Referenced by analyze().

edm::EDGetTokenT<CaloMETCollection> HLTJetMETValidation::CaloMETColl

private

Definition at line 84 of file HLTJetMETValidation.h.

Referenced by analyze().

int HLTJetMETValidation::evtCnt

private

Definition at line 139 of file HLTJetMETValidation.h.

Referenced by analyze(), and HLTJetMETValidation().

edm::EDGetTokenT<GenJetCollection> HLTJetMETValidation::GenJetAlgorithm

private

Definition at line 83 of file HLTJetMETValidation.h.

Referenced by analyze().

edm::EDGetTokenT<GenMETCollection> HLTJetMETValidation::GenMETColl

private

Definition at line 85 of file HLTJetMETValidation.h.

Referenced by analyze().

HLTConfigProvider HLTJetMETValidation::hltConfig_

private

Definition at line 141 of file HLTJetMETValidation.h.

Referenced by beginRun().

bool HLTJetMETValidation::HLTinit_

private

Definition at line 152 of file HLTJetMETValidation.h.

Referenced by getHLTResults().

edm::EDGetTokenT<edm::TriggerResults> HLTJetMETValidation::HLTriggerResults

private

Definition at line 86 of file HLTJetMETValidation.h.

Referenced by analyze().

std::vector<std::string> HLTJetMETValidation::hltTrgJet

private

Definition at line 142 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<std::string> HLTJetMETValidation::hltTrgJetLow

private

Definition at line 143 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<std::string> HLTJetMETValidation::hltTrgMet

private

Definition at line 144 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::vector<std::string> HLTJetMETValidation::hltTrgMetLow

private

Definition at line 145 of file HLTJetMETValidation.h.

Referenced by analyze(), and beginRun().

std::map<std::string,bool> HLTJetMETValidation::hltTriggerMap

private

Definition at line 149 of file HLTJetMETValidation.h.

Referenced by analyze(), and getHLTResults().

std::vector<bool> HLTJetMETValidation::hlttrigs

private

Definition at line 148 of file HLTJetMETValidation.h.

std::string HLTJetMETValidation::MyTrigger

private

Definition at line 89 of file HLTJetMETValidation.h.

std::string HLTJetMETValidation::outFile_

private

Definition at line 94 of file HLTJetMETValidation.h.

Referenced by endJob().

std::string HLTJetMETValidation::patternJetTrg_

private

Definition at line 89 of file HLTJetMETValidation.h.

Referenced by beginRun().

std::string HLTJetMETValidation::patternMetTrg_

private

Definition at line 89 of file HLTJetMETValidation.h.

Referenced by beginRun().

std::string HLTJetMETValidation::patternMuTrg_

private

Definition at line 89 of file HLTJetMETValidation.h.

Referenced by beginRun().

DQMStore* HLTJetMETValidation::store

private

Definition at line 96 of file HLTJetMETValidation.h.

Referenced by beginRun(), and HLTJetMETValidation().

std::map<std::string,bool>::iterator HLTJetMETValidation::trig_iter

private

Definition at line 150 of file HLTJetMETValidation.h.

Referenced by analyze(), and getHLTResults().

edm::EDGetTokenT<TriggerEventWithRefs> HLTJetMETValidation::triggerEventObject_

private

InputTag of TriggerEventWithRefs to analyze.

Definition at line 81 of file HLTJetMETValidation.h.

Referenced by analyze().

std::string HLTJetMETValidation::triggerTag_

private

Definition at line 89 of file HLTJetMETValidation.h.

Referenced by beginRun().

bool HLTJetMETValidation::writeFile_

private

Definition at line 155 of file HLTJetMETValidation.h.

Referenced by endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation