#include <HLTJetMETValidation.h>

Inheritance diagram for HLTJetMETValidation:

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

	~HLTJetMETValidation ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &iRun, edm::EventSetup const &iSetup) override

virtual void	dqmBeginRun (edm::Run const &iRun, edm::EventSetup const &iSetup) override

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 * >	_meHLTJetEta

std::vector< MonitorElement * >	_meHLTJetEtaTrg

std::vector< MonitorElement * >	_meHLTJetEtaTrgLow

std::vector< MonitorElement * >	_meHLTJetEtaTrgMC

std::vector< MonitorElement * >	_meHLTJetPhi

std::vector< MonitorElement * >	_meHLTJetPhiTrg

std::vector< MonitorElement * >	_meHLTJetPhiTrgLow

std::vector< MonitorElement * >	_meHLTJetPhiTrgMC

std::vector< MonitorElement * >	_meHLTJetPt

std::vector< MonitorElement * >	_meHLTJetPtTrg

std::vector< MonitorElement * >	_meHLTJetPtTrgLow

std::vector< MonitorElement * >	_meHLTJetPtTrgMC

std::vector< MonitorElement * >	_meHLTMET

std::vector< MonitorElement * >	_meHLTMETTrg

std::vector< MonitorElement * >	_meHLTMETTrgLow

std::vector< MonitorElement * >	_meHLTMETTrgMC

MonitorElement *	_triggerResults

edm::EDGetTokenT< reco::CaloMETCollection >	CaloMETColl

int	evtCnt

edm::EDGetTokenT< reco::GenJetCollection >	GenJetAlgorithm

edm::EDGetTokenT< reco::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	patternJetTrg_

std::string	patternMetTrg_

std::string	patternMuTrg_

edm::EDGetTokenT< reco::PFJetCollection >	PFJetAlgorithm

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

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

std::string	triggerTag_

bool	writeFile_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 57 of file HLTJetMETValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file HLTJetMETValidation.cc.

References evtCnt.

                                                                   : 
   triggerEventObject_(consumes<TriggerEventWithRefs>(ps.getUntrackedParameter<edm::InputTag>("triggerEventObject"))),
   PFJetAlgorithm( consumes<PFJetCollection>(ps.getUntrackedParameter<edm::InputTag>( "PFJetAlgorithm" ) )),
   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","")),
   HLTinit_(false)
 {
   evtCnt=0;
 }

HLTJetMETValidation::~HLTJetMETValidation ( )

Definition at line 30 of file HLTJetMETValidation.cc.

31 {

32 }

Member Function Documentation

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

overrideprivatevirtual

Definition at line 186 of file HLTJetMETValidation.cc.

