SimplePhotonAnalyzer Class Reference

#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
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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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)

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

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

Definition at line 36 of file SimplePhotonAnalyzer.h.

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

Definition at line 37 of file SimplePhotonAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file SimplePhotonAnalyzer.cc.

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

{
  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");
}

SimplePhotonAnalyzer::~SimplePhotonAnalyzer ( )

override

Definition at line 52 of file SimplePhotonAnalyzer.cc.

{}

Member Function Documentation

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

override

Match reconstructed photon candidates with the nearest generated photonPho;

Plot kinematic disctributions for matched photons

Definition at line 137 of file SimplePhotonAnalyzer.cc.

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(), 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_.

                                                                               {
  //========================================================================

  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
  }
}

void SimplePhotonAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 57 of file SimplePhotonAnalyzer.cc.

References dqm::dqmstoreimpl::DQMStore::book1D(), dbe_, HLT_2018_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, Utilities::operator, and sample_.

                                    {
  //========================================================================

  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.);
}

void SimplePhotonAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 335 of file SimplePhotonAnalyzer.cc.

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

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

private

Definition at line 297 of file SimplePhotonAnalyzer.cc.

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

Referenced by analyze().

                                                                              {
  //---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
}

Member Data Documentation

edm::InputTag SimplePhotonAnalyzer::barrelEcalHits_

private

Definition at line 55 of file SimplePhotonAnalyzer.h.

Referenced by SimplePhotonAnalyzer().

DQMStore* SimplePhotonAnalyzer::dbe_

private

Definition at line 63 of file SimplePhotonAnalyzer.h.

Referenced by beginJob().

edm::InputTag SimplePhotonAnalyzer::endcapEcalHits_

private

Definition at line 56 of file SimplePhotonAnalyzer.h.

Referenced by SimplePhotonAnalyzer().

MonitorElement* SimplePhotonAnalyzer::h1_deltaEta_

private

Definition at line 86 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_deltaEtaSC_

private

Definition at line 66 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_deltaPhi_

private

Definition at line 87 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_E_

private

Definition at line 67 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_ecalIsoBarrel_

private

Definition at line 77 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_ecalIsoEndcap_

private

Definition at line 78 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_Et_

private

Definition at line 68 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_Eta_

private

Definition at line 69 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_hcalIsoBarrel_

private

Definition at line 79 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_hcalIsoEndcap_

private

Definition at line 80 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_hOverEBarrel_

private

Definition at line 75 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_hOverEEndcap_

private

Definition at line 76 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_Phi_

private

Definition at line 70 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_R9Barrel_

private

Definition at line 71 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_R9Endcap_

private

Definition at line 72 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_sigmaIetaIetaBarrel_

private

Definition at line 73 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_sigmaIetaIetaEndcap_

private

Definition at line 74 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_trkIsoBarrel_

private

Definition at line 81 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_pho_trkIsoEndcap_

private

Definition at line 82 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_recEoverTrueEBarrel_

private

Definition at line 84 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_recEoverTrueEEndcap_

private

Definition at line 85 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SimplePhotonAnalyzer::h1_scEta_

private

Definition at line 65 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and beginJob().

std::string SimplePhotonAnalyzer::mcCollection_

private

Definition at line 50 of file SimplePhotonAnalyzer.h.

std::string SimplePhotonAnalyzer::mcProducer_

private

Definition at line 49 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

edm::InputTag SimplePhotonAnalyzer::pfEgammaCandidates_

private

Definition at line 54 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

std::string SimplePhotonAnalyzer::photonCollection_

private

Definition at line 52 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

std::string SimplePhotonAnalyzer::photonCollectionProducer_

private

Definition at line 51 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

float SimplePhotonAnalyzer::sample_

private

Definition at line 61 of file SimplePhotonAnalyzer.h.

Referenced by beginJob(), and SimplePhotonAnalyzer().

edm::ESHandle<CaloTopology> SimplePhotonAnalyzer::theCaloTopo_

private

Definition at line 58 of file SimplePhotonAnalyzer.h.

Referenced by analyze().

std::string SimplePhotonAnalyzer::valueMapPFCandPhoton_

private

Definition at line 53 of file SimplePhotonAnalyzer.h.

Referenced by analyze(), and SimplePhotonAnalyzer().

std::string SimplePhotonAnalyzer::vertexProducer_

private

Definition at line 60 of file SimplePhotonAnalyzer.h.

Referenced by SimplePhotonAnalyzer().