#include <HTMHTAnalyzer.h>

Inheritance diagram for HTMHTAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &, const edm::TriggerResults &)
	Get the analysis. More...

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

	HTMHTAnalyzer (const edm::ParameterSet &)
	Constructor. More...

	~HTMHTAnalyzer () override
	Destructor. More...

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

Public Attributes
int	evtCounter

Private Attributes
double	_ptThreshold

std::string	_source

int	_trig_JetMB

int	_verbose

MonitorElement *	hHT

std::vector< std::string >	HLTPathsJetMBByName_

MonitorElement *	hMHT

MonitorElement *	hMHTPhi

MonitorElement *	hMHx

MonitorElement *	hMHy

MonitorElement *	hNevents

MonitorElement *	hNJets

MonitorElement *	jetME

std::string	metname

edm::ParameterSet	parameters

edm::InputTag	theJetCollectionForHTMHTLabel

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM monitoring source for HTMHT

Author: K. Hatakeyama, Rockefeller University

Definition at line 31 of file HTMHTAnalyzer.h.

Constructor & Destructor Documentation

HTMHTAnalyzer::HTMHTAnalyzer ( const edm::ParameterSet & pSet )

Constructor.

Definition at line 29 of file HTMHTAnalyzer.cc.

                                                         {
 
   parameters = pSet;
   _ptThreshold = 30.;
 
 }

HTMHTAnalyzer::~HTMHTAnalyzer ( )

override

Destructor.

Definition at line 37 of file HTMHTAnalyzer.cc.

37 { }

Member Function Documentation

void HTMHTAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::TriggerResults &	triggerResults
	)

Get the analysis.

Definition at line 75 of file HTMHTAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), gather_cfg::cout, edm::Event::getByLabel(), mps_fire::i, edm::HandleBase::isValid(), LogTrace, metname, nanoDQM_cfi::MHT, njet, edm::HLTGlobalStatus::size(), edm::TriggerNames::triggerIndex(), and edm::Event::triggerNames().

                                                            {
 
   LogTrace(metname)<<"[HTMHTAnalyzer] Analyze HT & MHT";
 
   jetME->Fill(4);
 
   // ==========================================================  
   // Trigger information 
   //
 
 
   //
   //
   // Check how many HLT triggers are in triggerResults 
   int ntrigs = triggerResults.size();
   if (_verbose) std::cout << "ntrigs=" << ntrigs << std::endl;
   
   //
   //
   // If index=ntrigs, this HLT trigger doesn't exist in the HLT table for this data.
   const edm::TriggerNames & triggerNames = iEvent.triggerNames(triggerResults);
   
   //
   //
   // count number of requested Jet or MB HLT paths which have fired
   for (unsigned int i=0; i!=HLTPathsJetMBByName_.size(); i++) {
     unsigned int triggerIndex = triggerNames.triggerIndex(HLTPathsJetMBByName_[i]);
     if (triggerIndex<triggerResults.size()) {
       if (triggerResults.accept(triggerIndex)) {
         _trig_JetMB++;
       }
     }
   }
   // for empty input vectors (n==0), take all HLT triggers!
   if (HLTPathsJetMBByName_.empty()) _trig_JetMB=triggerResults.size()-1;
 
   if (_trig_JetMB==0) return;
    
   // ==========================================================
 
   // **** Get the Calo Jet container
   edm::Handle<reco::CaloJetCollection> jetcoll;
 
   // **** Get the SISCone Jet container
   iEvent.getByLabel(theJetCollectionForHTMHTLabel, jetcoll);
 
   if(!jetcoll.isValid()) return;
 
   // ==========================================================
   // Reconstructed HT & MHT Information
 
   int    njet=0;
   double MHx=0.;
   double MHy=0.;
   double MHT=0.;
   double MHTPhi=0.;
   double HT=0.;
 
   for (reco::CaloJetCollection::const_iterator calojet = jetcoll->begin(); calojet!=jetcoll->end(); ++calojet){
     if (calojet->pt()>_ptThreshold){
       njet++;
       MHx += -1.*calojet->px();
       MHy += -1.*calojet->py();
       HT  += calojet->pt();
     }
   }
 
   TVector2 *tv2 = new TVector2(MHx,MHy);
   MHT   =tv2->Mod();
   MHTPhi=tv2->Phi();
 
   //std::cout << "HTMHT " << MHT << std::endl;
 
   hNevents->Fill(1.);
   hNJets->Fill(njet);
   if (njet>0){
     hMHx->Fill(MHx);
     hMHy->Fill(MHy);
     hMHT->Fill(MHT);
     hMHTPhi->Fill(MHTPhi);
     hHT->Fill(HT);
   }
 
 }

