#include <HTMonitor.h>

Inheritance diagram for HTMonitor:

Classes
struct	HTME

struct	MEHTbinning

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

	~HTMonitor () override

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

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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)

static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)

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
void	analyze (edm::Event const &iEvent, edm::EventSetup const &iSetup) override

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

void	bookME (DQMStore::IBooker &, HTME &me, const std::string &histname, const std::string &histtitle, int nbins, double xmin, double xmax)

void	bookME (DQMStore::IBooker &, HTME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX)

void	bookME (DQMStore::IBooker &, HTME &me, const std::string &histname, const std::string &histtitle, int nbinsX, double xmin, double xmax, double ymin, double ymax)

void	bookME (DQMStore::IBooker &, HTME &me, const std::string &histname, const std::string &histtitle, int nbinsX, double xmin, double xmax, int nbinsY, double ymin, double ymax)

void	bookME (DQMStore::IBooker &, HTME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX, const std::vector< double > &binningY)

void	setHTitle (HTME &me, const std::string &titleX, const std::string &titleY)

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)

Static Private Member Functions
static MEHTbinning	getHistoLSPSet (const edm::ParameterSet &pset)

static MEHTbinning	getHistoPSet (const edm::ParameterSet &pset)

Private Attributes
HTME	deltaphij1j2ME_

HTME	deltaphimetj1ME_

std::unique_ptr< GenericTriggerEventFlag >	den_genTriggerEventFlag_

double	dEtaCut_

edm::InputTag	eleInputTag_

StringCutObjectSelector< reco::GsfElectron, true >	eleSelection_

edm::EDGetTokenT< reco::GsfElectronCollection >	eleToken_

std::string	folderName_

std::string	histoSuffix_

MEHTbinning	ht_binning_

std::vector< double >	ht_variable_binning_

HTME	htVsLS_

edm::InputTag	jetInputTag_

StringCutObjectSelector< reco::Jet, true >	jetSelection_

StringCutObjectSelector< reco::Jet, true >	jetSelection_HT_

edm::EDGetTokenT< reco::JetView >	jetToken_

MEHTbinning	ls_binning_

edm::InputTag	metInputTag_

StringCutObjectSelector< reco::MET, true >	metSelection_

edm::EDGetTokenT< reco::PFMETCollection >	metToken_

edm::InputTag	muoInputTag_

StringCutObjectSelector< reco::Muon, true >	muoSelection_

edm::EDGetTokenT< reco::MuonCollection >	muoToken_

unsigned	nelectrons_

unsigned	njets_

unsigned	nmuons_

std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_

HTME	qME_variableBinning_

quant	quantity_

edm::InputTag	vtxInputTag_

edm::EDGetTokenT< reco::VertexCollection >	vtxToken_

std::vector< bool >	warningPrinted4token_

Static Private Attributes
static double	MAXedge_PHI = 3.2

static int	Nbin_PHI = 64

static MEHTbinning	phi_binning_

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

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 50 of file HTMonitor.h.

Constructor & Destructor Documentation

HTMonitor::HTMonitor ( const edm::ParameterSet & iConfig )

Definition at line 12 of file HTMonitor.cc.

