#include <SimplePhotonAnalyzer.h>

Inheritance diagram for SimplePhotonAnalyzer:

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

void	beginJob () override

void	endJob () override

	SimplePhotonAnalyzer (const edm::ParameterSet &)

	~SimplePhotonAnalyzer () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > 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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
float	etaTransformation (float a, float b)

Private Attributes
edm::InputTag	barrelEcalHits_

DQMStore *	dbe_

edm::InputTag	endcapEcalHits_

MonitorElement *	h1_deltaEta_

MonitorElement *	h1_deltaEtaSC_

MonitorElement *	h1_deltaPhi_

MonitorElement *	h1_pho_E_

MonitorElement *	h1_pho_ecalIsoBarrel_

MonitorElement *	h1_pho_ecalIsoEndcap_

MonitorElement *	h1_pho_Et_

MonitorElement *	h1_pho_Eta_

MonitorElement *	h1_pho_hcalIsoBarrel_

MonitorElement *	h1_pho_hcalIsoEndcap_

MonitorElement *	h1_pho_hOverEBarrel_

MonitorElement *	h1_pho_hOverEEndcap_

MonitorElement *	h1_pho_Phi_

MonitorElement *	h1_pho_R9Barrel_

MonitorElement *	h1_pho_R9Endcap_

MonitorElement *	h1_pho_sigmaIetaIetaBarrel_

MonitorElement *	h1_pho_sigmaIetaIetaEndcap_

MonitorElement *	h1_pho_trkIsoBarrel_

MonitorElement *	h1_pho_trkIsoEndcap_

MonitorElement *	h1_recEoverTrueEBarrel_

MonitorElement *	h1_recEoverTrueEEndcap_

MonitorElement *	h1_scEta_

std::string	mcCollection_

std::string	mcProducer_

edm::InputTag	pfEgammaCandidates_

std::string	photonCollection_

std::string	photonCollectionProducer_

float	sample_

edm::ESHandle< CaloTopology >	theCaloTopo_

std::string	valueMapPFCandPhoton_

std::string	vertexProducer_

Additional Inherited Members
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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

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

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

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

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

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: Get Photon collection from the event and make very basic histos

Author: Nancy Marinelli, U. of Notre Dame, US

Definition at line 34 of file SimplePhotonAnalyzer.h.

Member Typedef Documentation

◆ DQMStore

typedef dqm::legacy::DQMStore SimplePhotonAnalyzer::DQMStore

Definition at line 36 of file SimplePhotonAnalyzer.h.

◆ MonitorElement

typedef dqm::legacy::MonitorElement SimplePhotonAnalyzer::MonitorElement

Definition at line 37 of file SimplePhotonAnalyzer.h.

Constructor & Destructor Documentation

◆ SimplePhotonAnalyzer()

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

explicit

Definition at line 33 of file SimplePhotonAnalyzer.cc.

 {
   photonCollectionProducer_ = ps.getParameter<std::string>("phoProducer");
   photonCollection_ = ps.getParameter<std::string>("photonCollection");
  
   barrelEcalHits_ = ps.getParameter<edm::InputTag>("barrelEcalHits");
   endcapEcalHits_ = ps.getParameter<edm::InputTag>("endcapEcalHits");
  
   pfEgammaCandidates_ = ps.getParameter<edm::InputTag>("pfEgammaCandidates");
   valueMapPFCandPhoton_ = ps.getParameter<std::string>("valueMapPhotons");
  
   mcProducer_ = ps.getParameter<std::string>("mcProducer");
   //mcCollection_ = ps.getParameter<std::string>("mcCollection");
   vertexProducer_ = ps.getParameter<std::string>("primaryVertexProducer");
   sample_ = ps.getParameter<int>("sample");
 }

References barrelEcalHits_, endcapEcalHits_, edm::ParameterSet::getParameter(), mcProducer_, pfEgammaCandidates_, photonCollection_, photonCollectionProducer_, sample_, AlCaHLTBitMon_QueryRunRegistry::string, valueMapPFCandPhoton_, and vertexProducer_.

◆ ~SimplePhotonAnalyzer()

SimplePhotonAnalyzer::~SimplePhotonAnalyzer ( )

