#include <ObjMonitor.h>

Inheritance diagram for ObjMonitor:

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

	~ObjMonitor () override

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

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

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

	~EDAnalyzerBase () 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

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

virtual	~EDConsumerBase () noexcept(false)

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

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

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

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)

Private Member Functions
bool	looseJetId (const double &abseta, const double &NHF, const double &NEMF, const double &CHF, const double &CEMF, const unsigned &NumNeutralParticles, const unsigned &CHM)

bool	tightJetId (const double &abseta, const double &NHF, const double &NEMF, const double &CHF, const double &CEMF, const unsigned &NumNeutralParticles, const unsigned &CHM)

Private Attributes
std::unique_ptr< GenericTriggerEventFlag >	den_genTriggerEventFlag_

bool	do_ht_

bool	do_jet_

bool	do_met_

StringCutObjectSelector< reco::GsfElectron, true >	eleSelection_

edm::EDGetTokenT< reco::GsfElectronCollection >	eleToken_

std::string	folderName_

std::string	histoSuffix_

HTDQM	htDQM_

StringCutObjectSelector< reco::PFJet, true >	htjetSelection_

JetDQM	jetDQM_

std::string	jetId_

StringCutObjectSelector< reco::PFJet, true >	jetSelection_

edm::EDGetTokenT< reco::PFJetCollection >	jetToken_

METDQM	metDQM_

StringCutObjectSelector< reco::MET, true >	metSelection_

edm::EDGetTokenT< reco::PFMETCollection >	metToken_

StringCutObjectSelector< reco::Muon, true >	muoSelection_

edm::EDGetTokenT< reco::MuonCollection >	muoToken_

int	nelectrons_

int	njets_

int	nmuons_

std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_

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

Detailed Description

Definition at line 53 of file ObjMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 18 of file ObjMonitor.cc.