References Exception, edm::ParameterSet::getParameter(), HT, MJJ, tauVariables_cff::quantity, quantity_, SOFTDROP, AlCaHLTBitMon_QueryRunRegistry::string, warningPrinted4token_, and ~HTMonitor().

                                                      : 
   folderName_             ( iConfig.getParameter<std::string>("FolderName") )
   , metInputTag_          ( iConfig.getParameter<edm::InputTag>    ("met")          )
   , jetInputTag_          ( iConfig.getParameter<edm::InputTag>    ("jets")         )
   , eleInputTag_          ( iConfig.getParameter<edm::InputTag>    ("electrons")    ) 
   , muoInputTag_          ( iConfig.getParameter<edm::InputTag>    ("muons")        ) 
   , vtxInputTag_          ( iConfig.getParameter<edm::InputTag>    ("vertices")     ) 
   , metToken_             ( consumes<reco::PFMETCollection>        ( metInputTag_ ) )
   , jetToken_             ( mayConsume<reco::JetView>              ( jetInputTag_ ) )
   , eleToken_             ( mayConsume<reco::GsfElectronCollection>( eleInputTag_ ) )
   , muoToken_             ( mayConsume<reco::MuonCollection>       ( muoInputTag_ ) )
   , vtxToken_             ( mayConsume<reco::VertexCollection>     ( vtxInputTag_ ) )
   , ht_variable_binning_  ( iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("htBinning") )
   , ht_binning_           ( getHistoPSet (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>   ("htPSet")    ) )
   , ls_binning_           ( getHistoPSet (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>   ("lsPSet")     ) )
   , num_genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"),consumesCollector(), *this))
   , den_genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"),consumesCollector(), *this))
   , metSelection_ ( iConfig.getParameter<std::string>("metSelection") )
   , jetSelection_ ( iConfig.getParameter<std::string>("jetSelection") )
   , eleSelection_ ( iConfig.getParameter<std::string>("eleSelection") )
   , muoSelection_ ( iConfig.getParameter<std::string>("muoSelection") )
   , jetSelection_HT_ ( iConfig.getParameter<std::string>("jetSelection_HT") )
   , njets_      ( iConfig.getParameter<unsigned>("njets" )      )
   , nelectrons_ ( iConfig.getParameter<unsigned>("nelectrons" ) )
   , nmuons_     ( iConfig.getParameter<unsigned>("nmuons" )     )
   , dEtaCut_    ( iConfig.getParameter<double>("dEtaCut")       )
 {
   /* mia: THIS CODE SHOULD BE DELETED !!!! */
     string quantity = iConfig.getParameter<std::string>("quantity");
     if(quantity == "HT")
     {
         quantity_ = HT;
     }
     else if(quantity == "Mjj")
     {
         quantity_ = MJJ;
     }
     else if(quantity == "softdrop")
     {
         quantity_ = SOFTDROP;
     }
     else
     {
         throw cms::Exception("quantity not defined") << "the quantity '" << quantity << "' is undefined. Please check your config!" << std::endl;
     }
 
     // this vector has to be alligned to the the number of Tokens accessed by this module
     warningPrinted4token_.push_back(false); // PFMETCollection
     warningPrinted4token_.push_back(false); // JetCollection
     warningPrinted4token_.push_back(false); // GsfElectronCollection
     warningPrinted4token_.push_back(false); // MuonCollection
     warningPrinted4token_.push_back(false); // VertexCollection
 }

HTMonitor::~HTMonitor ( )

overridedefault

Referenced by HTMonitor().

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 190 of file HTMonitor.cc.