override

Definition at line 52 of file SimplePhotonAnalyzer.cc.

54 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Match reconstructed photon candidates with the nearest generated photonPho;

Plot kinematic disctributions for matched photons

Implements edm::one::EDAnalyzerBase.

Definition at line 137 of file SimplePhotonAnalyzer.cc.

                                                                                {
   //========================================================================
  
   using namespace edm;  // needed for all fwk related classes
   edm::LogInfo("PhotonAnalyzer") << "Analyzing event number: " << evt.id() << "\n";
  
   // get the  calo topology  from the event setup:
   edm::ESHandle<CaloTopology> pTopology;
   es.get<CaloTopologyRecord>().get(theCaloTopo_);
  
   // Get the  corrected  photon collection (set in the configuration) which also contains infos about conversions
  
   Handle<reco::PhotonCollection> photonHandle;
   evt.getByLabel(photonCollectionProducer_, photonCollection_, photonHandle);
   const reco::PhotonCollection photonCollection = *(photonHandle.product());
  
   Handle<HepMCProduct> hepProd;
   evt.getByLabel(mcProducer_, hepProd);
   const HepMC::GenEvent* myGenEvent = hepProd->GetEvent();
  
   // Get the  PF refined cluster  collection
   Handle<reco::PFCandidateCollection> pfCandidateHandle;
   evt.getByLabel(pfEgammaCandidates_, pfCandidateHandle);
   if (!pfCandidateHandle.isValid()) {
     edm::LogError("SimplePhotonAnalyzer") << "Error! Can't get the product " << pfEgammaCandidates_.label();
   }
  
   edm::Handle<edm::ValueMap<reco::PhotonRef> > pfCandToPhotonMapHandle;
   edm::ValueMap<reco::PhotonRef> pfCandToPhotonMap;
   evt.getByLabel("gedPhotons", valueMapPFCandPhoton_, pfCandToPhotonMapHandle);
   if (!pfCandToPhotonMapHandle.isValid()) {
     edm::LogInfo("SimplePhotonAnalyzer") << "Error! Can't get the product: valueMapPhotons " << std::endl;
   }
   pfCandToPhotonMap = *(pfCandToPhotonMapHandle.product());
  
   std::cout << " SimplePhotonAnalyzer  valueMap size" << pfCandToPhotonMap.size() << std::endl;
   unsigned nObj = pfCandidateHandle->size();
   for (unsigned int lCand = 0; lCand < nObj; lCand++) {
     reco::PFCandidateRef pfCandRef(reco::PFCandidateRef(pfCandidateHandle, lCand));
     if (pfCandRef->particleId() != reco::PFCandidate::gamma)
       continue;
     reco::PhotonRef myPho = (pfCandToPhotonMap)[pfCandRef];
     if (myPho.isNonnull())
       std::cout << " PF SC " << pfCandRef->superClusterRef()->energy() << " Photon SC "
                 << myPho->superCluster()->energy() << std::endl;
   }
  
   for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
        ++p) {
     if (!((*p)->pdg_id() == 22 && (*p)->status() == 1))
       continue;
  
     // single primary photons or photons from Higgs or RS Graviton
     HepMC::GenParticle* mother = nullptr;
     if ((*p)->production_vertex()) {
       if ((*p)->production_vertex()->particles_begin(HepMC::parents) !=
           (*p)->production_vertex()->particles_end(HepMC::parents))
         mother = *((*p)->production_vertex()->particles_begin(HepMC::parents));
     }
     if (((mother == nullptr) || ((mother != nullptr) && (mother->pdg_id() == 25)) ||
          ((mother != nullptr) && (mother->pdg_id() == 22)))) {
       float minDelta = 10000.;
       std::vector<reco::Photon> localPhotons;
       int index = 0;
       int iMatch = -1;
  
       float phiPho = (*p)->momentum().phi();
       float etaPho = (*p)->momentum().eta();
       etaPho = etaTransformation(etaPho, (*p)->production_vertex()->position().z() / 10.);
  
       bool matched = false;
       // loop  Photon candidates
       for (reco::PhotonCollection::const_iterator iPho = photonCollection.begin(); iPho != photonCollection.end();
            iPho++) {
         reco::Photon localPho = reco::Photon(*iPho);
         localPhotons.push_back(localPho);
  
         float phiClu = localPho.phi();
         float etaClu = localPho.eta();
         float deltaPhi = phiClu - phiPho;
         float deltaEta = etaClu - etaPho;
  
         if (deltaPhi > pi)
           deltaPhi -= twopi;
         if (deltaPhi < -pi)
           deltaPhi += twopi;
         deltaPhi = std::pow(deltaPhi, 2);
         deltaEta = std::pow(deltaEta, 2);
         float delta = sqrt(deltaPhi + deltaEta);
         if (delta < 0.1 && delta < minDelta) {
           minDelta = delta;
           iMatch = index;
         }
         index++;
       }  // End loop over photons
  
       if (iMatch > -1)
         matched = true;
  
       if (matched) {
         reco::Photon matchingPho = localPhotons[iMatch];
  
         bool phoIsInBarrel = false;
         if (fabs(matchingPho.superCluster()->position().eta()) < 1.479) {
           phoIsInBarrel = true;
         }
         edm::Handle<EcalRecHitCollection> ecalRecHitHandle;
  
         h1_scEta_->Fill(matchingPho.superCluster()->position().eta());
         float trueEta = (*p)->momentum().eta();
         trueEta = etaTransformation(trueEta, (*p)->production_vertex()->position().z() / 10.);
         h1_deltaEtaSC_->Fill(localPhotons[iMatch].superCluster()->eta() - trueEta);
  
         float photonE = matchingPho.energy();
         float photonEt = matchingPho.et();
         float photonEta = matchingPho.eta();
         float photonPhi = matchingPho.phi();
  
         float r9 = matchingPho.r9();
         float sigmaIetaIeta = matchingPho.sigmaIetaIeta();
         float hOverE = matchingPho.hadronicOverEm();
         float ecalIso = matchingPho.ecalRecHitSumEtConeDR04();
         float hcalIso = matchingPho.hcalTowerSumEtConeDR04();
         float trkIso = matchingPho.trkSumPtSolidConeDR04();
  
         h1_pho_E_->Fill(photonE);
         h1_pho_Et_->Fill(photonEt);
         h1_pho_Eta_->Fill(photonEta);
         h1_pho_Phi_->Fill(photonPhi);
  
         h1_deltaEta_->Fill(photonEta - (*p)->momentum().eta());
         h1_deltaPhi_->Fill(photonPhi - (*p)->momentum().phi());
  
         if (phoIsInBarrel) {
           h1_recEoverTrueEBarrel_->Fill(photonE / (*p)->momentum().e());
           h1_pho_R9Barrel_->Fill(r9);
           h1_pho_sigmaIetaIetaBarrel_->Fill(sigmaIetaIeta);
           h1_pho_hOverEBarrel_->Fill(hOverE);
           h1_pho_ecalIsoBarrel_->Fill(ecalIso);
           h1_pho_hcalIsoBarrel_->Fill(hcalIso);
           h1_pho_trkIsoBarrel_->Fill(trkIso);
  
         } else {
           h1_recEoverTrueEEndcap_->Fill(photonE / (*p)->momentum().e());
           h1_pho_R9Endcap_->Fill(r9);
           h1_pho_sigmaIetaIetaEndcap_->Fill(sigmaIetaIeta);
           h1_pho_hOverEEndcap_->Fill(hOverE);
           h1_pho_ecalIsoEndcap_->Fill(ecalIso);
           h1_pho_hcalIsoEndcap_->Fill(hcalIso);
           h1_pho_trkIsoEndcap_->Fill(trkIso);
         }
  
       }  //  reco photon matching MC truth
  
     }  // End loop over MC particles
   }
 }

