Reimplemented from DQMEDAnalyzer.
Definition at line 193 of file PhotonValidatorMiniAOD.cc.
References pat::Photon::chargedHadronIso(), HGC3DClusterGenMatchSelector_cfi::dR, run3scouting_cff::ecalIso, reco::Photon::ecalRecHitSumEtConeDR04(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), dqm::impl::MonitorElement::Fill(), reco::Photon::full5x5_r9(), reco::Photon::full5x5_sigmaIetaIeta(), AJJGenJetFilter_cfi::genParticles, genpartToken_, reco::Photon::getCandidateP4type(), reco::Photon::getCorrectedEnergyError(), h_chHadIso_, h_ecalRecHitSumEtConeDR04_, h_full5x5_r9_, h_full5x5_sigmaIetaIeta_, h_hcalTowerBcSumEtConeDR04_, h_hcalTowerSumEtConeDR04_, h_hOverE_, h_isoTrkSolidConeDR04_, h_newhOverE_, h_nHadIso_, h_nTrkSolidConeDR04_, h_phoE_, h_phoERes_, h_phoEt_, h_phoIso_, h_phoSigmaEoE_, h_r1_, h_r2_, h_r9_, h_scEta_, h_scPhi_, h_sigmaIetaIeta_, reco::Photon::hadronicOverEm(), reco::Photon::hadTowOverEm(), run3scouting_cff::hcalIso, reco::Photon::hcalTowerSumEtBcConeDR04(), reco::Photon::hcalTowerSumEtConeDR04(), EgHLTOffHistBins_cfi::hOverE, iEvent, muonTagProbeFilters_cff::matched, pat::Photon::neutralHadronIso(), reco::Photon::nTrkSolidConeDR04(), reco::LeafCandidate::phi(), gedPhotonSequence_cff::photonIso, pat::Photon::photonIso(), BPHMonitor_cfi::photons, photonToken_, pi, conifer::pow(), diffTwoXMLs::r1, reco::Photon::r1x5(), diffTwoXMLs::r2, reco::Photon::r2x5(), reco::Photon::r9(), electrons_cff::r9, electrons_cff::sieie, reco::Photon::sigmaIetaIeta(), mathSSE::sqrt(), pat::Photon::superCluster(), and reco::Photon::trkSumPtSolidConeDR04().
200 const auto &
photons = *photonsHandle;
207 if (!(mcIter->pdgId() == 22))
209 if (mcIter->mother() !=
nullptr && !(mcIter->mother()->pdgId() == 25))
211 if (fabs(mcIter->eta()) > 2.5)
216 float mcPhi = mcIter->phi();
217 float mcEta = mcIter->eta();
219 float mcEnergy = mcIter->energy();
221 double dR = 9999999.;
222 float minDr = 10000.;
225 for (
size_t ipho = 0; ipho <
photons.size(); ipho++) {
228 double dphi = pho->
phi() - mcPhi;
230 dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
232 double deta = pho->
superCluster()->position().eta() - mcEta;
235 if (
dR < 0.1 &&
dR < minDr) {
248 bool phoIsInBarrel =
false;
249 bool phoIsInEndcap =
false;
251 float phoEta = matchingPho->
superCluster()->position().eta();
252 if (fabs(phoEta) < 1.479) {
253 phoIsInBarrel =
true;
255 phoIsInEndcap =
true;
258 float photonE = matchingPho->
energy();
260 float photonEt = matchingPho->
energy() / cosh(matchingPho->
eta());
263 float r9 = matchingPho->
r9();
265 float r1 = matchingPho->
r1x5();
266 float r2 = matchingPho->
r2x5();
edm::EDGetTokenT< reco::GenParticleCollection > genpartToken_
float neutralHadronIso() const
float ecalRecHitSumEtConeDR04() const
MonitorElement * h_full5x5_sigmaIetaIeta_[3][3]
MonitorElement * h_r1_[3][3]
MonitorElement * h_full5x5_r9_[3][3]
MonitorElement * h_phoSigmaEoE_[3][3]
float chargedHadronIso() const
MonitorElement * h_nHadIso_[3]
MonitorElement * h_ecalRecHitSumEtConeDR04_[3][3]
P4type getCandidateP4type() const
MonitorElement * h_sigmaIetaIeta_[3][3]
float trkSumPtSolidConeDR04() const
MonitorElement * h_chHadIso_[3]
MonitorElement * h_hcalTowerBcSumEtConeDR04_[3][3]
float hcalTowerSumEtConeDR04(int depth=0) const
MonitorElement * h_newhOverE_[3][3]
float sigmaIetaIeta() const
MonitorElement * h_isoTrkSolidConeDR04_[3][3]
MonitorElement * h_phoIso_[3]
float full5x5_sigmaIetaIeta() const
float getCorrectedEnergyError(P4type type) const
MonitorElement * h_nTrkSolidConeDR04_[3][3]
float hcalTowerSumEtBcConeDR04(int depth=0) const
float hadTowOverEm(int depth=0) const
MonitorElement * h_hOverE_[3][3]
MonitorElement * h_scPhi_[2]
MonitorElement * h_phoERes_[3][3]
MonitorElement * h_phoE_[2][3]
MonitorElement * h_r9_[3][3]
int nTrkSolidConeDR04() const
float hadronicOverEm(int depth=0) const
reco::SuperClusterRef superCluster() const override
override the superCluster method from CaloJet, to access the internal storage of the supercluster ...
edm::EDGetTokenT< edm::View< pat::Photon > > photonToken_
MonitorElement * h_phoEt_[2][3]
MonitorElement * h_hcalTowerSumEtConeDR04_[3][3]
MonitorElement * h_r2_[3][3]
double phi() const final
momentum azimuthal angle
MonitorElement * h_scEta_[2]
double energy() const final
energy
double eta() const final
momentum pseudorapidity