25 constexpr
auto size = 100;
29 -5.191, -4.889, -4.716, -4.538, -4.363, -4.191, -4.013, -3.839, -3.664, -3.489, -3.314, -3.139, -2.964, -2.853,
30 -2.65, -2.5, -2.322, -2.172, -2.043, -1.93, -1.83, -1.74, -1.653, -1.566, -1.479, -1.392, -1.305, -1.218,
31 -1.131, -1.044, -0.957, -0.879, -0.783, -0.696, -0.609, -0.522, -0.435, -0.348, -0.261, -0.174, -0.087, 0,
32 0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609, 0.696, 0.783, 0.879, 0.957, 1.044, 1.131, 1.218,
33 1.305, 1.392, 1.479, 1.566, 1.653, 1.74, 1.83, 1.93, 2.043, 2.172, 2.322, 2.5, 2.65, 2.853,
34 2.964, 3.139, 3.314, 3.489, 3.664, 3.839, 4.013, 4.191, 4.363, 4.538, 4.716, 4.889, 5.191};
55 strncpy(
histo,
"emean_vs_eta_EHFO",
size);
60 strncpy(
histo,
"Ratio_Esummed_ECAL_0",
size);
62 strncpy(
histo,
"Ratio_Esummed_HCAL_0",
size);
64 strncpy(
histo,
"Ratio_Esummed_HO_0",
size);
66 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_0",
size);
68 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_HO_0",
size);
71 strncpy(
histo,
"Ratio_Esummed_ECAL_1",
size);
73 strncpy(
histo,
"Ratio_Esummed_HCAL_1",
size);
75 strncpy(
histo,
"Ratio_Esummed_HO_1",
size);
77 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_1",
size);
79 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_HO_1",
size);
82 strncpy(
histo,
"Ratio_Esummed_ECAL_2",
size);
84 strncpy(
histo,
"Ratio_Esummed_HCAL_2",
size);
86 strncpy(
histo,
"Ratio_Esummed_HO_2",
size);
88 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_2",
size);
90 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_HO_2",
size);
93 strncpy(
histo,
"Ratio_Esummed_ECAL_3",
size);
95 strncpy(
histo,
"Ratio_Esummed_HCAL_3",
size);
97 strncpy(
histo,
"Ratio_Esummed_HO_3",
size);
99 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_3",
size);
101 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_HO_3",
size);
104 strncpy(
histo,
"Ratio_Esummed_ECAL_4",
size);
106 strncpy(
histo,
"Ratio_Esummed_HCAL_4",
size);
108 strncpy(
histo,
"Ratio_Esummed_HO_4",
size);
110 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_4",
size);
112 strncpy(
histo,
"Ratio_Esummed_ECAL_HCAL_HO_4",
size);
115 strncpy(
histo,
"Ratio_Esummed_HF_5",
size);
117 strncpy(
histo,
"Ratio_Esummed_HF_6",
size);
127 double energy_MC = 0.;
131 if (not
hepMC.isValid()) {
137 double maxPt = -99999.;
139 for (
auto p = myGenEvent->particles_begin();
p != myGenEvent->particles_end(); ++
p) {
140 double phip = (*p)->momentum().phi();
141 double etap = (*p)->momentum().eta();
142 double pt = (*p)->momentum().perp();
143 double energy = (*p)->momentum().e();
167 const double Econe =
sumEnergy(pfClusterECAL, eta_MC, phi_MC);
168 const double Hcone =
sumEnergy(pfClusterHCAL, eta_MC, phi_MC);
169 const double HOcone =
sumEnergy(pfClusterHO, eta_MC, phi_MC);
170 const double HFcone =
sumEnergy(pfClusterHF, eta_MC, phi_MC);
172 if (energy_MC > 0.) {
226 for (
auto pf = pfClusters->begin();
pf != pfClusters->end(); ++
pf) {