References gather_cfg::cout, dumpMFGeometry_cfg::delta, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::Photon::ecalRecHitSumEtConeDR04(), reco::LeafCandidate::energy(), reco::LeafCandidate::et(), PVValHelper::eta, reco::LeafCandidate::eta(), etaTransformation(), dqm::impl::MonitorElement::Fill(), reco::PFCandidate::gamma, GenParticle::GenParticle, edm::EventSetup::get(), get, edm::Event::getByLabel(), edm::HepMCProduct::GetEvent(), h1_deltaEta_, h1_deltaEtaSC_, h1_deltaPhi_, h1_pho_E_, h1_pho_ecalIsoBarrel_, h1_pho_ecalIsoEndcap_, h1_pho_Et_, h1_pho_Eta_, h1_pho_hcalIsoBarrel_, h1_pho_hcalIsoEndcap_, h1_pho_hOverEBarrel_, h1_pho_hOverEEndcap_, h1_pho_Phi_, h1_pho_R9Barrel_, h1_pho_R9Endcap_, h1_pho_sigmaIetaIetaBarrel_, h1_pho_sigmaIetaIetaEndcap_, h1_pho_trkIsoBarrel_, h1_pho_trkIsoEndcap_, h1_recEoverTrueEBarrel_, h1_recEoverTrueEEndcap_, h1_scEta_, reco::Photon::hadronicOverEm(), reco::Photon::hcalTowerSumEtConeDR04(), EgHLTOffHistBins_cfi::hOverE, edm::EventBase::id(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), edm::InputTag::label(), muonTagProbeFilters_cff::matched, mcProducer_, AlCaHLTBitMon_ParallelJobs::p, parents, pfEgammaCandidates_, reco::LeafCandidate::phi(), nanoDQM_cff::Photon, ExoticaDQM_cfi::photonCollection, photonCollection_, photonCollectionProducer_, HLTObjectMonitor_cfi::photonEta, HLTObjectMonitor_cfi::photonPhi, pi, funct::pow(), edm::Handle< T >::product(), reco::Photon::r9(), electrons_cff::r9, reco::Photon::sigmaIetaIeta(), edm::ValueMap< T >::size(), mathSSE::sqrt(), reco::Photon::superCluster(), theCaloTopo_, ZMuMuCategoriesSequences_cff::trkIso, reco::Photon::trkSumPtSolidConeDR04(), and valueMapPFCandPhoton_.

