#include <PhotonMonitor.h>

Inheritance diagram for PhotonMonitor:

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

	~PhotonMonitor () 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 () 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 &, PhotonME &me, const std::string &histname, const std::string &histtitle, unsigned int nbins, double xmin, double xmax)

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

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

void	bookME (DQMStore::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)

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

void	setTitle (PhotonME &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 MEbinning	getHistoLSPSet (edm::ParameterSet const &pset)

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

Private Attributes
std::unique_ptr< GenericTriggerEventFlag >	den_genTriggerEventFlag_

std::vector< double >	diphoton_mass_binning_

PhotonME	diphotonMassME_

StringCutObjectSelector< reco::GsfElectron, true >	eleSelection_

edm::EDGetTokenT< reco::GsfElectronCollection >	eleToken_

const MEbinning	eta_binning_

std::string	folderName_

std::string	histoSuffix_

const MEbinning	hoe_binning_

StringCutObjectSelector< reco::PFJet, true >	jetSelection_

edm::EDGetTokenT< reco::PFJetCollection >	jetToken_

MEbinning	ls_binning_

double	MAX_ETA = 1.4442

double	MAX_hoe = 0.02

double	MAX_PHI1 = 3.2

double	MAX_r9 = 1

StringCutObjectSelector< reco::MET, true >	metSelection_

edm::EDGetTokenT< reco::PFMETCollection >	metToken_

double	MIN_hoe = 0

double	MIN_r9 = 0

unsigned int	N_ETA = 34

unsigned int	N_PHI1 = 64

unsigned int	N_r9 = 50

unsigned int	nelectrons_

unsigned int	njets_

unsigned int	nphotons_

std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_

const MEbinning	phi_binning_1

MEbinning	photon_binning_

std::vector< double >	photon_variable_binning_

PhotonME	photonEtaME_

PhotonME	photonEtaPhiME_

PhotonME	photonHoverEME_

PhotonME	photonME_

PhotonME	photonME_variableBinning_

PhotonME	photonPhiME_

PhotonME	photonr9ME_

StringCutObjectSelector< reco::Photon, true >	photonSelection_

edm::EDGetTokenT< reco::PhotonCollection >	photonToken_

PhotonME	photonVsLS_

const MEbinning	r9_binning_

PhotonME	subphotonEtaME_

PhotonME	subphotonEtaPhiME_

PhotonME	subphotonHoverEME_

PhotonME	subphotonME_

PhotonME	subphotonME_variableBinning_

PhotonME	subphotonPhiME_

PhotonME	subphotonr9ME_

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 55 of file PhotonMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 13 of file PhotonMonitor.cc.

References PhotonME::denominator, diphotonMassME_, PhotonME::numerator, photonEtaME_, photonEtaPhiME_, photonHoverEME_, photonME_, photonME_variableBinning_, photonPhiME_, photonr9ME_, photonVsLS_, subphotonEtaME_, subphotonEtaPhiME_, subphotonHoverEME_, subphotonME_, subphotonME_variableBinning_, subphotonPhiME_, subphotonr9ME_, and ~PhotonMonitor().

                                                              : 
   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") )
   , diphoton_mass_binning_ ( iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("massBinning") )
   , photon_binning_          ( getHistoPSet   (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>   ("photonPSet")    ) )
   , 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") )
   , 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;
 
   diphotonMassME_.numerator   = nullptr;
   diphotonMassME_.denominator = nullptr;
 
   subphotonME_.numerator   = nullptr;
   subphotonME_.denominator = nullptr;
   subphotonME_variableBinning_.numerator   = nullptr;
   subphotonME_variableBinning_.denominator = nullptr;
   subphotonEtaME_.numerator   = nullptr;
   subphotonEtaME_.denominator = nullptr;
   subphotonPhiME_.numerator   = nullptr;
   subphotonPhiME_.denominator = nullptr;
   subphotonEtaPhiME_.numerator   = nullptr;       
   subphotonEtaPhiME_.denominator = nullptr;
   subphotonr9ME_.numerator   = nullptr;       
   subphotonr9ME_.denominator = nullptr;
   subphotonHoverEME_.numerator   = nullptr;       
   subphotonHoverEME_.denominator = nullptr;
   
 }

