|
|
EgammaCoreTools.
More...
#include <PhotonAnalyzer.h>
|
void | book2DHistoVector (DQMStore::IBooker &, std::vector< std::vector< MonitorElement * > > &toFill, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2) |
|
void | book3DHistoVector (DQMStore::IBooker &, std::vector< std::vector< std::vector< MonitorElement * > > > &toFill, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2) |
|
MonitorElement * | bookHisto (DQMStore::IBooker &, std::string histoName, std::string title, int bin, double min, double max) |
|
void | bookHistogramsConversions (DQMStore::IBooker &) |
|
void | bookHistogramsEfficiency (DQMStore::IBooker &) |
|
void | bookHistogramsForHistogramCounts (DQMStore::IBooker &) |
|
void | bookHistogramsInvMass (DQMStore::IBooker &) |
|
void | bookHistogramsPhotons (DQMStore::IBooker &) |
|
void | fill2DHistoVector (std::vector< std::vector< MonitorElement * > > &histoVector, double x, double y, int cut, int type) |
|
void | fill2DHistoVector (std::vector< std::vector< MonitorElement * > > &histoVector, double x, int cut, int type) |
|
void | fill3DHistoVector (std::vector< std::vector< std::vector< MonitorElement * > > > &histoVector, double x, double y, int cut, int type, int part) |
|
void | fill3DHistoVector (std::vector< std::vector< std::vector< MonitorElement * > > > &histoVector, double x, int cut, int type, int part) |
|
void | fillHistogramsForHistogramCounts (DQMStore::IBooker &) |
|
float | phiNormalization (float &a) |
|
bool | photonSelection (const reco::Photon *p) |
|
bool | photonSelectionSlimmed (const reco::Photon *p) |
|
EgammaCoreTools.
$Id: PhotonAnalyzer authors: Nancy Marinelli, U. of Notre Dame, US Jamie Antonelli, U. of Notre Dame, US
$Id: PhotonAnalyzer authors: Nancy Marinelli, U. of Notre Dame, US
Jamie Antonelli, U. of Notre Dame, US
Definition at line 97 of file PhotonAnalyzer.h.
◆ PhotonAnalyzer()
◆ ~PhotonAnalyzer()
PhotonAnalyzer::~PhotonAnalyzer |
( |
| ) |
|
|
override |
◆ analyze()
End loop over Reco photons
Reimplemented from DQMEDAnalyzer.
Definition at line 1415 of file PhotonAnalyzer.cc.
1416 using namespace edm;
1421 LogInfo(
fName_) <<
"PhotonAnalyzer Analyzing event number: " <<
e.id() <<
" Global Counter " <<
nEvt_ <<
"\n";
1424 bool validTriggerEvent =
true;
1428 if (!triggerEventHandle.
isValid()) {
1430 validTriggerEvent =
false;
1438 if (!photonHandle.
isValid()) {
1448 if (!loosePhotonFlag.
isValid()) {
1449 edm::LogInfo(
fName_) <<
"Error! Can't get the product: PhotonIDLoose_token_" << endl;
1458 if (!tightPhotonFlag.
isValid()) {
1459 edm::LogInfo(
fName_) <<
"Error! Can't get the product: PhotonIDTight_token_" << endl;
1472 int nPho[100][3][3];
1482 int nConv[100][3][3];
1498 if (
label.find(
"Photon") != string::npos) {
1499 for (
uint filterKeyIndex = 0; filterKeyIndex <
triggerEvent.filterKeys(filterIndex).size();
1511 if (
i != (
Keys.size() - 1)) {
1525 for (
unsigned int iPho = 0; iPho < photonHandle->size(); iPho++) {
1538 for (vector<int>::const_iterator objectKey =
Keys.begin(); objectKey !=
Keys.end();
1560 bool isLoosePhoton(
false), isTightPhoton(
false);
1562 isLoosePhoton =
true;
1565 bool phoIsInBarrel =
false;
1566 bool phoIsInEndcap =
false;
1568 if (fabs(etaPho) < 1.479)
1569 phoIsInBarrel =
true;
1571 phoIsInEndcap =
true;
1581 bool isIsolated =
false;
1583 isIsolated = isLoosePhoton;
1585 isIsolated = isTightPhoton;
1596 bool validEcalRecHits =
true;
1599 if (phoIsInBarrel) {
1602 if (!ecalRecHitHandle.
isValid()) {
1604 validEcalRecHits =
false;
1606 }
else if (phoIsInEndcap) {
1609 if (!ecalRecHitHandle.
isValid()) {
1611 validEcalRecHits =
false;
1614 if (validEcalRecHits)
1615 ecalRecHitCollection = *(ecalRecHitHandle.
product());
1620 if (isLoosePhoton) {
1624 if (isTightPhoton) {
1630 double Et = aPho->
et();
1631 bool passesCuts =
false;
1748 float(vtxH->size()),
1815 bool atLeastOneDeadChannel =
false;
1819 for (
vector<pair<DetId, float> >::const_iterator rhIt = (*bcIt)->hitsAndFractions().begin();
1820 rhIt != (*bcIt)->hitsAndFractions().end();
1824 it != ecalRecHitCollection.
end();
1826 if (rhIt->first == (*it).id()) {
1827 if ((*it).recoFlag() == 9) {
1828 atLeastOneDeadChannel =
true;
1835 if (atLeastOneDeadChannel) {
1852 for (
unsigned int iConv = 0; iConv <
conversions.size(); iConv++) {
1855 if (aConv->nTracks() < 2)
1862 if (!(aConv->conversionVertex().isValid()))
1865 float chi2Prob =
ChiSquaredProbability(aConv->conversionVertex().chi2(), aConv->conversionVertex().ndof());
1867 if (chi2Prob < 0.0005)
1876 if (
cut == 0 && isLoosePhoton) {
1880 if (
cut == 0 && isTightPhoton) {
1893 double convR =
sqrt(aConv->conversionVertex().position().perp2());
1894 double scalar = aConv->conversionVertex().position().x() * aConv->refittedPairMomentum().x() +
1895 aConv->conversionVertex().position().y() * aConv->refittedPairMomentum().y();
1900 aConv->conversionVertex().position().z(),
1906 if (fabs(aPho->
eta()) > 1.5) {
1908 }
else if (fabs(aPho->
eta()) < 1) {
1911 aConv->conversionVertex().position().x(),
1912 aConv->conversionVertex().position().y(),
1917 const std::vector<edm::RefToBase<reco::Track> >
tracks = aConv->tracks();
1919 for (
unsigned int i = 0;
i <
tracks.size();
i++) {
1928 float DPhiTracksAtVtx = -99;
1929 float dPhiTracksAtEcal = -99;
1930 float dEtaTracksAtEcal = -99;
1932 float phiTk1 = aConv->tracksPin()[0].phi();
1933 float phiTk2 = aConv->tracksPin()[1].phi();
1934 DPhiTracksAtVtx = phiTk1 - phiTk2;
1937 if (!aConv->bcMatchingWithTracks().empty() && aConv->bcMatchingWithTracks()[0].
isNonnull() &&
1938 aConv->bcMatchingWithTracks()[1].
isNonnull()) {
1939 float recoPhi1 = aConv->ecalImpactPosition()[0].phi();
1940 float recoPhi2 = aConv->ecalImpactPosition()[1].phi();
1941 float recoEta1 = aConv->ecalImpactPosition()[0].eta();
1942 float recoEta2 = aConv->ecalImpactPosition()[1].eta();
1947 dPhiTracksAtEcal = recoPhi1 - recoPhi2;
1949 dEtaTracksAtEcal = recoEta1 - recoEta2;
1965 for (
unsigned int iPho2 = iPho + 1; iPho2 < photonHandle->size(); iPho2++) {
1973 bool isTightPhoton2(
false), isLoosePhoton2(
false);
1975 isLoosePhoton2 =
true;
1982 bool isIsolated2 =
false;
1984 isIsolated2 = isLoosePhoton2;
1986 isIsolated2 = isTightPhoton2;
1995 float gamgamMass2 = p12.Dot(p12);
1998 if (aPho->
isEB() && aPho2->
isEB()) {
2000 }
else if (aPho->
isEE() && aPho2->
isEE()) {
2007 if (!conversions2.
empty() && conversions2[0]->nTracks() >= 2)
2011 }
else if (!conversions2.
empty() && conversions2[0]->nTracks() >= 2)
References edm::SortedCollection< T, SORT >::begin(), reco::Photon::chargedHadronIso(), ChiSquaredProbability(), pwdgSkimBPark_cfi::conversions, reco::Photon::conversions(), TkAlMuonSelectors_cfi::cut, reco::deltaR(), PbPb_ZMuSkimMuonDPG_cff::deltaR, deltaRMax, PFJetFilter_cfi::deltaRMin, MillePedeFileConverter_cfg::e, reco::Photon::e1x5(), reco::Photon::e2x5(), reco::Photon::e3x3(), reco::Photon::ecalRecHitSumEtConeDR04(), edm::RefVector< C, T, F >::empty(), edm::SortedCollection< T, SORT >::end(), reco::LeafCandidate::energy(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), dqmMemoryStats::float, reco::Photon::getCandidateP4type(), reco::Photon::getCorrectedEnergyError(), reco::Photon::hadronicDepth1OverEm(), reco::Photon::hadronicDepth2OverEm(), reco::Photon::hadronicOverEm(), reco::Photon::hasConversionTracks(), reco::Photon::hcalTowerSumEtConeDR04(), mps_fire::i, reco::Photon::isEB(), reco::Photon::isEE(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), label, reco::Photon::maxEnergyXtal(), reco::Photon::neutralHadronIso(), TrackCollections2monitor_cff::normalizedChi2, reco::Photon::nTrkHollowConeDR04(), reco::Photon::nTrkSolidConeDR04(), TrackCollections2monitor_cff::numberOfValidHits, reco::Photon::p4(), reco::LeafCandidate::phi(), ExoticaDQM_cfi::photonCollection, reco::Photon::photonIso(), PhotonMonitor_cff::photonSelection, edm::Handle< T >::product(), reco::Photon::r1x5(), reco::Photon::r2x5(), reco::Photon::r9(), hitfit::scalar(), reco::Photon::sigmaIetaIeta(), mathSSE::sqrt(), reco::Photon::superCluster(), PDWG_EXOHSCP_cff::tracks, PDWG_DiPhoton_SD_cff::triggerEvent, reco::Photon::trkSumPtHollowConeDR04(), reco::Photon::trkSumPtSolidConeDR04(), parallelization::uint, and trackerHitRTTI::vector.
◆ book2DHistoVector()
void PhotonAnalyzer::book2DHistoVector |
( |
DQMStore::IBooker & |
, |
|
|
std::vector< std::vector< MonitorElement * > > & |
toFill, |
|
|
std::string |
histoType, |
|
|
std::string |
histoName, |
|
|
std::string |
title, |
|
|
int |
xbin, |
|
|
double |
xmin, |
|
|
double |
xmax, |
|
|
int |
ybin = 1 , |
|
|
double |
ymin = 1 , |
|
|
double |
ymax = 2 |
|
) |
| |
|
private |
Definition at line 1253 of file PhotonAnalyzer.cc.
1264 int histo_index = 0;
1266 vector<MonitorElement*> temp1DVector;
1269 bool conversionPlot =
false;
1270 if (iBooker.pwd().find(
"Conversions") != string::npos)
1271 conversionPlot =
true;
1272 bool TwoDPlot =
false;
1273 if (
histoName.find(
"2D") != string::npos)
1276 if (conversionPlot) {
1284 stringstream histo_number_stream;
1285 histo_number_stream <<
"h_";
1286 if (histo_index < 10)
1287 histo_number_stream <<
"0";
1288 histo_number_stream << histo_index <<
"_";
1302 kind =
" Conversions: ";
1304 kind =
" Photons: ";
1306 if (histoType ==
"1D")
1307 temp1DVector.push_back(
1309 else if (histoType ==
"2D") {
1310 if ((TwoDPlot &&
type == 0) || !TwoDPlot) {
1311 temp1DVector.push_back(iBooker.book2D(
1314 }
else if (histoType ==
"Profile")
1315 temp1DVector.push_back(iBooker.bookProfile(
1318 cout <<
"bad histoType\n";
1321 temp2DVector.push_back(temp1DVector);
1322 temp1DVector.clear();
References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), gather_cfg::cout, TkAlMuonSelectors_cfi::cut, HltBtagPostValidation_cff::histoName, CalibrationSummaryClient_cfi::kind, dqm::implementation::NavigatorBase::pwd(), dqm::implementation::NavigatorBase::setCurrentFolder(), runGCPTkAlMap::title, parallelization::uint, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.
◆ book3DHistoVector()
void PhotonAnalyzer::book3DHistoVector |
( |
DQMStore::IBooker & |
, |
|
|
std::vector< std::vector< std::vector< MonitorElement * > > > & |
toFill, |
|
|
std::string |
histoType, |
|
|
std::string |
histoName, |
|
|
std::string |
title, |
|
|
int |
xbin, |
|
|
double |
xmin, |
|
|
double |
xmax, |
|
|
int |
ybin = 1 , |
|
|
double |
ymin = 1 , |
|
|
double |
ymax = 2 |
|
) |
| |
|
private |
Definition at line 1326 of file PhotonAnalyzer.cc.
1337 int histo_index = 0;
1339 vector<MonitorElement*> temp1DVector;
1340 vector<vector<MonitorElement*> > temp2DVector;
1343 bool conversionPlot =
false;
1344 if (iBooker.pwd().find(
"Conversions") != string::npos)
1345 conversionPlot =
true;
1347 if (conversionPlot) {
1355 stringstream histo_number_stream;
1356 histo_number_stream <<
"h_";
1357 if (histo_index < 10)
1358 histo_number_stream <<
"0";
1359 histo_number_stream << histo_index <<
"_";
1373 kind =
" Conversions: ";
1375 kind =
" Photons: ";
1377 if (histoType ==
"1D")
1378 temp1DVector.push_back(iBooker.book1D(histo_number_stream.str() +
histoName +
parts_[
part],
1383 else if (histoType ==
"2D")
1384 temp1DVector.push_back(iBooker.book2D(histo_number_stream.str() +
histoName +
parts_[
part],
1392 else if (histoType ==
"Profile")
1393 temp1DVector.push_back(iBooker.bookProfile(histo_number_stream.str() +
histoName +
parts_[
part],
1403 cout <<
"bad histoType\n";
1405 temp2DVector.push_back(temp1DVector);
1406 temp1DVector.clear();
1408 temp3DVector.push_back(temp2DVector);
1409 temp2DVector.clear();
References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), gather_cfg::cout, TkAlMuonSelectors_cfi::cut, HltBtagPostValidation_cff::histoName, CalibrationSummaryClient_cfi::kind, dqm::implementation::NavigatorBase::pwd(), dqm::implementation::NavigatorBase::setCurrentFolder(), runGCPTkAlMap::title, parallelization::uint, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, L1TOccupancyClient_cfi::ymax, and L1TOccupancyClient_cfi::ymin.
◆ bookHisto()
◆ bookHistograms()
◆ bookHistogramsConversions()
Definition at line 1095 of file PhotonAnalyzer.cc.
1112 "Number Of Conversions per Event ;# conversions",
1127 "#Delta#phi of Tracks at Vertex;#Delta#phi",
1135 "Abs(#Delta#phi) of Tracks at Ecal;#Delta#phi",
1143 "#Delta#eta of Tracks at Ecal;#Delta#eta",
1151 "#Deltacot(#theta) of Tracks;#Deltacot(#theta)",
1159 "Avg #Deltacot(#theta) of Tracks vs #eta;#eta;#Deltacot(#theta)",
1170 "Avg Number of Hits per Track vs #eta;#eta;# hits",
1183 "Avg #chi^{2} of Track Fitting vs #eta;#eta;#chi^{2}",
1196 "Vertex Position;Z (cm);R (cm)",
1211 "Vertex Position: #eta < 1;X (cm);Y (cm)",
1222 "#chi^{2} Probability of Vertex Fitting;#chi^{2}",
References dqm::implementation::NavigatorBase::setCurrentFolder().
◆ bookHistogramsEfficiency()
◆ bookHistogramsForHistogramCounts()
◆ bookHistogramsInvMass()
Definition at line 224 of file PhotonAnalyzer.cc.
229 "invMassAllIsolatedPhotons",
230 "Two photon invariant mass: All isolated photons;M (GeV)",
235 "invMassIsoPhotonsEBarrel",
236 "Two photon invariant mass: isolated photons in barrel; M (GeV)",
241 "invMassIsoPhotonsEEndcap",
242 "Two photon invariant mass: isolated photons in endcap; M (GeV)",
247 "invMassIsoPhotonsEEndcapEBarrel",
248 "Two photon invariant mass: isolated photons in endcap-barrel; M (GeV)",
254 iBooker,
"invMassZeroWithTracks",
"Two photon invariant mass: Neither has tracks;M (GeV)",
etBin_,
etMin_,
etMax_);
256 iBooker,
"invMassOneWithTracks",
"Two photon invariant mass: Only one has tracks;M (GeV)",
etBin_,
etMin_,
etMax_);
258 iBooker,
"invMassTwoWithTracks",
"Two photon invariant mass: Both have tracks;M (GeV)",
etBin_,
etMin_,
etMax_);
References dqm::implementation::NavigatorBase::setCurrentFolder().
◆ bookHistogramsPhotons()
Definition at line 263 of file PhotonAnalyzer.cc.
273 "phoSigmaEoverEvsNVtx",
274 "#sigma_{E}/E vs NVtx; N_{vtx}; #sigma_{E}/E",
304 "R9 vs E_{T};E_{T} (GeV);R9",
316 "Avg R9 vs E_{T};E_{T} (GeV);R9",
328 "R9 vs #eta;#eta;R9",
340 "Avg R9 vs #eta;#eta;R9",
353 "#sigma_{i#etai#eta};#sigma_{i#etai#eta}",
361 "sigmaIetaIetaVsEta2D",
362 "#sigma_{i#etai#eta} vs #eta;#eta;#sigma_{i#etai#eta}",
373 "sigmaIetaIetaVsEta",
374 "Avg #sigma_{i#etai#eta} vs #eta;#eta;#sigma_{i#etai#eta}",
388 "E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",
400 "Avg E1x5 vs E_{T};E_{T} (GeV);E1X5 (GeV)",
412 "E1x5 vs #eta;#eta;E1X5 (GeV)",
424 "Avg E1x5 vs #eta;#eta;E1X5 (GeV)",
438 "E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",
450 "Avg E2x5 vs E_{T};E_{T} (GeV);E2X5 (GeV)",
462 "E2x5 vs #eta;#eta;E2X5 (GeV)",
474 "Avg E2x5 vs #eta;#eta;E2X5 (GeV)",
488 "R1x5 vs E_{T};E_{T} (GeV);R1X5",
500 "Avg R1x5 vs E_{T};E_{T} (GeV);R1X5",
512 "R1x5 vs #eta;#eta;R1X5",
524 "Avg R1x5 vs #eta;#eta;R1X5",
538 "R2x5 vs E_{T};E_{T} (GeV);R2X5",
550 "Avg R2x5 vs E_{T};E_{T} (GeV);R2X5",
562 "R2x5 vs #eta;#eta;R2X5",
574 "Avg R2x5 vs #eta;#eta;R2X5",
587 "maxEXtalOver3x3VsEt2D",
588 "(Max Xtal E)/E3x3 vs E_{T};E_{T} (GeV);(Max Xtal E)/E3x3",
599 "maxEXtalOver3x3VsEt",
600 "Avg (Max Xtal E)/E3x3 vs E_{T};E_{T} (GeV);(Max Xtal E)/E3x3",
611 "maxEXtalOver3x3VsEta2D",
612 "(Max Xtal E)/E3x3 vs #eta;#eta;(Max Xtal E)/E3x3",
623 "maxEXtalOver3x3VsEta",
624 "Avg (Max Xtal E)/E3x3 vs #eta;#eta;(Max Xtal E)/E3x3",
638 "Number Of Tracks in the Solid Iso Cone;# tracks",
646 "nIsoTracksSolidVsEt2D",
647 "Number Of Tracks in the Solid Iso Cone vs E_{T};E_{T};# tracks",
658 "nIsoTracksSolidVsEt",
659 "Avg Number Of Tracks in the Solid Iso Cone vs E_{T};E_{T};# tracks",
670 "nIsoTracksSolidVsEta2D",
671 "Number Of Tracks in the Solid Iso Cone vs #eta;#eta;# tracks",
682 "nIsoTracksSolidVsEta",
683 "Avg Number Of Tracks in the Solid Iso Cone vs #eta;#eta;# tracks",
696 "Number Of Tracks in the Hollow Iso Cone;# tracks",
704 "nIsoTracksHollowVsEt2D",
705 "Number Of Tracks in the Hollow Iso Cone vs E_{T};E_{T};# tracks",
716 "nIsoTracksHollowVsEt",
717 "Avg Number Of Tracks in the Hollow Iso Cone vs E_{T};E_{T};# tracks",
728 "nIsoTracksHollowVsEta2D",
729 "Number Of Tracks in the Hollow Iso Cone vs #eta;#eta;# tracks",
740 "nIsoTracksHollowVsEta",
741 "Avg Number Of Tracks in the Hollow Iso Cone vs #eta;#eta;# tracks",
754 "Track P_{T} Sum in the Solid Iso Cone;P_{T} (GeV)",
762 "isoPtSumSolidVsEt2D",
763 "Track P_{T} Sum in the Solid Iso Cone;E_{T} (GeV);P_{T} (GeV)",
775 "Avg Track P_{T} Sum in the Solid Iso Cone vs E_{T};E_{T} (GeV);P_{T} (GeV)",
786 "isoPtSumSolidVsEta2D",
787 "Track P_{T} Sum in the Solid Iso Cone;#eta;P_{T} (GeV)",
798 "isoPtSumSolidVsEta",
799 "Avg Track P_{T} Sum in the Solid Iso Cone vs #eta;#eta;P_{T} (GeV)",
812 "Track P_{T} Sum in the Hollow Iso Cone;P_{T} (GeV)",
820 "isoPtSumHollowVsEt2D",
821 "Track P_{T} Sum in the Hollow Iso Cone;E_{T} (GeV);P_{T} (GeV)",
832 "isoPtSumHollowVsEt",
833 "Avg Track P_{T} Sum in the Hollow Iso Cone vs E_{T};E_{T} (GeV);P_{T} (GeV)",
844 "isoPtSumHollowVsEta2D",
845 "Track P_{T} Sum in the Hollow Iso Cone;#eta;P_{T} (GeV)",
856 "isoPtSumHollowVsEta",
857 "Avg Track P_{T} Sum in the Hollow Iso Cone vs #eta;#eta;P_{T} (GeV)",
873 "Ecal Sum in the IsoCone for Barrel;E (GeV)",
881 "Ecal Sum in the IsoCone for Endcap;E (GeV)",
890 "Ecal Sum in the Iso Cone;E_{T} (GeV);E (GeV)",
902 "Avg Ecal Sum in the Iso Cone vs E_{T};E_{T} (GeV);E (GeV)",
914 "Ecal Sum in the Iso Cone;#eta;E (GeV)",
926 "Avg Ecal Sum in the Iso Cone vs #eta;#eta;E (GeV)",
941 "Hcal Sum in the IsoCone for Barrel;E (GeV)",
949 "Hcal Sum in the IsoCone for Endcap;E (GeV)",
958 "Hcal Sum in the Iso Cone;E_{T} (GeV);E (GeV)",
970 "Avg Hcal Sum in the Iso Cone vs E_{T};E_{T} (GeV);E (GeV)",
982 "Hcal Sum in the Iso Cone;#eta;E (GeV)",
994 "Avg Hcal Sum in the Iso Cone vs #eta;#eta;E (GeV)",
1008 "Avg H/E vs Et;E_{T} (GeV);H/E",
1019 "Avg H/E vs #eta;#eta;H/E",
1038 "PF charged Had iso Barrel;E (GeV)",
1046 "PF charged Had iso Endcap;E (GeV)",
1053 "neutralHadIsoBarrel",
1054 "PF neutral Had iso Barrel;E (GeV)",
1061 "neutralHadIsoEndcap",
1062 "PF neutral Had iso Endcap;E (GeV)",
1073 "Fraction Containing Bad Channels: E_{T};E_{T} (GeV)",
1080 "phoEtaBadChannels",
1081 "Fraction Containing Bad Channels: #eta;#eta",
1088 "phoPhiBadChannels",
1089 "Fraction Containing Bad Channels: #phi;#phi",
◆ fill2DHistoVector() [1/2]
void PhotonAnalyzer::fill2DHistoVector |
( |
std::vector< std::vector< MonitorElement * > > & |
histoVector, |
|
|
double |
x, |
|
|
double |
y, |
|
|
int |
cut, |
|
|
int |
type |
|
) |
| |
|
private |
◆ fill2DHistoVector() [2/2]
void PhotonAnalyzer::fill2DHistoVector |
( |
std::vector< std::vector< MonitorElement * > > & |
histoVector, |
|
|
double |
x, |
|
|
int |
cut, |
|
|
int |
type |
|
) |
| |
|
private |
◆ fill3DHistoVector() [1/2]
void PhotonAnalyzer::fill3DHistoVector |
( |
std::vector< std::vector< std::vector< MonitorElement * > > > & |
histoVector, |
|
|
double |
x, |
|
|
double |
y, |
|
|
int |
cut, |
|
|
int |
type, |
|
|
int |
part |
|
) |
| |
|
private |
◆ fill3DHistoVector() [2/2]
void PhotonAnalyzer::fill3DHistoVector |
( |
std::vector< std::vector< std::vector< MonitorElement * > > > & |
histoVector, |
|
|
double |
x, |
|
|
int |
cut, |
|
|
int |
type, |
|
|
int |
part |
|
) |
| |
|
private |
◆ fillHistogramsForHistogramCounts()
◆ phiNormalization()
float PhotonAnalyzer::phiNormalization |
( |
float & |
a | ) |
|
|
private |
◆ photonSelection()
bool PhotonAnalyzer::photonSelection |
( |
const reco::Photon * |
p | ) |
|
|
private |
remove after moriond if (EtCorrEcalIso>4.0) result=false;
remove after moriond if (EtCorrEcalIso>50.0) result=false;
Definition at line 2075 of file PhotonAnalyzer.cc.
2081 if (pho->isEBEEGap())
2084 double EtCorrHcalIso = pho->hcalTowerSumEtConeDR03() - 0.005 * pho->pt();
2085 double EtCorrTrkIso = pho->trkSumPtHollowConeDR03() - 0.002 * pho->pt();
2087 if (pho->r9() <= 0.9) {
2088 if (pho->isEB() && (pho->hadTowOverEm() > 0.075 || pho->sigmaIetaIeta() > 0.014))
2090 if (pho->isEE() && (pho->hadTowOverEm() > 0.075 || pho->sigmaIetaIeta() > 0.034))
2093 if (EtCorrHcalIso > 4.0)
2095 if (EtCorrTrkIso > 4.0)
2097 if (pho->chargedHadronIso() > 4)
2100 if (pho->isEB() && (pho->hadTowOverEm() > 0.082 || pho->sigmaIetaIeta() > 0.014))
2102 if (pho->isEE() && (pho->hadTowOverEm() > 0.075 || pho->sigmaIetaIeta() > 0.034))
2105 if (EtCorrHcalIso > 50.0)
2107 if (EtCorrTrkIso > 50.0)
2109 if (pho->chargedHadronIso() > 4)
References reco::Photon::chargedHadronIso(), reco::LeafCandidate::eta(), reco::Photon::hadTowOverEm(), reco::Photon::hcalTowerSumEtConeDR03(), reco::Photon::isEB(), reco::Photon::isEBEEGap(), reco::Photon::isEE(), reco::LeafCandidate::pt(), reco::Photon::r9(), mps_fire::result, reco::Photon::sigmaIetaIeta(), and reco::Photon::trkSumPtHollowConeDR03().
◆ photonSelectionSlimmed()
bool PhotonAnalyzer::photonSelectionSlimmed |
( |
const reco::Photon * |
p | ) |
|
|
private |
◆ barrelRecHit_token_
◆ chi2Bin_
int PhotonAnalyzer::chi2Bin_ |
|
private |
◆ chi2Max_
double PhotonAnalyzer::chi2Max_ |
|
private |
◆ chi2Min_
double PhotonAnalyzer::chi2Min_ |
|
private |
◆ currentFolder_
std::stringstream PhotonAnalyzer::currentFolder_ |
|
private |
◆ cutStep_
double PhotonAnalyzer::cutStep_ |
|
private |
◆ dEtaTracksBin_
int PhotonAnalyzer::dEtaTracksBin_ |
|
private |
◆ dEtaTracksMax_
double PhotonAnalyzer::dEtaTracksMax_ |
|
private |
◆ dEtaTracksMin_
double PhotonAnalyzer::dEtaTracksMin_ |
|
private |
◆ dPhiTracksBin_
int PhotonAnalyzer::dPhiTracksBin_ |
|
private |
◆ dPhiTracksMax_
double PhotonAnalyzer::dPhiTracksMax_ |
|
private |
◆ dPhiTracksMin_
double PhotonAnalyzer::dPhiTracksMin_ |
|
private |
◆ eBin_
int PhotonAnalyzer::eBin_ |
|
private |
◆ eMax_
double PhotonAnalyzer::eMax_ |
|
private |
◆ eMin_
double PhotonAnalyzer::eMin_ |
|
private |
◆ endcapRecHit_token_
◆ eOverPBin_
int PhotonAnalyzer::eOverPBin_ |
|
private |
◆ eOverPMax_
double PhotonAnalyzer::eOverPMax_ |
|
private |
◆ eOverPMin_
double PhotonAnalyzer::eOverPMin_ |
|
private |
◆ etaBin_
int PhotonAnalyzer::etaBin_ |
|
private |
◆ etaMax_
double PhotonAnalyzer::etaMax_ |
|
private |
◆ etaMin_
double PhotonAnalyzer::etaMin_ |
|
private |
◆ etBin_
int PhotonAnalyzer::etBin_ |
|
private |
◆ etMax_
double PhotonAnalyzer::etMax_ |
|
private |
◆ etMin_
double PhotonAnalyzer::etMin_ |
|
private |
◆ excludeBkgHistos_
bool PhotonAnalyzer::excludeBkgHistos_ |
|
private |
◆ fName_
std::string PhotonAnalyzer::fName_ |
|
private |
◆ h_chHadIsoBarrel_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_chHadIsoBarrel_ |
|
private |
◆ h_chHadIsoEndcap_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_chHadIsoEndcap_ |
|
private |
◆ h_convEt_Loose_
◆ h_convEt_Tight_
◆ h_convEta_Loose_
◆ h_convEta_Tight_
◆ h_convVtxR_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_convVtxR_ |
|
private |
◆ h_convVtxRvsZ_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_convVtxRvsZ_ |
|
private |
◆ h_convVtxYvsX_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_convVtxYvsX_ |
|
private |
◆ h_convVtxZ_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_convVtxZ_ |
|
private |
◆ h_convVtxZEndcap_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_convVtxZEndcap_ |
|
private |
◆ h_dCotTracks_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_dCotTracks_ |
|
private |
◆ h_dEtaTracksAtEcal_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_dEtaTracksAtEcal_ |
|
private |
◆ h_dPhiTracksAtEcal_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_dPhiTracksAtEcal_ |
|
private |
◆ h_dPhiTracksAtVtx_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_dPhiTracksAtVtx_ |
|
private |
◆ h_e1x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_e1x5VsEt_ |
|
private |
◆ h_e1x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_e1x5VsEta_ |
|
private |
◆ h_e2x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_e2x5VsEt_ |
|
private |
◆ h_e2x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_e2x5VsEta_ |
|
private |
◆ h_ecalSum_
◆ h_ecalSumEBarrel_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_ecalSumEBarrel_ |
|
private |
◆ h_ecalSumEEndcap_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_ecalSumEEndcap_ |
|
private |
◆ h_ecalSumVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_ecalSumVsEt_ |
|
private |
◆ h_ecalSumVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_ecalSumVsEta_ |
|
private |
◆ h_eOverPTracks_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_eOverPTracks_ |
|
private |
◆ h_h1OverE_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_h1OverE_ |
|
private |
◆ h_h2OverE_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_h2OverE_ |
|
private |
◆ h_hcalSum_
◆ h_hcalSumEBarrel_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_hcalSumEBarrel_ |
|
private |
◆ h_hcalSumEEndcap_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_hcalSumEEndcap_ |
|
private |
◆ h_hcalSumVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_hcalSumVsEt_ |
|
private |
◆ h_hcalSumVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_hcalSumVsEta_ |
|
private |
◆ h_hOverE_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_hOverE_ |
|
private |
◆ h_invMassAllPhotons_
◆ h_invMassOneWithTracks_
◆ h_invMassPhotonsEBarrel_
◆ h_invMassPhotonsEEndcap_
◆ h_invMassPhotonsEEndcapEBarrel_
◆ h_invMassTwoWithTracks_
◆ h_invMassZeroWithTracks_
◆ h_maxEXtalOver3x3VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_maxEXtalOver3x3VsEt_ |
|
private |
◆ h_maxEXtalOver3x3VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_maxEXtalOver3x3VsEta_ |
|
private |
◆ h_nConv_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_nConv_ |
|
private |
◆ h_nHadIsoBarrel_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nHadIsoBarrel_ |
|
private |
◆ h_nHadIsoEndcap_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nHadIsoEndcap_ |
|
private |
◆ h_nPho_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_nPho_ |
|
private |
◆ h_nRecoVtx_
◆ h_nTrackIsolHollow_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nTrackIsolHollow_ |
|
private |
◆ h_nTrackIsolHollowVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nTrackIsolHollowVsEt_ |
|
private |
◆ h_nTrackIsolHollowVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nTrackIsolHollowVsEta_ |
|
private |
◆ h_nTrackIsolSolid_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nTrackIsolSolid_ |
|
private |
◆ h_nTrackIsolSolidVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nTrackIsolSolidVsEt_ |
|
private |
◆ h_nTrackIsolSolidVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_nTrackIsolSolidVsEta_ |
|
private |
◆ h_phoConvE_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoConvE_ |
|
private |
◆ h_phoConvEt_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoConvEt_ |
|
private |
◆ h_phoConvEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoConvEta_ |
|
private |
◆ h_phoConvEtaForEfficiency_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoConvEtaForEfficiency_ |
|
private |
◆ h_phoConvPhi_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoConvPhi_ |
|
private |
◆ h_phoConvPhiForEfficiency_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoConvPhiForEfficiency_ |
|
private |
◆ h_phoConvR9_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoConvR9_ |
|
private |
◆ h_phoE_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoE_ |
|
private |
◆ h_phoEt_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoEt_ |
|
private |
◆ h_phoEt_BadChannels_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoEt_BadChannels_ |
|
private |
◆ h_phoEt_Loose_
◆ h_phoEt_postHLT_
◆ h_phoEt_preHLT_
◆ h_phoEt_Tight_
◆ h_phoEta_
◆ h_phoEta_BadChannels_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoEta_BadChannels_ |
|
private |
◆ h_phoEta_Loose_
◆ h_phoEta_postHLT_
◆ h_phoEta_preHLT_
◆ h_phoEta_Tight_
◆ h_phoEta_Vertex_
◆ h_phoIsoBarrel_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoIsoBarrel_ |
|
private |
◆ h_phoIsoEndcap_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoIsoEndcap_ |
|
private |
◆ h_phoPhi_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoPhi_ |
|
private |
◆ h_phoPhi_BadChannels_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_phoPhi_BadChannels_ |
|
private |
◆ h_phoSigmaEoverE_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoSigmaEoverE_ |
|
private |
◆ h_phoSigmaIetaIeta_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_phoSigmaIetaIeta_ |
|
private |
◆ h_pOverETracks_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_pOverETracks_ |
|
private |
◆ h_r1x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_r1x5VsEt_ |
|
private |
◆ h_r1x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_r1x5VsEta_ |
|
private |
◆ h_r2x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_r2x5VsEt_ |
|
private |
◆ h_r2x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_r2x5VsEta_ |
|
private |
◆ h_r9_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_r9_ |
|
private |
◆ h_r9VsEt_
◆ h_r9VsEta_
◆ h_scEta_
◆ h_scPhi_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::h_scPhi_ |
|
private |
◆ h_sigmaIetaIetaVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_sigmaIetaIetaVsEta_ |
|
private |
◆ h_tkChi2_
◆ h_trackPtSumHollow_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_trackPtSumHollow_ |
|
private |
◆ h_trackPtSumHollowVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_trackPtSumHollowVsEt_ |
|
private |
◆ h_trackPtSumHollowVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_trackPtSumHollowVsEta_ |
|
private |
◆ h_trackPtSumSolid_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_trackPtSumSolid_ |
|
private |
◆ h_trackPtSumSolidVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_trackPtSumSolidVsEt_ |
|
private |
◆ h_trackPtSumSolidVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_trackPtSumSolidVsEta_ |
|
private |
◆ h_vertexChi2Prob_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::h_vertexChi2Prob_ |
|
private |
◆ histo_index_conversions_
int PhotonAnalyzer::histo_index_conversions_ |
|
private |
◆ histo_index_efficiency_
int PhotonAnalyzer::histo_index_efficiency_ |
|
private |
◆ histo_index_invMass_
int PhotonAnalyzer::histo_index_invMass_ |
|
private |
◆ histo_index_photons_
int PhotonAnalyzer::histo_index_photons_ |
|
private |
◆ hOverEBin_
int PhotonAnalyzer::hOverEBin_ |
|
private |
◆ hOverEMax_
double PhotonAnalyzer::hOverEMax_ |
|
private |
◆ hOverEMin_
double PhotonAnalyzer::hOverEMin_ |
|
private |
◆ invMassEtCut_
double PhotonAnalyzer::invMassEtCut_ |
|
private |
◆ isHeavyIon_
bool PhotonAnalyzer::isHeavyIon_ |
|
private |
◆ isolationStrength_
int PhotonAnalyzer::isolationStrength_ |
|
private |
◆ minimalSetOfHistos_
bool PhotonAnalyzer::minimalSetOfHistos_ |
|
private |
◆ minPhoEtCut_
double PhotonAnalyzer::minPhoEtCut_ |
|
private |
◆ nEvt_
int PhotonAnalyzer::nEvt_ |
|
private |
◆ numberBin_
int PhotonAnalyzer::numberBin_ |
|
private |
◆ numberMax_
double PhotonAnalyzer::numberMax_ |
|
private |
◆ numberMin_
double PhotonAnalyzer::numberMin_ |
|
private |
◆ numberOfSteps_
int PhotonAnalyzer::numberOfSteps_ |
|
private |
◆ offline_pvToken_
◆ p_dCotTracksVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_dCotTracksVsEta_ |
|
private |
◆ p_e1x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_e1x5VsEt_ |
|
private |
◆ p_e1x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_e1x5VsEta_ |
|
private |
◆ p_e2x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_e2x5VsEt_ |
|
private |
◆ p_e2x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_e2x5VsEta_ |
|
private |
◆ p_ecalSumVsEt_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::p_ecalSumVsEt_ |
|
private |
◆ p_ecalSumVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_ecalSumVsEta_ |
|
private |
◆ p_hcalSumVsEt_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::p_hcalSumVsEt_ |
|
private |
◆ p_hcalSumVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_hcalSumVsEta_ |
|
private |
◆ p_hOverEVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_hOverEVsEt_ |
|
private |
◆ p_hOverEVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_hOverEVsEta_ |
|
private |
◆ p_maxEXtalOver3x3VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_maxEXtalOver3x3VsEt_ |
|
private |
◆ p_maxEXtalOver3x3VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_maxEXtalOver3x3VsEta_ |
|
private |
◆ p_nHitsVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_nHitsVsEta_ |
|
private |
◆ p_nTrackIsolHollowVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_nTrackIsolHollowVsEt_ |
|
private |
◆ p_nTrackIsolHollowVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_nTrackIsolHollowVsEta_ |
|
private |
◆ p_nTrackIsolSolidVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_nTrackIsolSolidVsEt_ |
|
private |
◆ p_nTrackIsolSolidVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_nTrackIsolSolidVsEta_ |
|
private |
◆ p_phoSigmaEoverEvsNVtx_
std::vector<std::vector<std::vector<MonitorElement*> > > PhotonAnalyzer::p_phoSigmaEoverEvsNVtx_ |
|
private |
◆ p_r1x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_r1x5VsEt_ |
|
private |
◆ p_r1x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_r1x5VsEta_ |
|
private |
◆ p_r2x5VsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_r2x5VsEt_ |
|
private |
◆ p_r2x5VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_r2x5VsEta_ |
|
private |
◆ p_r9VsEt_
◆ p_r9VsEta_
◆ p_sigmaIetaIetaVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_sigmaIetaIetaVsEta_ |
|
private |
◆ p_tkChi2VsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_tkChi2VsEta_ |
|
private |
◆ p_trackPtSumHollowVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_trackPtSumHollowVsEt_ |
|
private |
◆ p_trackPtSumHollowVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_trackPtSumHollowVsEta_ |
|
private |
◆ p_trackPtSumSolidVsEt_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_trackPtSumSolidVsEt_ |
|
private |
◆ p_trackPtSumSolidVsEta_
std::vector<std::vector<MonitorElement*> > PhotonAnalyzer::p_trackPtSumSolidVsEta_ |
|
private |
◆ parts_
std::vector<std::string> PhotonAnalyzer::parts_ |
|
private |
◆ phiBin_
int PhotonAnalyzer::phiBin_ |
|
private |
◆ phiMax_
double PhotonAnalyzer::phiMax_ |
|
private |
◆ phiMin_
double PhotonAnalyzer::phiMin_ |
|
private |
◆ photon_token_
◆ PhotonIDLoose_token_
◆ PhotonIDTight_token_
◆ photonMaxEta_
double PhotonAnalyzer::photonMaxEta_ |
|
private |
◆ prescaleFactor_
unsigned int PhotonAnalyzer::prescaleFactor_ |
|
private |
◆ r9Bin_
int PhotonAnalyzer::r9Bin_ |
|
private |
◆ r9Max_
double PhotonAnalyzer::r9Max_ |
|
private |
◆ r9Min_
double PhotonAnalyzer::r9Min_ |
|
private |
◆ rBin_
int PhotonAnalyzer::rBin_ |
|
private |
◆ reducedEtaBin_
int PhotonAnalyzer::reducedEtaBin_ |
|
private |
◆ reducedEtBin_
int PhotonAnalyzer::reducedEtBin_ |
|
private |
◆ reducedR9Bin_
int PhotonAnalyzer::reducedR9Bin_ |
|
private |
◆ reducedSumBin_
int PhotonAnalyzer::reducedSumBin_ |
|
private |
◆ rMax_
double PhotonAnalyzer::rMax_ |
|
private |
◆ rMin_
double PhotonAnalyzer::rMin_ |
|
private |
◆ sigmaIetaBin_
int PhotonAnalyzer::sigmaIetaBin_ |
|
private |
◆ sigmaIetaMax_
double PhotonAnalyzer::sigmaIetaMax_ |
|
private |
◆ sigmaIetaMin_
double PhotonAnalyzer::sigmaIetaMin_ |
|
private |
◆ standAlone_
bool PhotonAnalyzer::standAlone_ |
|
private |
◆ sumBin_
int PhotonAnalyzer::sumBin_ |
|
private |
◆ sumMax_
double PhotonAnalyzer::sumMax_ |
|
private |
◆ sumMin_
double PhotonAnalyzer::sumMin_ |
|
private |
◆ totalNumberOfHistos_conversionsFolder
◆ totalNumberOfHistos_efficiencyFolder
◆ totalNumberOfHistos_invMassFolder
◆ totalNumberOfHistos_photonsFolder
◆ triggerEvent_token_
◆ types_
std::vector<std::string> PhotonAnalyzer::types_ |
|
private |
◆ useBinning_
bool PhotonAnalyzer::useBinning_ |
|
private |
◆ useTriggerFiltering_
bool PhotonAnalyzer::useTriggerFiltering_ |
|
private |
◆ xBin_
int PhotonAnalyzer::xBin_ |
|
private |
◆ xMax_
double PhotonAnalyzer::xMax_ |
|
private |
◆ xMin_
double PhotonAnalyzer::xMin_ |
|
private |
◆ yBin_
int PhotonAnalyzer::yBin_ |
|
private |
◆ yMax_
double PhotonAnalyzer::yMax_ |
|
private |
◆ yMin_
double PhotonAnalyzer::yMin_ |
|
private |
◆ zBin_
int PhotonAnalyzer::zBin_ |
|
private |
◆ zMax_
double PhotonAnalyzer::zMax_ |
|
private |
◆ zMin_
double PhotonAnalyzer::zMin_ |
|
private |
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoConvE_
edm::EDGetTokenT< std::vector< reco::Photon > > photon_token_
std::vector< std::vector< MonitorElement * > > p_tkChi2VsEta_
std::vector< std::vector< std::vector< MonitorElement * > > > h_dEtaTracksAtEcal_
std::vector< std::vector< MonitorElement * > > h_r2x5VsEta_
std::vector< std::vector< MonitorElement * > > h_tkChi2_
std::vector< std::vector< MonitorElement * > > h_hcalSumEEndcap_
std::vector< EcalRecHit >::const_iterator const_iterator
int histo_index_efficiency_
std::vector< std::vector< MonitorElement * > > h_trackPtSumHollow_
reco::SuperClusterRef superCluster() const override
Ref to SuperCluster.
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoSigmaIetaIeta_
std::vector< std::vector< MonitorElement * > > h_trackPtSumSolid_
MonitorElement * h_invMassPhotonsEEndcapEBarrel_
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoEt_
MonitorElement * h_invMassTwoWithTracks_
T const * product() const
MonitorElement * h_phoEta_Tight_
MonitorElement * h_nRecoVtx_
std::vector< std::vector< MonitorElement * > > h_phoEta_
float sigmaIetaIeta() const
std::vector< std::vector< MonitorElement * > > p_r2x5VsEta_
std::vector< std::vector< MonitorElement * > > h_trackPtSumHollowVsEta_
std::vector< std::vector< MonitorElement * > > p_trackPtSumSolidVsEta_
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoE_
std::vector< std::vector< MonitorElement * > > h_hcalSumVsEta_
std::stringstream currentFolder_
std::vector< std::vector< MonitorElement * > > p_nTrackIsolSolidVsEt_
std::vector< std::vector< MonitorElement * > > p_r9VsEta_
std::vector< std::vector< MonitorElement * > > p_hOverEVsEta_
float ChiSquaredProbability(double chiSquared, double nrDOF)
void fill3DHistoVector(std::vector< std::vector< std::vector< MonitorElement * > > > &histoVector, double x, int cut, int type, int part)
std::vector< std::vector< MonitorElement * > > h_r2x5VsEt_
std::vector< std::vector< MonitorElement * > > h_convVtxZ_
std::vector< std::vector< MonitorElement * > > h_maxEXtalOver3x3VsEt_
void book2DHistoVector(DQMStore::IBooker &, std::vector< std::vector< MonitorElement * > > &toFill, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)
MonitorElement * h_convEt_Tight_
std::vector< std::vector< MonitorElement * > > h_trackPtSumSolidVsEt_
edm::EDGetTokenT< edm::ValueMap< bool > > PhotonIDLoose_token_
std::vector< std::vector< MonitorElement * > > h_phoIsoBarrel_
const LorentzVector & p4(P4type type) const
virtual void setCurrentFolder(std::string const &fullpath)
std::vector< std::vector< MonitorElement * > > h_convVtxZEndcap_
edm::EDGetTokenT< edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit > > > endcapRecHit_token_
std::vector< std::vector< MonitorElement * > > p_e1x5VsEta_
std::vector< std::vector< MonitorElement * > > h_scEta_
std::vector< std::vector< MonitorElement * > > h_convVtxYvsX_
std::vector< std::vector< MonitorElement * > > h_r1x5VsEt_
Log< level::Info, false > LogInfo
std::vector< std::vector< MonitorElement * > > p_trackPtSumHollowVsEta_
bool useTriggerFiltering_
std::vector< std::vector< MonitorElement * > > h_nTrackIsolSolidVsEta_
std::vector< std::vector< MonitorElement * > > h_phoEt_BadChannels_
float trkSumPtHollowConeDR04() const
int histo_index_conversions_
std::vector< std::vector< MonitorElement * > > p_nTrackIsolSolidVsEta_
float neutralHadronIso() const
std::vector< std::vector< MonitorElement * > > p_hOverEVsEt_
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoConvPhiForEfficiency_
std::vector< std::vector< std::vector< MonitorElement * > > > h_r9_
std::vector< std::vector< MonitorElement * > > p_nTrackIsolHollowVsEta_
std::vector< std::vector< MonitorElement * > > h_phoEta_BadChannels_
edm::EDGetTokenT< trigger::TriggerEvent > triggerEvent_token_
std::vector< std::vector< MonitorElement * > > h_nTrackIsolSolid_
MonitorElement * h_phoEt_postHLT_
bool photonSelection(const reco::Photon *p)
bool hasConversionTracks() const
Bool flagging photons with a vector of refereces to conversions with size >0.
std::vector< std::vector< MonitorElement * > > h_hcalSumVsEt_
MonitorElement * h_phoEt_Loose_
std::vector< std::vector< MonitorElement * > > h_nTrackIsolHollowVsEta_
bool empty() const
Is the RefVector empty.
edm::EDGetTokenT< edm::ValueMap< bool > > PhotonIDTight_token_
std::vector< std::vector< std::vector< MonitorElement * > > > h_nConv_
MonitorElement * h_invMassPhotonsEEndcap_
const_iterator begin() const
std::vector< std::vector< std::vector< MonitorElement * > > > p_hcalSumVsEt_
std::vector< std::vector< std::vector< MonitorElement * > > > h_hOverE_
int nTrkSolidConeDR04() const
void fill2DHistoVector(std::vector< std::vector< MonitorElement * > > &histoVector, double x, int cut, int type)
void bookHistogramsInvMass(DQMStore::IBooker &)
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoConvR9_
std::vector< std::vector< std::vector< MonitorElement * > > > h_nPho_
int nTrkHollowConeDR04() const
double scalar(const CLHEP::HepGenMatrix &m)
Return the matrix as a scalar. Raise an assertion if the matris is not .
float chargedHadronIso() const
Accessors for Particle Flow Isolation variables.
MonitorElement * h_phoEta_preHLT_
std::vector< std::vector< MonitorElement * > > p_e2x5VsEt_
std::vector< std::vector< MonitorElement * > > h_chHadIsoEndcap_
void bookHistogramsPhotons(DQMStore::IBooker &)
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoPhi_
std::vector< std::vector< MonitorElement * > > h_e2x5VsEta_
MonitorElement * totalNumberOfHistos_invMassFolder
MonitorElement * h_phoEta_postHLT_
std::vector< std::vector< std::vector< MonitorElement * > > > p_ecalSumVsEt_
std::vector< std::vector< MonitorElement * > > h_phoConvEtaForEfficiency_
std::vector< std::vector< MonitorElement * > > h_r9VsEt_
std::vector< std::vector< std::vector< MonitorElement * > > > h_h1OverE_
std::vector< std::vector< MonitorElement * > > p_e2x5VsEta_
std::vector< std::vector< std::vector< MonitorElement * > > > p_phoSigmaEoverEvsNVtx_
std::vector< std::vector< MonitorElement * > > p_dCotTracksVsEta_
std::vector< std::vector< MonitorElement * > > p_nTrackIsolHollowVsEt_
void fillHistogramsForHistogramCounts(DQMStore::IBooker &)
std::vector< std::vector< std::vector< MonitorElement * > > > h_dPhiTracksAtEcal_
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoConvPhi_
std::vector< std::vector< MonitorElement * > > p_trackPtSumSolidVsEt_
std::vector< std::vector< std::vector< MonitorElement * > > > h_eOverPTracks_
MonitorElement * h_phoEt_Tight_
MonitorElement * h_convEta_Tight_
std::vector< std::vector< std::vector< MonitorElement * > > > h_scPhi_
void bookHistogramsEfficiency(DQMStore::IBooker &)
double eta() const final
momentum pseudorapidity
std::vector< std::vector< MonitorElement * > > p_maxEXtalOver3x3VsEt_
std::vector< std::vector< MonitorElement * > > p_sigmaIetaIetaVsEta_
std::vector< std::string > parts_
std::vector< std::vector< MonitorElement * > > p_trackPtSumHollowVsEt_
float getCorrectedEnergyError(P4type type) const
std::vector< std::vector< MonitorElement * > > h_vertexChi2Prob_
const_iterator end() const
std::vector< std::vector< MonitorElement * > > h_e1x5VsEt_
MonitorElement * totalNumberOfHistos_conversionsFolder
bool isNonnull() const
Checks for non-null.
edm::EDGetTokenT< reco::VertexCollection > offline_pvToken_
float phiNormalization(float &a)
MonitorElement * h_invMassAllPhotons_
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoConvEt_
std::vector< std::vector< std::vector< MonitorElement * > > > h_dPhiTracksAtVtx_
MonitorElement * bookInt(TString const &name, FUNC onbooking=NOOP())
MonitorElement * h_convEta_Loose_
std::vector< std::vector< MonitorElement * > > p_r1x5VsEta_
std::vector< std::vector< MonitorElement * > > h_trackPtSumSolidVsEta_
std::vector< size_type > Keys
std::vector< std::vector< std::vector< MonitorElement * > > > h_pOverETracks_
std::vector< std::vector< MonitorElement * > > h_nTrackIsolHollow_
void book3DHistoVector(DQMStore::IBooker &, std::vector< std::vector< std::vector< MonitorElement * > > > &toFill, std::string histoType, std::string histoName, std::string title, int xbin, double xmin, double xmax, int ybin=1, double ymin=1, double ymax=2)
float hadronicDepth1OverEm() const
the hadronic release in depth1 over electromagnetic fraction
std::vector< std::vector< MonitorElement * > > h_chHadIsoBarrel_
The single EDProduct to be saved for each event (AOD case)
Log< level::Error, false > LogError
std::vector< std::vector< MonitorElement * > > h_sigmaIetaIetaVsEta_
bool photonSelectionSlimmed(const reco::Photon *p)
std::vector< std::vector< MonitorElement * > > h_ecalSumEEndcap_
MonitorElement * h_phoEt_preHLT_
std::vector< std::vector< MonitorElement * > > h_nTrackIsolSolidVsEt_
MonitorElement * totalNumberOfHistos_photonsFolder
std::vector< std::vector< MonitorElement * > > h_phoIsoEndcap_
double et() const final
transverse energy
MonitorElement * totalNumberOfHistos_efficiencyFolder
std::vector< std::vector< MonitorElement * > > h_nTrackIsolHollowVsEt_
std::vector< std::vector< MonitorElement * > > h_ecalSumVsEt_
std::vector< std::vector< MonitorElement * > > h_e1x5VsEta_
float hadronicDepth2OverEm() const
the hadronic release in depth2 over electromagnetic fraction
float hadronicOverEm() const
the total hadronic over electromagnetic fraction
double phi() const final
momentum azimuthal angle
std::vector< std::vector< MonitorElement * > > h_convVtxRvsZ_
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
float trkSumPtSolidConeDR04() const
float e1x5() const
Shower shape variables.
std::vector< std::vector< MonitorElement * > > h_ecalSumEBarrel_
edm::EDGetTokenT< edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit > > > barrelRecHit_token_
reco::ConversionRefVector conversions() const
vector of references to Conversion's
unsigned int prescaleFactor_
std::vector< std::vector< MonitorElement * > > h_phoPhi_BadChannels_
std::vector< std::string > types_
std::vector< std::vector< MonitorElement * > > p_r2x5VsEt_
std::vector< std::vector< std::vector< MonitorElement * > > > h_dCotTracks_
std::vector< std::vector< MonitorElement * > > p_e1x5VsEt_
P4type getCandidateP4type() const
std::vector< std::vector< MonitorElement * > > p_r1x5VsEt_
std::vector< std::vector< MonitorElement * > > p_hcalSumVsEta_
float ecalRecHitSumEtConeDR04() const
std::vector< std::vector< MonitorElement * > > h_r1x5VsEta_
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
MonitorElement * h_convEt_Loose_
double energy() const final
energy
std::vector< std::vector< MonitorElement * > > h_maxEXtalOver3x3VsEta_
std::vector< std::vector< MonitorElement * > > p_nHitsVsEta_
std::vector< Photon > PhotonCollection
collectin of Photon objects
MonitorElement * h_invMassPhotonsEBarrel_
MonitorElement * h_phoEta_Vertex_
std::vector< std::vector< MonitorElement * > > h_e2x5VsEt_
MonitorElement * h_invMassZeroWithTracks_
std::vector< std::vector< MonitorElement * > > h_phoConvEta_
void bookHistogramsConversions(DQMStore::IBooker &)
std::vector< std::vector< MonitorElement * > > h_nHadIsoBarrel_
std::vector< std::vector< MonitorElement * > > h_ecalSumVsEta_
MonitorElement * h_invMassOneWithTracks_
std::vector< std::vector< std::vector< MonitorElement * > > > h_phoSigmaEoverE_
std::vector< std::vector< MonitorElement * > > h_ecalSum_
void bookHistogramsForHistogramCounts(DQMStore::IBooker &)
std::vector< std::vector< MonitorElement * > > h_convVtxR_
std::vector< std::vector< MonitorElement * > > h_hcalSum_
std::vector< std::vector< MonitorElement * > > h_nHadIsoEndcap_
std::vector< std::vector< MonitorElement * > > h_r9VsEta_
std::vector< std::vector< MonitorElement * > > p_maxEXtalOver3x3VsEta_
float hcalTowerSumEtConeDR04() const
Hcal isolation sum.
MonitorElement * h_phoEta_Loose_
std::vector< std::vector< std::vector< MonitorElement * > > > h_h2OverE_
std::vector< std::vector< MonitorElement * > > h_hcalSumEBarrel_
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
float maxEnergyXtal() const
std::vector< std::vector< MonitorElement * > > p_r9VsEt_
MonitorElement * bookHisto(DQMStore::IBooker &, std::string histoName, std::string title, int bin, double min, double max)
std::vector< std::vector< MonitorElement * > > h_trackPtSumHollowVsEt_
std::vector< std::vector< MonitorElement * > > p_ecalSumVsEta_