#include <EgammaObjects.h>

Inheritance diagram for EgammaObjects:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

	EgammaObjects (const edm::ParameterSet &)

void	endJob () override

	~EgammaObjects () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyzeElectrons (const edm::Event &, const edm::EventSetup &)

void	analyzePhotons (const edm::Event &, const edm::EventSetup &)

void	applyLabels ()

void	createBookedHistoObjects ()

void	createTempHistoObjects ()

float	ecalEta (float EtaParticle, float Zvertex, float plane_Radius)

double	findRecoMass (const reco::Photon &pOne, const reco::Photon &pTwo)

double	findRecoMass (const reco::GsfElectron &eOne, const reco::GsfElectron &eTwo)

void	fitHistos ()

void	getDeltaResHistosViaSlicing ()

void	getEfficiencyHistosViaDividing ()

void	loadCMSSWObjects (const edm::ParameterSet &ps)

void	loadHistoParameters (const edm::ParameterSet &ps)

void	saveHistos ()

void	setDrawOptions ()

Private Attributes
TH2D *	_TEMP_scatterPlot_deltaEtaVsE_

TH2D *	_TEMP_scatterPlot_deltaEtaVsEt_

TH2D *	_TEMP_scatterPlot_deltaEtaVsEta_

TH2D *	_TEMP_scatterPlot_deltaEtaVsPhi_

TH2D *	_TEMP_scatterPlot_deltaPhiVsE_

TH2D *	_TEMP_scatterPlot_deltaPhiVsEt_

TH2D *	_TEMP_scatterPlot_deltaPhiVsEta_

TH2D *	_TEMP_scatterPlot_deltaPhiVsPhi_

TH2D *	_TEMP_scatterPlot_EOverTruthVsE_

TH2D *	_TEMP_scatterPlot_EOverTruthVsEt_

TH2D *	_TEMP_scatterPlot_EOverTruthVsEta_

TH2D *	_TEMP_scatterPlot_EOverTruthVsPhi_

TH2D *	_TEMP_scatterPlot_EtOverTruthVsE_

TH2D *	_TEMP_scatterPlot_EtOverTruthVsEt_

TH2D *	_TEMP_scatterPlot_EtOverTruthVsEta_

TH2D *	_TEMP_scatterPlot_EtOverTruthVsPhi_

int	EtCut

TH1D *	hist_All_recoMass_

TH1D *	hist_BarrelOnly_recoMass_

int	hist_bins_deltaEta_

int	hist_bins_deltaPhi_

int	hist_bins_E_

int	hist_bins_EOverTruth_

int	hist_bins_Et_

int	hist_bins_Eta_

int	hist_bins_EtaOverTruth_

int	hist_bins_EtOverTruth_

int	hist_bins_Phi_

int	hist_bins_PhiOverTruth_

int	hist_bins_recoMass_

TH1D *	hist_deltaEtaVsE_

TH1D *	hist_deltaEtaVsEt_

TH1D *	hist_deltaEtaVsEta_

TH1D *	hist_deltaEtaVsPhi_

TH1D *	hist_deltaPhiVsE_

TH1D *	hist_deltaPhiVsEt_

TH1D *	hist_deltaPhiVsEta_

TH1D *	hist_deltaPhiVsPhi_

TH1D *	hist_E_

TH1D *	hist_EEfficiency_

TH1D *	hist_EndcapOnly_recoMass_

TH1D *	hist_ENumRecoOverNumTrue_

TH1D *	hist_EOverTruth_

TH1D *	hist_EOverTruthVsE_

TH1D *	hist_EOverTruthVsEt_

TH1D *	hist_EOverTruthVsEta_

TH1D *	hist_EOverTruthVsPhi_

TH1D *	hist_Et_

TH1D *	hist_Eta_

TH1D *	hist_EtaEfficiency_

TH1D *	hist_EtaNumRecoOverNumTrue_

TH1D *	hist_EtaOverTruth_

TH1D *	hist_EtEfficiency_

TH1D *	hist_EtNumRecoOverNumTrue_

TH1D *	hist_EtOverTruth_

TH1D *	hist_EtOverTruthVsE_

TH1D *	hist_EtOverTruthVsEt_

TH1D *	hist_EtOverTruthVsEta_

TH1D *	hist_EtOverTruthVsPhi_

double	hist_max_deltaEta_

double	hist_max_deltaPhi_

double	hist_max_E_

double	hist_max_EOverTruth_

double	hist_max_Et_

double	hist_max_Eta_

double	hist_max_EtaOverTruth_

double	hist_max_EtOverTruth_

double	hist_max_Phi_

double	hist_max_PhiOverTruth_

double	hist_max_recoMass_

double	hist_min_deltaEta_

double	hist_min_deltaPhi_

double	hist_min_E_

double	hist_min_EOverTruth_

double	hist_min_Et_

double	hist_min_Eta_

double	hist_min_EtaOverTruth_

double	hist_min_EtOverTruth_

double	hist_min_Phi_

double	hist_min_PhiOverTruth_

double	hist_min_recoMass_

TH1D *	hist_Mixed_recoMass_

TH1D *	hist_Phi_

TH1D *	hist_PhiEfficiency_

TH1D *	hist_PhiNumRecoOverNumTrue_

TH1D *	hist_PhiOverTruth_

TH1D *	hist_recoMass_withBackgroud_10EtCut_

TH1D *	hist_recoMass_withBackgroud_20EtCut_

TH1D *	hist_recoMass_withBackgroud_5EtCut_

TH1D *	hist_recoMass_withBackgroud_NoEtCut_

TH1D *	hist_resolutionEtaVsE_

TH1D *	hist_resolutionEtaVsEt_

TH1D *	hist_resolutionEtaVsEta_

TH1D *	hist_resolutionEtaVsPhi_

TH1D *	hist_resolutionEtVsE_

TH1D *	hist_resolutionEtVsEt_

TH1D *	hist_resolutionEtVsEta_

TH1D *	hist_resolutionEtVsPhi_

TH1D *	hist_resolutionEVsE_

TH1D *	hist_resolutionEVsEt_

TH1D *	hist_resolutionEVsEta_

TH1D *	hist_resolutionEVsPhi_

TH1D *	hist_resolutionPhiVsE_

TH1D *	hist_resolutionPhiVsEt_

TH1D *	hist_resolutionPhiVsEta_

TH1D *	hist_resolutionPhiVsPhi_

edm::EDGetTokenT< edm::HepMCProduct >	MCTruthCollectionT_

int	particleID

std::string	particleString

edm::EDGetTokenT< reco::GsfElectronCollection >	RecoCollectionT_

TFile *	rootFile_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 25 of file EgammaObjects.h.

Constructor & Destructor Documentation

◆ EgammaObjects()

EgammaObjects::EgammaObjects ( const edm::ParameterSet & ps )

explicit

Definition at line 19 of file EgammaObjects.cc.

                                                       {
   particleID = ps.getParameter<int>("particleID");
   EtCut = ps.getParameter<int>("EtCut");
 
   if (particleID == 22)
     particleString = "Photon";
   else if (particleID == 11)
     particleString = "Electron";
   else
     throw(std::runtime_error(
         "\n\nEgammaObjects: Only 11 = Photon and 22 = Electron are acceptable parictleIDs! Exiting...\n\n"));
 
   loadCMSSWObjects(ps);
   loadHistoParameters(ps);
 
   rootFile_ = TFile::Open(ps.getParameter<std::string>("outputFile").c_str(), "RECREATE");
 
   hist_EtaEfficiency_ = nullptr;
   hist_EtaNumRecoOverNumTrue_ = nullptr;
   hist_deltaEtaVsEt_ = nullptr;
   hist_deltaEtaVsE_ = nullptr;
   hist_deltaEtaVsEta_ = nullptr;
   hist_deltaEtaVsPhi_ = nullptr;
   hist_resolutionEtaVsEt_ = nullptr;
   hist_resolutionEtaVsE_ = nullptr;
   hist_resolutionEtaVsEta_ = nullptr;
   hist_resolutionEtaVsPhi_ = nullptr;
 
   hist_Phi_ = nullptr;
   hist_PhiOverTruth_ = nullptr;
   hist_PhiEfficiency_ = nullptr;
   hist_PhiNumRecoOverNumTrue_ = nullptr;
   hist_deltaPhiVsEt_ = nullptr;
   hist_deltaPhiVsE_ = nullptr;
   hist_deltaPhiVsEta_ = nullptr;
   hist_deltaPhiVsPhi_ = nullptr;
   hist_resolutionPhiVsEt_ = nullptr;
   hist_resolutionPhiVsE_ = nullptr;
   hist_resolutionPhiVsEta_ = nullptr;
   hist_resolutionPhiVsPhi_ = nullptr;
 
   hist_All_recoMass_ = nullptr;
   hist_BarrelOnly_recoMass_ = nullptr;
   hist_EndcapOnly_recoMass_ = nullptr;
   hist_Mixed_recoMass_ = nullptr;
 
   hist_recoMass_withBackgroud_NoEtCut_ = nullptr;
   hist_recoMass_withBackgroud_5EtCut_ = nullptr;
   hist_recoMass_withBackgroud_10EtCut_ = nullptr;
   hist_recoMass_withBackgroud_20EtCut_ = nullptr;
 
   _TEMP_scatterPlot_EtOverTruthVsEt_ = nullptr;
   _TEMP_scatterPlot_EtOverTruthVsE_ = nullptr;
   _TEMP_scatterPlot_EtOverTruthVsEta_ = nullptr;
   _TEMP_scatterPlot_EtOverTruthVsPhi_ = nullptr;
 
   _TEMP_scatterPlot_EOverTruthVsEt_ = nullptr;
   _TEMP_scatterPlot_EOverTruthVsE_ = nullptr;
   _TEMP_scatterPlot_EOverTruthVsEta_ = nullptr;
   _TEMP_scatterPlot_EOverTruthVsPhi_ = nullptr;
 
   _TEMP_scatterPlot_deltaEtaVsEt_ = nullptr;
   _TEMP_scatterPlot_deltaEtaVsE_ = nullptr;
   _TEMP_scatterPlot_deltaEtaVsEta_ = nullptr;
   _TEMP_scatterPlot_deltaEtaVsPhi_ = nullptr;
 
   _TEMP_scatterPlot_deltaPhiVsEt_ = nullptr;
   _TEMP_scatterPlot_deltaPhiVsE_ = nullptr;
   _TEMP_scatterPlot_deltaPhiVsEta_ = nullptr;
   _TEMP_scatterPlot_deltaPhiVsPhi_ = nullptr;
 }

◆ ~EgammaObjects()

EgammaObjects::~EgammaObjects ( )

override

Definition at line 142 of file EgammaObjects.cc.

References rootFile_.

142 { delete rootFile_; }

EgammaObjects::rootFile_

TFile * rootFile_

Definition: EgammaObjects.h:35

Member Function Documentation

◆ analyze()

void EgammaObjects::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 424 of file EgammaObjects.cc.

References analyzeElectrons(), analyzePhotons(), and particleID.

                                                                         {
   if (particleID == 22)
     analyzePhotons(evt, es);
   else if (particleID == 11)
     analyzeElectrons(evt, es);
 }

◆ analyzeElectrons()

void EgammaObjects::analyzeElectrons	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

private

Definition at line 594 of file EgammaObjects.cc.

