53 h_time =
new TH2F(
"h_time",
"Detector # vs. time", 20, 0., 20., 1000, 0., 750.);
55 h_timeWeighted =
new TH2F(
"h_timeWeigthed",
"Detector # vs. time weighted by the energy deposited", 20, 0., 20., 1000, 0., 750.);
57 h_timeThres =
new TH2F(
"h_timeThres",
"Detector # vs. time", 20, 0., 20., 750, 0., 750.);
59 h_timeWeightedThres =
new TH2F(
"h_timeWeigthedThres",
"Detector # vs. time weighted by the energy deposited", 20, 0., 20., 750, 0.,
62 h_edep =
new TH2F(
"h_edep",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
64 h_edep1 =
new TH2F(
"h_edep1",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
66 h_edep2 =
new TH2F(
"h_edep2",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
68 h_edep3 =
new TH2F(
"h_edep3",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
71 h_edepThres =
new TH2F(
"h_edepThres",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
73 h_edepThres1 =
new TH2F(
"h_edepThres1",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
75 h_edepThres2 =
new TH2F(
"h_edepThres2",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
77 h_edepThres3 =
new TH2F(
"h_edepThres3",
"Time bin # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
79 for (
int i = 0; i < 3; i++) {
80 h_zvedep[i] =
new TH1F(TString::Format(
"h_zvedep_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25.);
83 h_xvzvedep[i] =
new TH2F(TString::Format(
"h_xvzvedep_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
86 h_yvzvedep[i] =
new TH2F(TString::Format(
"h_yvzvedep_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
89 h_rvzvedep[i] =
new TH2F(TString::Format(
"h_rvzvedep_%d", i),
"edep [MeV] vs. z [cm]", 2000, 0., 25., 2000, -25., 25.);
92 h_xvyvedep[i] =
new TH2F(TString::Format(
"h_xvyvedep_%d", i),
"edep [MeV] vs. y [cm]", 2000, -25., 25., 2000, -25., 25.);
95 h_zvedepW[i] =
new TH1F(TString::Format(
"h_zvedepW_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25.);
98 h_xvzvedepW[i] =
new TH2F(TString::Format(
"h_xvzvedepW_%d", i),
"edep [MeV] vs. x vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
101 h_yvzvedepW[i] =
new TH2F(TString::Format(
"h_yvzvedepW_%d", i),
"edep [MeV] vs. y vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
104 h_xvyvedepW[i] =
new TH2F(TString::Format(
"h_xvyvedepW_%d", i),
"edep [MeV] vs. x vs. y [cm]", 2000, -25., 25., 2000, -25., 25.);
107 h_rvzvedepW[i] =
new TH2F(TString::Format(
"h_rvzvedepW_%d", i),
"edep [MeV] vs. z [cm]", 2000, 0., 25., 2000, -25., 25.);
111 h_zvedepT[i] =
new TH1F(TString::Format(
"h_zvedepT_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25.);
114 h_xvzvedepT[i] =
new TH2F(TString::Format(
"h_xvzvedepT_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
117 h_yvzvedepT[i] =
new TH2F(TString::Format(
"h_yvzvedepT_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
120 h_rvzvedepT[i] =
new TH2F(TString::Format(
"h_rvzvedepT_%d", i),
"edep [MeV] vs. z [cm]", 2000, 0., 25., 2000, -25., 25.);
123 h_xvyvedepT[i] =
new TH2F(TString::Format(
"h_xvyvedepT_%d", i),
"edep [MeV] vs. y [cm]", 2000, -25., 25., 2000, -25., 25.);
126 h_zvedepWT[i] =
new TH1F(TString::Format(
"h_zvedepWT_%d", i),
"edep [MeV] vs. z [cm]", 2000, -25., 25.);
129 h_xvzvedepWT[i] =
new TH2F(TString::Format(
"h_xvzvedepWT_%d", i),
"edep [MeV] vs. x vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
132 h_yvzvedepWT[i] =
new TH2F(TString::Format(
"h_yvzvedepWT_%d", i),
"edep [MeV] vs. y vs. z [cm]", 2000, -25., 25., 2000, -25., 25.);
135 h_xvyvedepWT[i] =
new TH2F(TString::Format(
"h_xvyvedepWT_%d", i),
"edep [MeV] vs. x vs. y [cm]", 2000, -25., 25., 2000, -25., 25.);
138 h_rvzvedepWT[i] =
new TH2F(TString::Format(
"h_rvzvedepWT_%d", i),
"edep [MeV] vs. z [cm]", 2000, 0., 25., 2000, -25., 25.);
141 h_Edep =
new TH2F(
"h_Edep",
"det # # vs. energy deposited", 20, 0., 20., 1000, 0., 10.);
142 h_pxNb =
new TH2F(
"h_pxNb",
"det # # vs. nb pixel", 20, 0., 20., 1000, 0., 1000.);
143 for (
int i = 0; i < 15; i++) {
144 h_cvr[i] =
new TH2F(TString::Format(
"cvr_%d", i),
" col v. row", 80, 0., 80., 336, 0., 336.);
191 if (olddetNb != detNb) {
193 h_Edep->Fill(detNb, esum);
194 h_pxNb->Fill(detNb, ipix);
198 for (
int j = 0; j < maxSIZE; j++) {
212 h_cvr[detNb]->Fill(col, row);
220 for (
int i = 0; i < nSimHits; i++) {
225 double timeBin = aHit->
getTime();
228 ROOT::Math::XYZVector position = aHit->
getPosEntry();
229 double r =
sqrt(position.X() * position.X() + position.Y() * position.Y());
230 int detNB = (lad - 1) * 5 + sen - 1;
233 h_time->Fill(detNB, timeBin);
234 h_edep->Fill(detNB, adep * 1e3);
238 if (adep > 50.*1
e-6) {
245 h_zvedep[lad - 1]->Fill(position.Z());
246 h_xvzvedep[lad - 1]->Fill(position.X(), position.Z());
247 h_yvzvedep[lad - 1]->Fill(position.Y(), position.Z());
248 h_xvyvedep[lad - 1]->Fill(position.X(), position.Y());
250 h_zvedepW[lad - 1]->Fill(position.Z(), adep * 1e3);
251 h_xvzvedepW[lad - 1]->Fill(position.X(), position.Z(), adep * 1e3);
252 h_yvzvedepW[lad - 1]->Fill(position.Y(), position.Z(), adep * 1e3);
253 h_xvyvedepW[lad - 1]->Fill(position.X(), position.Y(), adep * 1e3);
254 h_rvzvedepW[lad - 1]->Fill(r, position.Z(), adep * 1e3);
255 if (adep > 50.*1
e-6) {
257 h_xvzvedepT[lad - 1]->Fill(position.X(), position.Z());
258 h_yvzvedepT[lad - 1]->Fill(position.Y(), position.Z());
259 h_xvyvedepT[lad - 1]->Fill(position.X(), position.Y());
261 h_zvedepWT[lad - 1]->Fill(position.Z(), adep * 1e3);
262 h_xvzvedepWT[lad - 1]->Fill(position.X(), position.Z(), adep * 1e3);
263 h_yvzvedepWT[lad - 1]->Fill(position.Y(), position.Z(), adep * 1e3);
264 h_xvyvedepWT[lad - 1]->Fill(position.X(), position.Y(), adep * 1e3);
265 h_rvzvedepWT[lad - 1]->Fill(r, position.Z(), adep * 1e3);