References hiPixelPairStep_cff::deltaPhi, deltaphij1j2ME_, deltaphimetj1ME_, den_genTriggerEventFlag_, HTMonitor::HTME::denominator, dEtaCut_, MillePedeFileConverter_cfg::e, electrons_cff::electrons, eleInputTag_, eleSelection_, eleToken_, reco::tau::disc::Eta(), HcalObjRepresent::Fill(), edm::Event::getByToken(), htVsLS_, edm::EventBase::id(), edm::HandleBase::isValid(), jetInputTag_, fwrapper::jets, jetSelection_, jetSelection_HT_, jetToken_, edm::InputTag::label(), eostools::ls(), edm::EventID::luminosityBlock(), funct::m, metInputTag_, metSelection_, metToken_, MJJ, muoInputTag_, extraflags_cff::muons, muoSelection_, muoToken_, nelectrons_, njets_, nmuons_, num_genTriggerEventFlag_, HTMonitor::HTME::numerator, p4, phi, reco::LeafCandidate::phi(), reco::Vertex::position(), qME_variableBinning_, quantity_, edm::View< T >::size(), SOFTDROP, findQualityFiles::v, badGlobalMuonTaggersAOD_cff::vtx, vtxInputTag_, vtxToken_, and warningPrinted4token_.

                                                                             {
 
   // Filter out events if Trigger Filtering is requested
   if (den_genTriggerEventFlag_->on() && ! den_genTriggerEventFlag_->accept( iEvent, iSetup) ) return;
 
   edm::Handle<reco::PFMETCollection> metHandle;
   iEvent.getByToken( metToken_, metHandle );
   if ( !metHandle.isValid() ) {
     if (!warningPrinted4token_[0]) {
       edm::LogWarning("HTMonitor") << "skipping events because the collection " << metInputTag_.label().c_str() << " is not available";
       warningPrinted4token_[0] = true;
     }
     return;
   }
   reco::PFMET pfmet = metHandle->front();
   if ( ! metSelection_( pfmet ) ) return;
 
   edm::Handle<reco::JetView> jetHandle; //add a configurable jet collection & jet pt selection
   iEvent.getByToken( jetToken_, jetHandle );
   if ( !jetHandle.isValid() ) {
     if (!warningPrinted4token_[1]) {
       edm::LogWarning("HTMonitor") << "skipping events because the collection " << jetInputTag_.label().c_str() << " is not available";
       warningPrinted4token_[1] = true;
     }
     return;
   }
   std::vector<reco::Jet> jets;
   if ( jetHandle->size() < njets_ ) return;
   for ( auto const & j : *jetHandle ) {
     if ( jetSelection_( j ) ) {
       jets.push_back(j);
     }
   }
 
   if ( jets.size() < njets_ ) return;
 
   float deltaPhi_met_j1 = 10.0;
   float deltaPhi_j1_j2 = 10.0;
 
   if (!jets.empty()) deltaPhi_met_j1 = fabs( deltaPhi( pfmet.phi(),  jets[0].phi() ));
   if (jets.size() >= 2) deltaPhi_j1_j2 = fabs( deltaPhi( jets[0].phi(),  jets[1].phi() ));
 
   std::vector<reco::GsfElectron> electrons;
   edm::Handle<reco::GsfElectronCollection> eleHandle;
   iEvent.getByToken( eleToken_, eleHandle );
   if ( eleHandle.isValid() ) {
     if ( eleHandle->size() < nelectrons_ ) return;
     for ( auto const & e : *eleHandle ) {
       if ( eleSelection_( e ) ) electrons.push_back(e);
     }
     if ( electrons.size() < nelectrons_ ) return;
   } else {
     if (!warningPrinted4token_[2]) {
       warningPrinted4token_[2] = true;
       if ( eleInputTag_.label().empty() ) 
     edm::LogWarning("HTMonitor") << "GsfElectronCollection not set";
       else
     edm::LogWarning("HTMonitor") << "skipping events because the collection " << eleInputTag_.label().c_str() << " is not available";
     }
     if ( !eleInputTag_.label().empty() ) return;
   }
   
 
   reco::Vertex vtx;
   edm::Handle<reco::VertexCollection> vtxHandle;
   iEvent.getByToken(vtxToken_, vtxHandle);
   if ( vtxHandle.isValid() ) {
     for (auto const & v : *vtxHandle) {
       bool isFake =  v.isFake() ;
       
       if (!isFake) {
     vtx = v;
     break;
       }
     }
   } else {
     if (!warningPrinted4token_[3]) {
       warningPrinted4token_[3] = true;
       if ( vtxInputTag_.label().empty() )
     edm::LogWarning("HTMonitor") << "VertexCollection not set";
       else
     edm::LogWarning("HTMonitor") << "skipping events because the collection " << vtxInputTag_.label().c_str() << " is not available";
     }
     if ( !vtxInputTag_.label().empty() ) return;
   }
 
   std::vector<reco::Muon> muons;
   edm::Handle<reco::MuonCollection> muoHandle;
   iEvent.getByToken( muoToken_, muoHandle );
   if ( muoHandle.isValid() ) {
     if ( muoHandle->size() < nmuons_ ) return;
     for ( auto const & m : *muoHandle ) {
       if ( muoSelection_( m ) && m.isGlobalMuon() && m.isPFMuon() && m.globalTrack()->normalizedChi2() < 10. && m.globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && m.numberOfMatchedStations() > 1 && fabs(m.muonBestTrack()->dxy(vtx.position())) < 0.2 && fabs(m.muonBestTrack()->dz(vtx.position())) < 0.5 && m.innerTrack()->hitPattern().numberOfValidPixelHits() > 0 && m.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 )  muons.push_back(m);
     }
     if ( muons.size() < nmuons_ ) return;
   } else {
     if (!warningPrinted4token_[4]) {
       warningPrinted4token_[4] = true;
       if ( muoInputTag_.label().empty() )
     edm::LogWarning("HTMonitor") << "MuonCollection not set";
       else
     edm::LogWarning("HTMonitor") << "skipping events because the collection " << muoInputTag_.label().c_str() << " is not available";
     }
     if ( !muoInputTag_.label().empty() ) return;
   }
 
   // fill histograms
   switch(quantity_)
   {
     case HT:
     {
         float ht = 0.0;
         for ( auto const & j : *jetHandle )
         {
             if ( jetSelection_HT_(j)) ht += j.pt();
         }
 
         // filling histograms (denominator)  
         qME_variableBinning_.denominator -> Fill(ht);
 
         deltaphimetj1ME_.denominator -> Fill(deltaPhi_met_j1);
         deltaphij1j2ME_.denominator -> Fill(deltaPhi_j1_j2);
 
         int ls = iEvent.id().luminosityBlock();
         htVsLS_.denominator -> Fill(ls, ht);
 
         // applying selection for numerator
         if (num_genTriggerEventFlag_->on() && ! num_genTriggerEventFlag_->accept( iEvent, iSetup) ) return;
 
         // filling histograms (num_genTriggerEventFlag_)  
         qME_variableBinning_.numerator -> Fill(ht);
 
         htVsLS_.numerator -> Fill(ls, ht);
         deltaphimetj1ME_.numerator  -> Fill(deltaPhi_met_j1); 
         deltaphij1j2ME_.numerator  -> Fill(deltaPhi_j1_j2);
         break;
     }
 
     case MJJ:
     {
         if (jets.size() < 2) return;
 
         // deltaEta cut
         if(fabs(jets[0].p4().Eta() - jets[1].p4().Eta()) >= dEtaCut_) return;
         float mjj = (jets[0].p4() + jets[1].p4()).M();
 
         qME_variableBinning_.denominator -> Fill(mjj);
 
         // applying selection for numerator
         if (num_genTriggerEventFlag_->on() && ! num_genTriggerEventFlag_->accept( iEvent, iSetup) ) return;
 
         qME_variableBinning_.numerator -> Fill(mjj);
         break;
     }
 
     case SOFTDROP:
     {
         if (jets.size() < 2) return;
 
         // deltaEta cut
         if(fabs(jets[0].p4().Eta() - jets[1].p4().Eta()) >= dEtaCut_) return;
 
         float softdrop = jets[0].p4().M();
 
         qME_variableBinning_.denominator -> Fill(softdrop);
 
         // applying selection for numerator
         if (num_genTriggerEventFlag_->on() && ! num_genTriggerEventFlag_->accept( iEvent, iSetup) ) return;
 
         qME_variableBinning_.numerator -> Fill(softdrop);
         break;
     }
   }
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 131 of file HTMonitor.cc.

References bookME(), deltaphij1j2ME_, deltaphimetj1ME_, den_genTriggerEventFlag_, folderName_, ht_binning_, ht_variable_binning_, htVsLS_, ls_binning_, MJJ, HTMonitor::MEHTbinning::nbins, num_genTriggerEventFlag_, phi_binning_, qME_variableBinning_, quantity_, DQMStore::IBooker::setCurrentFolder(), setHTitle(), SOFTDROP, AlCaHLTBitMon_QueryRunRegistry::string, HTMonitor::MEHTbinning::xmax, and HTMonitor::MEHTbinning::xmin.

 {
   std::string histname, histtitle;
 
   std::string currentFolder = folderName_ ;
   ibooker.setCurrentFolder(currentFolder);
 
   switch(quantity_)
   {
     case HT:
     {
         histname = "ht_variable"; histtitle = "HT";
         bookME(ibooker,qME_variableBinning_,histname,histtitle,ht_variable_binning_);
         setHTitle(qME_variableBinning_,"HT [GeV]","events / [GeV]");
 
         histname = "htVsLS"; histtitle = "HT vs LS";
         bookME(ibooker,htVsLS_,histname,histtitle,ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax,ht_binning_.xmin, ht_binning_.xmax);
         setHTitle(htVsLS_,"LS","HT [GeV]");
 
         histname = "deltaphi_metjet1"; histtitle = "DPHI_METJ1";
         bookME(ibooker,deltaphimetj1ME_,histname,histtitle,phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax);
         setHTitle(deltaphimetj1ME_,"delta phi (met, j1)","events / 0.1 rad");
 
         histname = "deltaphi_jet1jet2"; histtitle = "DPHI_J1J2";
         bookME(ibooker,deltaphij1j2ME_,histname,histtitle,phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax);
         setHTitle(deltaphij1j2ME_,"delta phi (j1, j2)","events / 0.1 rad");
         break;
     }
 
 
     case MJJ:
     {
         histname = "mjj_variable"; histtitle = "Mjj";
         bookME(ibooker,qME_variableBinning_,histname,histtitle, ht_variable_binning_);
         setHTitle(qME_variableBinning_,"Mjj [GeV]","events / [GeV]");
         break;
     }
 
 
     case SOFTDROP:
     {
         histname = "softdrop_variable"; histtitle = "softdropmass";
         bookME(ibooker,qME_variableBinning_,histname,histtitle, ht_variable_binning_);
         setHTitle(qME_variableBinning_,"leading jet softdropmass [GeV]","events / [GeV]");
         break;
     }
   }
 
   // Initialize the GenericTriggerEventFlag
   if ( num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() ) num_genTriggerEventFlag_->initRun( iRun, iSetup );
   if ( den_genTriggerEventFlag_ && den_genTriggerEventFlag_->on() ) den_genTriggerEventFlag_->initRun( iRun, iSetup );
 }

void HTMonitor::bookME	(	DQMStore::IBooker &	ibooker,
		HTME &	me,
		const std::string &	histname,
		const std::string &	histtitle,
		int	nbins,
		double	xmin,
		double	xmax
	)

protected

Definition at line 95 of file HTMonitor.cc.

References DQMStore::IBooker::book1D(), HTMonitor::HTME::denominator, and HTMonitor::HTME::numerator.

Referenced by bookHistograms().

 {
   me.numerator   = ibooker.book1D(histname+"_numerator",   histtitle+" (numerator)",   nbins, min, max);
   me.denominator = ibooker.book1D(histname+"_denominator", histtitle+" (denominator)", nbins, min, max);
 }

void HTMonitor::bookME	(	DQMStore::IBooker &	ibooker,
		HTME &	me,
		const std::string &	histname,
		const std::string &	histtitle,
		const std::vector< double > &	binningX
	)

protected

Definition at line 100 of file HTMonitor.cc.

References DQMStore::IBooker::book1D(), HTMonitor::HTME::denominator, pileupCalc::nbins, and HTMonitor::HTME::numerator.

 {
   int nbins = binning.size()-1;
   std::vector<float> fbinning(binning.begin(),binning.end());
   float* arr = &fbinning[0];
   me.numerator   = ibooker.book1D(histname+"_numerator",   histtitle+" (numerator)",   nbins, arr);
   me.denominator = ibooker.book1D(histname+"_denominator", histtitle+" (denominator)", nbins, arr);
 }

void HTMonitor::bookME	(	DQMStore::IBooker &	ibooker,
		HTME &	me,
		const std::string &	histname,
		const std::string &	histtitle,
		int	nbinsX,
		double	xmin,
		double	xmax,
		double	ymin,
		double	ymax
	)

protected

Definition at line 108 of file HTMonitor.cc.

References DQMStore::IBooker::bookProfile(), HTMonitor::HTME::denominator, and HTMonitor::HTME::numerator.

 {
   me.numerator   = ibooker.bookProfile(histname+"_numerator",   histtitle+" (numerator)",   nbinsX, xmin, xmax, ymin, ymax);
   me.denominator = ibooker.bookProfile(histname+"_denominator", histtitle+" (denominator)", nbinsX, xmin, xmax, ymin, ymax);
 }

void HTMonitor::bookME	(	DQMStore::IBooker &	ibooker,
		HTME &	me,
		const std::string &	histname,
		const std::string &	histtitle,
		int	nbinsX,
		double	xmin,
		double	xmax,
		int	nbinsY,
		double	ymin,
		double	ymax
	)

protected

Definition at line 113 of file HTMonitor.cc.

References DQMStore::IBooker::book2D(), HTMonitor::HTME::denominator, and HTMonitor::HTME::numerator.

 {
   me.numerator   = ibooker.book2D(histname+"_numerator",   histtitle+" (numerator)",   nbinsX, xmin, xmax, nbinsY, ymin, ymax);
   me.denominator = ibooker.book2D(histname+"_denominator", histtitle+" (denominator)", nbinsX, xmin, xmax, nbinsY, ymin, ymax);
 }

void HTMonitor::bookME	(	DQMStore::IBooker &	ibooker,
		HTME &	me,
		const std::string &	histname,
		const std::string &	histtitle,
		const std::vector< double > &	binningX,
		const std::vector< double > &	binningY
	)

protected

Definition at line 118 of file HTMonitor.cc.

References DQMStore::IBooker::book2D(), HTMonitor::HTME::denominator, and HTMonitor::HTME::numerator.

 {
   int nbinsX = binningX.size()-1;
   std::vector<float> fbinningX(binningX.begin(),binningX.end());
   float* arrX = &fbinningX[0];
   int nbinsY = binningY.size()-1;
   std::vector<float> fbinningY(binningY.begin(),binningY.end());
   float* arrY = &fbinningY[0];
 
   me.numerator   = ibooker.book2D(histname+"_numerator",   histtitle+" (numerator)",   nbinsX, arrX, nbinsY, arrY);
   me.denominator = ibooker.book2D(histname+"_denominator", histtitle+" (denominator)", nbinsX, arrX, nbinsY, arrY);
 }

void HTMonitor::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 379 of file HTMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), create_public_pileup_plots::bins, DEFINE_FWK_MODULE, fillHistoLSPSetDescription(), fillHistoPSetDescription(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   edm::ParameterSetDescription desc;
   desc.add<std::string>  ( "FolderName", "HLT/HT" );
   desc.add<std::string>  ( "quantity", "HT" );
 
 
   desc.add<edm::InputTag>( "met",      edm::InputTag("pfMet") );
   desc.add<edm::InputTag>( "jets",     edm::InputTag("ak4PFJetsCHS") );
   desc.add<edm::InputTag>( "electrons",edm::InputTag("gedGsfElectrons") );
   desc.add<edm::InputTag>( "muons",    edm::InputTag("muons") );
   desc.add<edm::InputTag>( "vertices",edm::InputTag("offlinePrimaryVertices") );
   desc.add<std::string>("metSelection", "pt > 0");
   desc.add<std::string>("jetSelection", "pt > 0");
   desc.add<std::string>("eleSelection", "pt > 0");
   desc.add<std::string>("muoSelection", "pt > 0");
   desc.add<std::string>("jetSelection_HT", "pt > 30 && eta < 2.5");
   desc.add<unsigned>("njets",      0);
   desc.add<unsigned>("nelectrons", 0);
   desc.add<unsigned>("nmuons",     0);
   desc.add<double>("dEtaCut",      1.3);
 
   edm::ParameterSetDescription genericTriggerEventPSet;
   genericTriggerEventPSet.add<bool>("andOr");
   genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi") );
   genericTriggerEventPSet.add<std::vector<int> >("dcsPartitions",{});
   genericTriggerEventPSet.add<bool>("andOrDcs", false);
   genericTriggerEventPSet.add<bool>("errorReplyDcs", true);
   genericTriggerEventPSet.add<std::string>("dbLabel","");
   genericTriggerEventPSet.add<bool>("andOrHlt", true);
   genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT") );
   genericTriggerEventPSet.add<std::vector<std::string> >("hltPaths",{});
   genericTriggerEventPSet.add<std::string>("hltDBKey","");
   genericTriggerEventPSet.add<bool>("errorReplyHlt",false);
   genericTriggerEventPSet.add<unsigned int>("verbosityLevel",1);
 
   desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
   desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
 
   edm::ParameterSetDescription histoPSet;
   edm::ParameterSetDescription htPSet;
   fillHistoPSetDescription(htPSet);
   histoPSet.add<edm::ParameterSetDescription>("htPSet", htPSet);
   std::vector<double> bins = {0.,20.,40.,60.,80.,90.,100.,110.,120.,130.,140.,150.,160.,170.,180.,190.,200.,220.,240.,260.,280.,300.,350.,400.,450.,500.,550.,600.,650.,700.,750.,800.,850.,900.,950.,1000.,1050.,1100.,1200.,1300.,1400.,1500.,2000.,2500.};
   histoPSet.add<std::vector<double> >("htBinning", bins);
 
   edm::ParameterSetDescription lsPSet;
   fillHistoLSPSetDescription(lsPSet);
   histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
 
   desc.add<edm::ParameterSetDescription>("histoPSet",histoPSet);
 
   descriptions.add("htMonitoring", desc);
 }