Referenced by analyze().

                                                                                  {
   edm::Handle<reco::GsfElectronCollection> pElectrons;
   evt.getByToken(RecoCollectionT_, pElectrons);
   if (!pElectrons.isValid()) {
     Labels l;
     labelsForToken(RecoCollectionT_, l);
     edm::LogError("DOEPlotsProducerElectrons") << "Error! can't get collection with label " << l.module;
   }
 
   const reco::GsfElectronCollection* electrons = pElectrons.product();
   std::vector<reco::GsfElectron> electronsMCMatched;
 
   for (reco::GsfElectronCollection::const_iterator aClus = electrons->begin(); aClus != electrons->end(); aClus++) {
     if (aClus->et() >= EtCut) {
       hist_Et_->Fill(aClus->et());
       hist_E_->Fill(aClus->energy());
       hist_Eta_->Fill(aClus->eta());
       hist_Phi_->Fill(aClus->phi());
     }
   }
 
   for (int firstElectron = 0, numElectrons = electrons->size(); firstElectron < numElectrons - 1; firstElectron++)
     for (int secondElectron = firstElectron + 1; secondElectron < numElectrons; secondElectron++) {
       reco::GsfElectron eOne = (*electrons)[firstElectron];
       reco::GsfElectron eTwo = (*electrons)[secondElectron];
 
       double recoMass = findRecoMass(eOne, eTwo);
 
       hist_recoMass_withBackgroud_NoEtCut_->Fill(recoMass);
 
       if (eOne.et() >= 5 && eTwo.et() >= 5)
         hist_recoMass_withBackgroud_5EtCut_->Fill(recoMass);
 
       if (eOne.et() >= 10 && eTwo.et() >= 10)
         hist_recoMass_withBackgroud_10EtCut_->Fill(recoMass);
 
       if (eOne.et() >= 20 && eTwo.et() >= 20)
         hist_recoMass_withBackgroud_20EtCut_->Fill(recoMass);
     }
 
   edm::Handle<edm::HepMCProduct> pMCTruth;
   evt.getByToken(MCTruthCollectionT_, pMCTruth);
   if (!pMCTruth.isValid()) {
     Labels l;
     labelsForToken(MCTruthCollectionT_, l);
     edm::LogError("DOEPlotsProducerElectrons") << "Error! can't get collection with label " << l.module;
   }
 
   const HepMC::GenEvent* genEvent = pMCTruth->GetEvent();
   for (HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
        currentParticle != genEvent->particles_end();
        currentParticle++) {
     if (abs((*currentParticle)->pdg_id()) == 11 && (*currentParticle)->status() == 1 &&
         (*currentParticle)->momentum().e() / cosh((*currentParticle)->momentum().eta()) >= EtCut) {
       double phiTrue = (*currentParticle)->momentum().phi();
       double etaTrue = (*currentParticle)->momentum().eta();
       double eTrue = (*currentParticle)->momentum().e();
       double etTrue = (*currentParticle)->momentum().e() / cosh(etaTrue);
 
       double etaCurrent, etaFound = -999;
       double phiCurrent, phiFound = -999;
       double etCurrent, etFound = -999;
       double eCurrent, eFound = -999;
 
       reco::GsfElectron bestMatchElectron;
 
       double closestParticleDistance = 999;
 
       for (reco::GsfElectronCollection::const_iterator aClus = electrons->begin(); aClus != electrons->end(); aClus++) {
         if (aClus->et() > EtCut) {
           etaCurrent = aClus->eta();
           phiCurrent = aClus->phi();
           etCurrent = aClus->et();
           eCurrent = aClus->energy();
 
           double deltaPhi = phiCurrent - phiTrue;
           if (deltaPhi > Geom::pi())
             deltaPhi -= 2. * Geom::pi();
           if (deltaPhi < -Geom::pi())
             deltaPhi += 2. * Geom::pi();
           double deltaR = std::sqrt(std::pow(etaCurrent - etaTrue, 2) + std::pow(deltaPhi, 2));
 
           if (deltaR < closestParticleDistance) {
             etFound = etCurrent;
             eFound = eCurrent;
             etaFound = etaCurrent;
             phiFound = phiCurrent;
             closestParticleDistance = deltaR;
             bestMatchElectron = *aClus;
           }
         }
       }
 
       if (closestParticleDistance < 0.05 && etFound / etTrue > .5 && etFound / etTrue < 1.5) {
         hist_EtOverTruth_->Fill(etFound / etTrue);
         hist_EOverTruth_->Fill(eFound / eTrue);
         hist_EtaOverTruth_->Fill(etaFound / etaTrue);
         hist_PhiOverTruth_->Fill(phiFound / phiTrue);
 
         hist_EtEfficiency_->Fill(etTrue);
         hist_EEfficiency_->Fill(eTrue);
         hist_EtaEfficiency_->Fill(etaTrue);
         hist_PhiEfficiency_->Fill(phiTrue);
 
         double deltaPhi = phiFound - phiTrue;
         if (deltaPhi > Geom::pi())
           deltaPhi -= 2. * Geom::pi();
         if (deltaPhi < -Geom::pi())
           deltaPhi += 2. * Geom::pi();
 
         _TEMP_scatterPlot_EtOverTruthVsEt_->Fill(etFound, etFound / etTrue);
         _TEMP_scatterPlot_EtOverTruthVsE_->Fill(eFound, etFound / etTrue);
         _TEMP_scatterPlot_EtOverTruthVsEta_->Fill(etaFound, etFound / etTrue);
         _TEMP_scatterPlot_EtOverTruthVsPhi_->Fill(phiFound, etFound / etTrue);
 
         _TEMP_scatterPlot_EOverTruthVsEt_->Fill(etFound, eFound / eTrue);
         _TEMP_scatterPlot_EOverTruthVsE_->Fill(eFound, eFound / eTrue);
         _TEMP_scatterPlot_EOverTruthVsEta_->Fill(etaFound, eFound / eTrue);
         _TEMP_scatterPlot_EOverTruthVsPhi_->Fill(phiFound, eFound / eTrue);
 
         _TEMP_scatterPlot_deltaEtaVsEt_->Fill(etFound, etaFound - etaTrue);
         _TEMP_scatterPlot_deltaEtaVsE_->Fill(eFound, etaFound - etaTrue);
         _TEMP_scatterPlot_deltaEtaVsEta_->Fill(etaFound, etaFound - etaTrue);
         _TEMP_scatterPlot_deltaEtaVsPhi_->Fill(phiFound, etaFound - etaTrue);
 
         _TEMP_scatterPlot_deltaPhiVsEt_->Fill(etFound, deltaPhi);
         _TEMP_scatterPlot_deltaPhiVsE_->Fill(eFound, deltaPhi);
         _TEMP_scatterPlot_deltaPhiVsEta_->Fill(etaFound, deltaPhi);
         _TEMP_scatterPlot_deltaPhiVsPhi_->Fill(phiFound, deltaPhi);
 
         electronsMCMatched.push_back(bestMatchElectron);
       }
 
       hist_EtNumRecoOverNumTrue_->Fill(etTrue);
       hist_ENumRecoOverNumTrue_->Fill(eTrue);
       hist_EtaNumRecoOverNumTrue_->Fill(etaTrue);
       hist_PhiNumRecoOverNumTrue_->Fill(phiTrue);
     }
   }
 
   if (electronsMCMatched.size() == 2) {
     reco::GsfElectron eOne = electronsMCMatched[0];
     reco::GsfElectron eTwo = electronsMCMatched[1];
 
     double recoMass = findRecoMass(eOne, eTwo);
 
     hist_All_recoMass_->Fill(recoMass);
 
     if (eOne.superCluster()->seed()->algo() == 1 && eTwo.superCluster()->seed()->algo() == 1)
       hist_BarrelOnly_recoMass_->Fill(recoMass);
     else if (eOne.superCluster()->seed()->algo() == 0 && eTwo.superCluster()->seed()->algo() == 0)
       hist_EndcapOnly_recoMass_->Fill(recoMass);
     else
       hist_Mixed_recoMass_->Fill(recoMass);
   }
 }

◆ analyzePhotons()

void EgammaObjects::analyzePhotons	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

private

Definition at line 431 of file EgammaObjects.cc.

Referenced by analyze().

                                                                                {
   edm::Handle<reco::PhotonCollection> pPhotons;
   evt.getByToken(RecoCollectionT_, pPhotons);
   if (!pPhotons.isValid()) {
     Labels l;
     labelsForToken(RecoCollectionT_, l);
     edm::LogError("EgammaObjects") << "Error! can't get collection with label " << l.module;
   }
 
   const reco::PhotonCollection* photons = pPhotons.product();
   std::vector<reco::Photon> photonsMCMatched;
 
   for (reco::PhotonCollection::const_iterator aClus = photons->begin(); aClus != photons->end(); aClus++) {
     if (aClus->et() >= EtCut) {
       hist_Et_->Fill(aClus->et());
       hist_E_->Fill(aClus->energy());
       hist_Eta_->Fill(aClus->eta());
       hist_Phi_->Fill(aClus->phi());
     }
   }
 
   for (int firstPhoton = 0, numPhotons = photons->size(); firstPhoton < numPhotons - 1; firstPhoton++)
     for (int secondPhoton = firstPhoton + 1; secondPhoton < numPhotons; secondPhoton++) {
       reco::Photon pOne = (*photons)[firstPhoton];
       reco::Photon pTwo = (*photons)[secondPhoton];
 
       double recoMass = findRecoMass(pOne, pTwo);
 
       hist_recoMass_withBackgroud_NoEtCut_->Fill(recoMass);
 
       if (pOne.et() >= 5 && pTwo.et() >= 5)
         hist_recoMass_withBackgroud_5EtCut_->Fill(recoMass);
 
       if (pOne.et() >= 10 && pTwo.et() >= 10)
         hist_recoMass_withBackgroud_10EtCut_->Fill(recoMass);
 
       if (pOne.et() >= 20 && pTwo.et() >= 20)
         hist_recoMass_withBackgroud_20EtCut_->Fill(recoMass);
     }
 
   edm::Handle<edm::HepMCProduct> pMCTruth;
   evt.getByToken(MCTruthCollectionT_, pMCTruth);
   if (!pMCTruth.isValid()) {
     Labels l;
     labelsForToken(MCTruthCollectionT_, l);
     edm::LogError("EgammaObjects") << "Error! can't get collection with label " << l.module;
   }
 
   const HepMC::GenEvent* genEvent = pMCTruth->GetEvent();
 
   for (HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
        currentParticle != genEvent->particles_end();
        currentParticle++) {
     if (abs((*currentParticle)->pdg_id()) == 22 && (*currentParticle)->status() == 1 &&
         (*currentParticle)->momentum().e() /
                 cosh(ecalEta((*currentParticle)->momentum().eta(),
                              (*currentParticle)->production_vertex()->position().z() / 10.,
                              (*currentParticle)->production_vertex()->position().perp() / 10.)) >=
             EtCut) {
       HepMC::FourVector vtx = (*currentParticle)->production_vertex()->position();
       double phiTrue = (*currentParticle)->momentum().phi();
       double etaTrue = ecalEta((*currentParticle)->momentum().eta(), vtx.z() / 10., vtx.perp() / 10.);
       double eTrue = (*currentParticle)->momentum().e();
       double etTrue = (*currentParticle)->momentum().e() / cosh(etaTrue);
 
       double etaCurrent, etaFound = -999;
       double phiCurrent, phiFound = -999;
       double etCurrent, etFound = -999;
       double eCurrent, eFound = -999;
 
       reco::Photon bestMatchPhoton;
 
       double closestParticleDistance = 999;
 
       for (reco::PhotonCollection::const_iterator aClus = photons->begin(); aClus != photons->end(); aClus++) {
         if (aClus->et() > EtCut) {
           etaCurrent = aClus->eta();
           phiCurrent = aClus->phi();
           etCurrent = aClus->et();
           eCurrent = aClus->energy();
 
           double deltaPhi = phiCurrent - phiTrue;
           if (deltaPhi > Geom::pi())
             deltaPhi -= 2. * Geom::pi();
           if (deltaPhi < -Geom::pi())
             deltaPhi += 2. * Geom::pi();
           double deltaR = std::sqrt(std::pow(etaCurrent - etaTrue, 2) + std::pow(deltaPhi, 2));
 
           if (deltaR < closestParticleDistance) {
             etFound = etCurrent;
             eFound = eCurrent;
             etaFound = etaCurrent;
             phiFound = phiCurrent;
             closestParticleDistance = deltaR;
             bestMatchPhoton = *aClus;
           }
         }
       }
 
       if (closestParticleDistance < 0.05 && etFound / etTrue > .5 && etFound / etTrue < 1.5) {
         hist_EtOverTruth_->Fill(etFound / etTrue);
         hist_EOverTruth_->Fill(eFound / eTrue);
         hist_EtaOverTruth_->Fill(etaFound / etaTrue);
         hist_PhiOverTruth_->Fill(phiFound / phiTrue);
 
         hist_EtEfficiency_->Fill(etTrue);
         hist_EEfficiency_->Fill(eTrue);
         hist_EtaEfficiency_->Fill(etaTrue);
         hist_PhiEfficiency_->Fill(phiTrue);
 
         double deltaPhi = phiFound - phiTrue;
         if (deltaPhi > Geom::pi())
           deltaPhi -= 2. * Geom::pi();
         if (deltaPhi < -Geom::pi())
           deltaPhi += 2. * Geom::pi();
 
         _TEMP_scatterPlot_EtOverTruthVsEt_->Fill(etFound, etFound / etTrue);
         _TEMP_scatterPlot_EtOverTruthVsE_->Fill(eFound, etFound / etTrue);
         _TEMP_scatterPlot_EtOverTruthVsEta_->Fill(etaFound, etFound / etTrue);
         _TEMP_scatterPlot_EtOverTruthVsPhi_->Fill(phiFound, etFound / etTrue);
 
         _TEMP_scatterPlot_EOverTruthVsEt_->Fill(etFound, eFound / eTrue);
         _TEMP_scatterPlot_EOverTruthVsE_->Fill(eFound, eFound / eTrue);
         _TEMP_scatterPlot_EOverTruthVsEta_->Fill(etaFound, eFound / eTrue);
         _TEMP_scatterPlot_EOverTruthVsPhi_->Fill(phiFound, eFound / eTrue);
 
         _TEMP_scatterPlot_deltaEtaVsEt_->Fill(etFound, etaFound - etaTrue);
         _TEMP_scatterPlot_deltaEtaVsE_->Fill(eFound, etaFound - etaTrue);
         _TEMP_scatterPlot_deltaEtaVsEta_->Fill(etaFound, etaFound - etaTrue);
         _TEMP_scatterPlot_deltaEtaVsPhi_->Fill(phiFound, etaFound - etaTrue);
 
         _TEMP_scatterPlot_deltaPhiVsEt_->Fill(etFound, deltaPhi);
         _TEMP_scatterPlot_deltaPhiVsE_->Fill(eFound, deltaPhi);
         _TEMP_scatterPlot_deltaPhiVsEta_->Fill(etaFound, deltaPhi);
         _TEMP_scatterPlot_deltaPhiVsPhi_->Fill(phiFound, deltaPhi);
 
         photonsMCMatched.push_back(bestMatchPhoton);
       }
 
       hist_EtNumRecoOverNumTrue_->Fill(etTrue);
       hist_ENumRecoOverNumTrue_->Fill(eTrue);
       hist_EtaNumRecoOverNumTrue_->Fill(etaTrue);
       hist_PhiNumRecoOverNumTrue_->Fill(phiTrue);
     }
   }
 
   if (photonsMCMatched.size() == 2) {
     reco::Photon pOne = photonsMCMatched[0];
     reco::Photon pTwo = photonsMCMatched[1];
 
     double recoMass = findRecoMass(pOne, pTwo);
 
     hist_All_recoMass_->Fill(recoMass);
 
     if (pOne.superCluster()->seed()->algo() == 1 && pTwo.superCluster()->seed()->algo() == 1)
       hist_BarrelOnly_recoMass_->Fill(recoMass);
     else if (pOne.superCluster()->seed()->algo() == 0 && pTwo.superCluster()->seed()->algo() == 0)
       hist_EndcapOnly_recoMass_->Fill(recoMass);
     else
       hist_Mixed_recoMass_->Fill(recoMass);
   }
 }

