#include <ObjMonitor.h>

Inheritance diagram for ObjMonitor:

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

	~ObjMonitor () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::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

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_hmg_

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_

HMesonGammaDQM	hmgDQM_

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_

unsigned	nelectrons_

unsigned	njets_

unsigned	nmesons_

unsigned	nmuons_

unsigned	nphotons_

std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_

StringCutObjectSelector< reco::Photon, true >	phoSelection_

edm::EDGetTokenT< reco::PhotonCollection >	phoToken_

StringCutObjectSelector< reco::Track, true >	trkSelection_

edm::EDGetTokenT< reco::TrackCollection >	trkToken_

Detailed Description

Definition at line 54 of file ObjMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 19 of file ObjMonitor.cc.

References do_hmg_, do_ht_, do_jet_, do_met_, hmgDQM_, htDQM_, METDQM::initialise(), HTDQM::initialise(), JetDQM::initialise(), HMesonGammaDQM::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")     ) )
   , phoToken_             ( mayConsume<reco::PhotonCollection>     (iConfig.getParameter<edm::InputTag>("photons")   ) )
   , trkToken_             ( mayConsume<reco::TrackCollection>     (iConfig.getParameter<edm::InputTag>("tracks")   ) )
   , do_met_ (iConfig.getParameter<bool>("doMETHistos") )
   , do_jet_ (iConfig.getParameter<bool>("doJetHistos") )
   , do_ht_  (iConfig.getParameter<bool>("doHTHistos")  )
   , do_hmg_  (iConfig.getParameter<bool>("doHMesonGammaHistos")  )
   , 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") )
   , phoSelection_ ( iConfig.getParameter<std::string>("phoSelection") )
   , trkSelection_ ( iConfig.getParameter<std::string>("trkSelection") )
   , njets_      ( iConfig.getParameter<int>("njets" )      )
   , nelectrons_ ( iConfig.getParameter<int>("nelectrons" ) )
   , nmuons_     ( iConfig.getParameter<int>("nmuons" )     )
   , nphotons_   ( iConfig.getParameter<int>("nphotons")    )
   , nmesons_    ( iConfig.getParameter<int>("nmesons")    )
 {
   if (do_met_){
     metDQM_.initialise(iConfig);
   }
   if (do_jet_){
     jetDQM_.initialise(iConfig);
   }
   if (do_ht_ ){
     htDQM_.initialise(iConfig);
   }
   if (do_hmg_ ){
     hmgDQM_.initialise(iConfig);
   }
 }

ObjMonitor::~ObjMonitor ( )

overridedefault

Referenced by ObjMonitor().

Member Function Documentation

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

overrideprotected

Definition at line 82 of file ObjMonitor.cc.

References funct::abs(), objects.autophobj::absIso, reco::TrackBase::charge(), boostedElectronIsolation_cff::deltaR, den_genTriggerEventFlag_, do_hmg_, do_ht_, do_jet_, do_met_, reco::TrackBase::dz(), MillePedeFileConverter_cfg::e, nano_cff::electrons, eleSelection_, eleToken_, reco::TrackBase::eta(), METDQM::fillHistograms(), HTDQM::fillHistograms(), JetDQM::fillHistograms(), HMesonGammaDQM::fillHistograms(), edm::Event::getByToken(), hmgDQM_, htDQM_, htjetSelection_, mps_fire::i, edm::EventBase::id(), jetDQM_, jetId_, fwrapper::jets, jetSelection_, jetToken_, gen::k, looseJetId(), eostools::ls(), edm::EventID::luminosityBlock(), funct::m, RazorAnalyzer::met, metDQM_, metSelection_, metToken_, extraflags_cff::muons, muoSelection_, muoToken_, nelectrons_, njets_, nmesons_, nmuons_, nphotons_, num_genTriggerEventFlag_, phi, reco::LeafCandidate::phi(), reco::TrackBase::phi(), phoSelection_, phoToken_, muons_cff::photons, reco::LeafCandidate::pt(), reco::TrackBase::pt(), tightJetId(), trkSelection_, and trkToken_.

                                                                              {
 
   // 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 ( 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 ( jets.size() < njets_ ) return;
   
   edm::Handle<reco::GsfElectronCollection> eleHandle;
   iEvent.getByToken( eleToken_, eleHandle );
   std::vector<reco::GsfElectron> electrons;
   if ( eleHandle->size() < nelectrons_ ) return;
   for ( auto const & e : *eleHandle ) {
     if ( eleSelection_( e ) ) electrons.push_back(e);
   }
   if ( electrons.size() < nelectrons_ ) return;
   
   edm::Handle<reco::MuonCollection> muoHandle;
   iEvent.getByToken( muoToken_, muoHandle );
   if ( muoHandle->size() < nmuons_ ) return;
   std::vector<reco::Muon> muons;
   for ( auto const & m : *muoHandle ) {
     if ( muoSelection_( m ) ) muons.push_back(m);
   }
   if ( muons.size() < nmuons_ ) return;
 
   edm::Handle<reco::PhotonCollection> phoHandle;
   iEvent.getByToken( phoToken_, phoHandle );
   if ( phoHandle->size() < nphotons_ ) return;
   std::vector<reco::Photon> photons;
   for ( auto const & m : *phoHandle ) {
     if ( phoSelection_( m ) ) photons.push_back(m);
   }
   if ( photons.size() < nphotons_ ) return;
 
   std::vector<TLorentzVector> passedMesons;
   if (do_hmg_ ) {
     edm::Handle<reco::TrackCollection> trkHandle;
     iEvent.getByToken( trkToken_, trkHandle );
     // find isolated mesons (phi or rho)
     TLorentzVector t1,t2;
     float hadronMassHyp[2] = {0.1396,0.4937};  // pi or K mass
     float loMassLim[2] = {0.5,0.9};  // rho or phi mass 
     float hiMassLim[2] = {1.0,1.11};  // rho or phi mass 
     
     for (size_t i = 0; i < trkHandle->size(); ++i){
       const reco::Track trk1 = trkHandle->at(i);
       if (!trkSelection_(trk1)) continue;
       for (size_t j = i+1; j < trkHandle->size(); ++j){
     const reco::Track trk2 = trkHandle->at(j);
     if (!trkSelection_(trk2)) continue;
     if (trk1.charge() * trk2.charge() != -1) continue;
     
     for (unsigned hyp = 0; hyp < 2; ++hyp){
       
       t1.SetPtEtaPhiM(trk1.pt(),trk1.eta(),trk1.phi(),hadronMassHyp[hyp]);
       t2.SetPtEtaPhiM(trk2.pt(),trk2.eta(),trk2.phi(),hadronMassHyp[hyp]);
       TLorentzVector mesCand = t1 + t2;
       
       // cuts
       if (mesCand.M() < loMassLim[hyp] || mesCand.M() > hiMassLim[hyp]) continue;   //mass
       if (mesCand.Pt() < 35. || fabs(mesCand.Rapidity()) > 2.1 ) continue;   //pT eta
       
       // isolation
       float absIso = 0.;
       for (size_t k = 0; k < trkHandle->size(); ++k){
         if (k == i || k == j) continue;
         const reco::Track trkN = trkHandle->at(k);
         if (trkN.charge() == 0 || trkN.pt() < 0.5 ||
         ( trkN.dz() > 0.1 ) || 
         deltaR(trkN.eta(),trkN.phi(),mesCand.Eta(),mesCand.Phi()) > 0.5) continue;
         absIso += trkN.pt();
       }
       if (absIso/mesCand.Pt() > 0.2) continue;
       passedMesons.push_back(mesCand);
     }
       }
     }
     if ( passedMesons.size() < nmesons_ ) 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);
   if (do_hmg_ ) hmgDQM_.fillHistograms(photons,passedMesons,ls,passNumCond); 
 
 }

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

overrideprotected

Definition at line 63 of file ObjMonitor.cc.

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

 {  
 
   std::string currentFolder = folderName_ ;
   ibooker.setCurrentFolder(currentFolder);
 
   if (do_met_) metDQM_.bookHistograms(ibooker);
   if (do_jet_) jetDQM_.bookHistograms(ibooker);
   if (do_ht_ ) htDQM_.bookHistograms(ibooker);
   if (do_hmg_ ) hmgDQM_.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 242 of file ObjMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, HMesonGammaDQM::fillHmgDescription(), 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<edm::InputTag>( "photons",    edm::InputTag("gedPhotons") );
   desc.add<edm::InputTag>( "tracks",    edm::InputTag("generalTracks") );
   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<std::string>("phoSelection", "pt > 0");
   desc.add<std::string>("trkSelection", "pt > 0");
   desc.add<int>("njets",      0);
   desc.add<int>("nelectrons", 0);
   desc.add<int>("nmuons",     0);
   desc.add<int>("nphotons",     0);
   desc.add<int>("nmesons",     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<bool>("doHMesonGammaHistos", true);
   HMesonGammaDQM::fillHmgDescription(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 202 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 222 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);
   }
 }