void HTMonitor::fillHistoLSPSetDescription ( edm::ParameterSetDescription & pset )

static

Definition at line 372 of file HTMonitor.cc.

References edm::ParameterSetDescription::add().

Referenced by fillDescriptions().

 {
   pset.add<unsigned int>   ( "nbins", 2500);
   pset.add<double>         ( "xmin",     0.);
   pset.add<double>         ( "xmax",  2500.);
 }

void HTMonitor::fillHistoPSetDescription ( edm::ParameterSetDescription & pset )

static

Definition at line 365 of file HTMonitor.cc.

References edm::ParameterSetDescription::add().

Referenced by fillDescriptions().

 {
   pset.add<unsigned int>   ( "nbins");
   pset.add<double>( "xmin" );
   pset.add<double>( "xmax" );
 }

HTMonitor::MEHTbinning HTMonitor::getHistoLSPSet ( const edm::ParameterSet & pset )

staticprivate

Definition at line 77 of file HTMonitor.cc.

References edm::ParameterSet::getParameter().

 {
   return HTMonitor::MEHTbinning{
     pset.getParameter<unsigned>("nbins"),
       0.,
       double(pset.getParameter<unsigned>("nbins"))
       };
 }

HTMonitor::MEHTbinning HTMonitor::getHistoPSet ( const edm::ParameterSet & pset )