◆ applyLabels()

void EgammaObjects::applyLabels ( )

private

Definition at line 1002 of file EgammaObjects.cc.

Referenced by endJob().

                                 {
   hist_Et_->GetXaxis()->SetTitle("Et (GeV)");
   hist_Et_->GetYaxis()->SetTitle("# per Et Bin");
   hist_EtOverTruth_->GetXaxis()->SetTitle("Reco Et/True Et");
   hist_EtOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
   hist_EtEfficiency_->GetXaxis()->SetTitle("Et (GeV)");
   hist_EtEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per Et Bin");
   hist_EtNumRecoOverNumTrue_->GetXaxis()->SetTitle("Et (GeV)");
   hist_EtNumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per Et Bin");
   hist_EtOverTruthVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_EtOverTruthVsEt_->GetYaxis()->SetTitle("Reco Et/True Et per Et Bin");
   hist_EtOverTruthVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_EtOverTruthVsE_->GetYaxis()->SetTitle("Reco Et/True Et per E Bin");
   hist_EtOverTruthVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_EtOverTruthVsEta_->GetYaxis()->SetTitle("Reco Et/True Et per Eta Bin");
   hist_EtOverTruthVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_EtOverTruthVsPhi_->GetYaxis()->SetTitle("Reco Et/True Et per Phi Bin");
   hist_resolutionEtVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_resolutionEtVsEt_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per Et Bin");
   hist_resolutionEtVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_resolutionEtVsE_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per E Bin");
   hist_resolutionEtVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_resolutionEtVsEta_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per Eta Bin");
   hist_resolutionEtVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_resolutionEtVsPhi_->GetYaxis()->SetTitle("#sigma of Reco Et/True Et per Phi Bin");
 
   hist_E_->GetXaxis()->SetTitle("E (GeV)");
   hist_E_->GetYaxis()->SetTitle("# per E Bin");
   hist_EOverTruth_->GetXaxis()->SetTitle("Reco E/True E");
   hist_EOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
   hist_EEfficiency_->GetXaxis()->SetTitle("E (GeV)");
   hist_EEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per E Bin");
   hist_ENumRecoOverNumTrue_->GetXaxis()->SetTitle("E (GeV)");
   hist_ENumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per E Bin");
   hist_EOverTruthVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_EOverTruthVsEt_->GetYaxis()->SetTitle("Reco E/True E per Et Bin");
   hist_EOverTruthVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_EOverTruthVsE_->GetYaxis()->SetTitle("Reco E/True E per E Bin");
   hist_EOverTruthVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_EOverTruthVsEta_->GetYaxis()->SetTitle("Reco E/True E per Eta Bin");
   hist_EOverTruthVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_EOverTruthVsPhi_->GetYaxis()->SetTitle("Reco E/True E per Phi Bin");
   hist_resolutionEVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_resolutionEVsEt_->GetYaxis()->SetTitle("#sigma of Reco E/True E per Et Bin");
   hist_resolutionEVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_resolutionEVsE_->GetYaxis()->SetTitle("#sigma of Reco E/True E per E Bin");
   hist_resolutionEVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_resolutionEVsEta_->GetYaxis()->SetTitle("#sigma of Reco E/True E per Eta Bin");
   hist_resolutionEVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_resolutionEVsPhi_->GetYaxis()->SetTitle("#sigma of Reco E/True E per Phi Bin");
 
   hist_Eta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_Eta_->GetYaxis()->SetTitle("# per Eta Bin");
   hist_EtaOverTruth_->GetXaxis()->SetTitle("Reco Eta/True Eta");
   hist_EtaOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
   hist_EtaEfficiency_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_EtaEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per Eta Bin");
   hist_EtaNumRecoOverNumTrue_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_EtaNumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per Eta Bin");
   hist_deltaEtaVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_deltaEtaVsEt_->GetYaxis()->SetTitle("Reco Eta - True Eta per Et Bin");
   hist_deltaEtaVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_deltaEtaVsE_->GetYaxis()->SetTitle("Reco Eta - True Eta per E Bin");
   hist_deltaEtaVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_deltaEtaVsEta_->GetYaxis()->SetTitle("Reco Eta - True Eta per Eta Bin");
   hist_deltaEtaVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_deltaEtaVsPhi_->GetYaxis()->SetTitle("Reco Eta - True Eta per Phi Bin");
   hist_resolutionEtaVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_resolutionEtaVsEt_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per Et Bin");
   hist_resolutionEtaVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_resolutionEtaVsE_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per E Bin");
   hist_resolutionEtaVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_resolutionEtaVsEta_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per Eta Bin");
   hist_resolutionEtaVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_resolutionEtaVsPhi_->GetYaxis()->SetTitle("#sigma of Reco Eta - True Eta per Phi Bin");
 
   hist_Phi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_Phi_->GetYaxis()->SetTitle("# per Phi Bin");
   hist_PhiOverTruth_->GetXaxis()->SetTitle("Reco Phi/True Phi");
   hist_PhiOverTruth_->GetYaxis()->SetTitle("# per Ratio Bin");
   hist_PhiEfficiency_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_PhiEfficiency_->GetYaxis()->SetTitle("# True Reconstructed/# True per Phi Bin");
   hist_PhiNumRecoOverNumTrue_->GetXaxis()->SetTitle("#Phi (Radians)");
   hist_PhiNumRecoOverNumTrue_->GetYaxis()->SetTitle("# Reco/# True per Phi Bin");
   hist_deltaPhiVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_deltaPhiVsEt_->GetYaxis()->SetTitle("Reco Phi - True Phi per Et Bin");
   hist_deltaPhiVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_deltaPhiVsE_->GetYaxis()->SetTitle("Reco Phi - True Phi per E Bin");
   hist_deltaPhiVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_deltaPhiVsEta_->GetYaxis()->SetTitle("Reco Phi - True Phi per Eta Bin");
   hist_deltaPhiVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_deltaPhiVsPhi_->GetYaxis()->SetTitle("Reco Phi - True Phi per Phi Bin");
   hist_resolutionPhiVsEt_->GetXaxis()->SetTitle("Et (GeV)");
   hist_resolutionPhiVsEt_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per Et Bin");
   hist_resolutionPhiVsE_->GetXaxis()->SetTitle("E (GeV)");
   hist_resolutionPhiVsE_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per E Bin");
   hist_resolutionPhiVsEta_->GetXaxis()->SetTitle("#eta (Radians)");
   hist_resolutionPhiVsEta_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per Eta Bin");
   hist_resolutionPhiVsPhi_->GetXaxis()->SetTitle("#phi (Radians)");
   hist_resolutionPhiVsPhi_->GetYaxis()->SetTitle("#sigma of Reco Phi - True Phi per Phi Bin");
 
   std::string recoParticleName;
 
   if (particleID == 22)
     recoParticleName = "Reco Higgs";
   else if (particleID == 11)
     recoParticleName = "Reco Z";
 
   hist_All_recoMass_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_All_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_BarrelOnly_recoMass_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_BarrelOnly_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_EndcapOnly_recoMass_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_EndcapOnly_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_Mixed_recoMass_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_Mixed_recoMass_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_recoMass_withBackgroud_NoEtCut_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_recoMass_withBackgroud_NoEtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_recoMass_withBackgroud_5EtCut_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_recoMass_withBackgroud_5EtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_recoMass_withBackgroud_10EtCut_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_recoMass_withBackgroud_10EtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
   hist_recoMass_withBackgroud_20EtCut_->GetXaxis()->SetTitle((recoParticleName + " Mass (GeV)").c_str());
   hist_recoMass_withBackgroud_20EtCut_->GetYaxis()->SetTitle("# of Reco Masses per Mass Bin");
 }

◆ beginJob()

void EgammaObjects::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 144 of file EgammaObjects.cc.

References createBookedHistoObjects(), and createTempHistoObjects().

                              {
   TH1::SetDefaultSumw2(true);
 
   createBookedHistoObjects();
   createTempHistoObjects();
 }

◆ createBookedHistoObjects()

void EgammaObjects::createBookedHistoObjects ( )

private

Definition at line 151 of file EgammaObjects.cc.