void HTMHTAnalyzer::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 41 of file HTMHTAnalyzer.cc.

References DQMStore::IBooker::book1D(), LogTrace, metname, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                              {
   evtCounter = 0;
   metname = "HTMHTAnalyzer";
 
   // PFMET information
   theJetCollectionForHTMHTLabel = parameters.getParameter<edm::InputTag>("JetCollectionForHTMHTLabel");
   _source                       = parameters.getParameter<std::string>("Source");
 
   LogTrace(metname)<<"[HTMHTAnalyzer] Parameters initialization";
   ibooker.setCurrentFolder("JetMET/MET/"+_source);
 
   HLTPathsJetMBByName_ = parameters.getParameter<std::vector<std::string > >("HLTPathsJetMB");
 
   // misc
   _verbose     = parameters.getParameter<int>("verbose");
   _ptThreshold = parameters.getParameter<double>("ptThreshold");
 
   jetME = ibooker.book1D("metReco", "metReco", 4, 1, 5);
   jetME->setBinLabel(4,"HTMHT",1);
 
   hNevents = ibooker.book1D("METTask_Nevents",   "METTask_Nevents",   1,0,1);
   hNJets   = ibooker.book1D("METTask_NJets",     "METTask_NJets",     100, 0, 100);
   hMHx     = ibooker.book1D("METTask_MHx",       "METTask_MHx",       500,-500,500);
   hMHy     = ibooker.book1D("METTask_MHy",       "METTask_MHy",       500,-500,500);
   hMHT     = ibooker.book1D("METTask_MHT",       "METTask_MHT",       500,0,1000);
   hMHTPhi  = ibooker.book1D("METTask_MhTPhi",    "METTask_MhTPhi",    80,-4,4);
   hHT      = ibooker.book1D("METTask_HT",        "METTask_HT",        500,0,2000);
 
 }

Member Data Documentation

double HTMHTAnalyzer::_ptThreshold

private

Definition at line 66 of file HTMHTAnalyzer.h.

std::string HTMHTAnalyzer::_source

private

Definition at line 56 of file HTMHTAnalyzer.h.

int HTMHTAnalyzer::_trig_JetMB

private

Definition at line 63 of file HTMHTAnalyzer.h.

int HTMHTAnalyzer::_verbose

private

Definition at line 52 of file HTMHTAnalyzer.h.

Referenced by MassReplace.MassSearchReplaceAnyInputTagVisitor::doIt(), MassReplace.MassSearchReplaceParamVisitor::enter(), cmsPerfSuite.PerfSuite::mkCandleDir(), cmsPerfSuite.PerfSuite::optionParse(), cmsPerfSuite.PerfSuite::printFlush(), cmsPerfSuite.PerfSuite::runcmd(), and cmsPerfSuite.PerfSuite::runCmdSet().

int HTMHTAnalyzer::evtCounter

Definition at line 45 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hHT

private

Definition at line 80 of file HTMHTAnalyzer.h.

std::vector<std::string > HTMHTAnalyzer::HLTPathsJetMBByName_

private

Definition at line 61 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHT

private

Definition at line 77 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHTPhi

private

Definition at line 78 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHx

private

Definition at line 75 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHy

private

Definition at line 76 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hNevents

private

Definition at line 71 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hNJets

private

Definition at line 73 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::jetME

private

Definition at line 69 of file HTMHTAnalyzer.h.

std::string HTMHTAnalyzer::metname

private

Definition at line 54 of file HTMHTAnalyzer.h.

edm::ParameterSet HTMHTAnalyzer::parameters

private

Definition at line 50 of file HTMHTAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

edm::InputTag HTMHTAnalyzer::theJetCollectionForHTMHTLabel

private

Definition at line 58 of file HTMHTAnalyzer.h.

Public Member Functions

Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation