14 TString
N = name_.c_str();
16 h4DSegmXPosInCham =
new TH1F(N +
"_h4DSegmXPosInCham",
"4D Segment x position (cm) in Chamber RF", 200, -200, 200);
17 h4DSegmYPosInCham =
new TH1F(N +
"_h4DSegmYPosInCham",
"4D Segment y position (cm) in Chamber RF", 200, -200, 200);
19 new TH1F(N +
"_h4DSegmPhiAngleInCham",
"4D Segment phi angle (rad) in Chamber RF", 180, -180, 180);
21 new TH1F(N +
"_h4DSegmThetaAngleInCham",
"4D Segment theta angle (rad) in Chamber RF", 180, -180, 180);
23 new TH1F(N +
"_h4DSegmImpactAngleInCham",
"4D Segment impact angle (rad) in Chamber RF", 180, -180, 180);
73 TString
N = name_.c_str();
75 h2DSegmPosInCham =
new TH1F(N +
"_h2DSegmPosInCham",
"2D Segment position (cm) in Chamber RF", 200, -200, 200);
76 h2DSegmAngleInCham =
new TH1F(N +
"_h2DSegmAngleInCham",
"2D Segment angle (rad) in Chamber RF", 200, -2, 2);
77 h2DSegmCosAngleInCham =
new TH1F(N +
"_h2DSegmCosAngleInCham",
"2D Segment cos(angle) in Chamber RF", 200, -2, 2);
82 h2DSegmPosInCham = (TH1F *)file->Get(
name +
"_h2DSegmPosInCham");
83 h2DSegmAngleInCham = (TH1F *)file->Get(
name +
"_h2DSegmAngleInCham");
84 h2DSegmCosAngleInCham = (TH1F *)file->Get(
name +
"_h2DSegmCosAngleInCham");
87 delete h2DSegmPosInCham;
88 delete h2DSegmAngleInCham;
89 delete h2DSegmCosAngleInCham;
92 h2DSegmPosInCham->Fill(pos);
93 h2DSegmAngleInCham->Fill(atan(localAngle));
94 h2DSegmCosAngleInCham->Fill(
cos(atan(localAngle)));
97 h2DSegmPosInCham->Write();
98 h2DSegmAngleInCham->Write();
99 h2DSegmCosAngleInCham->Write();
117 hTmax123 =
new TH1F(
name +
"_Tmax123",
"Tmax123 value", 2000, -1000., 1000.);
118 hTmax124s72 =
new TH1F(
name +
"_Tmax124_s72",
"Tmax124 sigma=sqrt(7/2) value", 2000, -1000., 1000.);
119 hTmax124s78 =
new TH1F(
name +
"_Tmax124_s78",
"Tmax124 sigma=sqrt(7/8) value", 2000, -1000., 1000.);
120 hTmax134s72 =
new TH1F(
name +
"_Tmax134_s72",
"Tmax134 sigma=sqrt(7/2) value", 2000, -1000., 1000.);
121 hTmax134s78 =
new TH1F(
name +
"_Tmax134_s78",
"Tmax134 sigma=sqrt(7/8) value", 2000, -1000., 1000.);
122 hTmax234 =
new TH1F(
name +
"_Tmax234",
"Tmax234 value", 2000, -1000., 1000.);
123 hTmax_3t0 =
new TH1F(
name +
"_3t0",
"Tmax+3*Delta(t0)", 2000, -1000., 1000.);
124 hTmax_3t0_0 =
new TH1F(
name +
"_3t0_0",
"Tmax+3*Delta(t0); 3 hits", 2000, -1000., 1000.);
125 hTmax_3t0_1 =
new TH1F(
name +
"_3t0_1",
"Tmax+3*Delta(t0); one t<5ns", 2000, -1000., 1000.);
126 hTmax_3t0_2 =
new TH1F(
name +
"_3t0_2",
"Tmax+3*Delta(t0); one t<10ns", 2000, -1000., 1000.);
127 hTmax_3t0_3 =
new TH1F(
name +
"_3t0_3",
"Tmax+3*Delta(t0); one t<20ns", 2000, -1000., 1000.);
128 hTmax_3t0_4 =
new TH1F(
name +
"_3t0_4",
"Tmax+3*Delta(t0); one t<50ns", 2000, -1000., 1000.);
129 hTmax_3t0_5 =
new TH1F(
name +
"_3t0_5",
"Tmax+3*Delta(t0); all t>50ns", 2000, -1000., 1000.);
130 hTmax_2t0 =
new TH1F(
name +
"_2t0",
"Tmax+2*Delta(t0)", 2000, -1000., 1000.);
131 hTmax_2t0_0 =
new TH1F(
name +
"_2t0_0",
"Tmax+2*Delta(t0); 3 hits", 2000, -1000., 1000.);
132 hTmax_2t0_1 =
new TH1F(
name +
"_2t0_1",
"Tmax+2*Delta(t0); one t<5ns", 2000, -1000., 1000.);
133 hTmax_2t0_2 =
new TH1F(
name +
"_2t0_2",
"Tmax+2*Delta(t0); one t<10ns", 2000, -1000., 1000.);
134 hTmax_2t0_3 =
new TH1F(
name +
"_2t0_3",
"Tmax+2*Delta(t0); one t<20ns", 2000, -1000., 1000.);
135 hTmax_2t0_4 =
new TH1F(
name +
"_2t0_4",
"Tmax+2*Delta(t0); one t<50ns", 2000, -1000., 1000.);
136 hTmax_2t0_5 =
new TH1F(
name +
"_2t0_5",
"Tmax+2*Delta(t0); all t>50ns", 2000, -1000., 1000.);
137 hTmax_t0 =
new TH1F(
name +
"_t0",
"Tmax+Delta(t0)", 2000, -1000., 1000.);
138 hTmax_t0_0 =
new TH1F(
name +
"_t0_0",
"Tmax+Delta(t0); 3 hits", 2000, -1000., 1000.);
139 hTmax_t0_1 =
new TH1F(
name +
"_t0_1",
"Tmax+Delta(t0); one t<5ns", 2000, -1000., 1000.);
140 hTmax_t0_2 =
new TH1F(
name +
"_t0_2",
"Tmax+Delta(t0); one t<10ns", 2000, -1000., 1000.);
141 hTmax_t0_3 =
new TH1F(
name +
"_t0_3",
"Tmax+Delta(t0); one t<20ns", 2000, -1000., 1000.);
142 hTmax_t0_4 =
new TH1F(
name +
"_t0_4",
"Tmax+Delta(t0); one t<50ns", 2000, -1000., 1000.);
143 hTmax_t0_5 =
new TH1F(
name +
"_t0_5",
"Tmax+Delta(t0); all t>50ns", 2000, -1000., 1000.);
144 hTmax_0 =
new TH1F(
name +
"_0",
"Tmax", 2000, -1000., 1000.);
149 hTmax123 = (TH1F *)file->Get(
name +
"_Tmax123");
150 hTmax124s72 = (TH1F *)file->Get(
name +
"_Tmax124_s72");
151 hTmax124s78 = (TH1F *)file->Get(
name +
"_Tmax124_s78");
152 hTmax134s72 = (TH1F *)file->Get(
name +
"_Tmax134_s72");
153 hTmax134s78 = (TH1F *)file->Get(
name +
"_Tmax134_s78");
154 hTmax234 = (TH1F *)file->Get(
name +
"_Tmax234");
155 hTmax_3t0 = (TH1F *)file->Get(
name +
"_3t0");
156 hTmax_3t0_0 = (TH1F *)file->Get(
name +
"_3t0_0");
157 hTmax_3t0_1 = (TH1F *)file->Get(
name +
"_3t0_1");
158 hTmax_3t0_2 = (TH1F *)file->Get(
name +
"_3t0_2");
159 hTmax_3t0_3 = (TH1F *)file->Get(
name +
"_3t0_3");
160 hTmax_3t0_4 = (TH1F *)file->Get(
name +
"_3t0_4");
161 hTmax_3t0_5 = (TH1F *)file->Get(
name +
"_3t0_5");
162 hTmax_2t0 = (TH1F *)file->Get(
name +
"_2t0");
163 hTmax_2t0_0 = (TH1F *)file->Get(
name +
"_2t0_0");
164 hTmax_2t0_1 = (TH1F *)file->Get(
name +
"_2t0_1");
165 hTmax_2t0_2 = (TH1F *)file->Get(
name +
"_2t0_2");
166 hTmax_2t0_3 = (TH1F *)file->Get(
name +
"_2t0_3");
167 hTmax_2t0_4 = (TH1F *)file->Get(
name +
"_2t0_4");
168 hTmax_2t0_5 = (TH1F *)file->Get(
name +
"_2t0_5");
169 hTmax_t0 = (TH1F *)file->Get(
name +
"_t0");
170 hTmax_t0_1 = (TH1F *)file->Get(
name +
"_t0_1");
171 hTmax_t0_2 = (TH1F *)file->Get(
name +
"_t0_2");
172 hTmax_t0_3 = (TH1F *)file->Get(
name +
"_t0_3");
173 hTmax_t0_4 = (TH1F *)file->Get(
name +
"_t0_4");
174 hTmax_t0_5 = (TH1F *)file->Get(
name +
"_t0_5");
175 hTmax_0 = (TH1F *)file->Get(
name +
"_0");
222 hTmax123->Fill(tmax123);
224 hTmax_t0->Fill(tmax123);
227 hTmax_t0_0->Fill(tmax123);
230 hTmax_t0_1->Fill(tmax123);
233 hTmax_t0_2->Fill(tmax123);
236 hTmax_t0_3->Fill(tmax123);
239 hTmax_t0_4->Fill(tmax123);
242 hTmax_t0_5->Fill(tmax123);
246 hTmax_2t0->Fill(tmax123);
249 hTmax_2t0_0->Fill(tmax123);
252 hTmax_2t0_1->Fill(tmax123);
255 hTmax_2t0_2->Fill(tmax123);
258 hTmax_2t0_3->Fill(tmax123);
261 hTmax_2t0_4->Fill(tmax123);
264 hTmax_2t0_5->Fill(tmax123);
270 (s124 ==
dttmaxenums::r72) ? hTmax124s72->Fill(tmax124) : hTmax124s78->Fill(tmax124);
272 hTmax_0->Fill(tmax124);
273 else if (t0_124 == 1) {
274 hTmax_t0->Fill(tmax124);
277 hTmax_t0_0->Fill(tmax124);
280 hTmax_t0_1->Fill(tmax124);
283 hTmax_t0_2->Fill(tmax124);
286 hTmax_t0_3->Fill(tmax124);
289 hTmax_t0_4->Fill(tmax124);
292 hTmax_t0_5->Fill(tmax124);
295 }
else if (t0_124 == 2) {
296 hTmax_2t0->Fill(tmax124);
299 hTmax_2t0_0->Fill(tmax124);
302 hTmax_2t0_1->Fill(tmax124);
305 hTmax_2t0_2->Fill(tmax124);
308 hTmax_2t0_3->Fill(tmax124);
311 hTmax_2t0_4->Fill(tmax124);
314 hTmax_2t0_5->Fill(tmax124);
317 }
else if (t0_124 == 3) {
318 hTmax_3t0->Fill(tmax124);
321 hTmax_3t0_0->Fill(tmax124);
324 hTmax_3t0_1->Fill(tmax124);
327 hTmax_3t0_2->Fill(tmax124);
330 hTmax_3t0_3->Fill(tmax124);
333 hTmax_3t0_4->Fill(tmax124);
336 hTmax_3t0_5->Fill(tmax124);
342 (s134 ==
dttmaxenums::r72) ? hTmax134s72->Fill(tmax134) : hTmax134s78->Fill(tmax134);
344 hTmax_0->Fill(tmax134);
345 else if (t0_134 == 1) {
346 hTmax_t0->Fill(tmax134);
349 hTmax_t0_0->Fill(tmax134);
352 hTmax_t0_1->Fill(tmax134);
355 hTmax_t0_2->Fill(tmax134);
358 hTmax_t0_3->Fill(tmax134);
361 hTmax_t0_4->Fill(tmax134);
364 hTmax_t0_5->Fill(tmax134);
367 }
else if (t0_134 == 2) {
368 hTmax_2t0->Fill(tmax134);
371 hTmax_2t0_0->Fill(tmax134);
374 hTmax_2t0_1->Fill(tmax134);
377 hTmax_2t0_2->Fill(tmax134);
380 hTmax_2t0_3->Fill(tmax134);
383 hTmax_2t0_4->Fill(tmax134);
386 hTmax_2t0_5->Fill(tmax134);
389 }
else if (t0_134 == 3) {
390 hTmax_3t0->Fill(tmax134);
393 hTmax_3t0_0->Fill(tmax134);
396 hTmax_3t0_1->Fill(tmax134);
399 hTmax_3t0_2->Fill(tmax134);
402 hTmax_3t0_3->Fill(tmax134);
405 hTmax_3t0_4->Fill(tmax134);
408 hTmax_3t0_5->Fill(tmax134);
414 hTmax234->Fill(tmax234);
416 hTmax_t0->Fill(tmax234);
419 hTmax_t0_0->Fill(tmax234);
422 hTmax_t0_1->Fill(tmax234);
425 hTmax_t0_2->Fill(tmax234);
428 hTmax_t0_3->Fill(tmax234);
431 hTmax_t0_4->Fill(tmax234);
434 hTmax_t0_5->Fill(tmax234);
438 hTmax_2t0->Fill(tmax234);
441 hTmax_2t0_0->Fill(tmax234);
444 hTmax_2t0_1->Fill(tmax234);
447 hTmax_2t0_2->Fill(tmax234);
450 hTmax_2t0_3->Fill(tmax234);
453 hTmax_2t0_4->Fill(tmax234);
456 hTmax_2t0_5->Fill(tmax234);
466 hTmax124s72->Write();
467 hTmax124s78->Write();
468 hTmax134s72->Write();
469 hTmax134s78->Write();
472 hTmax_3t0_0->Write();
473 hTmax_3t0_1->Write();
474 hTmax_3t0_2->Write();
475 hTmax_3t0_3->Write();
476 hTmax_3t0_4->Write();
477 hTmax_3t0_5->Write();
479 hTmax_2t0_0->Write();
480 hTmax_2t0_1->Write();
481 hTmax_2t0_2->Write();
482 hTmax_2t0_3->Write();
483 hTmax_2t0_4->Write();
484 hTmax_2t0_5->Write();
498 if (hist == hTmax_3t0)
500 else if (hist == hTmax_2t0)
502 else if (hist == hTmax_t0)
504 else if (hist == hTmax_0)
TH1F * h4DSegmThetaAngleInCham
h2DSegm(const TString &name_, TFile *file)
hTMaxCell(const TString &name_, TFile *file)
int GetT0Factor(TH1F *hist)
void Fill(float x, float phi)
Geom::Theta< T > theta() const
void Fill(float x, float y, float phi, float theta, float impact)
void Fill(float pos, float localAngle)
TH1F * h2DSegmAngleInCham
h4DSegm(const TString &name_, TFile *file)
TH1F * h4DSegmImpactAngleInCham
Cos< T >::type cos(const T &t)
void Fill(float tmax123, float tmax124, float tmax134, float tmax234, dttmaxenums::SigmaFactor s124, dttmaxenums::SigmaFactor s134, unsigned t0_123, unsigned t0_124, unsigned t0_134, unsigned t0_234, unsigned hSubGroup)
hTMaxCell(const TString &name_)
TH1F * h2DSegmCosAngleInCham
h2DSegm(std::string name_)
TH1F * h4DSegmPhiAngleInCham
h4DSegm(std::string name_)