References hist_All_recoMass_, hist_BarrelOnly_recoMass_, hist_bins_E_, hist_bins_EOverTruth_, hist_bins_Et_, hist_bins_Eta_, hist_bins_EtaOverTruth_, hist_bins_EtOverTruth_, hist_bins_Phi_, hist_bins_PhiOverTruth_, hist_bins_recoMass_, hist_E_, hist_EEfficiency_, hist_EndcapOnly_recoMass_, hist_ENumRecoOverNumTrue_, hist_EOverTruth_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtaOverTruth_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_EtOverTruth_, hist_max_E_, hist_max_EOverTruth_, hist_max_Et_, hist_max_Eta_, hist_max_EtaOverTruth_, hist_max_EtOverTruth_, hist_max_Phi_, hist_max_PhiOverTruth_, hist_max_recoMass_, hist_min_E_, hist_min_EOverTruth_, hist_min_Et_, hist_min_Eta_, hist_min_EtaOverTruth_, hist_min_EtOverTruth_, hist_min_Phi_, hist_min_PhiOverTruth_, hist_min_recoMass_, hist_Mixed_recoMass_, hist_Phi_, hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, hist_PhiOverTruth_, hist_recoMass_withBackgroud_10EtCut_, hist_recoMass_withBackgroud_20EtCut_, hist_recoMass_withBackgroud_5EtCut_, hist_recoMass_withBackgroud_NoEtCut_, particleID, particleString, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by beginJob().

                                              {
   hist_Et_ =
       new TH1D("hist_Et_", ("Et Distribution of " + particleString).c_str(), hist_bins_Et_, hist_min_Et_, hist_max_Et_);
   hist_EtOverTruth_ = new TH1D("hist_EtOverTruth_",
                                ("Reco Et over True Et of " + particleString).c_str(),
                                hist_bins_EtOverTruth_,
                                hist_min_EtOverTruth_,
                                hist_max_EtOverTruth_);
   hist_EtEfficiency_ = new TH1D("hist_EtEfficiency_",
                                 ("# of True " + particleString + " Reconstructed over # of True " + particleString +
                                  " VS Et of " + particleString)
                                     .c_str(),
                                 hist_bins_Et_,
                                 hist_min_Et_,
                                 hist_max_Et_);
   hist_EtNumRecoOverNumTrue_ = new TH1D(
       "hist_EtNumRecoOverNumTrue_",
       ("# of Reco " + particleString + " over # of True " + particleString + " VS Et of " + particleString).c_str(),
       hist_bins_Et_,
       hist_min_Et_,
       hist_max_Et_);
 
   hist_E_ =
       new TH1D("hist_E_", ("E Distribution of " + particleString).c_str(), hist_bins_E_, hist_min_E_, hist_max_E_);
   hist_EOverTruth_ = new TH1D("hist_EOverTruth_",
                               ("Reco E over True E of " + particleString).c_str(),
                               hist_bins_EOverTruth_,
                               hist_min_EOverTruth_,
                               hist_max_EOverTruth_);
   hist_EEfficiency_ = new TH1D(
       "hist_EEfficiency_",
       ("# of True " + particleString + " Reconstructed over # of True " + particleString + " VS E of " + particleString)
           .c_str(),
       hist_bins_E_,
       hist_min_E_,
       hist_max_E_);
   hist_ENumRecoOverNumTrue_ = new TH1D(
       "hist_ENumRecoOverNumTrue_",
       ("# of Reco " + particleString + " over # of True " + particleString + " VS E of " + particleString).c_str(),
       hist_bins_E_,
       hist_min_E_,
       hist_max_E_);
 
   hist_Eta_ = new TH1D(
       "hist_Eta_", ("Eta Distribution of " + particleString).c_str(), hist_bins_Eta_, hist_min_Eta_, hist_max_Eta_);
   hist_EtaOverTruth_ = new TH1D("hist_EtaOverTruth_",
                                 ("Reco Eta over True Eta of " + particleString).c_str(),
                                 hist_bins_EtaOverTruth_,
                                 hist_min_EtaOverTruth_,
                                 hist_max_EtaOverTruth_);
   hist_EtaEfficiency_ = new TH1D("hist_EtaEfficiency_",
                                  ("# of True " + particleString + " Reconstructed over # of True " + particleString +
                                   " VS Eta of " + particleString)
                                      .c_str(),
                                  hist_bins_Eta_,
                                  hist_min_Eta_,
                                  hist_max_Eta_);
   hist_EtaNumRecoOverNumTrue_ = new TH1D(
       "hist_EtaNumRecoOverNumTrue_",
       ("# of Reco " + particleString + " over # of True " + particleString + " VS Eta of " + particleString).c_str(),
       hist_bins_Eta_,
       hist_min_Eta_,
       hist_max_Eta_);
 
   hist_Phi_ = new TH1D(
       "hist_Phi_", ("Phi Distribution of " + particleString).c_str(), hist_bins_Phi_, hist_min_Phi_, hist_max_Phi_);
   hist_PhiOverTruth_ = new TH1D("hist_PhiOverTruth_",
                                 ("Reco Phi over True Phi of " + particleString).c_str(),
                                 hist_bins_PhiOverTruth_,
                                 hist_min_PhiOverTruth_,
                                 hist_max_PhiOverTruth_);
   hist_PhiEfficiency_ = new TH1D("hist_PhiEfficiency_",
                                  ("# of True " + particleString + " Reconstructed over # of True " + particleString +
                                   " VS Phi of " + particleString)
                                      .c_str(),
                                  hist_bins_Phi_,
                                  hist_min_Phi_,
                                  hist_max_Phi_);
   hist_PhiNumRecoOverNumTrue_ = new TH1D(
       "hist_PhiNumRecoOverNumTrue_",
       ("# of Reco " + particleString + " over # of True " + particleString + " VS Phi of " + particleString).c_str(),
       hist_bins_Phi_,
       hist_min_Phi_,
       hist_max_Phi_);
 
   std::string recoParticleName;
 
   if (particleID == 22)
     recoParticleName = "Higgs";
   else if (particleID == 11)
     recoParticleName = "Z";
 
   hist_All_recoMass_ = new TH1D("hist_All_recoMass_",
                                 (recoParticleName + " Mass from " + particleString + " in All Regions").c_str(),
                                 hist_bins_recoMass_,
                                 hist_min_recoMass_,
                                 hist_max_recoMass_);
   hist_BarrelOnly_recoMass_ = new TH1D("hist_BarrelOnly_recoMass_",
                                        (recoParticleName + " Mass from " + particleString + " in Barrel").c_str(),
                                        hist_bins_recoMass_,
                                        hist_min_recoMass_,
                                        hist_max_recoMass_);
   hist_EndcapOnly_recoMass_ = new TH1D("hist_EndcapOnly_recoMass_",
                                        (recoParticleName + " Mass from " + particleString + " in EndCap").c_str(),
                                        hist_bins_recoMass_,
                                        hist_min_recoMass_,
                                        hist_max_recoMass_);
   hist_Mixed_recoMass_ = new TH1D("hist_Mixed_recoMass_",
                                   (recoParticleName + " Mass from " + particleString + " in Split Detectors").c_str(),
                                   hist_bins_recoMass_,
                                   hist_min_recoMass_,
                                   hist_max_recoMass_);
 
   hist_recoMass_withBackgroud_NoEtCut_ =
       new TH1D("hist_recoMass_withBackgroud_NoEtCut_",
                (recoParticleName + " Mass from " + particleString + " with Background, No Et Cut").c_str(),
                hist_bins_recoMass_,
                hist_min_recoMass_,
                hist_max_recoMass_);
   hist_recoMass_withBackgroud_5EtCut_ =
       new TH1D("hist_recoMass_withBackgroud_5EtCut_",
                (recoParticleName + " Mass from " + particleString + " with Background, 5 Et Cut").c_str(),
                hist_bins_recoMass_,
                hist_min_recoMass_,
                hist_max_recoMass_);
   hist_recoMass_withBackgroud_10EtCut_ =
       new TH1D("hist_recoMass_withBackgroud_10EtCut_",
                (recoParticleName + " Mass from " + particleString + " with Background, 10 Et Cut").c_str(),
                hist_bins_recoMass_,
                hist_min_recoMass_,
                hist_max_recoMass_);
   hist_recoMass_withBackgroud_20EtCut_ =
       new TH1D("hist_recoMass_withBackgroud_20EtCut_",
                (recoParticleName + " Mass from " + particleString + " with Background, 20 Et Cut").c_str(),
                hist_bins_recoMass_,
                hist_min_recoMass_,
                hist_max_recoMass_);
 }

◆ createTempHistoObjects()

void EgammaObjects::createTempHistoObjects ( )

private

Definition at line 290 of file EgammaObjects.cc.

Referenced by beginJob().

                                            {
   _TEMP_scatterPlot_EtOverTruthVsEt_ = new TH2D("_TEMP_scatterPlot_EtOverTruthVsEt_",
                                                 "_TEMP_scatterPlot_EtOverTruthVsEt_",
                                                 hist_bins_Et_,
                                                 hist_min_Et_,
                                                 hist_max_Et_,
                                                 hist_bins_EtOverTruth_,
                                                 hist_min_EtOverTruth_,
                                                 hist_max_EtOverTruth_);
   _TEMP_scatterPlot_EtOverTruthVsE_ = new TH2D("_TEMP_scatterPlot_EtOverTruthVsE_",
                                                "_TEMP_scatterPlot_EtOverTruthVsE_",
                                                hist_bins_E_,
                                                hist_min_E_,
                                                hist_max_E_,
                                                hist_bins_EtOverTruth_,
                                                hist_min_EtOverTruth_,
                                                hist_max_EtOverTruth_);
   _TEMP_scatterPlot_EtOverTruthVsEta_ = new TH2D("_TEMP_scatterPlot_EtOverTruthVsEta_",
                                                  "_TEMP_scatterPlot_EtOverTruthVsEta_",
                                                  hist_bins_Eta_,
                                                  hist_min_Eta_,
                                                  hist_max_Eta_,
                                                  hist_bins_EtOverTruth_,
                                                  hist_min_EtOverTruth_,
                                                  hist_max_EtOverTruth_);
   _TEMP_scatterPlot_EtOverTruthVsPhi_ = new TH2D("_TEMP_scatterPlot_EtOverTruthVsPhi_",
                                                  "_TEMP_scatterPlot_EtOverTruthVsPhi_",
                                                  hist_bins_Phi_,
                                                  hist_min_Phi_,
                                                  hist_max_Phi_,
                                                  hist_bins_EtOverTruth_,
                                                  hist_min_EtOverTruth_,
                                                  hist_max_EtOverTruth_);
 
   _TEMP_scatterPlot_EOverTruthVsEt_ = new TH2D("_TEMP_scatterPlot_EOverTruthVsEt_",
                                                "_TEMP_scatterPlot_EOverTruthVsEt_",
                                                hist_bins_Et_,
                                                hist_min_Et_,
                                                hist_max_Et_,
                                                hist_bins_EOverTruth_,
                                                hist_min_EOverTruth_,
                                                hist_max_EOverTruth_);
   _TEMP_scatterPlot_EOverTruthVsE_ = new TH2D("_TEMP_scatterPlot_EOverTruthVsE_",
                                               "_TEMP_scatterPlot_EOverTruthVsE_",
                                               hist_bins_E_,
                                               hist_min_E_,
                                               hist_max_E_,
                                               hist_bins_EOverTruth_,
                                               hist_min_EOverTruth_,
                                               hist_max_EOverTruth_);
   _TEMP_scatterPlot_EOverTruthVsEta_ = new TH2D("_TEMP_scatterPlot_EOverTruthVsEta_",
                                                 "_TEMP_scatterPlot_EOverTruthVsEta_",
                                                 hist_bins_Eta_,
                                                 hist_min_Eta_,
                                                 hist_max_Eta_,
                                                 hist_bins_EOverTruth_,
                                                 hist_min_EOverTruth_,
                                                 hist_max_EOverTruth_);
   _TEMP_scatterPlot_EOverTruthVsPhi_ = new TH2D("_TEMP_scatterPlot_EOverTruthVsPhi_",
                                                 "_TEMP_scatterPlot_EOverTruthVsPhi_",
                                                 hist_bins_Phi_,
                                                 hist_min_Phi_,
                                                 hist_max_Phi_,
                                                 hist_bins_EOverTruth_,
                                                 hist_min_EOverTruth_,
                                                 hist_max_EOverTruth_);
 
   _TEMP_scatterPlot_deltaEtaVsEt_ = new TH2D("_TEMP_scatterPlot_deltaEtaVsEt_",
                                              "_TEMP_scatterPlot_deltaEtaVsEt_",
                                              hist_bins_Et_,
                                              hist_min_Et_,
                                              hist_max_Et_,
                                              hist_bins_deltaEta_,
                                              hist_min_deltaEta_,
                                              hist_max_deltaEta_);
   _TEMP_scatterPlot_deltaEtaVsE_ = new TH2D("_TEMP_scatterPlot_deltaEtaVsE_",
                                             "_TEMP_scatterPlot_deltaEtaVsE_",
                                             hist_bins_E_,
                                             hist_min_E_,
                                             hist_max_E_,
                                             hist_bins_deltaEta_,
                                             hist_min_deltaEta_,
                                             hist_max_deltaEta_);
   _TEMP_scatterPlot_deltaEtaVsEta_ = new TH2D("_TEMP_scatterPlot_deltaEtaVsEta_",
                                               "_TEMP_scatterPlot_deltaEtaVsEta_",
                                               hist_bins_Eta_,
                                               hist_min_Eta_,
                                               hist_max_Eta_,
                                               hist_bins_deltaEta_,
                                               hist_min_deltaEta_,
                                               hist_max_deltaEta_);
   _TEMP_scatterPlot_deltaEtaVsPhi_ = new TH2D("_TEMP_scatterPlot_deltaEtaVsPhi_",
                                               "_TEMP_scatterPlot_deltaEtaVsPhi_",
                                               hist_bins_Phi_,
                                               hist_min_Phi_,
                                               hist_max_Phi_,
                                               hist_bins_deltaEta_,
                                               hist_min_deltaEta_,
                                               hist_max_deltaEta_);
 
   _TEMP_scatterPlot_deltaPhiVsEt_ = new TH2D("_TEMP_scatterPlot_deltaPhiVsEt_",
                                              "_TEMP_scatterPlot_deltaPhiVsEt_",
                                              hist_bins_Et_,
                                              hist_min_Et_,
                                              hist_max_Et_,
                                              hist_bins_deltaPhi_,
                                              hist_min_deltaPhi_,
                                              hist_max_deltaPhi_);
   _TEMP_scatterPlot_deltaPhiVsE_ = new TH2D("_TEMP_scatterPlot_deltaPhiVsE_",
                                             "_TEMP_scatterPlot_deltaPhiVsE_",
                                             hist_bins_E_,
                                             hist_min_E_,
                                             hist_max_E_,
                                             hist_bins_deltaPhi_,
                                             hist_min_deltaPhi_,
                                             hist_max_deltaPhi_);
   _TEMP_scatterPlot_deltaPhiVsEta_ = new TH2D("_TEMP_scatterPlot_deltaPhiVsEta_",
                                               "_TEMP_scatterPlot_deltaPhiVsEta_",
                                               hist_bins_Eta_,
                                               hist_min_Eta_,
                                               hist_max_Eta_,
                                               hist_bins_deltaPhi_,
                                               hist_min_deltaPhi_,
                                               hist_max_deltaPhi_);
   _TEMP_scatterPlot_deltaPhiVsPhi_ = new TH2D("_TEMP_scatterPlot_deltaPhiVsPhi_",
                                               "_TEMP_scatterPlot_deltaPhiVsPhi_",
                                               hist_bins_Phi_,
                                               hist_min_Phi_,
                                               hist_max_Phi_,
                                               hist_bins_deltaPhi_,
                                               hist_min_deltaPhi_,
                                               hist_max_deltaPhi_);
 }