References _meGenJetEta, _meGenJetEtaTrg, _meGenJetEtaTrgLow, _meGenJetEtaTrgMC, _meGenJetPhi, _meGenJetPhiTrg, _meGenJetPhiTrgLow, _meGenJetPhiTrgMC, _meGenJetPt, _meGenJetPtTrg, _meGenJetPtTrgLow, _meGenJetPtTrgMC, _meGenMET, _meGenMETTrg, _meGenMETTrgLow, _meGenMETTrgMC, _meHLTJetEta, _meHLTJetEtaTrg, _meHLTJetEtaTrgLow, _meHLTJetEtaTrgMC, _meHLTJetPhi, _meHLTJetPhiTrg, _meHLTJetPhiTrgLow, _meHLTJetPhiTrgMC, _meHLTJetPt, _meHLTJetPtTrg, _meHLTJetPtTrgLow, _meHLTJetPtTrgMC, _meHLTMET, _meHLTMETTrg, _meHLTMETTrgLow, _meHLTMETTrgMC, CaloMETColl, evtCnt, HcalObjRepresent::Fill(), relval_steps::gen(), GenJetAlgorithm, patPFMETCorrections_cff::genJets, HepMCValidationHelper::genMet(), RazorAnalyzer::genmet, GenMETColl, edm::Event::getByToken(), getHLTResults(), HLTriggerResults, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, hltTriggerMap, mps_fire::i, edm::HandleBase::isValid(), packedPFCandidateRefMixer_cfi::pf, PFJetAlgorithm, pfJetBenchmark_cfi::pfJets, 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;
   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;
   iEvent.getByToken(HLTriggerResults,hltresults);
   if (! hltresults.isValid() ) { 
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No HLTRESULTS";    
     gotHLT=false;
   }
 
   if (gotHLT) {
     const edm::TriggerNames & triggerNames = iEvent.triggerNames(*hltresults);
     getHLTResults(*hltresults, triggerNames);
 
     //---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;
     }
   }
 
   // --- Fill histos for PFJet paths ---
   // HLT jets namely hltAK4PFJets
   Handle<PFJetCollection> pfJets;
   iEvent.getByToken( PFJetAlgorithm, pfJets );
   double pfJetPt=-1.;
   double pfJetEta=-999.;
   double pfJetPhi=-999.;
 
   if (pfJets.isValid()) { 
     //Loop over the PFJets and fill some histograms
     int jetInd = 0;
     for( PFJetCollection::const_iterator pf = pfJets->begin(); pf != pfJets->end(); ++ pf ) {
       //std::cout << "PF JET #" << jetInd << std::endl << pf->print() << std::endl;
       if (jetInd == 0){
     pfJetPt=pf->pt();
     pfJetEta=pf->eta();
     pfJetPhi=pf->phi();
     for (size_t it=0;it<hltTrgJet.size();it++) {
       _meHLTJetPt[it]->Fill( pfJetPt );
       _meHLTJetEta[it]->Fill( pfJetEta );
       _meHLTJetPhi[it]->Fill( pfJetPhi );
       if (myTrigJ[it]) _meHLTJetPtTrgMC[it]->Fill( pfJetPt );
       if (myTrigJ[it]) _meHLTJetEtaTrgMC[it]->Fill( pfJetEta );
       if (myTrigJ[it]) _meHLTJetPhiTrgMC[it]->Fill( pfJetPhi );
       if (myTrigJ[it] && myTrigJLow[it]) _meHLTJetPtTrg[it]->Fill( pfJetPt );
       if (myTrigJ[it] && myTrigJLow[it]) _meHLTJetEtaTrg[it]->Fill( pfJetEta );
       if (myTrigJ[it] && myTrigJLow[it]) _meHLTJetPhiTrg[it]->Fill( pfJetPhi );
       if (myTrigJLow[it]) _meHLTJetPtTrgLow[it]->Fill( pfJetPt );
       if (myTrigJLow[it]) _meHLTJetEtaTrgLow[it]->Fill( pfJetEta );
       if (myTrigJLow[it]) _meHLTJetPhiTrgLow[it]->Fill( pfJetPhi );
     }
     jetInd++;
       }
     }//loop over pfjets
   } 
   else{
     //std::cout << "  -- No PFJets" << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No PFJets"; 
   }
 
   //GenJets
   Handle<GenJetCollection> genJets;
   iEvent.getByToken( GenJetAlgorithm, genJets );
   double genJetPt=-1.;
   double genJetEta=-999.;
   double genJetPhi=-999.;
 
   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 ) {
       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++;
       }
     }
   } 
   else{
     //std::cout << "  -- No GenJets" << std::endl;
     //if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No GenJets"; 
   }
 
   // --- Fill histos for PFMET paths ---
   // HLT MET namely hltmet
   edm::Handle<CaloMETCollection> recmet;
   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++) {
     _meHLTMET[it] -> Fill(calMet);
     if (myTrigM.size() > it && myTrigM[it]) _meHLTMETTrgMC[it] -> Fill(calMet);
     if (myTrigM.size() > it && myTrigMLow.size() > it && myTrigM[it] && myTrigMLow[it]) _meHLTMETTrg[it] -> Fill(calMet);
     if (myTrigMLow.size() > it && myTrigMLow[it]) _meHLTMETTrgLow[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;
   iEvent.getByToken(GenMETColl, genmet);
 
   double genMet=-1;
   if (genmet.isValid()) { 
     typedef GenMETCollection::const_iterator cmiter;
     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.size() > it && myTrigM[it]) _meGenMETTrgMC[it] -> Fill(genMet);
     if (myTrigM.size() > it && myTrigMLow.size() > it && myTrigM[it] && myTrigMLow[it]) _meGenMETTrg[it] -> Fill(genMet);
     if (myTrigMLow.size() > it && 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::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 126 of file HLTJetMETValidation.cc.

References _meGenJetEta, _meGenJetEtaTrg, _meGenJetEtaTrgLow, _meGenJetEtaTrgMC, _meGenJetPhi, _meGenJetPhiTrg, _meGenJetPhiTrgLow, _meGenJetPhiTrgMC, _meGenJetPt, _meGenJetPtTrg, _meGenJetPtTrgLow, _meGenJetPtTrgMC, _meGenMET, _meGenMETTrg, _meGenMETTrgLow, _meGenMETTrgMC, _meHLTJetEta, _meHLTJetEtaTrg, _meHLTJetEtaTrgLow, _meHLTJetEtaTrgMC, _meHLTJetPhi, _meHLTJetPhiTrg, _meHLTJetPhiTrgLow, _meHLTJetPhiTrgMC, _meHLTJetPt, _meHLTJetPtTrg, _meHLTJetPtTrgLow, _meHLTJetPtTrgMC, _meHLTMET, _meHLTMETTrg, _meHLTMETTrgLow, _meHLTMETTrgMC, DQMStore::IBooker::book1D(), hltTrgJet, hltTrgMet, HLTConfigProvider::removeVersion(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and triggerTag_.

 {
 
   //----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;
     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;
     iBooker.setCurrentFolder(trgPathName);
     _meHLTJetPt.push_back(iBooker.book1D("_meHLTJetPt","Single HLT Jet Pt",100,0,500));
     _meHLTJetPtTrgMC.push_back(iBooker.book1D("_meHLTJetPtTrgMC","Single HLT Jet Pt - HLT Triggered",100,0,500));
     _meHLTJetPtTrg.push_back(iBooker.book1D("_meHLTJetPtTrg","Single HLT Jet Pt - HLT Triggered",100,0,500));
     _meHLTJetPtTrgLow.push_back(iBooker.book1D("_meHLTJetPtTrgLow","Single HLT Jet Pt - HLT Triggered Low",100,0,500));
     
     _meHLTJetEta.push_back(iBooker.book1D("_meHLTJetEta","Single HLT Jet Eta",100,-10,10));
     _meHLTJetEtaTrgMC.push_back(iBooker.book1D("_meHLTJetEtaTrgMC","Single HLT Jet Eta - HLT Triggered",100,-10,10));
     _meHLTJetEtaTrg.push_back(iBooker.book1D("_meHLTJetEtaTrg","Single HLT Jet Eta - HLT Triggered",100,-10,10));
     _meHLTJetEtaTrgLow.push_back(iBooker.book1D("_meHLTJetEtaTrgLow","Single HLT Jet Eta - HLT Triggered Low",100,-10,10));
     
     _meHLTJetPhi.push_back(iBooker.book1D("_meHLTJetPhi","Single HLT Jet Phi",100,-4.,4.));
     _meHLTJetPhiTrgMC.push_back(iBooker.book1D("_meHLTJetPhiTrgMC","Single HLT Jet Phi - HLT Triggered",100,-4.,4.));
     _meHLTJetPhiTrg.push_back(iBooker.book1D("_meHLTJetPhiTrg","Single HLT Jet Phi - HLT Triggered",100,-4.,4.));
     _meHLTJetPhiTrgLow.push_back(iBooker.book1D("_meHLTJetPhiTrgLow","Single HLT Jet Phi - HLT Triggered Low",100,-4.,4.));
     
     _meGenJetPt.push_back(iBooker.book1D("_meGenJetPt","Single Generated Jet Pt",100,0,500));
     _meGenJetPtTrgMC.push_back(iBooker.book1D("_meGenJetPtTrgMC","Single Generated Jet Pt - HLT Triggered",100,0,500));
     _meGenJetPtTrg.push_back(iBooker.book1D("_meGenJetPtTrg","Single Generated Jet Pt - HLT Triggered",100,0,500));
     _meGenJetPtTrgLow.push_back(iBooker.book1D("_meGenJetPtTrgLow","Single Generated Jet Pt - HLT Triggered Low",100,0,500));
     
     _meGenJetEta.push_back(iBooker.book1D("_meGenJetEta","Single Generated Jet Eta",100,-10,10));
     _meGenJetEtaTrgMC.push_back(iBooker.book1D("_meGenJetEtaTrgMC","Single Generated Jet Eta - HLT Triggered",100,-10,10));
     _meGenJetEtaTrg.push_back(iBooker.book1D("_meGenJetEtaTrg","Single Generated Jet Eta - HLT Triggered",100,-10,10));
     _meGenJetEtaTrgLow.push_back(iBooker.book1D("_meGenJetEtaTrgLow","Single Generated Jet Eta - HLT Triggered Low",100,-10,10));
     
     _meGenJetPhi.push_back(iBooker.book1D("_meGenJetPhi","Single Generated Jet Phi",100,-4.,4.));
     _meGenJetPhiTrgMC.push_back(iBooker.book1D("_meGenJetPhiTrgMC","Single Generated Jet Phi - HLT Triggered",100,-4.,4.));
     _meGenJetPhiTrg.push_back(iBooker.book1D("_meGenJetPhiTrg","Single Generated Jet Phi - HLT Triggered",100,-4.,4.));
     _meGenJetPhiTrgLow.push_back(iBooker.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;
     std::string trgPathName = HLTConfigProvider::removeVersion(triggerTag_+hltTrgMet[it].c_str());
     iBooker.setCurrentFolder(trgPathName);
     _meHLTMET.push_back(iBooker.book1D("_meHLTMET","HLT Missing ET",100,0,500));
     _meHLTMETTrgMC.push_back(iBooker.book1D("_meHLTMETTrgMC","HLT Missing ET - HLT Triggered",100,0,500));
     _meHLTMETTrg.push_back(iBooker.book1D("_meHLTMETTrg","HLT Missing ET - HLT Triggered",100,0,500));
     _meHLTMETTrgLow.push_back(iBooker.book1D("_meHLTMETTrgLow","HLT Missing ET - HLT Triggered Low",100,0,500));
     
     _meGenMET.push_back(iBooker.book1D("_meGenMET","Generated Missing ET",100,0,500));
     _meGenMETTrgMC.push_back(iBooker.book1D("_meGenMETTrgMC","Generated Missing ET - HLT Triggered",100,0,500));
     _meGenMETTrg.push_back(iBooker.book1D("_meGenMETTrg","Generated Missing ET - HLT Triggered",100,0,500));
     _meGenMETTrgLow.push_back(iBooker.book1D("_meGenMETTrgLow","Generated Missing ET - HLT Triggered Low",100,0,500));
   }
 
 }

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 40 of file HLTJetMETValidation.cc.

References hltConfig_, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, HLTConfigProvider::init(), patternJetTrg_, patternMetTrg_, patternMuTrg_, AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerNames().

                                                                                 {
 
   bool foundMuTrg = false;
   std::string trgMuNm;
   bool changedConfig = true;
   
   //--define search patterns
   TPRegexp patternJet(patternJetTrg_);
   TPRegexp patternMet(patternMetTrg_);
   TPRegexp patternMu(patternMuTrg_);
 
   if (!hltConfig_.init(iRun, iSetup, "HLT", changedConfig)) {
     edm::LogError("HLTJetMETValidation") << "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;
   }
 
 }

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

private

Definition at line 392 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;
     
     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;
   }
 }