110 hadleakvset =
Book2D(
"hadleakvset",
"Hadronic leakage vs E_{T} of " +
m_sParticleType+
"; E_{hadleak} (GeV) ; E_{T} (GeV) ", 100, -0.07, 0.13, 200, 0., 200);
159 float eta2 = fabs (
egamma.caloCluster()->etaBE (2));
160 float et37 =
egamma.caloCluster()->
e() / cosh (eta2);
166 raphad = et37 > 0. ? ethad / et37 : 0.;
169 raphad1 = et37 > 0. ? ethad1 / et37 : 0.;
171 float hadrleak = (eta2 >= 0.8 && eta2 < 1.37) ? raphad : raphad1;
172 hadleak->Fill(hadrleak, weight);
187 float rateta37 = e277!= 0. ? e237/e277:0.;
188 float rateta33 = e277!= 0. ? e233/e277:0.;
189 reta37->Fill(rateta37, weight);
190 reta33->Fill(rateta33, weight);
201 weta1->Fill(shweta1, weight);
202 weta2->Fill(shweta2, weight);
216 float deltemax1 = fabs (emax + emax2) > 0. ? (emax - emax2) / (emax + emax2) : 0.;
217 demax1->Fill(deltemax1, weight);
218 demax2->Fill(emax2 / (1. + 0.009 * et37/
GeV)/
GeV, weight);
219 de->Fill( (emax2 - emin)/
GeV, weight);
241 fracs1->Fill(shfracs1, weight);
247 wtots1->Fill(shwtots1, weight);
253 f1->Fill(fracf1, weight);
259 f3->Fill(fracf3, weight);
264 Eratio->Fill(eRatio, weight);
270 pos7->Fill(difpos7, weight);
276 Rhad->Fill(rhad, weight);
282 Rhad1->Fill(rhad1, weight);
288 Reta->Fill(reta, weight);
294 Rphi->Fill(rphi, weight);
299 float topoEtconeIso(0);
303 topoEtconeIso = et37!=0. ? shetcone20/et37 : 0.;
304 clusiso->Fill(topoEtconeIso, weight);
310 const double e0 =
egamma.caloCluster()->energyBE(0);
311 const double e1 =
egamma.caloCluster()->energyBE(1);
312 const double e2 =
egamma.caloCluster()->energyBE(2);
313 const double e3 =
egamma.caloCluster()->energyBE(3);
315 E0->Fill(e0 /
GeV, weight);
316 E1->Fill(e1 /
GeV, weight);
317 E2->Fill(e2 /
GeV, weight);
318 E3->Fill(e3 /
GeV, weight);