◆ ecalEta()

float EgammaObjects::ecalEta	(	float	EtaParticle,
		float	Zvertex,
		float	plane_Radius
	)

private

Definition at line 771 of file EgammaObjects.cc.

References funct::abs(), ETA, etaBarrelEndcap, CrabHelper::log, Geom::pi(), R_ECAL, funct::tan(), and Z_Endcap.

Referenced by analyzePhotons().

                                                                                  {
   const float R_ECAL = 136.5;
   const float Z_Endcap = 328.0;
   const float etaBarrelEndcap = 1.479;
 
   if (EtaParticle != 0.) {
     float Theta = 0.0;
     float ZEcal = (R_ECAL - plane_Radius) * sinh(EtaParticle) + Zvertex;
 
     if (ZEcal != 0.0)
       Theta = atan(R_ECAL / ZEcal);
     if (Theta < 0.0)
       Theta = Theta + Geom::pi();
 
     float ETA = -log(tan(0.5 * Theta));
 
     if (std::abs(ETA) > etaBarrelEndcap) {
       float Zend = Z_Endcap;
       if (EtaParticle < 0.0)
         Zend = -Zend;
       float Zlen = Zend - Zvertex;
       float RR = Zlen / sinh(EtaParticle);
       Theta = atan((RR + plane_Radius) / Zend);
       if (Theta < 0.0)
         Theta = Theta + Geom::pi();
       ETA = -log(tan(0.5 * Theta));
     }
 
     return ETA;
   } else {
     edm::LogWarning("") << "[EgammaObjects::ecalEta] Warning: Eta equals to zero, not correcting";
     return EtaParticle;
   }
 }

◆ endJob()

void EgammaObjects::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 806 of file EgammaObjects.cc.

References applyLabels(), fitHistos(), getDeltaResHistosViaSlicing(), getEfficiencyHistosViaDividing(), particleString, rootFile_, saveHistos(), and setDrawOptions().

                            {
   rootFile_->cd();
   rootFile_->mkdir(particleString.c_str());
 
   getDeltaResHistosViaSlicing();
   getEfficiencyHistosViaDividing();
   fitHistos();
 
   applyLabels();
   setDrawOptions();
   saveHistos();
   rootFile_->Close();
 }

◆ findRecoMass() [1/2]

double EgammaObjects::findRecoMass	(	const reco::Photon &	pOne,
		const reco::Photon &	pTwo
	)

private

Definition at line 751 of file EgammaObjects.cc.

References funct::cos(), HBHEDarkening_cff::energy, mathSSE::sqrt(), and reco::Photon::superCluster().

Referenced by analyzeElectrons(), and analyzePhotons().

                                                                                  {
   double cosTheta = (cos(pOne.superCluster()->phi() - pTwo.superCluster()->phi()) +
                      sinh(pOne.superCluster()->eta()) * sinh(pTwo.superCluster()->eta())) /
                     (cosh(pOne.superCluster()->eta()) * cosh(pTwo.superCluster()->eta()));
 
   double recoMass = sqrt(2 * (pOne.superCluster())->energy() * (pTwo.superCluster())->energy() * (1 - cosTheta));
 
   return recoMass;
 }

◆ findRecoMass() [2/2]

double EgammaObjects::findRecoMass	(	const reco::GsfElectron &	eOne,
		const reco::GsfElectron &	eTwo
	)

private

Definition at line 761 of file EgammaObjects.cc.

References reco::GsfElectron::caloEnergy(), reco::GsfElectron::caloPosition(), funct::cos(), and mathSSE::sqrt().

                                                                                            {
   double cosTheta = (cos(eOne.caloPosition().phi() - eTwo.caloPosition().phi()) +
                      sinh(eOne.caloPosition().eta()) * sinh(eTwo.caloPosition().eta())) /
                     (cosh(eOne.caloPosition().eta()) * cosh(eTwo.caloPosition().eta()));
 
   double recoMass = sqrt(2 * eOne.caloEnergy() * eTwo.caloEnergy() * (1 - cosTheta));
 
   return recoMass;
 }

◆ fitHistos()

void EgammaObjects::fitHistos ( )

private

Definition at line 956 of file EgammaObjects.cc.

References hist_EOverTruth_, hist_EtaOverTruth_, hist_EtOverTruth_, and hist_PhiOverTruth_.

Referenced by endJob().

                               {
   //Use our own copy for thread safety
   TF1 gaus("mygaus", "gaus");
   hist_EtOverTruth_->Fit(&gaus, "QEM");
   //  hist_EtNumRecoOverNumTrue_->Fit("pol1","QEM");
 
   hist_EOverTruth_->Fit(&gaus, "QEM");
   //  hist_ENumRecoOverNumTrue_->Fit("pol1","QEM");
 
   hist_EtaOverTruth_->Fit(&gaus, "QEM");
   //  hist_EtaNumRecoOverNumTrue_->Fit("pol1","QEM");
 
   hist_PhiOverTruth_->Fit(&gaus, "QEM");
   //  hist_PhiNumRecoOverNumTrue_->Fit("pol1","QEM");
 
   /*
   hist_EtOverTruthVsEt_->Fit("pol1","QEM");
   hist_EtOverTruthVsEta_->Fit("pol1","QEM");
   hist_EtOverTruthVsPhi_->Fit("pol1","QEM");
   hist_resolutionEtVsEt_->Fit("pol1","QEM");
   hist_resolutionEtVsEta_->Fit("pol1","QEM");
   hist_resolutionEtVsPhi_->Fit("pol1","QEM");
 
   hist_EOverTruthVsEt_->Fit("pol1","QEM");
   hist_EOverTruthVsEta_->Fit("pol1","QEM");
   hist_EOverTruthVsPhi_->Fit("pol1","QEM");
   hist_resolutionEVsEt_->Fit("pol1","QEM");
   hist_resolutionEVsEta_->Fit("pol1","QEM");
   hist_resolutionEVsPhi_->Fit("pol1","QEM");
 
   hist_deltaEtaVsEt_->Fit("pol1","QEM");
   hist_deltaEtaVsEta_->Fit("pol1","QEM");
   hist_deltaEtaVsPhi_->Fit("pol1","QEM");
   hist_resolutionEtaVsEt_->Fit("pol1","QEM");
   hist_resolutionEtaVsEta_->Fit("pol1","QEM");
   hist_resolutionEtaVsPhi_->Fit("pol1","QEM");
 
   hist_deltaPhiVsEt_->Fit("pol1","QEM");
   hist_deltaPhiVsEta_->Fit("pol1","QEM");
   hist_deltaPhiVsPhi_->Fit("pol1","QEM");
   hist_resolutionPhiVsEt_->Fit("pol1","QEM");
   hist_resolutionPhiVsEta_->Fit("pol1","QEM");
   hist_resolutionPhiVsPhi_->Fit("pol1","QEM");
   */
 }

◆ getDeltaResHistosViaSlicing()

void EgammaObjects::getDeltaResHistosViaSlicing ( )

private

Definition at line 820 of file EgammaObjects.cc.

Referenced by endJob().

                                                 {
   _TEMP_scatterPlot_EtOverTruthVsEt_->FitSlicesY(nullptr, 1, hist_bins_Et_, 10, "QRG3");
   _TEMP_scatterPlot_EtOverTruthVsE_->FitSlicesY(nullptr, 1, hist_bins_E_, 10, "QRG3");
   _TEMP_scatterPlot_EtOverTruthVsEta_->FitSlicesY(nullptr, 1, hist_bins_Eta_, 10, "QRG2");
   _TEMP_scatterPlot_EtOverTruthVsPhi_->FitSlicesY(nullptr, 1, hist_bins_Phi_, 10, "QRG2");
 
   _TEMP_scatterPlot_EOverTruthVsEt_->FitSlicesY(nullptr, 1, hist_bins_Et_, 10, "QRG3");
   _TEMP_scatterPlot_EOverTruthVsE_->FitSlicesY(nullptr, 1, hist_bins_E_, 10, "QRG3");
   _TEMP_scatterPlot_EOverTruthVsEta_->FitSlicesY(nullptr, 1, hist_bins_Eta_, 10, "QRG2");
   _TEMP_scatterPlot_EOverTruthVsPhi_->FitSlicesY(nullptr, 1, hist_bins_Phi_, 10, "QRG2");
 
   _TEMP_scatterPlot_deltaEtaVsEt_->FitSlicesY(nullptr, 1, hist_bins_Et_, 10, "QRG3");
   _TEMP_scatterPlot_deltaEtaVsE_->FitSlicesY(nullptr, 1, hist_bins_E_, 10, "QRG3");
   _TEMP_scatterPlot_deltaEtaVsEta_->FitSlicesY(nullptr, 1, hist_bins_Eta_, 10, "QRG2");
   _TEMP_scatterPlot_deltaEtaVsPhi_->FitSlicesY(nullptr, 1, hist_bins_Phi_, 10, "QRG2");
 
   _TEMP_scatterPlot_deltaPhiVsEt_->FitSlicesY(nullptr, 1, hist_bins_Et_, 10, "QRG3");
   _TEMP_scatterPlot_deltaPhiVsE_->FitSlicesY(nullptr, 1, hist_bins_E_, 10, "QRG3");
   _TEMP_scatterPlot_deltaPhiVsEta_->FitSlicesY(nullptr, 1, hist_bins_Eta_, 10, "QRG2");
   _TEMP_scatterPlot_deltaPhiVsPhi_->FitSlicesY(nullptr, 1, hist_bins_Phi_, 10, "QRG2");
 
   hist_EtOverTruthVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEt__1");
   hist_EtOverTruthVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsE__1");
   hist_EtOverTruthVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEta__1");
   hist_EtOverTruthVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsPhi__1");
 
   hist_EOverTruthVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEt__1");
   hist_EOverTruthVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsE__1");
   hist_EOverTruthVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEta__1");
   hist_EOverTruthVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsPhi__1");
 
   hist_deltaEtaVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEt__1");
   hist_deltaEtaVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsE__1");
   hist_deltaEtaVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEta__1");
   hist_deltaEtaVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsPhi__1");
 
   hist_deltaPhiVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEt__1");
   hist_deltaPhiVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsE__1");
   hist_deltaPhiVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEta__1");
   hist_deltaPhiVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsPhi__1");
 
   hist_EtOverTruthVsEt_->SetNameTitle("hist_EtOverTruthVsEt_",
                                       ("Reco Et over True Et VS Et of " + particleString).c_str());
   hist_EtOverTruthVsE_->SetNameTitle("hist_EtOverTruthVsE_",
                                      ("Reco Et over True Et VS E of " + particleString).c_str());
   hist_EtOverTruthVsEta_->SetNameTitle("hist_EtOverTruthVsEta_",
                                        ("Reco Et over True Et VS Eta of " + particleString).c_str());
   hist_EtOverTruthVsPhi_->SetNameTitle("hist_EtOverTruthVsPhi_",
                                        ("Reco Et over True Et VS Phi of " + particleString).c_str());
 
   hist_EOverTruthVsEt_->SetNameTitle("hist_EOverTruthVsEt_", ("Reco E over True E VS Et of " + particleString).c_str());
   hist_EOverTruthVsE_->SetNameTitle("hist_EOverTruthVsE_", ("Reco E over True E VS E of " + particleString).c_str());
   hist_EOverTruthVsEta_->SetNameTitle("hist_EOverTruthVsEta_",
                                       ("Reco E over True E VS Eta of " + particleString).c_str());
   hist_EOverTruthVsPhi_->SetNameTitle("hist_EOverTruthVsPhi_",
                                       ("Reco E over True E VS Phi of " + particleString).c_str());
 
   hist_deltaEtaVsEt_->SetNameTitle("hist_deltaEtaVsEt_", ("delta Eta VS Et of " + particleString).c_str());
   hist_deltaEtaVsE_->SetNameTitle("hist_deltaEtaVsE_", ("delta Eta VS E of " + particleString).c_str());
   hist_deltaEtaVsEta_->SetNameTitle("hist_deltaEtaVsEta_", ("delta Eta VS Eta of " + particleString).c_str());
   hist_deltaEtaVsPhi_->SetNameTitle("hist_deltaEtaVsPhi_", ("delta Eta VS Phi of " + particleString).c_str());
 
   hist_deltaPhiVsEt_->SetNameTitle("hist_deltaPhiVsEt_", ("delta Phi VS Et of " + particleString).c_str());
   hist_deltaPhiVsE_->SetNameTitle("hist_deltaPhiVsE_", ("delta Phi VS E of " + particleString).c_str());
   hist_deltaPhiVsEta_->SetNameTitle("hist_deltaPhiVsEta_", ("delta Phi VS Eta of " + particleString).c_str());
   hist_deltaPhiVsPhi_->SetNameTitle("hist_deltaPhiVsPhi_", ("delta Phi VS Phi of " + particleString).c_str());
 
   hist_resolutionEtVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEt__2");
   hist_resolutionEtVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsE__2");
   hist_resolutionEtVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsEta__2");
   hist_resolutionEtVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EtOverTruthVsPhi__2");
 
   hist_resolutionEVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEt__2");
   hist_resolutionEVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsE__2");
   hist_resolutionEVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsEta__2");
   hist_resolutionEVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_EOverTruthVsPhi__2");
 
   hist_resolutionEtaVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEt__2");
   hist_resolutionEtaVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsE__2");
   hist_resolutionEtaVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsEta__2");
   hist_resolutionEtaVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaEtaVsPhi__2");
 
   hist_resolutionPhiVsEt_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEt__2");
   hist_resolutionPhiVsE_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsE__2");
   hist_resolutionPhiVsEta_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsEta__2");
   hist_resolutionPhiVsPhi_ = (TH1D*)gDirectory->Get("_TEMP_scatterPlot_deltaPhiVsPhi__2");
 
   hist_resolutionEtVsEt_->SetNameTitle("hist_resolutionEtVsEt_",
                                        ("#sigma of Reco Et over True Et VS Et of " + particleString).c_str());
   hist_resolutionEtVsE_->SetNameTitle("hist_resolutionEtVsE_",
                                       ("#sigma of Reco Et over True Et VS E of " + particleString).c_str());
   hist_resolutionEtVsEta_->SetNameTitle("hist_resolutionEtVsEta_",
                                         ("#sigma of Reco Et over True Et VS Eta of " + particleString).c_str());
   hist_resolutionEtVsPhi_->SetNameTitle("hist_resolutionEtVsPhi_",
                                         ("#sigma of Reco Et over True Et VS Phi of " + particleString).c_str());
 
   hist_resolutionEVsEt_->SetNameTitle("hist_resolutionEVsEt_",
                                       ("#sigma of Reco E over True E VS Et of " + particleString).c_str());
   hist_resolutionEVsE_->SetNameTitle("hist_resolutionEVsE_",
                                      ("#sigma of Reco E over True E VS E of " + particleString).c_str());
   hist_resolutionEVsEta_->SetNameTitle("hist_resolutionEVsEta_",
                                        ("#sigma of Reco E over True E VS Eta of " + particleString).c_str());
   hist_resolutionEVsPhi_->SetNameTitle("hist_resolutionEVsPhi_",
                                        ("#sigma of Reco E over True E VS Phi of " + particleString).c_str());
 
   hist_resolutionEtaVsEt_->SetNameTitle("hist_resolutionEtaVsEt_",
                                         ("#sigma of delta Eta VS Et of " + particleString).c_str());
   hist_resolutionEtaVsE_->SetNameTitle("hist_resolutionEtaVsE_",
                                        ("#sigma of delta Eta VS E of " + particleString).c_str());
   hist_resolutionEtaVsEta_->SetNameTitle("hist_resolutionEtaVsEta_",
                                          ("#sigma of delta Eta VS Eta of " + particleString).c_str());
   hist_resolutionEtaVsPhi_->SetNameTitle("hist_resolutionEtaVsPhi_",
                                          ("#sigma of delta Eta VS Phi of " + particleString).c_str());
 
   hist_resolutionPhiVsEt_->SetNameTitle("hist_resolutionPhiVsEt_",
                                         ("#sigma of delta Phi VS Et of " + particleString).c_str());
   hist_resolutionPhiVsE_->SetNameTitle("hist_resolutionPhiVsE_",
                                        ("#sigma of delta Phi VS E of " + particleString).c_str());
   hist_resolutionPhiVsEta_->SetNameTitle("hist_resolutionPhiVsEta_",
                                          ("#sigma of delta Phi VS Eta of " + particleString).c_str());
   hist_resolutionPhiVsPhi_->SetNameTitle("hist_resolutionPhiVsPhi_",
                                          ("#sigma of delta Phi VS Phi of " + particleString).c_str());
 }