staticprivate

Definition at line 68 of file HTMonitor.cc.

References edm::ParameterSet::getParameter().

 {
   return HTMonitor::MEHTbinning{
     pset.getParameter<unsigned>("nbins"),
       pset.getParameter<double>("xmin"),
       pset.getParameter<double>("xmax"),
       };
 }

void HTMonitor::setHTitle	(	HTME &	me,
		const std::string &	titleX,
		const std::string &	titleY
	)

protected

Definition at line 86 of file HTMonitor.cc.

References HTMonitor::HTME::denominator, HTMonitor::HTME::numerator, and MonitorElement::setAxisTitle().

Referenced by bookHistograms().

 {
   me.numerator->setAxisTitle(titleX,1);
   me.numerator->setAxisTitle(titleY,2);
   me.denominator->setAxisTitle(titleX,1);
   me.denominator->setAxisTitle(titleY,2);
 
 }

Member Data Documentation

HTME HTMonitor::deltaphij1j2ME_

private

Definition at line 112 of file HTMonitor.h.

Referenced by analyze(), and bookHistograms().

HTME HTMonitor::deltaphimetj1ME_

private

Definition at line 111 of file HTMonitor.h.

Referenced by analyze(), and bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> HTMonitor::den_genTriggerEventFlag_

private

