14 TString
N = name_.c_str();
17 "4D Segment x position (cm) in Chamber RF", 200, -200, 200);
19 "4D Segment y position (cm) in Chamber RF", 200, -200, 200);
21 "4D Segment phi angle (rad) in Chamber RF", 180, -180, 180);
23 "4D Segment theta angle (rad) in Chamber RF", 180, -180, 180);
25 "4D Segment impact angle (rad) in Chamber RF", 180, -180, 180);
76 TString
N = name_.c_str();
78 h2DSegmPosInCham =
new TH1F(N+
"_h2DSegmPosInCham",
79 "2D Segment position (cm) in Chamber RF", 200, -200, 200);
80 h2DSegmAngleInCham =
new TH1F(N+
"_h2DSegmAngleInCham",
81 "2D Segment angle (rad) in Chamber RF", 200, -2, 2);
82 h2DSegmCosAngleInCham =
new TH1F(N+
"_h2DSegmCosAngleInCham",
83 "2D Segment cos(angle) in Chamber RF", 200, -2, 2);
88 h2DSegmPosInCham = (TH1F *) file->Get(
name+
"_h2DSegmPosInCham");
89 h2DSegmAngleInCham = (TH1F *) file->Get(
name+
"_h2DSegmAngleInCham");
90 h2DSegmCosAngleInCham = (TH1F *) file->Get(
name+
"_h2DSegmCosAngleInCham");
93 delete h2DSegmPosInCham;
94 delete h2DSegmAngleInCham;
95 delete h2DSegmCosAngleInCham;
99 h2DSegmPosInCham->Fill(pos);
100 h2DSegmAngleInCham->Fill(atan(localAngle));
101 h2DSegmCosAngleInCham->Fill(
cos(atan(localAngle)));
104 h2DSegmPosInCham->Write();
105 h2DSegmAngleInCham->Write();
106 h2DSegmCosAngleInCham->Write();
124 hTmax123 =
new TH1F (
name+
"_Tmax123",
"Tmax123 value", 2000, -1000., 1000.);
125 hTmax124s72 =
new TH1F (
name+
"_Tmax124_s72",
"Tmax124 sigma=sqrt(7/2) value", 2000, -1000., 1000.);
126 hTmax124s78 =
new TH1F (
name+
"_Tmax124_s78",
"Tmax124 sigma=sqrt(7/8) value", 2000, -1000., 1000.);
127 hTmax134s72 =
new TH1F (
name+
"_Tmax134_s72",
"Tmax134 sigma=sqrt(7/2) value", 2000, -1000., 1000.);
128 hTmax134s78 =
new TH1F (
name+
"_Tmax134_s78",
"Tmax134 sigma=sqrt(7/8) value", 2000, -1000., 1000.);
129 hTmax234 =
new TH1F (
name+
"_Tmax234",
"Tmax234 value", 2000, -1000., 1000.);
130 hTmax_3t0 =
new TH1F (
name+
"_3t0",
"Tmax+3*Delta(t0)", 2000, -1000., 1000.);
131 hTmax_3t0_0 =
new TH1F (
name+
"_3t0_0",
"Tmax+3*Delta(t0); 3 hits", 2000, -1000., 1000.);
132 hTmax_3t0_1 =
new TH1F (
name+
"_3t0_1",
"Tmax+3*Delta(t0); one t<5ns", 2000, -1000., 1000.);
133 hTmax_3t0_2 =
new TH1F (
name+
"_3t0_2",
"Tmax+3*Delta(t0); one t<10ns", 2000, -1000., 1000.);
134 hTmax_3t0_3 =
new TH1F (
name+
"_3t0_3",
"Tmax+3*Delta(t0); one t<20ns", 2000, -1000., 1000.);
135 hTmax_3t0_4 =
new TH1F (
name+
"_3t0_4",
"Tmax+3*Delta(t0); one t<50ns", 2000, -1000., 1000.);
136 hTmax_3t0_5 =
new TH1F (
name+
"_3t0_5",
"Tmax+3*Delta(t0); all t>50ns", 2000, -1000., 1000.);
137 hTmax_2t0 =
new TH1F (
name+
"_2t0",
"Tmax+2*Delta(t0)", 2000, -1000., 1000.);
138 hTmax_2t0_0 =
new TH1F (
name+
"_2t0_0",
"Tmax+2*Delta(t0); 3 hits", 2000, -1000., 1000.);
139 hTmax_2t0_1 =
new TH1F (
name+
"_2t0_1",
"Tmax+2*Delta(t0); one t<5ns", 2000, -1000., 1000.);
140 hTmax_2t0_2 =
new TH1F (
name+
"_2t0_2",
"Tmax+2*Delta(t0); one t<10ns", 2000, -1000., 1000.);
141 hTmax_2t0_3 =
new TH1F (
name+
"_2t0_3",
"Tmax+2*Delta(t0); one t<20ns", 2000, -1000., 1000.);
142 hTmax_2t0_4 =
new TH1F (
name+
"_2t0_4",
"Tmax+2*Delta(t0); one t<50ns", 2000, -1000., 1000.);
143 hTmax_2t0_5 =
new TH1F (
name+
"_2t0_5",
"Tmax+2*Delta(t0); all t>50ns", 2000, -1000., 1000.);
144 hTmax_t0 =
new TH1F (
name+
"_t0",
"Tmax+Delta(t0)", 2000, -1000., 1000.);
145 hTmax_t0_0 =
new TH1F (
name+
"_t0_0",
"Tmax+Delta(t0); 3 hits", 2000, -1000., 1000.);
146 hTmax_t0_1 =
new TH1F (
name+
"_t0_1",
"Tmax+Delta(t0); one t<5ns", 2000, -1000., 1000.);
147 hTmax_t0_2 =
new TH1F (
name+
"_t0_2",
"Tmax+Delta(t0); one t<10ns", 2000, -1000., 1000.);
148 hTmax_t0_3 =
new TH1F (
name+
"_t0_3",
"Tmax+Delta(t0); one t<20ns", 2000, -1000., 1000.);
149 hTmax_t0_4 =
new TH1F (
name+
"_t0_4",
"Tmax+Delta(t0); one t<50ns", 2000, -1000., 1000.);
150 hTmax_t0_5 =
new TH1F (
name+
"_t0_5",
"Tmax+Delta(t0); all t>50ns", 2000, -1000., 1000.);
151 hTmax_0 =
new TH1F (
name+
"_0",
"Tmax", 2000, -1000., 1000.);
157 hTmax123 = (TH1F *) file->Get(
name+
"_Tmax123");
158 hTmax124s72 = (TH1F *) file->Get(
name+
"_Tmax124_s72");
159 hTmax124s78 = (TH1F *) file->Get(
name+
"_Tmax124_s78");
160 hTmax134s72 = (TH1F *) file->Get(
name+
"_Tmax134_s72");
161 hTmax134s78 = (TH1F *) file->Get(
name+
"_Tmax134_s78");
162 hTmax234 = (TH1F *) file->Get(
name+
"_Tmax234");
163 hTmax_3t0 = (TH1F *) file->Get(
name+
"_3t0");
164 hTmax_3t0_0 = (TH1F *) file->Get(
name+
"_3t0_0");
165 hTmax_3t0_1 = (TH1F *) file->Get(
name+
"_3t0_1");
166 hTmax_3t0_2 = (TH1F *) file->Get(
name+
"_3t0_2");
167 hTmax_3t0_3 = (TH1F *) file->Get(
name+
"_3t0_3");
168 hTmax_3t0_4 = (TH1F *) file->Get(
name+
"_3t0_4");
169 hTmax_3t0_5 = (TH1F *) file->Get(
name+
"_3t0_5");
170 hTmax_2t0 = (TH1F *) file->Get(
name+
"_2t0");
171 hTmax_2t0_0 = (TH1F *) file->Get(
name+
"_2t0_0");
172 hTmax_2t0_1 = (TH1F *) file->Get(
name+
"_2t0_1");
173 hTmax_2t0_2 = (TH1F *) file->Get(
name+
"_2t0_2");
174 hTmax_2t0_3 = (TH1F *) file->Get(
name+
"_2t0_3");
175 hTmax_2t0_4 = (TH1F *) file->Get(
name+
"_2t0_4");
176 hTmax_2t0_5 = (TH1F *) file->Get(
name+
"_2t0_5");
177 hTmax_t0 = (TH1F *) file->Get(
name+
"_t0");
178 hTmax_t0_1 = (TH1F *) file->Get(
name+
"_t0_1");
179 hTmax_t0_2 = (TH1F *) file->Get(
name+
"_t0_2");
180 hTmax_t0_3 = (TH1F *) file->Get(
name+
"_t0_3");
181 hTmax_t0_4 = (TH1F *) file->Get(
name+
"_t0_4");
182 hTmax_t0_5 = (TH1F *) file->Get(
name+
"_t0_5");
183 hTmax_0 = (TH1F *) file->Get(
name+
"_0");
220 void Fill(
float tmax123,
float tmax124,
float tmax134,
float tmax234,
231 hTmax123->Fill(tmax123);
233 hTmax_t0->Fill(tmax123);
235 case 0: hTmax_t0_0->Fill(tmax123);
break;
236 case 1: hTmax_t0_1->Fill(tmax123);
break;
237 case 2: hTmax_t0_2->Fill(tmax123);
break;
238 case 3: hTmax_t0_3->Fill(tmax123);
break;
239 case 4: hTmax_t0_4->Fill(tmax123);
break;
240 case 99: hTmax_t0_5->Fill(tmax123);
break;
244 hTmax_2t0->Fill(tmax123);
246 case 0: hTmax_2t0_0->Fill(tmax123);
break;
247 case 1: hTmax_2t0_1->Fill(tmax123);
break;
248 case 2: hTmax_2t0_2->Fill(tmax123);
break;
249 case 3: hTmax_2t0_3->Fill(tmax123);
break;
250 case 4: hTmax_2t0_4->Fill(tmax123);
break;
251 case 99: hTmax_2t0_5->Fill(tmax123);
break;
256 (s124==
dttmaxenums::r72)? hTmax124s72->Fill(tmax124):hTmax124s78->Fill(tmax124);
258 hTmax_0->Fill(tmax124);
260 hTmax_t0->Fill(tmax124);
262 case 0: hTmax_t0_0->Fill(tmax124);
break;
263 case 1: hTmax_t0_1->Fill(tmax124);
break;
264 case 2: hTmax_t0_2->Fill(tmax124);
break;
265 case 3: hTmax_t0_3->Fill(tmax124);
break;
266 case 4: hTmax_t0_4->Fill(tmax124);
break;
267 case 99: hTmax_t0_5->Fill(tmax124);
break;
270 else if(t0_124== 2) {
271 hTmax_2t0->Fill(tmax124);
273 case 0: hTmax_2t0_0->Fill(tmax124);
break;
274 case 1: hTmax_2t0_1->Fill(tmax124);
break;
275 case 2: hTmax_2t0_2->Fill(tmax124);
break;
276 case 3: hTmax_2t0_3->Fill(tmax124);
break;
277 case 4: hTmax_2t0_4->Fill(tmax124);
break;
278 case 99: hTmax_2t0_5->Fill(tmax124);
break;
282 hTmax_3t0->Fill(tmax124);
284 case 0: hTmax_3t0_0->Fill(tmax124);
break;
285 case 1: hTmax_3t0_1->Fill(tmax124);
break;
286 case 2: hTmax_3t0_2->Fill(tmax124);
break;
287 case 3: hTmax_3t0_3->Fill(tmax124);
break;
288 case 4: hTmax_3t0_4->Fill(tmax124);
break;
289 case 99: hTmax_3t0_5->Fill(tmax124);
break;
294 (s134==
dttmaxenums::r72)? hTmax134s72->Fill(tmax134):hTmax134s78->Fill(tmax134);
296 hTmax_0->Fill(tmax134);
298 hTmax_t0->Fill(tmax134);
300 case 0: hTmax_t0_0->Fill(tmax134);
break;
301 case 1: hTmax_t0_1->Fill(tmax134);
break;
302 case 2: hTmax_t0_2->Fill(tmax134);
break;
303 case 3: hTmax_t0_3->Fill(tmax134);
break;
304 case 4: hTmax_t0_4->Fill(tmax134);
break;
305 case 99: hTmax_t0_5->Fill(tmax134);
break;
308 else if(t0_134== 2) {
309 hTmax_2t0->Fill(tmax134);
311 case 0: hTmax_2t0_0->Fill(tmax134);
break;
312 case 1: hTmax_2t0_1->Fill(tmax134);
break;
313 case 2: hTmax_2t0_2->Fill(tmax134);
break;
314 case 3: hTmax_2t0_3->Fill(tmax134);
break;
315 case 4: hTmax_2t0_4->Fill(tmax134);
break;
316 case 99: hTmax_2t0_5->Fill(tmax134);
break;
320 hTmax_3t0->Fill(tmax134);
322 case 0: hTmax_3t0_0->Fill(tmax134);
break;
323 case 1: hTmax_3t0_1->Fill(tmax134);
break;
324 case 2: hTmax_3t0_2->Fill(tmax134);
break;
325 case 3: hTmax_3t0_3->Fill(tmax134);
break;
326 case 4: hTmax_3t0_4->Fill(tmax134);
break;
327 case 99: hTmax_3t0_5->Fill(tmax134);
break;
332 hTmax234->Fill(tmax234);
334 hTmax_t0->Fill(tmax234);
336 case 0: hTmax_t0_0->Fill(tmax234);
break;
337 case 1: hTmax_t0_1->Fill(tmax234);
break;
338 case 2: hTmax_t0_2->Fill(tmax234);
break;
339 case 3: hTmax_t0_3->Fill(tmax234);
break;
340 case 4: hTmax_t0_4->Fill(tmax234);
break;
341 case 99: hTmax_t0_5->Fill(tmax234);
break;
345 hTmax_2t0->Fill(tmax234);
347 case 0: hTmax_2t0_0->Fill(tmax234);
break;
348 case 1: hTmax_2t0_1->Fill(tmax234);
break;
349 case 2: hTmax_2t0_2->Fill(tmax234);
break;
350 case 3: hTmax_2t0_3->Fill(tmax234);
break;
351 case 4: hTmax_2t0_4->Fill(tmax234);
break;
352 case 99: hTmax_2t0_5->Fill(tmax234);
break;
361 hTmax124s72->Write();
362 hTmax124s78->Write();
363 hTmax134s72->Write();
364 hTmax134s78->Write();
367 hTmax_3t0_0->Write();
368 hTmax_3t0_1->Write();
369 hTmax_3t0_2->Write();
370 hTmax_3t0_3->Write();
371 hTmax_3t0_4->Write();
372 hTmax_3t0_5->Write();
374 hTmax_2t0_0->Write();
375 hTmax_2t0_1->Write();
376 hTmax_2t0_2->Write();
377 hTmax_2t0_3->Write();
378 hTmax_2t0_4->Write();
379 hTmax_2t0_5->Write();
394 if (hist == hTmax_3t0)
396 else if(hist == hTmax_2t0)
398 else if(hist == hTmax_t0)
400 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_)