◆ getEfficiencyHistosViaDividing()

void EgammaObjects::getEfficiencyHistosViaDividing ( )

private

Definition at line 944 of file EgammaObjects.cc.

References hist_E_, hist_EEfficiency_, hist_ENumRecoOverNumTrue_, hist_Et_, hist_Eta_, hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, hist_Phi_, hist_PhiEfficiency_, and hist_PhiNumRecoOverNumTrue_.

Referenced by endJob().

                                                    {
   hist_EtEfficiency_->Divide(hist_EtEfficiency_, hist_EtNumRecoOverNumTrue_, 1, 1);
   hist_EEfficiency_->Divide(hist_EEfficiency_, hist_ENumRecoOverNumTrue_, 1, 1);
   hist_EtaEfficiency_->Divide(hist_EtaEfficiency_, hist_EtaNumRecoOverNumTrue_, 1, 1);
   hist_PhiEfficiency_->Divide(hist_PhiEfficiency_, hist_PhiNumRecoOverNumTrue_, 1, 1);
 
   hist_EtNumRecoOverNumTrue_->Divide(hist_Et_, hist_EtNumRecoOverNumTrue_, 1, 1);
   hist_ENumRecoOverNumTrue_->Divide(hist_E_, hist_ENumRecoOverNumTrue_, 1, 1);
   hist_EtaNumRecoOverNumTrue_->Divide(hist_Eta_, hist_EtaNumRecoOverNumTrue_, 1, 1);
   hist_PhiNumRecoOverNumTrue_->Divide(hist_Phi_, hist_PhiNumRecoOverNumTrue_, 1, 1);
 }

◆ loadCMSSWObjects()

void EgammaObjects::loadCMSSWObjects ( const edm::ParameterSet & ps )

private

Definition at line 91 of file EgammaObjects.cc.

References edm::ParameterSet::getParameter(), MCTruthCollectionT_, and RecoCollectionT_.

Referenced by EgammaObjects().

                                                               {
   MCTruthCollectionT_ = consumes<edm::HepMCProduct>(ps.getParameter<edm::InputTag>("MCTruthCollection"));
   RecoCollectionT_ = consumes<reco::GsfElectronCollection>(ps.getParameter<edm::InputTag>("RecoCollection"));
 }

◆ loadHistoParameters()

void EgammaObjects::loadHistoParameters ( const edm::ParameterSet & ps )

private

Definition at line 96 of file EgammaObjects.cc.

References edm::ParameterSet::getParameter(), hist_bins_deltaEta_, hist_bins_deltaPhi_, hist_bins_E_, hist_bins_EOverTruth_, hist_bins_Et_, hist_bins_Eta_, hist_bins_EtaOverTruth_, hist_bins_EtOverTruth_, hist_bins_Phi_, hist_bins_PhiOverTruth_, hist_bins_recoMass_, hist_max_deltaEta_, hist_max_deltaPhi_, hist_max_E_, hist_max_EOverTruth_, hist_max_Et_, hist_max_Eta_, hist_max_EtaOverTruth_, hist_max_EtOverTruth_, hist_max_Phi_, hist_max_PhiOverTruth_, hist_max_recoMass_, hist_min_deltaEta_, hist_min_deltaPhi_, hist_min_E_, hist_min_EOverTruth_, hist_min_Et_, hist_min_Eta_, hist_min_EtaOverTruth_, hist_min_EtOverTruth_, hist_min_Phi_, hist_min_PhiOverTruth_, and hist_min_recoMass_.

Referenced by EgammaObjects().

                                                                  {
   hist_min_Et_ = ps.getParameter<double>("hist_min_Et");
   hist_max_Et_ = ps.getParameter<double>("hist_max_Et");
   hist_bins_Et_ = ps.getParameter<int>("hist_bins_Et");
 
   hist_min_E_ = ps.getParameter<double>("hist_min_E");
   hist_max_E_ = ps.getParameter<double>("hist_max_E");
   hist_bins_E_ = ps.getParameter<int>("hist_bins_E");
 
   hist_min_Eta_ = ps.getParameter<double>("hist_min_Eta");
   hist_max_Eta_ = ps.getParameter<double>("hist_max_Eta");
   hist_bins_Eta_ = ps.getParameter<int>("hist_bins_Eta");
 
   hist_min_Phi_ = ps.getParameter<double>("hist_min_Phi");
   hist_max_Phi_ = ps.getParameter<double>("hist_max_Phi");
   hist_bins_Phi_ = ps.getParameter<int>("hist_bins_Phi");
 
   hist_min_EtOverTruth_ = ps.getParameter<double>("hist_min_EtOverTruth");
   hist_max_EtOverTruth_ = ps.getParameter<double>("hist_max_EtOverTruth");
   hist_bins_EtOverTruth_ = ps.getParameter<int>("hist_bins_EtOverTruth");
 
   hist_min_EOverTruth_ = ps.getParameter<double>("hist_min_EOverTruth");
   hist_max_EOverTruth_ = ps.getParameter<double>("hist_max_EOverTruth");
   hist_bins_EOverTruth_ = ps.getParameter<int>("hist_bins_EOverTruth");
 
   hist_min_EtaOverTruth_ = ps.getParameter<double>("hist_min_EtaOverTruth");
   hist_max_EtaOverTruth_ = ps.getParameter<double>("hist_max_EtaOverTruth");
   hist_bins_EtaOverTruth_ = ps.getParameter<int>("hist_bins_EtaOverTruth");
 
   hist_min_PhiOverTruth_ = ps.getParameter<double>("hist_min_PhiOverTruth");
   hist_max_PhiOverTruth_ = ps.getParameter<double>("hist_max_PhiOverTruth");
   hist_bins_PhiOverTruth_ = ps.getParameter<int>("hist_bins_PhiOverTruth");
 
   hist_min_deltaEta_ = ps.getParameter<double>("hist_min_deltaEta");
   hist_max_deltaEta_ = ps.getParameter<double>("hist_max_deltaEta");
   hist_bins_deltaEta_ = ps.getParameter<int>("hist_bins_deltaEta");
 
   hist_min_deltaPhi_ = ps.getParameter<double>("hist_min_deltaPhi");
   hist_max_deltaPhi_ = ps.getParameter<double>("hist_max_deltaPhi");
   hist_bins_deltaPhi_ = ps.getParameter<int>("hist_bins_deltaPhi");
 
   hist_min_recoMass_ = ps.getParameter<double>("hist_min_recoMass");
   hist_max_recoMass_ = ps.getParameter<double>("hist_max_recoMass");
   hist_bins_recoMass_ = ps.getParameter<int>("hist_bins_recoMass");
 }

◆ saveHistos()

void EgammaObjects::saveHistos ( )

private

Definition at line 1191 of file EgammaObjects.cc.