PhotonMonitor::~PhotonMonitor ( )

overridedefault

Referenced by PhotonMonitor().

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 221 of file PhotonMonitor.cc.

References funct::cos(), den_genTriggerEventFlag_, PhotonME::denominator, diphotonMassME_, MillePedeFileConverter_cfg::e, electrons_cff::electrons, eleSelection_, eleToken_, PVValHelper::eta, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::Event::getByToken(), edm::EventBase::id(), fwrapper::jets, jetSelection_, jetToken_, eostools::ls(), edm::EventID::luminosityBlock(), metSelection_, metToken_, nelectrons_, njets_, nphotons_, num_genTriggerEventFlag_, PhotonME::numerator, AlCaHLTBitMon_ParallelJobs::p, phi, photonEtaME_, photonEtaPhiME_, photonHoverEME_, photonME_, photonME_variableBinning_, photonPhiME_, photonr9ME_, nano_cff::photons, photonSelection_, photonToken_, photonVsLS_, EnergyCorrector::pt, electrons_cff::r9, mathSSE::sqrt(), subphotonEtaME_, subphotonEtaPhiME_, subphotonHoverEME_, subphotonME_, subphotonPhiME_, and subphotonr9ME_.

                                                                                 {
 
   // 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());
     }
   
   if(nphotons_>1) 
     //filling diphoton histograms
     {
       subphotonME_.denominator -> Fill(photons[1].pt());
       subphotonEtaME_.denominator -> Fill(photons[1].eta());
       subphotonPhiME_.denominator -> Fill(photons[1].phi());
       subphotonEtaPhiME_.denominator -> Fill(photons[1].eta(), photons[1].phi());
       subphotonr9ME_.denominator->Fill(photons[1].r9());
       subphotonHoverEME_.denominator->Fill(photons[1].hadTowOverEm());
       diphotonMassME_.denominator -> Fill(sqrt(2*photons[0].pt()*photons[1].pt()*(cosh(photons[0].eta()-photons[1].eta())-cos(photons[0].phi()-photons[1].phi()))));
     }
 
   // 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());
     }
   if(nphotons_>1) 
     //filling diphoton histograms
     {
       subphotonME_.numerator -> Fill(photons[1].pt());
       subphotonEtaME_.numerator -> Fill(photons[1].eta());
       subphotonPhiME_.numerator -> Fill(photons[1].phi());
       subphotonEtaPhiME_.numerator -> Fill(photons[1].eta(), photons[1].phi());
       subphotonr9ME_.numerator->Fill(photons[1].r9());
       subphotonHoverEME_.numerator->Fill(photons[1].hadTowOverEm());
       diphotonMassME_.numerator -> Fill(sqrt(2*photons[0].pt()*photons[1].pt()*(cosh(photons[0].eta()-photons[1].eta())-cos(photons[0].phi()-photons[1].phi()))));
     }
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 136 of file PhotonMonitor.cc.