◆ beginJob()

void SimplePhotonAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 57 of file SimplePhotonAnalyzer.cc.

                                     {
   //========================================================================
  
   dbe_ = nullptr;
   dbe_ = edm::Service<DQMStore>().operator->();
  
   float hiE = 0;
   float loE = 0;
   float hiEt = 0;
   float loEt = 0;
   float dPhi = 0;
   float loRes = 0;
   float hiRes = 0;
   if (sample_ == 1) {
     loE = 0.;
     hiE = 30.;
     loEt = 0.;
     hiEt = 30.;
     dPhi = 0.2;
     loRes = 0.;
     hiRes = 1.2;
   } else if (sample_ == 2) {
     loE = 0.;
     hiE = 200.;
     loEt = 0.;
     hiEt = 50.;
     dPhi = 0.05;
     loRes = 0.7;
     hiRes = 1.2;
   } else if (sample_ == 3) {
     loE = 0.;
     hiE = 500.;
     loEt = 0.;
     hiEt = 500.;
     dPhi = 0.05;
     loRes = 0.7;
     hiRes = 1.2;
   } else if (sample_ == 4) {
     loE = 0.;
     hiE = 6000.;
     loEt = 0.;
     hiEt = 1200.;
     dPhi = 0.05;
     loRes = 0.7;
     hiRes = 1.2;
   }
  
   h1_deltaEta_ = dbe_->book1D("deltaEta", " Reco photon Eta minus Generated photon Eta  ", 100, -0.2, 0.2);
   h1_deltaPhi_ = dbe_->book1D("deltaPhi", "Reco photon Phi minus Generated photon Phi ", 100, -dPhi, dPhi);
   h1_pho_Eta_ = dbe_->book1D("phoEta", "Photon  Eta ", 40, -3., 3.);
   h1_pho_Phi_ = dbe_->book1D("phoPhi", "Photon  Phi ", 40, -3.14, 3.14);
   h1_pho_E_ = dbe_->book1D("phoE", "Photon Energy ", 100, loE, hiE);
   h1_pho_Et_ = dbe_->book1D("phoEt", "Photon Et ", 100, loEt, hiEt);
  
   h1_scEta_ = dbe_->book1D("scEta", " SC Eta ", 40, -3., 3.);
   h1_deltaEtaSC_ = dbe_->book1D("deltaEtaSC", " SC Eta minus Generated photon Eta  ", 100, -0.02, 0.02);
  
   //
   h1_recEoverTrueEBarrel_ = dbe_->book1D(
       "recEoverTrueEBarrel", " Reco photon Energy over Generated photon Energy: Barrel ", 100, loRes, hiRes);
   h1_recEoverTrueEEndcap_ = dbe_->book1D(
       "recEoverTrueEEndcap", " Reco photon Energy over Generated photon Energy: Endcap ", 100, loRes, hiRes);
  
   //
  
   h1_pho_R9Barrel_ = dbe_->book1D("phoR9Barrel", "Photon  3x3 energy / SuperCluster energy : Barrel ", 100, 0., 1.2);
   h1_pho_R9Endcap_ = dbe_->book1D("phoR9Endcap", "Photon  3x3 energy / SuperCluster energy : Endcap ", 100, 0., 1.2);
   h1_pho_sigmaIetaIetaBarrel_ = dbe_->book1D("sigmaIetaIetaBarrel", "sigmaIetaIeta: Barrel", 100, 0., 0.05);
   h1_pho_sigmaIetaIetaEndcap_ = dbe_->book1D("sigmaIetaIetaEndcap", "sigmaIetaIeta: Endcap", 100, 0., 0.1);
   h1_pho_hOverEBarrel_ = dbe_->book1D("hOverEBarrel", "H/E: Barrel", 100, 0., 0.1);
   h1_pho_hOverEEndcap_ = dbe_->book1D("hOverEEndcap", "H/E: Endcap", 100, 0., 0.1);
   h1_pho_ecalIsoBarrel_ = dbe_->book1D("ecalIsolBarrel", "isolation et sum in Ecal: Barrel", 100, 0., 100.);
   h1_pho_ecalIsoEndcap_ = dbe_->book1D("ecalIsolEndcap", "isolation et sum in Ecal: Endcap", 100, 0., 100.);
   h1_pho_hcalIsoBarrel_ = dbe_->book1D("hcalIsolBarrel", "isolation et sum in Hcal: Barrel", 100, 0., 100.);
   h1_pho_hcalIsoEndcap_ = dbe_->book1D("hcalIsolEndcap", "isolation et sum in Hcal: Endcap", 100, 0., 100.);
   h1_pho_trkIsoBarrel_ = dbe_->book1D("trkIsolBarrel", "isolation pt sum in the tracker: Barrel", 100, 0., 100.);
   h1_pho_trkIsoEndcap_ = dbe_->book1D("trkIsolEndcap", "isolation pt sum in the tracker: Endcap", 100, 0., 100.);
 }