Referenced by endJob().

                                {
   rootFile_->cd();
   rootFile_->GetDirectory(particleString.c_str())->mkdir("ET");
   rootFile_->cd(("/" + particleString + "/ET").c_str());
 
   hist_Et_->Write();
   hist_EtOverTruth_->Write();
   hist_EtEfficiency_->Write();
   hist_EtNumRecoOverNumTrue_->Write();
   hist_EtOverTruthVsEt_->Write();
   hist_EtOverTruthVsE_->Write();
   hist_EtOverTruthVsEta_->Write();
   hist_EtOverTruthVsPhi_->Write();
   hist_resolutionEtVsEt_->Write();
   hist_resolutionEtVsE_->Write();
   hist_resolutionEtVsEta_->Write();
   hist_resolutionEtVsPhi_->Write();
 
   rootFile_->cd();
   rootFile_->GetDirectory(particleString.c_str())->mkdir("E");
   rootFile_->cd(("/" + particleString + "/E").c_str());
 
   hist_E_->Write();
   hist_EOverTruth_->Write();
   hist_EEfficiency_->Write();
   hist_ENumRecoOverNumTrue_->Write();
   hist_EOverTruthVsEt_->Write();
   hist_EOverTruthVsE_->Write();
   hist_EOverTruthVsEta_->Write();
   hist_EOverTruthVsPhi_->Write();
   hist_resolutionEVsEt_->Write();
   hist_resolutionEVsE_->Write();
   hist_resolutionEVsEta_->Write();
   hist_resolutionEVsPhi_->Write();
 
   rootFile_->cd();
   rootFile_->GetDirectory(particleString.c_str())->mkdir("Eta");
   rootFile_->cd(("/" + particleString + "/Eta").c_str());
 
   hist_Eta_->Write();
   hist_EtaOverTruth_->Write();
   hist_EtaEfficiency_->Write();
   hist_EtaNumRecoOverNumTrue_->Write();
   hist_deltaEtaVsEt_->Write();
   hist_deltaEtaVsE_->Write();
   hist_deltaEtaVsEta_->Write();
   hist_deltaEtaVsPhi_->Write();
   hist_resolutionEtaVsEt_->Write();
   hist_resolutionEtaVsE_->Write();
   hist_resolutionEtaVsEta_->Write();
   hist_resolutionEtaVsPhi_->Write();
 
   rootFile_->cd();
   rootFile_->GetDirectory(particleString.c_str())->mkdir("Phi");
   rootFile_->cd(("/" + particleString + "/Phi").c_str());
 
   hist_Phi_->Write();
   hist_PhiOverTruth_->Write();
   hist_PhiEfficiency_->Write();
   hist_PhiNumRecoOverNumTrue_->Write();
   hist_deltaPhiVsEt_->Write();
   hist_deltaPhiVsE_->Write();
   hist_deltaPhiVsEta_->Write();
   hist_deltaPhiVsPhi_->Write();
   hist_resolutionPhiVsEt_->Write();
   hist_resolutionPhiVsE_->Write();
   hist_resolutionPhiVsEta_->Write();
   hist_resolutionPhiVsPhi_->Write();
 
   std::string recoParticleName;
 
   if (particleID == 22)
     recoParticleName = "HiggsRecoMass";
   else if (particleID == 11)
     recoParticleName = "ZRecoMass";
 
   rootFile_->cd();
   rootFile_->GetDirectory(particleString.c_str())->mkdir(recoParticleName.c_str());
   rootFile_->cd(("/" + particleString + "/" + recoParticleName).c_str());
 
   hist_All_recoMass_->Write();
   hist_BarrelOnly_recoMass_->Write();
   hist_EndcapOnly_recoMass_->Write();
   hist_Mixed_recoMass_->Write();
   hist_recoMass_withBackgroud_NoEtCut_->Write();
   hist_recoMass_withBackgroud_5EtCut_->Write();
   hist_recoMass_withBackgroud_10EtCut_->Write();
   hist_recoMass_withBackgroud_20EtCut_->Write();
 
   rootFile_->cd();
   rootFile_->GetDirectory(particleString.c_str())->mkdir("_TempScatterPlots");
   rootFile_->cd(("/" + particleString + "/_TempScatterPlots").c_str());
 
   _TEMP_scatterPlot_EtOverTruthVsEt_->Write();
   _TEMP_scatterPlot_EtOverTruthVsE_->Write();
   _TEMP_scatterPlot_EtOverTruthVsEta_->Write();
   _TEMP_scatterPlot_EtOverTruthVsPhi_->Write();
 
   _TEMP_scatterPlot_EOverTruthVsEt_->Write();
   _TEMP_scatterPlot_EOverTruthVsE_->Write();
   _TEMP_scatterPlot_EOverTruthVsEta_->Write();
   _TEMP_scatterPlot_EOverTruthVsPhi_->Write();
 
   _TEMP_scatterPlot_deltaEtaVsEt_->Write();
   _TEMP_scatterPlot_deltaEtaVsE_->Write();
   _TEMP_scatterPlot_deltaEtaVsEta_->Write();
   _TEMP_scatterPlot_deltaEtaVsPhi_->Write();
 
   _TEMP_scatterPlot_deltaPhiVsEt_->Write();
   _TEMP_scatterPlot_deltaPhiVsE_->Write();
   _TEMP_scatterPlot_deltaPhiVsEta_->Write();
   _TEMP_scatterPlot_deltaPhiVsPhi_->Write();
 
   rootFile_->cd();
 }

◆ setDrawOptions()

void EgammaObjects::setDrawOptions ( )

private

Definition at line 1128 of file EgammaObjects.cc.

Referenced by endJob().

                                    {
   hist_Et_->SetOption("e");
   hist_EtOverTruth_->SetOption("e");
   hist_EtEfficiency_->SetOption("e");
   hist_EtNumRecoOverNumTrue_->SetOption("e");
   hist_EtOverTruthVsEt_->SetOption("e");
   hist_EtOverTruthVsE_->SetOption("e");
   hist_EtOverTruthVsEta_->SetOption("e");
   hist_EtOverTruthVsPhi_->SetOption("e");
   hist_resolutionEtVsEt_->SetOption("e");
   hist_resolutionEtVsE_->SetOption("e");
   hist_resolutionEtVsEta_->SetOption("e");
   hist_resolutionEtVsPhi_->SetOption("e");
 
   hist_E_->SetOption("e");
   hist_EOverTruth_->SetOption("e");
   hist_EEfficiency_->SetOption("e");
   hist_ENumRecoOverNumTrue_->SetOption("e");
   hist_EOverTruthVsEt_->SetOption("e");
   hist_EOverTruthVsE_->SetOption("e");
   hist_EOverTruthVsEta_->SetOption("e");
   hist_EOverTruthVsPhi_->SetOption("e");
   hist_resolutionEVsEt_->SetOption("e");
   hist_resolutionEVsE_->SetOption("e");
   hist_resolutionEVsEta_->SetOption("e");
   hist_resolutionEVsPhi_->SetOption("e");
 
   hist_Eta_->SetOption("e");
   hist_EtaOverTruth_->SetOption("e");
   hist_EtaEfficiency_->SetOption("e");
   hist_EtaNumRecoOverNumTrue_->SetOption("e");
   hist_deltaEtaVsEt_->SetOption("e");
   hist_deltaEtaVsE_->SetOption("e");
   hist_deltaEtaVsEta_->SetOption("e");
   hist_deltaEtaVsPhi_->SetOption("e");
   hist_resolutionEtaVsEt_->SetOption("e");
   hist_resolutionEtaVsE_->SetOption("e");
   hist_resolutionEtaVsEta_->SetOption("e");
   hist_resolutionEtaVsPhi_->SetOption("e");
 
   hist_Phi_->SetOption("e");
   hist_PhiOverTruth_->SetOption("e");
   hist_PhiEfficiency_->SetOption("e");
   hist_PhiNumRecoOverNumTrue_->SetOption("e");
   hist_deltaPhiVsEt_->SetOption("e");
   hist_deltaPhiVsE_->SetOption("e");
   hist_deltaPhiVsEta_->SetOption("e");
   hist_deltaPhiVsPhi_->SetOption("e");
   hist_resolutionPhiVsEt_->SetOption("e");
   hist_resolutionPhiVsE_->SetOption("e");
   hist_resolutionPhiVsEta_->SetOption("e");
   hist_resolutionPhiVsPhi_->SetOption("e");
 
   hist_All_recoMass_->SetOption("e");
   hist_BarrelOnly_recoMass_->SetOption("e");
   hist_EndcapOnly_recoMass_->SetOption("e");
   hist_Mixed_recoMass_->SetOption("e");
   hist_recoMass_withBackgroud_NoEtCut_->SetOption("e");
   hist_recoMass_withBackgroud_5EtCut_->SetOption("e");
   hist_recoMass_withBackgroud_10EtCut_->SetOption("e");
   hist_recoMass_withBackgroud_20EtCut_->SetOption("e");
 }

Member Data Documentation

◆ _TEMP_scatterPlot_deltaEtaVsE_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaEtaVsE_

private

Definition at line 162 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaEtaVsEt_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaEtaVsEt_

private

Definition at line 161 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaEtaVsEta_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaEtaVsEta_

private

Definition at line 163 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaEtaVsPhi_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaEtaVsPhi_

private

Definition at line 164 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaPhiVsE_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaPhiVsE_

private

Definition at line 167 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaPhiVsEt_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaPhiVsEt_

private

Definition at line 166 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaPhiVsEta_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaPhiVsEta_

private

Definition at line 168 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_deltaPhiVsPhi_

TH2D* EgammaObjects::_TEMP_scatterPlot_deltaPhiVsPhi_

private

Definition at line 169 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EOverTruthVsE_

TH2D* EgammaObjects::_TEMP_scatterPlot_EOverTruthVsE_

private

Definition at line 157 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EOverTruthVsEt_

TH2D* EgammaObjects::_TEMP_scatterPlot_EOverTruthVsEt_

private

Definition at line 156 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EOverTruthVsEta_

TH2D* EgammaObjects::_TEMP_scatterPlot_EOverTruthVsEta_

private

Definition at line 158 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EOverTruthVsPhi_

TH2D* EgammaObjects::_TEMP_scatterPlot_EOverTruthVsPhi_

private

Definition at line 159 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EtOverTruthVsE_

TH2D* EgammaObjects::_TEMP_scatterPlot_EtOverTruthVsE_

private

Definition at line 152 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EtOverTruthVsEt_

TH2D* EgammaObjects::_TEMP_scatterPlot_EtOverTruthVsEt_

private

Definition at line 151 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EtOverTruthVsEta_

TH2D* EgammaObjects::_TEMP_scatterPlot_EtOverTruthVsEta_

private

Definition at line 153 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ _TEMP_scatterPlot_EtOverTruthVsPhi_

TH2D* EgammaObjects::_TEMP_scatterPlot_EtOverTruthVsPhi_

private

Definition at line 154 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), createTempHistoObjects(), EgammaObjects(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ EtCut

int EgammaObjects::EtCut

private

Definition at line 43 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), and EgammaObjects().

◆ hist_All_recoMass_

TH1D* EgammaObjects::hist_All_recoMass_

private

Definition at line 141 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_BarrelOnly_recoMass_

TH1D* EgammaObjects::hist_BarrelOnly_recoMass_

private

Definition at line 142 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_bins_deltaEta_

int EgammaObjects::hist_bins_deltaEta_

private

Definition at line 79 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

◆ hist_bins_deltaPhi_

int EgammaObjects::hist_bins_deltaPhi_

private

Definition at line 83 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

◆ hist_bins_E_

int EgammaObjects::hist_bins_E_

private

Definition at line 51 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), getDeltaResHistosViaSlicing(), and loadHistoParameters().

◆ hist_bins_EOverTruth_

int EgammaObjects::hist_bins_EOverTruth_

private

Definition at line 67 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_bins_Et_

int EgammaObjects::hist_bins_Et_

private

Definition at line 47 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), getDeltaResHistosViaSlicing(), and loadHistoParameters().

◆ hist_bins_Eta_

int EgammaObjects::hist_bins_Eta_

private

Definition at line 55 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), getDeltaResHistosViaSlicing(), and loadHistoParameters().

◆ hist_bins_EtaOverTruth_

int EgammaObjects::hist_bins_EtaOverTruth_

private

Definition at line 71 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_bins_EtOverTruth_

int EgammaObjects::hist_bins_EtOverTruth_

private

Definition at line 63 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_bins_Phi_

int EgammaObjects::hist_bins_Phi_

private

Definition at line 59 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), getDeltaResHistosViaSlicing(), and loadHistoParameters().

◆ hist_bins_PhiOverTruth_

int EgammaObjects::hist_bins_PhiOverTruth_

private

Definition at line 75 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_bins_recoMass_

int EgammaObjects::hist_bins_recoMass_

private

Definition at line 87 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_deltaEtaVsE_

TH1D* EgammaObjects::hist_deltaEtaVsE_

private

Definition at line 120 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaEtaVsEt_

TH1D* EgammaObjects::hist_deltaEtaVsEt_

private

Definition at line 119 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaEtaVsEta_

TH1D* EgammaObjects::hist_deltaEtaVsEta_

private

Definition at line 121 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaEtaVsPhi_

TH1D* EgammaObjects::hist_deltaEtaVsPhi_

private

Definition at line 122 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaPhiVsE_

TH1D* EgammaObjects::hist_deltaPhiVsE_

private

Definition at line 133 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaPhiVsEt_

TH1D* EgammaObjects::hist_deltaPhiVsEt_

private

Definition at line 132 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaPhiVsEta_

TH1D* EgammaObjects::hist_deltaPhiVsEta_

private

Definition at line 134 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_deltaPhiVsPhi_

TH1D* EgammaObjects::hist_deltaPhiVsPhi_

private

Definition at line 135 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_E_

TH1D* EgammaObjects::hist_E_

private

Definition at line 102 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EEfficiency_

TH1D* EgammaObjects::hist_EEfficiency_

private

Definition at line 104 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EndcapOnly_recoMass_

TH1D* EgammaObjects::hist_EndcapOnly_recoMass_