References bookME(), den_genTriggerEventFlag_, diphoton_mass_binning_, diphotonMassME_, eta_binning_, folderName_, hoe_binning_, ls_binning_, nphotons_, num_genTriggerEventFlag_, phi_binning_1, photon_binning_, photon_variable_binning_, photonEtaME_, photonEtaPhiME_, photonHoverEME_, photonME_, photonME_variableBinning_, photonPhiME_, photonr9ME_, photonVsLS_, r9_binning_, DQMStore::IBooker::setCurrentFolder(), setTitle(), AlCaHLTBitMon_QueryRunRegistry::string, subphotonEtaME_, subphotonEtaPhiME_, subphotonHoverEME_, subphotonME_, subphotonPhiME_, and subphotonr9ME_.

 {  
   
   std::string histname, histtitle;
 
   std::string currentFolder = folderName_ ;
   ibooker.setCurrentFolder(currentFolder);
 
   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"); 
 
   //for diphotons
   if(nphotons_>1)
     {
       histname = "diphoton_mass"; histtitle = "Diphoton mass";
       bookME(ibooker,diphotonMassME_,histname,histtitle, diphoton_mass_binning_);
       setTitle(diphotonMassME_,"Diphoton mass","events / 0.1");
       
       histname = "subphoton_pt"; histtitle = "subphoton PT";
       bookME(ibooker,subphotonME_,histname,histtitle,photon_binning_.nbins,photon_binning_.xmin, photon_binning_.xmax);
       setTitle(subphotonME_,"subPhoton pT [GeV]","events / [GeV]");
       
       histname = "subphoton_eta"; histtitle = "subPhoton eta";
       bookME(ibooker,subphotonEtaME_,histname,histtitle, eta_binning_.nbins, eta_binning_.xmin,eta_binning_.xmax);
       setTitle(subphotonEtaME_,"subPhoton #eta","events / 0.1");
       
       histname = "subphoton_phi"; histtitle = "subPhoton phi";
       bookME(ibooker,subphotonPhiME_,histname,histtitle, phi_binning_1.nbins, phi_binning_1.xmin, phi_binning_1.xmax);
       setTitle(subphotonPhiME_,"subPhoton #phi","events / 0.1 rad");
       
       histname = "subphoton_r9"; histtitle = "subPhoton r9";
       bookME(ibooker,subphotonr9ME_,histname,histtitle, r9_binning_.nbins, r9_binning_.xmin, r9_binning_.xmax);
       setTitle(subphotonr9ME_,"subPhoton r9","events");
       
       histname = "subphoton_hoE"; histtitle = "subPhoton hoverE";
       bookME(ibooker,subphotonHoverEME_,histname,histtitle, hoe_binning_.nbins, hoe_binning_.xmin, hoe_binning_.xmax);
       setTitle(subphotonHoverEME_,"subPhoton hoE","events");
       
       histname = "subphoton_etaphi"; histtitle = "subPhoton eta-phi"; 
       bookME(ibooker,subphotonEtaPhiME_,histname,histtitle, eta_binning_.nbins, eta_binning_.xmin, eta_binning_.xmax,phi_binning_1.nbins, phi_binning_1.xmin, phi_binning_1.xmax);
       setTitle(subphotonEtaPhiME_,"#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 );
 
 }

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

protected

Definition at line 99 of file PhotonMonitor.cc.

References DQMStore::IBooker::book1D(), PhotonME::denominator, and PhotonME::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 PhotonMonitor::bookME	(	DQMStore::IBooker &	ibooker,
		PhotonME &	me,
		const std::string &	histname,
		const std::string &	histtitle,
		const std::vector< double > &	binningX
	)

protected

Definition at line 104 of file PhotonMonitor.cc.

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

 {
   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
	)

protected

Definition at line 113 of file PhotonMonitor.cc.

References DQMStore::IBooker::bookProfile(), PhotonME::denominator, and PhotonME::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 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
	)

protected

Definition at line 118 of file PhotonMonitor.cc.

References DQMStore::IBooker::book2D(), PhotonME::denominator, and PhotonME::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 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
	)

protected

Definition at line 123 of file PhotonMonitor.cc.

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

 {
   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::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 334 of file PhotonMonitor.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/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);
   std::vector<double> massbins = {90.,91.,92.,93.,94.,95.,96.,97.,98.,99.,100.,101.,102.,103.,104.,105.,106.,107.,108.,109.,110.,115.,120.,130.,150.,200.};
   histoPSet.add<std::vector<double> >("massBinning", massbins);
   edm::ParameterSetDescription lsPSet;
   fillHistoLSPSetDescription(lsPSet);
   histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
 
   desc.add<edm::ParameterSetDescription>("histoPSet",histoPSet);
 
   descriptions.add("photonMonitoring", desc);
 }

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