References dqm::implementation::IBooker::book1D(), dbe_, HLT_FULL_cff::dPhi, h1_deltaEta_, h1_deltaEtaSC_, h1_deltaPhi_, h1_pho_E_, h1_pho_ecalIsoBarrel_, h1_pho_ecalIsoEndcap_, h1_pho_Et_, h1_pho_Eta_, h1_pho_hcalIsoBarrel_, h1_pho_hcalIsoEndcap_, h1_pho_hOverEBarrel_, h1_pho_hOverEEndcap_, h1_pho_Phi_, h1_pho_R9Barrel_, h1_pho_R9Endcap_, h1_pho_sigmaIetaIetaBarrel_, h1_pho_sigmaIetaIetaEndcap_, h1_pho_trkIsoBarrel_, h1_pho_trkIsoEndcap_, h1_recEoverTrueEBarrel_, h1_recEoverTrueEEndcap_, h1_scEta_, SiPixelOfflineDQM_source_cff::hiRes, and sample_.

◆ endJob()

void SimplePhotonAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 335 of file SimplePhotonAnalyzer.cc.

                                   {
   //========================================================================
 }

◆ etaTransformation()

float SimplePhotonAnalyzer::etaTransformation	(	float	a,
		float	b
	)

private

Definition at line 297 of file SimplePhotonAnalyzer.cc.

                                                                               {
   //---Definitions
   const float PI = 3.1415927;
   //UNUSED const float TWOPI = 2.0*PI;
  
   //---Definitions for ECAL
   const float R_ECAL = 136.5;
   const float Z_Endcap = 328.0;
   const float etaBarrelEndcap = 1.479;
  
   //---ETA correction
  
   float Theta = 0.0;
   float ZEcal = R_ECAL * sinh(EtaParticle) + Zvertex;
  
   if (ZEcal != 0.0)
     Theta = atan(R_ECAL / ZEcal);
   if (Theta < 0.0)
     Theta = Theta + PI;
   float ETA = -log(tan(0.5 * Theta));
  
   if (fabs(ETA) > etaBarrelEndcap) {
     float Zend = Z_Endcap;
     if (EtaParticle < 0.0)
       Zend = -Zend;
     float Zlen = Zend - Zvertex;
     float RR = Zlen / sinh(EtaParticle);
     Theta = atan(RR / Zend);
     if (Theta < 0.0)
       Theta = Theta + PI;
     ETA = -log(tan(0.5 * Theta));
   }
   //---Return the result
   return ETA;
   //---end
 }