private

Definition at line 143 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_ENumRecoOverNumTrue_

TH1D* EgammaObjects::hist_ENumRecoOverNumTrue_

private

Definition at line 105 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EOverTruth_

TH1D* EgammaObjects::hist_EOverTruth_

private

Definition at line 103 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), fitHistos(), saveHistos(), and setDrawOptions().

◆ hist_EOverTruthVsE_

TH1D* EgammaObjects::hist_EOverTruthVsE_

private

Definition at line 107 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_EOverTruthVsEt_

TH1D* EgammaObjects::hist_EOverTruthVsEt_

private

Definition at line 106 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_EOverTruthVsEta_

TH1D* EgammaObjects::hist_EOverTruthVsEta_

private

Definition at line 108 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_EOverTruthVsPhi_

TH1D* EgammaObjects::hist_EOverTruthVsPhi_

private

Definition at line 109 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_Et_

TH1D* EgammaObjects::hist_Et_

private

Definition at line 89 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_Eta_

TH1D* EgammaObjects::hist_Eta_

private

Definition at line 115 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EtaEfficiency_

TH1D* EgammaObjects::hist_EtaEfficiency_

private

Definition at line 117 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EtaNumRecoOverNumTrue_

TH1D* EgammaObjects::hist_EtaNumRecoOverNumTrue_

private

Definition at line 118 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EtaOverTruth_

TH1D* EgammaObjects::hist_EtaOverTruth_

private

Definition at line 116 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), fitHistos(), saveHistos(), and setDrawOptions().

◆ hist_EtEfficiency_

TH1D* EgammaObjects::hist_EtEfficiency_

private

Definition at line 91 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EtNumRecoOverNumTrue_

TH1D* EgammaObjects::hist_EtNumRecoOverNumTrue_

private

Definition at line 92 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_EtOverTruth_

TH1D* EgammaObjects::hist_EtOverTruth_

private

Definition at line 90 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), fitHistos(), saveHistos(), and setDrawOptions().

◆ hist_EtOverTruthVsE_

TH1D* EgammaObjects::hist_EtOverTruthVsE_

private

Definition at line 94 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_EtOverTruthVsEt_

TH1D* EgammaObjects::hist_EtOverTruthVsEt_

private

Definition at line 93 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_EtOverTruthVsEta_

TH1D* EgammaObjects::hist_EtOverTruthVsEta_

private

Definition at line 95 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_EtOverTruthVsPhi_

TH1D* EgammaObjects::hist_EtOverTruthVsPhi_

private

Definition at line 96 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_max_deltaEta_

double EgammaObjects::hist_max_deltaEta_

private

Definition at line 78 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_deltaPhi_

double EgammaObjects::hist_max_deltaPhi_

private

Definition at line 82 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_E_

double EgammaObjects::hist_max_E_

private

Definition at line 50 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_EOverTruth_

double EgammaObjects::hist_max_EOverTruth_

private

Definition at line 66 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_Et_

double EgammaObjects::hist_max_Et_

private

Definition at line 46 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_Eta_

double EgammaObjects::hist_max_Eta_

private

Definition at line 54 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_EtaOverTruth_

double EgammaObjects::hist_max_EtaOverTruth_

private

Definition at line 70 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_max_EtOverTruth_

double EgammaObjects::hist_max_EtOverTruth_

private

Definition at line 62 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_Phi_

double EgammaObjects::hist_max_Phi_

private

Definition at line 58 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_max_PhiOverTruth_

double EgammaObjects::hist_max_PhiOverTruth_

private

Definition at line 74 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_max_recoMass_

double EgammaObjects::hist_max_recoMass_

private

Definition at line 86 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_min_deltaEta_

double EgammaObjects::hist_min_deltaEta_

private

Definition at line 77 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_deltaPhi_

double EgammaObjects::hist_min_deltaPhi_

private

Definition at line 81 of file EgammaObjects.h.

Referenced by createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_E_

double EgammaObjects::hist_min_E_

private

Definition at line 49 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_EOverTruth_

double EgammaObjects::hist_min_EOverTruth_

private

Definition at line 65 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_Et_

double EgammaObjects::hist_min_Et_

private

Definition at line 45 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_Eta_

double EgammaObjects::hist_min_Eta_

private

Definition at line 53 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_EtaOverTruth_

double EgammaObjects::hist_min_EtaOverTruth_

private

Definition at line 69 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_min_EtOverTruth_

double EgammaObjects::hist_min_EtOverTruth_

private

Definition at line 61 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_Phi_

double EgammaObjects::hist_min_Phi_

private

Definition at line 57 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), createTempHistoObjects(), and loadHistoParameters().

◆ hist_min_PhiOverTruth_

double EgammaObjects::hist_min_PhiOverTruth_

private

Definition at line 73 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_min_recoMass_

double EgammaObjects::hist_min_recoMass_

private

Definition at line 85 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), and loadHistoParameters().

◆ hist_Mixed_recoMass_

TH1D* EgammaObjects::hist_Mixed_recoMass_

private

Definition at line 144 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_Phi_

TH1D* EgammaObjects::hist_Phi_

private

Definition at line 128 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_PhiEfficiency_

TH1D* EgammaObjects::hist_PhiEfficiency_

private

Definition at line 130 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_PhiNumRecoOverNumTrue_

TH1D* EgammaObjects::hist_PhiNumRecoOverNumTrue_

private

Definition at line 131 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), getEfficiencyHistosViaDividing(), saveHistos(), and setDrawOptions().

◆ hist_PhiOverTruth_

TH1D* EgammaObjects::hist_PhiOverTruth_

private

Definition at line 129 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), fitHistos(), saveHistos(), and setDrawOptions().

◆ hist_recoMass_withBackgroud_10EtCut_

TH1D* EgammaObjects::hist_recoMass_withBackgroud_10EtCut_

private

Definition at line 148 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_recoMass_withBackgroud_20EtCut_

TH1D* EgammaObjects::hist_recoMass_withBackgroud_20EtCut_

private

Definition at line 149 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_recoMass_withBackgroud_5EtCut_

TH1D* EgammaObjects::hist_recoMass_withBackgroud_5EtCut_

private

Definition at line 147 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_recoMass_withBackgroud_NoEtCut_

TH1D* EgammaObjects::hist_recoMass_withBackgroud_NoEtCut_

private

Definition at line 146 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtaVsE_

TH1D* EgammaObjects::hist_resolutionEtaVsE_

private

Definition at line 124 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtaVsEt_

TH1D* EgammaObjects::hist_resolutionEtaVsEt_

private

Definition at line 123 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtaVsEta_

TH1D* EgammaObjects::hist_resolutionEtaVsEta_

private

Definition at line 125 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtaVsPhi_

TH1D* EgammaObjects::hist_resolutionEtaVsPhi_

private

Definition at line 126 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtVsE_

TH1D* EgammaObjects::hist_resolutionEtVsE_

private

Definition at line 98 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtVsEt_

TH1D* EgammaObjects::hist_resolutionEtVsEt_

private

Definition at line 97 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtVsEta_

TH1D* EgammaObjects::hist_resolutionEtVsEta_

private

Definition at line 99 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEtVsPhi_

TH1D* EgammaObjects::hist_resolutionEtVsPhi_

private

Definition at line 100 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEVsE_

TH1D* EgammaObjects::hist_resolutionEVsE_

private

Definition at line 111 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEVsEt_

TH1D* EgammaObjects::hist_resolutionEVsEt_

private

Definition at line 110 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEVsEta_

TH1D* EgammaObjects::hist_resolutionEVsEta_

private

Definition at line 112 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionEVsPhi_

TH1D* EgammaObjects::hist_resolutionEVsPhi_

private

Definition at line 113 of file EgammaObjects.h.

Referenced by applyLabels(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionPhiVsE_

TH1D* EgammaObjects::hist_resolutionPhiVsE_

private

Definition at line 137 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionPhiVsEt_

TH1D* EgammaObjects::hist_resolutionPhiVsEt_

private

Definition at line 136 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionPhiVsEta_

TH1D* EgammaObjects::hist_resolutionPhiVsEta_

private

Definition at line 138 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ hist_resolutionPhiVsPhi_

TH1D* EgammaObjects::hist_resolutionPhiVsPhi_

private

Definition at line 139 of file EgammaObjects.h.

Referenced by applyLabels(), EgammaObjects(), getDeltaResHistosViaSlicing(), saveHistos(), and setDrawOptions().

◆ MCTruthCollectionT_

edm::EDGetTokenT<edm::HepMCProduct> EgammaObjects::MCTruthCollectionT_

private

Definition at line 37 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), and loadCMSSWObjects().

◆ particleID

int EgammaObjects::particleID

private

Definition at line 40 of file EgammaObjects.h.

Referenced by analyze(), applyLabels(), createBookedHistoObjects(), EgammaObjects(), and saveHistos().

◆ particleString

std::string EgammaObjects::particleString

private

Definition at line 41 of file EgammaObjects.h.

Referenced by createBookedHistoObjects(), EgammaObjects(), endJob(), getDeltaResHistosViaSlicing(), and saveHistos().

◆ RecoCollectionT_

edm::EDGetTokenT<reco::GsfElectronCollection> EgammaObjects::RecoCollectionT_

private

Definition at line 38 of file EgammaObjects.h.

Referenced by analyzeElectrons(), analyzePhotons(), and loadCMSSWObjects().

◆ rootFile_

TFile* EgammaObjects::rootFile_

private

Definition at line 35 of file EgammaObjects.h.

Referenced by EgammaObjects(), endJob(), saveHistos(), and ~EgammaObjects().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EgammaObjects()

◆ ~EgammaObjects()

Member Function Documentation

◆ analyze()

◆ analyzeElectrons()

◆ analyzePhotons()

◆ applyLabels()

◆ beginJob()

◆ createBookedHistoObjects()

◆ createTempHistoObjects()

◆ ecalEta()

◆ endJob()

◆ findRecoMass() [1/2]

◆ findRecoMass() [2/2]

◆ fitHistos()

◆ getDeltaResHistosViaSlicing()

◆ getEfficiencyHistosViaDividing()

◆ loadCMSSWObjects()

◆ loadHistoParameters()

◆ saveHistos()

◆ setDrawOptions()

Member Data Documentation

◆ _TEMP_scatterPlot_deltaEtaVsE_

◆ _TEMP_scatterPlot_deltaEtaVsEt_

◆ _TEMP_scatterPlot_deltaEtaVsEta_

◆ _TEMP_scatterPlot_deltaEtaVsPhi_

◆ _TEMP_scatterPlot_deltaPhiVsE_

◆ _TEMP_scatterPlot_deltaPhiVsEt_

◆ _TEMP_scatterPlot_deltaPhiVsEta_

◆ _TEMP_scatterPlot_deltaPhiVsPhi_

◆ _TEMP_scatterPlot_EOverTruthVsE_

◆ _TEMP_scatterPlot_EOverTruthVsEt_

◆ _TEMP_scatterPlot_EOverTruthVsEta_

◆ _TEMP_scatterPlot_EOverTruthVsPhi_

◆ _TEMP_scatterPlot_EtOverTruthVsE_

◆ _TEMP_scatterPlot_EtOverTruthVsEt_

◆ _TEMP_scatterPlot_EtOverTruthVsEta_

◆ _TEMP_scatterPlot_EtOverTruthVsPhi_

◆ EtCut

◆ hist_All_recoMass_

◆ hist_BarrelOnly_recoMass_

◆ hist_bins_deltaEta_

◆ hist_bins_deltaPhi_

◆ hist_bins_E_

◆ hist_bins_EOverTruth_

◆ hist_bins_Et_

◆ hist_bins_Eta_

◆ hist_bins_EtaOverTruth_

◆ hist_bins_EtOverTruth_

◆ hist_bins_Phi_

◆ hist_bins_PhiOverTruth_

◆ hist_bins_recoMass_

◆ hist_deltaEtaVsE_

◆ hist_deltaEtaVsEt_

◆ hist_deltaEtaVsEta_

◆ hist_deltaEtaVsPhi_

◆ hist_deltaPhiVsE_

◆ hist_deltaPhiVsEt_

◆ hist_deltaPhiVsEta_

◆ hist_deltaPhiVsPhi_

◆ hist_E_

◆ hist_EEfficiency_

◆ hist_EndcapOnly_recoMass_

◆ hist_ENumRecoOverNumTrue_

◆ hist_EOverTruth_

◆ hist_EOverTruthVsE_

◆ hist_EOverTruthVsEt_

◆ hist_EOverTruthVsEta_

◆ hist_EOverTruthVsPhi_

◆ hist_Et_

◆ hist_Eta_

◆ hist_EtaEfficiency_

◆ hist_EtaNumRecoOverNumTrue_

◆ hist_EtaOverTruth_