static

Definition at line 327 of file PhotonMonitor.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 PhotonMonitor::fillHistoPSetDescription ( edm::ParameterSetDescription & pset )

static

Definition at line 320 of file PhotonMonitor.cc.

References edm::ParameterSetDescription::add().

Referenced by fillDescriptions().

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

MEbinning PhotonMonitor::getHistoLSPSet ( edm::ParameterSet const & pset )

staticprivate

Definition at line 81 of file PhotonMonitor.cc.

References edm::ParameterSet::getParameter().

 {
   return MEbinning{
     pset.getParameter<unsigned int>("nbins"),
       0.,
       double(pset.getParameter<unsigned int>("nbins"))
       };
 }

MEbinning PhotonMonitor::getHistoPSet ( edm::ParameterSet const & pset )

staticprivate

Definition at line 72 of file PhotonMonitor.cc.

References edm::ParameterSet::getParameter().

 {
   return MEbinning{
     pset.getParameter<unsigned int>("nbins"),
       pset.getParameter<double>("xmin"),
       pset.getParameter<double>("xmax"),
       };
 }

void PhotonMonitor::setTitle	(	PhotonME &	me,
		const std::string &	titleX,
		const std::string &	titleY
	)

protected

Definition at line 90 of file PhotonMonitor.cc.

References PhotonME::denominator, PhotonME::numerator, and MonitorElement::setAxisTitle().

Referenced by bookHistograms(), Vispa.Gui.PortWidget.PortWidget::setName(), and Vispa.Views.LineDecayView.LineDecayContainer::setPxlObject().

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

Member Data Documentation

std::unique_ptr<GenericTriggerEventFlag> PhotonMonitor::den_genTriggerEventFlag_

private

Definition at line 152 of file PhotonMonitor.h.

Referenced by analyze(), and bookHistograms().

std::vector<double> PhotonMonitor::diphoton_mass_binning_

private

Definition at line 91 of file PhotonMonitor.h.

Referenced by bookHistograms().

PhotonME PhotonMonitor::diphotonMassME_

private

Definition at line 103 of file PhotonMonitor.h.

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

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

private

Definition at line 156 of file PhotonMonitor.h.

Referenced by analyze().

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

private

Definition at line 84 of file PhotonMonitor.h.

Referenced by analyze().

const MEbinning PhotonMonitor::eta_binning_

private

Initial value:

{
    N_ETA, -MAX_ETA, MAX_ETA
      }

Definition at line 123 of file PhotonMonitor.h.

Referenced by bookHistograms().

std::string PhotonMonitor::folderName_

private

Definition at line 79 of file PhotonMonitor.h.

Referenced by bookHistograms().

std::string PhotonMonitor::histoSuffix_

private

Definition at line 80 of file PhotonMonitor.h.

const MEbinning PhotonMonitor::hoe_binning_

private

Initial value:

{
    N_r9, MIN_hoe, MAX_hoe
      }

Definition at line 140 of file PhotonMonitor.h.

Referenced by bookHistograms().

StringCutObjectSelector<reco::PFJet,true > PhotonMonitor::jetSelection_

private

Definition at line 155 of file PhotonMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::PFJetCollection> PhotonMonitor::jetToken_

private

Definition at line 83 of file PhotonMonitor.h.

Referenced by analyze().

MEbinning PhotonMonitor::ls_binning_

private

Definition at line 94 of file PhotonMonitor.h.

Referenced by bookHistograms().

double PhotonMonitor::MAX_ETA = 1.4442

private

Definition at line 121 of file PhotonMonitor.h.

double PhotonMonitor::MAX_hoe = 0.02

private