References ETA, etaBarrelEndcap, dqm-mbProfile::log, PI, R_ECAL, funct::tan(), Z_Endcap, and ZEcal.

Referenced by analyze().

Member Data Documentation

◆ barrelEcalHits_

edm::InputTag SimplePhotonAnalyzer::barrelEcalHits_

private

Definition at line 55 of file SimplePhotonAnalyzer.h.

Referenced by SimplePhotonAnalyzer().

◆ dbe_

DQMStore* SimplePhotonAnalyzer::dbe_

private

Definition at line 63 of file SimplePhotonAnalyzer.h.

Referenced by beginJob().

◆ endcapEcalHits_

edm::InputTag SimplePhotonAnalyzer::endcapEcalHits_

private

Definition at line 56 of file SimplePhotonAnalyzer.h.

Referenced by SimplePhotonAnalyzer().

◆ h1_deltaEta_

MonitorElement* SimplePhotonAnalyzer::h1_deltaEta_

private

Definition at line 86 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_deltaEtaSC_

MonitorElement* SimplePhotonAnalyzer::h1_deltaEtaSC_

private

Definition at line 66 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_deltaPhi_

MonitorElement* SimplePhotonAnalyzer::h1_deltaPhi_

private

Definition at line 87 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_E_

MonitorElement* SimplePhotonAnalyzer::h1_pho_E_

private

Definition at line 67 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_ecalIsoBarrel_

MonitorElement* SimplePhotonAnalyzer::h1_pho_ecalIsoBarrel_

private

Definition at line 77 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_ecalIsoEndcap_

MonitorElement* SimplePhotonAnalyzer::h1_pho_ecalIsoEndcap_

private

Definition at line 78 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_Et_

MonitorElement* SimplePhotonAnalyzer::h1_pho_Et_

private

Definition at line 68 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_Eta_

MonitorElement* SimplePhotonAnalyzer::h1_pho_Eta_

private

Definition at line 69 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_hcalIsoBarrel_

MonitorElement* SimplePhotonAnalyzer::h1_pho_hcalIsoBarrel_

private

Definition at line 79 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_hcalIsoEndcap_

MonitorElement* SimplePhotonAnalyzer::h1_pho_hcalIsoEndcap_

private

Definition at line 80 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_hOverEBarrel_

MonitorElement* SimplePhotonAnalyzer::h1_pho_hOverEBarrel_

private

Definition at line 75 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_hOverEEndcap_

MonitorElement* SimplePhotonAnalyzer::h1_pho_hOverEEndcap_

private

Definition at line 76 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_Phi_

MonitorElement* SimplePhotonAnalyzer::h1_pho_Phi_

private

Definition at line 70 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_R9Barrel_

MonitorElement* SimplePhotonAnalyzer::h1_pho_R9Barrel_

private

Definition at line 71 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_R9Endcap_

MonitorElement* SimplePhotonAnalyzer::h1_pho_R9Endcap_

private