Definition at line 115 of file HTMonitor.h.

Referenced by analyze(), and bookHistograms().

double HTMonitor::dEtaCut_

private

Definition at line 125 of file HTMonitor.h.

Referenced by analyze().

edm::InputTag HTMonitor::eleInputTag_

private

Definition at line 95 of file HTMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::GsfElectron,true> HTMonitor::eleSelection_

private

Definition at line 119 of file HTMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::GsfElectronCollection> HTMonitor::eleToken_

private

Definition at line 101 of file HTMonitor.h.

Referenced by analyze().

std::string HTMonitor::folderName_

private

Definition at line 89 of file HTMonitor.h.

Referenced by bookHistograms().

std::string HTMonitor::histoSuffix_

private

Definition at line 90 of file HTMonitor.h.

MEHTbinning HTMonitor::ht_binning_

private

Definition at line 106 of file HTMonitor.h.

Referenced by bookHistograms().

std::vector<double> HTMonitor::ht_variable_binning_

private

Definition at line 105 of file HTMonitor.h.

Referenced by bookHistograms().

HTME HTMonitor::htVsLS_

private

Definition at line 110 of file HTMonitor.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag HTMonitor::jetInputTag_

private

Definition at line 94 of file HTMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Jet,true > HTMonitor::jetSelection_