Definition at line 138 of file PhotonMonitor.h.

double PhotonMonitor::MAX_PHI1 = 3.2

private

Definition at line 114 of file PhotonMonitor.h.

double PhotonMonitor::MAX_r9 = 1

private

Definition at line 129 of file PhotonMonitor.h.

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

private

Definition at line 154 of file PhotonMonitor.h.

Referenced by analyze().

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

private

Definition at line 82 of file PhotonMonitor.h.

Referenced by analyze().

double PhotonMonitor::MIN_hoe = 0

private

Definition at line 139 of file PhotonMonitor.h.

double PhotonMonitor::MIN_r9 = 0

private

Definition at line 130 of file PhotonMonitor.h.

unsigned int PhotonMonitor::N_ETA = 34

private

Definition at line 122 of file PhotonMonitor.h.

unsigned int PhotonMonitor::N_PHI1 = 64

private

Definition at line 115 of file PhotonMonitor.h.

unsigned int PhotonMonitor::N_r9 = 50

private

Definition at line 131 of file PhotonMonitor.h.

unsigned int PhotonMonitor::nelectrons_

private

Definition at line 160 of file PhotonMonitor.h.

Referenced by analyze().

unsigned int PhotonMonitor::njets_

private

Definition at line 158 of file PhotonMonitor.h.

Referenced by analyze().

unsigned int PhotonMonitor::nphotons_

private

Definition at line 159 of file PhotonMonitor.h.

Referenced by analyze(), and bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> PhotonMonitor::num_genTriggerEventFlag_

private

Definition at line 151 of file PhotonMonitor.h.

Referenced by analyze(), and bookHistograms().

const MEbinning PhotonMonitor::phi_binning_1

private

Initial value:

{
    N_PHI1, -MAX_PHI1, MAX_PHI1
      }

Definition at line 116 of file PhotonMonitor.h.

Referenced by bookHistograms().

MEbinning PhotonMonitor::photon_binning_

private

Definition at line 93 of file PhotonMonitor.h.

Referenced by bookHistograms().

std::vector<double> PhotonMonitor::photon_variable_binning_

private

Definition at line 90 of file PhotonMonitor.h.

Referenced by bookHistograms().

PhotonME PhotonMonitor::photonEtaME_

private

Definition at line 105 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::photonEtaPhiME_

private

Definition at line 110 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::photonHoverEME_

private

Definition at line 112 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::photonME_

private

Definition at line 106 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::photonME_variableBinning_

private

Definition at line 108 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::photonPhiME_

private

Definition at line 107 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::photonr9ME_

private

Definition at line 111 of file PhotonMonitor.h.

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

StringCutObjectSelector<reco::Photon,true> PhotonMonitor::photonSelection_

private

Definition at line 157 of file PhotonMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::PhotonCollection> PhotonMonitor::photonToken_

private

Definition at line 85 of file PhotonMonitor.h.

Referenced by analyze().

PhotonME PhotonMonitor::photonVsLS_

private

Definition at line 109 of file PhotonMonitor.h.

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

const MEbinning PhotonMonitor::r9_binning_

private

Initial value:

{
    N_r9, MIN_r9, MAX_r9
      }

Definition at line 132 of file PhotonMonitor.h.

Referenced by bookHistograms().

PhotonME PhotonMonitor::subphotonEtaME_

private

Definition at line 96 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::subphotonEtaPhiME_

private

Definition at line 100 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::subphotonHoverEME_

private

Definition at line 102 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::subphotonME_

private

Definition at line 97 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::subphotonME_variableBinning_

private

Definition at line 99 of file PhotonMonitor.h.

Referenced by PhotonMonitor().

PhotonME PhotonMonitor::subphotonPhiME_

private

Definition at line 98 of file PhotonMonitor.h.

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

PhotonME PhotonMonitor::subphotonr9ME_

private

Definition at line 101 of file PhotonMonitor.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation