Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #18971 from gomber/EXOTrigger_promptMonitoring92X
Exo trigger prompt monitoring92 x
- Loading branch information
Showing
7 changed files
with
615 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,324 @@ | ||
| #include "DQMOffline/Trigger/plugins/PhotonMonitor.h" | ||
|
|
||
| #include "FWCore/MessageLogger/interface/MessageLogger.h" | ||
|
|
||
| #include "DQM/TrackingMonitor/interface/GetLumi.h" | ||
|
|
||
| #include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h" | ||
|
|
||
| // ----------------------------- | ||
| // constructors and destructor | ||
| // ----------------------------- | ||
|
|
||
| PhotonMonitor::PhotonMonitor( const edm::ParameterSet& iConfig ) : | ||
| 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") ) ) | ||
| , photonToken_ ( mayConsume<reco::PhotonCollection> (iConfig.getParameter<edm::InputTag>("photons") ) ) | ||
| , photon_variable_binning_ ( iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("photonBinning") ) | ||
| , photon_binning_ ( getHistoPSet (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet> ("photonPSet") ) ) | ||
| , ls_binning_ ( getHistoLSPSet (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") ) | ||
| , photonSelection_ ( iConfig.getParameter<std::string>("photonSelection") ) | ||
| , njets_ ( iConfig.getParameter<unsigned int>("njets" ) ) | ||
| , nphotons_ ( iConfig.getParameter<unsigned int>("nphotons" ) ) | ||
| , nelectrons_ ( iConfig.getParameter<unsigned int>("nelectrons" ) ) | ||
| { | ||
|
|
||
| photonME_.numerator = nullptr; | ||
| photonME_.denominator = nullptr; | ||
| photonME_variableBinning_.numerator = nullptr; | ||
| photonME_variableBinning_.denominator = nullptr; | ||
| photonVsLS_.numerator = nullptr; | ||
| photonVsLS_.denominator = nullptr; | ||
| photonEtaME_.numerator = nullptr; | ||
| photonEtaME_.denominator = nullptr; | ||
| photonPhiME_.numerator = nullptr; | ||
| photonPhiME_.denominator = nullptr; | ||
| photonEtaPhiME_.numerator = nullptr; | ||
| photonEtaPhiME_.denominator = nullptr; | ||
| photonr9ME_.numerator = nullptr; | ||
| photonr9ME_.denominator = nullptr; | ||
| photonHoverEME_.numerator = nullptr; | ||
| photonHoverEME_.denominator = nullptr; | ||
|
|
||
|
|
||
|
|
||
|
|
||
| } | ||
|
|
||
|
|
||
| PhotonMonitor::~PhotonMonitor() | ||
| { | ||
|
|
||
| } | ||
|
|
||
| MEbinning PhotonMonitor::getHistoPSet(edm::ParameterSet const& pset) | ||
| { | ||
| return MEbinning{ | ||
| pset.getParameter<unsigned int>("nbins"), | ||
| pset.getParameter<double>("xmin"), | ||
| pset.getParameter<double>("xmax"), | ||
| }; | ||
| } | ||
|
|
||
| MEbinning PhotonMonitor::getHistoLSPSet(edm::ParameterSet const& pset) | ||
| { | ||
| return MEbinning{ | ||
| pset.getParameter<unsigned int>("nbins"), | ||
| 0., | ||
| double(pset.getParameter<unsigned int>("nbins")) | ||
| }; | ||
| } | ||
|
|
||
| void PhotonMonitor::setTitle(PhotonME& me, const std::string& titleX, const std::string& titleY) | ||
| { | ||
| me.numerator->setAxisTitle(titleX,1); | ||
| me.numerator->setAxisTitle(titleY,2); | ||
| me.denominator->setAxisTitle(titleX,1); | ||
| me.denominator->setAxisTitle(titleY,2); | ||
|
|
||
| } | ||
|
|
||
| void PhotonMonitor::bookME(DQMStore::IBooker &ibooker, PhotonME& me, const std::string& histname, const std::string& histtitle, unsigned int nbins, double min, double max) | ||
| { | ||
| me.numerator = ibooker.book1D(histname+"_numerator", histtitle+" (numerator)", nbins, min, max); | ||
| me.denominator = ibooker.book1D(histname+"_denominator", histtitle+" (denominator)", nbins, min, max); | ||
| } | ||
| void PhotonMonitor::bookME(DQMStore::IBooker &ibooker, PhotonME& me, const std::string& histname, const std::string& histtitle, const std::vector<double>& binning) | ||
| { | ||
| unsigned int nbins = binning.size()-1; | ||
| std::vector<float> fbinning(binning.begin(),binning.end()); | ||
| // float* arr = &fbinning[0]; | ||
| float* arr = fbinning.data(); | ||
| me.numerator = ibooker.book1D(histname+"_numerator", histtitle+" (numerator)", nbins, arr); | ||
| me.denominator = ibooker.book1D(histname+"_denominator", histtitle+" (denominator)", nbins, arr); | ||
| } | ||
| void PhotonMonitor::bookME(DQMStore::IBooker &ibooker, PhotonME& me, const std::string& histname, const std::string& histtitle, unsigned int nbinsX, double xmin, double xmax, double ymin, double ymax) | ||
| { | ||
| 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 PhotonMonitor::bookME(DQMStore::IBooker &ibooker, PhotonME& me, const std::string& histname, const std::string& histtitle, unsigned int nbinsX, double xmin, double xmax, unsigned int nbinsY, double ymin, double ymax) | ||
| { | ||
| 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 PhotonMonitor::bookME(DQMStore::IBooker &ibooker, PhotonME& me, const std::string& histname, const std::string& histtitle, const std::vector<double>& binningX, const std::vector<double>& binningY) | ||
| { | ||
| unsigned int nbinsX = binningX.size()-1; | ||
| std::vector<float> fbinningX(binningX.begin(),binningX.end()); | ||
| float* arrX = &fbinningX[0]; | ||
| unsigned 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 PhotonMonitor::bookHistograms(DQMStore::IBooker & ibooker, | ||
| edm::Run const & iRun, | ||
| edm::EventSetup const & iSetup) | ||
| { | ||
|
|
||
| std::string histname, histtitle; | ||
|
|
||
| std::string currentFolder = folderName_ ; | ||
| ibooker.setCurrentFolder(currentFolder.c_str()); | ||
|
|
||
| histname = "photon_pt"; histtitle = "photon PT"; | ||
| bookME(ibooker,photonME_,histname,histtitle,photon_binning_.nbins,photon_binning_.xmin, photon_binning_.xmax); | ||
| setTitle(photonME_,"Photon pT [GeV]","events / [GeV]"); | ||
|
|
||
| histname = "photon_pt_variable"; histtitle = "photon PT"; | ||
| bookME(ibooker,photonME_variableBinning_,histname,histtitle,photon_variable_binning_); | ||
| setTitle(photonME_variableBinning_,"Photon pT [GeV]","events / [GeV]"); | ||
|
|
||
| histname = "photonVsLS"; histtitle = "photon pt vs LS"; | ||
| bookME(ibooker,photonVsLS_,histname,histtitle,ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax,photon_binning_.xmin, photon_binning_.xmax); | ||
| setTitle(photonVsLS_,"LS","Photon pT [GeV]"); | ||
|
|
||
| histname = "photon_phi"; histtitle = "Photon phi"; | ||
| bookME(ibooker,photonPhiME_,histname,histtitle, phi_binning_1.nbins, phi_binning_1.xmin, phi_binning_1.xmax); | ||
| setTitle(photonPhiME_,"Photon #phi","events / 0.1 rad"); | ||
|
|
||
|
|
||
| histname = "photon_eta"; histtitle = "Photon eta"; | ||
| bookME(ibooker,photonEtaME_,histname,histtitle, eta_binning_.nbins, eta_binning_.xmin,eta_binning_.xmax); | ||
| setTitle(photonEtaME_,"Photon #eta","events"); | ||
|
|
||
|
|
||
| histname = "photon_r9"; histtitle = "Photon r9"; | ||
| bookME(ibooker,photonr9ME_,histname,histtitle, r9_binning_.nbins, r9_binning_.xmin, r9_binning_.xmax); | ||
| setTitle(photonr9ME_,"Photon r9","events"); | ||
|
|
||
|
|
||
| histname = "photon_hoE"; histtitle = "Photon hoverE"; | ||
| bookME(ibooker,photonHoverEME_,histname,histtitle, hoe_binning_.nbins, hoe_binning_.xmin, hoe_binning_.xmax); | ||
| setTitle(photonHoverEME_,"Photon hoE","events"); | ||
|
|
||
| histname = "photon_etaphi"; histtitle = "Photon eta-phi"; | ||
| bookME(ibooker,photonEtaPhiME_,histname,histtitle, eta_binning_.nbins, eta_binning_.xmin, eta_binning_.xmax,phi_binning_1.nbins, phi_binning_1.xmin, phi_binning_1.xmax); | ||
| setTitle(photonEtaPhiME_,"#eta","#phi"); | ||
|
|
||
| // 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 ); | ||
|
|
||
| } | ||
|
|
||
| #include "FWCore/Framework/interface/ESHandle.h" | ||
| #include "FWCore/Framework/interface/EventSetup.h" | ||
| #include "DataFormats/TrackerCommon/interface/TrackerTopology.h" | ||
| #include "Geometry/Records/interface/TrackerTopologyRcd.h" | ||
| void PhotonMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) { | ||
|
|
||
| // 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; | ||
| jets.clear(); | ||
| if ( jetHandle->size() < njets_ ) return; | ||
| for ( auto const & j : *jetHandle ) { | ||
| if ( jetSelection_( j ) ) jets.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::PhotonCollection> photonHandle; | ||
| iEvent.getByToken( photonToken_, photonHandle ); | ||
| std::vector<reco::Photon> photons; | ||
| photons.clear(); | ||
| if ( photonHandle->size() < nphotons_ ) return; | ||
| for ( auto const & p : *photonHandle ) { | ||
| if ( photonSelection_( p ) ) photons.push_back(p); | ||
| } | ||
| if ( photons.size() < nphotons_ ) return; | ||
|
|
||
|
|
||
| // filling histograms (denominator) | ||
| int ls = iEvent.id().luminosityBlock(); | ||
| if(!(photons.empty())) | ||
|
|
||
| { | ||
| photonME_.denominator -> Fill(photons[0].pt()); | ||
| photonME_variableBinning_.denominator -> Fill(photons[0].pt()); | ||
| photonPhiME_.denominator->Fill(photons[0].phi()); | ||
| photonEtaME_.denominator->Fill(photons[0].eta()); | ||
| photonVsLS_.denominator -> Fill(ls, photons[0].pt()); | ||
| photonEtaPhiME_.denominator -> Fill(photons[0].eta(), photons[0].phi()); | ||
| photonr9ME_.denominator->Fill(photons[0].r9()); | ||
| photonHoverEME_.denominator->Fill(photons[0].hadTowOverEm()); | ||
| } | ||
|
|
||
| // applying selection for numerator | ||
| if (num_genTriggerEventFlag_->on() && ! num_genTriggerEventFlag_->accept( iEvent, iSetup) ) return; | ||
|
|
||
| // filling histograms (num_genTriggerEventFlag_) | ||
| if(!(photons.empty())) | ||
| { | ||
| photonME_.numerator -> Fill(photons[0].pt()); | ||
| photonME_variableBinning_.numerator -> Fill(photons[0].pt()); | ||
| photonPhiME_.numerator->Fill(photons[0].phi()); | ||
| photonEtaME_.numerator->Fill(photons[0].eta()); | ||
| photonVsLS_.numerator -> Fill(ls, photons[0].pt()); | ||
| photonEtaPhiME_.numerator -> Fill(photons[0].eta(), photons[0].phi()); | ||
| photonr9ME_.numerator->Fill(photons[0].r9()); | ||
| photonHoverEME_.numerator->Fill(photons[0].hadTowOverEm()); | ||
|
|
||
|
|
||
| } | ||
|
|
||
| } | ||
|
|
||
| void PhotonMonitor::fillHistoPSetDescription(edm::ParameterSetDescription & pset) | ||
| { | ||
| pset.add<unsigned int> ( "nbins"); | ||
| pset.add<double>( "xmin" ); | ||
| pset.add<double>( "xmax" ); | ||
| } | ||
|
|
||
| void PhotonMonitor::fillHistoLSPSetDescription(edm::ParameterSetDescription & pset) | ||
| { | ||
| pset.add<unsigned int> ( "nbins", 2500); | ||
| } | ||
|
|
||
| void PhotonMonitor::fillDescriptions(edm::ConfigurationDescriptions & descriptions) | ||
| { | ||
| edm::ParameterSetDescription desc; | ||
| desc.add<std::string> ( "FolderName", "HLT/Photon" ); | ||
| 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>( "photons",edm::InputTag("gedPhotons") ); | ||
| 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>("photonSelection", "pt > 145 && eta<1.4442 && hadTowOverEm<0.0597 && full5x5_sigmaIetaIeta()<0.01031 && chargedHadronIso<1.295"); | ||
| //desc.add<std::string>("photonSelection", "pt > 145"); | ||
| desc.add<unsigned int>("njets", 0); | ||
| desc.add<unsigned int>("nelectrons", 0); | ||
| desc.add<unsigned int>("nphotons", 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); | ||
|
|
||
| edm::ParameterSetDescription histoPSet; | ||
| edm::ParameterSetDescription metPSet; | ||
| fillHistoPSetDescription(metPSet); | ||
| histoPSet.add<edm::ParameterSetDescription>("photonPSet", metPSet); | ||
| 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.,1000.}; | ||
| histoPSet.add<std::vector<double> >("photonBinning", bins); | ||
|
|
||
| edm::ParameterSetDescription lsPSet; | ||
| fillHistoLSPSetDescription(lsPSet); | ||
| histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet); | ||
|
|
||
| desc.add<edm::ParameterSetDescription>("histoPSet",histoPSet); | ||
|
|
||
| descriptions.add("photonMonitoring", desc); | ||
| } | ||
|
|
||
| // Define this as a plug-in | ||
| #include "FWCore/Framework/interface/MakerMacros.h" | ||
| DEFINE_FWK_MODULE(PhotonMonitor); |
Oops, something went wrong.