private

Definition at line 118 of file HTMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Jet,true > HTMonitor::jetSelection_HT_

private

Definition at line 121 of file HTMonitor.h.

Referenced by analyze().

edm::EDGetTokenT< reco::JetView > HTMonitor::jetToken_

private

Definition at line 100 of file HTMonitor.h.

Referenced by analyze().

MEHTbinning HTMonitor::ls_binning_

private

Definition at line 107 of file HTMonitor.h.

Referenced by bookHistograms().

double HTMonitor::MAXedge_PHI = 3.2

staticprivate

Definition at line 127 of file HTMonitor.h.

edm::InputTag HTMonitor::metInputTag_

private

Definition at line 93 of file HTMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::MET,true> HTMonitor::metSelection_

private

Definition at line 117 of file HTMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::PFMETCollection> HTMonitor::metToken_

private

Definition at line 99 of file HTMonitor.h.

Referenced by analyze().

edm::InputTag HTMonitor::muoInputTag_

private

Definition at line 96 of file HTMonitor.h.

Referenced by analyze().

StringCutObjectSelector<reco::Muon,true> HTMonitor::muoSelection_

private

Definition at line 120 of file HTMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::MuonCollection> HTMonitor::muoToken_

private

Definition at line 102 of file HTMonitor.h.