Definition at line 72 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_sigmaIetaIetaBarrel_

MonitorElement* SimplePhotonAnalyzer::h1_pho_sigmaIetaIetaBarrel_

private

Definition at line 73 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_sigmaIetaIetaEndcap_

MonitorElement* SimplePhotonAnalyzer::h1_pho_sigmaIetaIetaEndcap_

private

Definition at line 74 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_trkIsoBarrel_

MonitorElement* SimplePhotonAnalyzer::h1_pho_trkIsoBarrel_

private

Definition at line 81 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_pho_trkIsoEndcap_

MonitorElement* SimplePhotonAnalyzer::h1_pho_trkIsoEndcap_

private

Definition at line 82 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_recEoverTrueEBarrel_

MonitorElement* SimplePhotonAnalyzer::h1_recEoverTrueEBarrel_

private

Definition at line 84 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_recEoverTrueEEndcap_

MonitorElement* SimplePhotonAnalyzer::h1_recEoverTrueEEndcap_

private

Definition at line 85 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h1_scEta_

MonitorElement* SimplePhotonAnalyzer::h1_scEta_

private

Definition at line 65 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ mcCollection_

std::string SimplePhotonAnalyzer::mcCollection_

private

Definition at line 50 of file SimplePhotonAnalyzer.h.

◆ mcProducer_

std::string SimplePhotonAnalyzer::mcProducer_

private

Definition at line 49 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

◆ pfEgammaCandidates_

edm::InputTag SimplePhotonAnalyzer::pfEgammaCandidates_

private

Definition at line 54 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

◆ photonCollection_

std::string SimplePhotonAnalyzer::photonCollection_

private

Definition at line 52 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

◆ photonCollectionProducer_

std::string SimplePhotonAnalyzer::photonCollectionProducer_

private

Definition at line 51 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

◆ sample_

float SimplePhotonAnalyzer::sample_

private

Definition at line 61 of file SimplePhotonAnalyzer.h.

Referenced by beginJob(), and SimplePhotonAnalyzer().

◆ theCaloTopo_

edm::ESHandle<CaloTopology> SimplePhotonAnalyzer::theCaloTopo_

private

Definition at line 58 of file SimplePhotonAnalyzer.h.

Referenced by analyze().

◆ valueMapPFCandPhoton_

std::string SimplePhotonAnalyzer::valueMapPFCandPhoton_

private

Definition at line 53 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

◆ vertexProducer_

std::string SimplePhotonAnalyzer::vertexProducer_

private

Definition at line 60 of file SimplePhotonAnalyzer.h.

Referenced by SimplePhotonAnalyzer().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DQMStore

◆ MonitorElement

Constructor & Destructor Documentation

◆ SimplePhotonAnalyzer()

◆ ~SimplePhotonAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ etaTransformation()

Member Data Documentation

◆ barrelEcalHits_

◆ dbe_

◆ endcapEcalHits_

◆ h1_deltaEta_

◆ h1_deltaEtaSC_

◆ h1_deltaPhi_

◆ h1_pho_E_

◆ h1_pho_ecalIsoBarrel_

◆ h1_pho_ecalIsoEndcap_

◆ h1_pho_Et_

◆ h1_pho_Eta_

◆ h1_pho_hcalIsoBarrel_

◆ h1_pho_hcalIsoEndcap_

◆ h1_pho_hOverEBarrel_

◆ h1_pho_hOverEEndcap_

◆ h1_pho_Phi_

◆ h1_pho_R9Barrel_

◆ h1_pho_R9Endcap_

◆ h1_pho_sigmaIetaIetaBarrel_

◆ h1_pho_sigmaIetaIetaEndcap_

◆ h1_pho_trkIsoBarrel_

◆ h1_pho_trkIsoEndcap_

◆ h1_recEoverTrueEBarrel_

◆ h1_recEoverTrueEEndcap_

◆ h1_scEta_

◆ mcCollection_

◆ mcProducer_

◆ pfEgammaCandidates_

◆ photonCollection_

◆ photonCollectionProducer_

◆ sample_

◆ theCaloTopo_

◆ valueMapPFCandPhoton_

◆ vertexProducer_