References do_ht_, do_jet_, do_met_, htDQM_, METDQM::initialise(), HTDQM::initialise(), JetDQM::initialise(), jetDQM_, metDQM_, and ~ObjMonitor().

                                                        : 
   folderName_             ( iConfig.getParameter<std::string>("FolderName") )
   , metToken_             ( consumes<reco::PFMETCollection>      (iConfig.getParameter<edm::InputTag>("met")       ) )   
   , jetToken_             ( mayConsume<reco::PFJetCollection>      (iConfig.getParameter<edm::InputTag>("jets")      ) )   
   , eleToken_             ( mayConsume<reco::GsfElectronCollection>(iConfig.getParameter<edm::InputTag>("electrons") ) )   
   , muoToken_             ( mayConsume<reco::MuonCollection>       (iConfig.getParameter<edm::InputTag>("muons")     ) )
   , do_met_ (iConfig.getParameter<bool>("doMETHistos") )
   , do_jet_ (iConfig.getParameter<bool>("doJetHistos") )
   , do_ht_  (iConfig.getParameter<bool>("doHTHistos")  )
   , num_genTriggerEventFlag_(std::make_unique<GenericTriggerEventFlag>(iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"),consumesCollector(), *this))
   , den_genTriggerEventFlag_(std::make_unique<GenericTriggerEventFlag>(iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"),consumesCollector(), *this))
   , metSelection_ ( iConfig.getParameter<std::string>("metSelection") )
   , jetSelection_ ( iConfig.getParameter<std::string>("jetSelection") )
   , jetId_        ( iConfig.getParameter<std::string>("jetId")        )
   , htjetSelection_ ( iConfig.getParameter<std::string>("htjetSelection"))
   , eleSelection_ ( iConfig.getParameter<std::string>("eleSelection") )
   , muoSelection_ ( iConfig.getParameter<std::string>("muoSelection") )
   , njets_      ( iConfig.getParameter<int>("njets" )      )
   , nelectrons_ ( iConfig.getParameter<int>("nelectrons" ) )
   , nmuons_     ( iConfig.getParameter<int>("nmuons" )     )
 {
   if (do_met_){
     metDQM_.initialise(iConfig);
   }
   if (do_jet_){
     jetDQM_.initialise(iConfig);
   }
   if (do_ht_ ){
     htDQM_.initialise(iConfig);
   }
 }

ObjMonitor::~ObjMonitor ( )

overridedefault

Referenced by ObjMonitor().

Member Function Documentation

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

overrideprotectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 70 of file ObjMonitor.cc.

References funct::abs(), den_genTriggerEventFlag_, do_ht_, do_jet_, do_met_, MillePedeFileConverter_cfg::e, electrons_cff::electrons, eleSelection_, eleToken_, METDQM::fillHistograms(), HTDQM::fillHistograms(), JetDQM::fillHistograms(), edm::Event::getByToken(), htDQM_, htjetSelection_, edm::EventBase::id(), jetDQM_, jetId_, fwrapper::jets, jetSelection_, jetToken_, looseJetId(), eostools::ls(), edm::EventID::luminosityBlock(), funct::m, RazorAnalyzer::met, metDQM_, metSelection_, metToken_, extraflags_cff::muons, muoSelection_, muoToken_, nelectrons_, njets_, nmuons_, num_genTriggerEventFlag_, phi, reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), and tightJetId().

                                                                              {
 
   // 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 );
   reco::PFMET pfmet = metHandle->front();
   if ( ! metSelection_( pfmet ) ) return;
   
   float met = pfmet.pt();
   float phi = pfmet.phi();
 
   edm::Handle<reco::PFJetCollection> jetHandle;
   iEvent.getByToken( jetToken_, jetHandle );
   std::vector<reco::PFJet> jets;
   std::vector<reco::PFJet> htjets;
   if ( int(jetHandle->size()) < njets_ ) return;
   for ( auto const & j : *jetHandle ) {
     if ( jetSelection_( j ) ) {
       if (jetId_=="loose" || jetId_ =="tight"){
     double abseta = abs(j.eta());
     double NHF  = j.neutralHadronEnergyFraction();
     double NEMF = j.neutralEmEnergyFraction();
     double CHF  = j.chargedHadronEnergyFraction();
     double CEMF = j.chargedEmEnergyFraction();
     unsigned NumNeutralParticles =j.neutralMultiplicity();
     unsigned CHM      = j.chargedMultiplicity();
     bool passId = (jetId_=="loose" && looseJetId(abseta,NHF,NEMF,CHF,CEMF,NumNeutralParticles,CHM)) || (jetId_=="tight" && tightJetId(abseta,NHF,NEMF,CHF,CEMF,NumNeutralParticles,CHM));
     if (passId) jets.push_back(j);
       }
       else jets.push_back(j);
     }
     if ( htjetSelection_( j ) ) htjets.push_back(j);
   }
   if ( int(jets.size()) < njets_ ) return;
   
   edm::Handle<reco::GsfElectronCollection> eleHandle;
   iEvent.getByToken( eleToken_, eleHandle );
   std::vector<reco::GsfElectron> electrons;
   if ( int(eleHandle->size()) < nelectrons_ ) return;
   for ( auto const & e : *eleHandle ) {
     if ( eleSelection_( e ) ) electrons.push_back(e);
   }
   if ( int(electrons.size()) < nelectrons_ ) return;
   
   edm::Handle<reco::MuonCollection> muoHandle;
   iEvent.getByToken( muoToken_, muoHandle );
   if ( int(muoHandle->size()) < nmuons_ ) return;
   std::vector<reco::Muon> muons;
   for ( auto const & m : *muoHandle ) {
     if ( muoSelection_( m ) ) muons.push_back(m);
   }
   if ( int(muons.size()) < nmuons_ ) return;
 
 
   bool passNumCond = num_genTriggerEventFlag_->off() || num_genTriggerEventFlag_->accept( iEvent, iSetup);
   int ls = iEvent.id().luminosityBlock();
 
   if (do_met_) metDQM_.fillHistograms(met,phi,ls,passNumCond);
   if (do_jet_) jetDQM_.fillHistograms(jets,pfmet,ls,passNumCond);
   if (do_ht_ ) htDQM_.fillHistograms(htjets,met,ls,passNumCond);
 
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 52 of file ObjMonitor.cc.

References METDQM::bookHistograms(), HTDQM::bookHistograms(), JetDQM::bookHistograms(), den_genTriggerEventFlag_, do_ht_, do_jet_, do_met_, folderName_, htDQM_, jetDQM_, metDQM_, num_genTriggerEventFlag_, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {  
 
   std::string currentFolder = folderName_ ;
   ibooker.setCurrentFolder(currentFolder.c_str());
 
   if (do_met_) metDQM_.bookHistograms(ibooker);
   if (do_jet_) jetDQM_.bookHistograms(ibooker);
   if (do_ht_ ) htDQM_.bookHistograms(ibooker);
 
   // 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 ObjMonitor::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 175 of file ObjMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, HTDQM::fillHtDescription(), JetDQM::fillJetDescription(), METDQM::fillMetDescription(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   edm::ParameterSetDescription desc;
   desc.add<std::string>  ( "FolderName", "HLT/OBJ" );
 
   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<std::string>("metSelection", "pt > 0");
   desc.add<std::string>("jetSelection", "pt > 0");
   desc.add<std::string>("jetId", "");
   desc.add<std::string>("htjetSelection", "pt > 30");
   desc.add<std::string>("eleSelection", "pt > 0");
   desc.add<std::string>("muoSelection", "pt > 0");
   desc.add<int>("njets",      0);
   desc.add<int>("nelectrons", 0);
   desc.add<int>("nmuons",     0);
 
   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);
 
   desc.add<bool>("doMETHistos", true);
   edm::ParameterSetDescription histoPSet;
   METDQM::fillMetDescription(histoPSet);
   desc.add<bool>("doJetHistos", true);
   JetDQM::fillJetDescription(histoPSet);
   desc.add<bool>("doHTHistos", true);
   HTDQM::fillHtDescription(histoPSet);
 
   desc.add<edm::ParameterSetDescription>("histoPSet",histoPSet);
 
   descriptions.add("objMonitoring", desc);
 }

bool ObjMonitor::looseJetId	(	const double &	abseta,
		const double &	NHF,
		const double &	NEMF,
		const double &	CHF,
		const double &	CEMF,
		const unsigned &	NumNeutralParticles,
		const unsigned &	CHM
	)

private

Definition at line 135 of file ObjMonitor.cc.

Referenced by analyze().

 {
   if (abseta<=2.7){
     unsigned NumConst = CHM+NumNeutralParticles;
 
     return ((NumConst>1 && NHF<0.99 && NEMF<0.99) && ((abseta<=2.4 && CHF>0 && CHM>0 && CEMF<0.99) || abseta>2.4));
   }
   else if (abseta<=3){
     return (NumNeutralParticles>2 && NEMF>0.01 && NHF<0.98);
   }
   else {
     return NumNeutralParticles>10 && NEMF<0.90;
   }
 }

bool ObjMonitor::tightJetId	(	const double &	abseta,
		const double &	NHF,
		const double &	NEMF,
		const double &	CHF,
		const double &	CEMF,
		const unsigned &	NumNeutralParticles,
		const unsigned &	CHM
	)

private

Definition at line 155 of file ObjMonitor.cc.

Referenced by analyze().

 {
   if (abseta<=2.7){
     unsigned NumConst = CHM+NumNeutralParticles;
     return (NumConst>1 && NHF<0.90 && NEMF<0.90 ) && ((abseta<=2.4 && CHF>0 && CHM>0 && CEMF<0.99) || abseta>2.4);
   }
   else if (abseta<=3){
     return (NHF<0.98 && NEMF>0.01 && NumNeutralParticles>2);
   }
   else {
     return (NEMF<0.90 && NumNeutralParticles>10);
   }
 }