Referenced by analyze().

int HTMonitor::Nbin_PHI = 64

staticprivate

Definition at line 128 of file HTMonitor.h.

unsigned HTMonitor::nelectrons_

private

Definition at line 123 of file HTMonitor.h.

Referenced by analyze().

unsigned HTMonitor::njets_

private

Definition at line 122 of file HTMonitor.h.

Referenced by analyze().

unsigned HTMonitor::nmuons_

private

Definition at line 124 of file HTMonitor.h.

Referenced by analyze().

std::unique_ptr<GenericTriggerEventFlag> HTMonitor::num_genTriggerEventFlag_

private

Definition at line 114 of file HTMonitor.h.

Referenced by analyze(), and bookHistograms().

MEHTbinning HTMonitor::phi_binning_

staticprivate

Initial value:

{
    Nbin_PHI, -MAXedge_PHI, MAXedge_PHI
  }

Definition at line 129 of file HTMonitor.h.

Referenced by bookHistograms().

HTME HTMonitor::qME_variableBinning_

private

Definition at line 109 of file HTMonitor.h.

Referenced by analyze(), and bookHistograms().

quant HTMonitor::quantity_

private

Definition at line 91 of file HTMonitor.h.

Referenced by analyze(), bookHistograms(), and HTMonitor().

edm::InputTag HTMonitor::vtxInputTag_

private

Definition at line 97 of file HTMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::VertexCollection> HTMonitor::vtxToken_

private

Definition at line 103 of file HTMonitor.h.

Referenced by analyze().

std::vector<bool> HTMonitor::warningPrinted4token_

private

Definition at line 133 of file HTMonitor.h.

Referenced by analyze(), and HTMonitor().

Classes

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation