35 |
// |
// |
36 |
} |
} |
37 |
|
|
38 |
|
void CaloEnergy::Delete(){ |
39 |
|
Clear(); |
40 |
|
delete this; |
41 |
|
} |
42 |
|
|
43 |
|
void CaloEnergy::UseLevel2(){ |
44 |
|
if ( clong ){ |
45 |
|
delete clong; |
46 |
|
clong = new CaloLong(L2); |
47 |
|
clong->SplitInto(0,22); |
48 |
|
// clong->HeavyTail(true); |
49 |
|
clong->HeavyTail(false); |
50 |
|
}; |
51 |
|
if ( cp ) delete cp; |
52 |
|
cp = NULL; |
53 |
|
} |
54 |
|
|
55 |
void CaloEnergy::UseCaloPreSampler(){ |
void CaloEnergy::UseCaloPreSampler(){ |
56 |
// |
// |
57 |
// use the presampler setting forcefitmode to 1000 means to force the DV routine to find the track inside the calorimeter using the "shower" approach developed for electrons |
// use the presampler setting forcefitmode to 1000 means to force the DV routine to find the track inside the calorimeter using the "shower" approach developed for electrons |
63 |
// cp->ForceCaloFit(); |
// cp->ForceCaloFit(); |
64 |
// cp->SetDebug(true); |
// cp->SetDebug(true); |
65 |
// cp->Process(); |
// cp->Process(); |
66 |
if ( clong ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer()); |
if ( clong ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer()); |
67 |
} |
} |
68 |
|
|
69 |
|
|
72 |
fLong = true; |
fLong = true; |
73 |
if ( !clong ){ |
if ( !clong ){ |
74 |
clong = new CaloLong(L2); |
clong = new CaloLong(L2); |
75 |
if ( cp ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer()); |
if ( cp ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer()); |
76 |
clong->SplitInto(0,22); |
clong->SplitInto(0,22); |
77 |
|
// clong->HeavyTail(true); |
78 |
|
clong->HeavyTail(false); |
79 |
}; |
}; |
80 |
// |
// |
81 |
} |
} |
98 |
// |
// |
99 |
indep = false; |
indep = false; |
100 |
// |
// |
101 |
|
fAllpl = true; |
102 |
fLong = false; |
fLong = false; |
103 |
fPl = 1; |
fPl = 1; |
104 |
fRad = -1; |
fRad = -1; |
106 |
clong = NULL; |
clong = NULL; |
107 |
x0max = -1.; |
x0max = -1.; |
108 |
// |
// |
109 |
|
multicol = false; |
110 |
|
// |
111 |
this->DefineGeometry(); |
this->DefineGeometry(); |
112 |
fXosel =true; |
fXosel =true; |
113 |
fXesel = true; |
fXesel = true; |
121 |
fXemin = 1000; |
fXemin = 1000; |
122 |
fYomin = 1000; |
fYomin = 1000; |
123 |
fYemin = 1000; |
fYemin = 1000; |
124 |
|
// |
125 |
|
this->UseCaloPreSampler(); // use it by default, to go back to "standard" mode use CaloEnergy::UseLevel2(). |
126 |
|
// |
127 |
} |
} |
|
|
|
|
void CaloEnergy::SetMinimumContainment(TString section, Int_t plane){ |
|
|
section.ToUpper(); |
|
|
if ( section.Contains("XO") ) fXomin = plane; |
|
|
if ( section.Contains("XE") ) fXemin = plane; |
|
|
if ( section.Contains("YO") ) fYomin = plane; |
|
|
if ( section.Contains("YE") ) fYemin = plane; |
|
|
} |
|
|
|
|
|
void CaloEnergy::SetMinimumContainment(Int_t plane){ |
|
|
this->SetMinimumContainment("XEXOYEYO",plane); |
|
|
} |
|
|
|
|
|
Int_t CaloEnergy::GetMinimumContainment(TString section){ |
|
|
section.ToUpper(); |
|
|
if ( section.Contains("XO") ) return(fXomin); |
|
|
if ( section.Contains("XE") ) return(fXemin); |
|
|
if ( section.Contains("YE") ) return(fYemin); |
|
|
if ( section.Contains("YO") ) return(fYomin); |
|
|
printf(" ERROR: section not recognized \n"); |
|
|
return(-1000); |
|
|
} |
|
|
|
|
|
void CaloEnergy::SetConversionFactor(TString section, Float_t conv){ |
|
|
section.ToUpper(); |
|
|
if ( section.Contains("XO") ) fConv_rxo = conv; |
|
|
if ( section.Contains("XE") ) fConv_rxe = conv; |
|
|
if ( section.Contains("YO") ) fConv_ryo = conv; |
|
|
if ( section.Contains("YE") ) fConv_rye = conv; |
|
|
} |
|
|
|
|
|
void CaloEnergy::SetConversionFactor(Float_t conv){ |
|
|
this->SetConversionFactor("XEXOYEYO",conv); |
|
|
} |
|
|
|
|
|
Float_t CaloEnergy::GetConversionFactor(TString section){ |
|
|
section.ToUpper(); |
|
|
if ( section.Contains("XO") ) return(fConv_rxo); |
|
|
if ( section.Contains("XE") ) return(fConv_rxe); |
|
|
if ( section.Contains("YO") ) return(fConv_ryo); |
|
|
if ( section.Contains("YE") ) return(fConv_rye); |
|
|
printf(" ERROR: section not recognized \n"); |
|
|
return(-1000.); |
|
|
} |
|
|
|
|
|
Int_t CaloEnergy::GetMaxplane(TString section){ |
|
|
section.ToUpper(); |
|
|
if ( section.Contains("XO") ) return fMax_planexo; |
|
|
if ( section.Contains("XE") ) return fMax_planexe; |
|
|
if ( section.Contains("YO") ) return fMax_planeyo; |
|
|
if ( section.Contains("YE") ) return fMax_planeye; |
|
|
return(-1); |
|
|
} |
|
|
|
|
|
Float_t CaloEnergy::GetEnergyAtMaxplane(TString section){ |
|
|
section.ToUpper(); |
|
|
if ( section.Contains("XO") ) return xomax_en; |
|
|
if ( section.Contains("XE") ) return xemax_en; |
|
|
if ( section.Contains("YO") ) return yomax_en; |
|
|
if ( section.Contains("YE") ) return yemax_en; |
|
|
return(-1); |
|
|
} |
|
|
|
|
|
Float_t CaloEnergy::GetMaxEnergy(TString section){ |
|
|
section.ToUpper(); |
|
|
if ( fLong ){ |
|
|
this->Process(section); |
|
|
return fXOen_maxplane; |
|
|
} else { |
|
|
if ( section.Contains("XO") ) return fXOen_maxplane; |
|
|
if ( section.Contains("XE") ) return fXEen_maxplane; |
|
|
if ( section.Contains("YO") ) return fYOen_maxplane; |
|
|
if ( section.Contains("YE") ) return fYEen_maxplane; |
|
|
}; |
|
|
return(-1); |
|
|
} |
|
|
|
|
|
Float_t CaloEnergy::GetMaxEnergy(){ |
|
|
if ( fLong ){ |
|
|
if ( debug ) printf(" oh! call process! with asntr %s and sntr %s \n",asntr.Data(),sntr.Data()); |
|
|
this->Process(asntr); |
|
|
}; |
|
|
return((fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane)); |
|
|
} |
|
|
|
|
128 |
void CaloEnergy::DefineGeometry(){ |
void CaloEnergy::DefineGeometry(){ |
129 |
// |
// |
130 |
// Use CaloStrip to determine once the position of border strips for each section |
// Use CaloStrip to determine once the position of border strips for each section |
131 |
// |
// |
132 |
// |
// fM = 2. + 0.096; // real position from cbar BUG the 0.096 is already taken into account in the border calculation made by Giovanna |
133 |
fM = 2. + 0.096; // real position from cbar |
fM = 2. ; // real position from cbar |
134 |
// fM1 = 2. - 0.122 - 0.096; // due to calculation of xe1 etc. |
// fM1 = 2. - 0.122 - 0.096; // due to calculation of xe1 etc. BUG! this way we count from the silicon border not from the silicon sensitive area |
135 |
fM1 = 2. - 0.122 + 0.096; // due to calculation of xe1 etc. |
fM1 = 2. - 0.122 - 0.096 + 0.096; // due to calculation of xe1 etc. |
136 |
if ( fM1 < 0. ) fM1 = 0.; |
if ( fM1 < 0. ) fM1 = 0.; |
137 |
// |
// |
138 |
CaloStrip *cs = new CaloStrip(fSimu); |
CaloStrip *cs = new CaloStrip(fSimu); |
140 |
// view y plane 0 strip 0 |
// view y plane 0 strip 0 |
141 |
cs->Set(1,0,0); |
cs->Set(1,0,0); |
142 |
xe1= cs->GetY(); |
xe1= cs->GetY(); |
|
// z1= cs->GetZ(); |
|
143 |
// view y plane 0 strip 31 |
// view y plane 0 strip 31 |
144 |
cs->Set(1,0,31); |
cs->Set(1,0,31); |
145 |
xe2= cs->GetY(); |
xe2= cs->GetY(); |
155 |
// view y plane 0 strip 95 |
// view y plane 0 strip 95 |
156 |
cs->Set(1,0,95); |
cs->Set(1,0,95); |
157 |
xe6= cs->GetY(); |
xe6= cs->GetY(); |
|
|
|
158 |
// view x plane 0 strip 0 |
// view x plane 0 strip 0 |
159 |
cs->Set(0,0,0); |
cs->Set(0,0,0); |
160 |
yo1= cs->GetX(); |
yo1= cs->GetX(); |
|
// z2= cs->GetZ(); |
|
161 |
// view x plane 0 strip 31 |
// view x plane 0 strip 31 |
162 |
cs->Set(0,0,31); |
cs->Set(0,0,31); |
163 |
yo2= cs->GetX(); |
yo2= cs->GetX(); |
173 |
// view x plane 0 strip 95 |
// view x plane 0 strip 95 |
174 |
cs->Set(0,0,95); |
cs->Set(0,0,95); |
175 |
yo6= cs->GetX(); |
yo6= cs->GetX(); |
|
|
|
|
|
|
176 |
// view y plane 1 strip 0 |
// view y plane 1 strip 0 |
177 |
cs->Set(1,1,0); |
cs->Set(1,1,0); |
178 |
xo1= cs->GetY(); |
xo1= cs->GetY(); |
|
// z3= cs->GetZ(); |
|
179 |
// view y plane 1 strip 31 |
// view y plane 1 strip 31 |
180 |
cs->Set(1,1,31); |
cs->Set(1,1,31); |
181 |
xo2= cs->GetY(); |
xo2= cs->GetY(); |
190 |
xo5= cs->GetY(); |
xo5= cs->GetY(); |
191 |
// view y plane 1 strip 95 |
// view y plane 1 strip 95 |
192 |
cs->Set(1,1,95); |
cs->Set(1,1,95); |
193 |
xo6= cs->GetY(); |
xo6= cs->GetY(); |
|
|
|
194 |
// view x plane 1 strip 0 |
// view x plane 1 strip 0 |
195 |
cs->Set(0,1,0); |
cs->Set(0,1,0); |
196 |
ye1= cs->GetX(); |
ye1= cs->GetX(); |
|
// z4= cs->GetZ(); |
|
197 |
// view x plane 1 strip 31 |
// view x plane 1 strip 31 |
198 |
cs->Set(0,1,31); |
cs->Set(0,1,31); |
199 |
ye2= cs->GetX(); |
ye2= cs->GetX(); |
209 |
// view x plane 1 strip 95 |
// view x plane 1 strip 95 |
210 |
cs->Set(0,1,95); |
cs->Set(0,1,95); |
211 |
ye6= cs->GetX(); |
ye6= cs->GetX(); |
|
|
|
212 |
// |
// |
213 |
for (Int_t p = 0; p<22; p ++){ |
for (Int_t p = 0; p<22; p ++){ |
214 |
for (Int_t v = 0; v<2; v++ ){ |
for (Int_t v = 0; v<2; v++ ){ |
247 |
yomax_en= 0.; |
yomax_en= 0.; |
248 |
yemax_en= 0.; |
yemax_en= 0.; |
249 |
// |
// |
250 |
memset(enstrip,0,2*22*96*(sizeof(Float_t))); |
memset(enstrip,0,2*22*96*(sizeof(Float_t))); |
251 |
en = 0.; |
en = 0.; |
252 |
view = 0; |
view = 0; |
253 |
plane = 0; |
plane = 0; |
269 |
memset(en_xop,0,11*sizeof(Float_t)); |
memset(en_xop,0,11*sizeof(Float_t)); |
270 |
memset(en_yop,0,11*sizeof(Float_t)); |
memset(en_yop,0,11*sizeof(Float_t)); |
271 |
// |
// |
272 |
|
fColumn = -1; |
273 |
|
fColXE = -1; |
274 |
|
fColXO = -1; |
275 |
|
fColYE = -1; |
276 |
|
fColYO = -1; |
277 |
|
memset(encol,0,2*3*sizeof(Float_t)); |
278 |
|
entot[0] = 0.; |
279 |
|
entot[1] = 0.; |
280 |
|
// |
281 |
|
X0pl = 0.76; |
282 |
|
// |
283 |
} |
} |
284 |
|
|
285 |
void CaloEnergy::Print(){ |
void CaloEnergy::Print(){ |
312 |
printf(" fM1 :.............. %f \n",fM1); |
printf(" fM1 :.............. %f \n",fM1); |
313 |
printf(" fRad :.............. %i \n",fRad); |
printf(" fRad :.............. %i \n",fRad); |
314 |
printf(" fPl :.............. %i \n",fPl); |
printf(" fPl :.............. %i \n",fPl); |
315 |
|
printf(" fColumn :.............. %i \n",fColumn); |
316 |
|
printf(" multicol:.............. %i \n",multicol); |
317 |
|
printf(" encol x :.............. %f \n",this->GetEncol(0)); |
318 |
|
printf(" encol y :.............. %f \n",this->GetEncol(1)); |
319 |
|
printf(" entot x :.............. %f \n",this->GetEntot(0)); |
320 |
|
printf(" entot y :.............. %f \n",this->GetEntot(1)); |
321 |
|
printf(" fColXE :.............. %i \n",fColXE); |
322 |
|
printf(" fColXO :.............. %i \n",fColXO); |
323 |
|
printf(" fColYE :.............. %i \n",fColYE); |
324 |
|
printf(" fColYO :.............. %i \n",fColYO); |
325 |
printf(" fConv_rxe ............. %f \n",fConv_rxe); |
printf(" fConv_rxe ............. %f \n",fConv_rxe); |
326 |
printf(" fConv_rxo ............. %f \n",fConv_rxo); |
printf(" fConv_rxo ............. %f \n",fConv_rxo); |
327 |
printf(" fConv_ryo ............. %f \n",fConv_ryo); |
printf(" fConv_ryo ............. %f \n",fConv_ryo); |
336 |
printf(" fYEen_maxplane:........ %f \n",fYEen_maxplane); |
printf(" fYEen_maxplane:........ %f \n",fYEen_maxplane); |
337 |
printf(" fYOen_maxplane:........ %f \n",fYOen_maxplane); |
printf(" fYOen_maxplane:........ %f \n",fYOen_maxplane); |
338 |
printf(" x0max :.............. %f \n",x0max); |
printf(" x0max :.............. %f \n",x0max); |
339 |
|
printf(" X0pl :.............. %f \n",X0pl); |
340 |
printf(" debug :.............. %i \n",debug); |
printf(" debug :.............. %i \n",debug); |
341 |
|
|
342 |
printf("========================================================================\n"); |
printf("========================================================================\n"); |
343 |
// |
// |
344 |
} |
} |
345 |
|
|
346 |
void CaloEnergy::Delete(){ |
void CaloEnergy::SetMinimumContainment(TString section, Int_t plane){ |
347 |
Clear(); |
section.ToUpper(); |
348 |
delete this; |
if ( section.Contains("XO") ) fXomin = plane; |
349 |
|
if ( section.Contains("XE") ) fXemin = plane; |
350 |
|
if ( section.Contains("YO") ) fYomin = plane; |
351 |
|
if ( section.Contains("YE") ) fYemin = plane; |
352 |
|
} |
353 |
|
|
354 |
|
void CaloEnergy::SetMinimumContainment(Int_t plane){ |
355 |
|
this->SetMinimumContainment("XEXOYEYO",plane); |
356 |
|
} |
357 |
|
|
358 |
|
void CaloEnergy::SetConversionFactor(TString section, Float_t conv){ |
359 |
|
section.ToUpper(); |
360 |
|
if ( section.Contains("XO") ) fConv_rxo = conv; |
361 |
|
if ( section.Contains("XE") ) fConv_rxe = conv; |
362 |
|
if ( section.Contains("YO") ) fConv_ryo = conv; |
363 |
|
if ( section.Contains("YE") ) fConv_rye = conv; |
364 |
|
} |
365 |
|
|
366 |
|
void CaloEnergy::SetConversionFactor(Float_t conv){ |
367 |
|
this->SetConversionFactor("XEXOYEYO",conv); |
368 |
|
} |
369 |
|
|
370 |
|
Int_t CaloEnergy::GetMinimumContainment(TString section){ |
371 |
|
section.ToUpper(); |
372 |
|
if ( section.Contains("XO") ) return(fXomin); |
373 |
|
if ( section.Contains("XE") ) return(fXemin); |
374 |
|
if ( section.Contains("YE") ) return(fYemin); |
375 |
|
if ( section.Contains("YO") ) return(fYomin); |
376 |
|
printf(" ERROR: section not recognized \n"); |
377 |
|
return(-1000); |
378 |
|
} |
379 |
|
|
380 |
|
Float_t CaloEnergy::GetConversionFactor(TString section){ |
381 |
|
section.ToUpper(); |
382 |
|
if ( section.Contains("XO") ) return(fConv_rxo); |
383 |
|
if ( section.Contains("XE") ) return(fConv_rxe); |
384 |
|
if ( section.Contains("YO") ) return(fConv_ryo); |
385 |
|
if ( section.Contains("YE") ) return(fConv_rye); |
386 |
|
printf(" ERROR: section not recognized \n"); |
387 |
|
return(-1000.); |
388 |
|
} |
389 |
|
|
390 |
|
Int_t CaloEnergy::GetMaxplane(TString section){ |
391 |
|
section.ToUpper(); |
392 |
|
if ( section.Contains("XO") ) return fMax_planexo; |
393 |
|
if ( section.Contains("XE") ) return fMax_planexe; |
394 |
|
if ( section.Contains("YO") ) return fMax_planeyo; |
395 |
|
if ( section.Contains("YE") ) return fMax_planeye; |
396 |
|
return(-1); |
397 |
|
} |
398 |
|
|
399 |
|
Int_t CaloEnergy::GetColumn(TString section){ |
400 |
|
section.ToUpper(); |
401 |
|
if ( section.Contains("XO") ) return fColXO; |
402 |
|
if ( section.Contains("XE") ) return fColXE; |
403 |
|
if ( section.Contains("YO") ) return fColYO; |
404 |
|
if ( section.Contains("YE") ) return fColYE; |
405 |
|
return(-1); |
406 |
|
} |
407 |
|
|
408 |
|
Float_t CaloEnergy::GetMipEnergyAtMaxplane(TString section){ |
409 |
|
printf(" WARNING: OBSOLETE METHOD, use GetMipEnergyAtMaxplane(TString) instead! \n"); |
410 |
|
return (this->GetEnergyAtMaxplane(section)); |
411 |
|
} |
412 |
|
|
413 |
|
Float_t CaloEnergy::GetEnergyAtMaxplane(TString section){ |
414 |
|
section.ToUpper(); |
415 |
|
if ( section.Contains("XO") ) return xomax_en; |
416 |
|
if ( section.Contains("XE") ) return xemax_en; |
417 |
|
if ( section.Contains("YO") ) return yomax_en; |
418 |
|
if ( section.Contains("YE") ) return yemax_en; |
419 |
|
return(-1); |
420 |
|
} |
421 |
|
|
422 |
|
Float_t CaloEnergy::GetMaxEnergy(TString section){ |
423 |
|
printf(" WARNING: OBSOLETE METHOD, use GetMipEnergy(TString) instead! \n"); |
424 |
|
return (this->GetMipEnergy(section)); |
425 |
|
} |
426 |
|
|
427 |
|
Float_t CaloEnergy::GetMipEnergy(TString section){ |
428 |
|
section.ToUpper(); |
429 |
|
if ( fLong ){ |
430 |
|
this->Process(section); |
431 |
|
return fXOen_maxplane; |
432 |
|
} else { |
433 |
|
if ( section.Contains("XO") ) return fXOen_maxplane; |
434 |
|
if ( section.Contains("XE") ) return fXEen_maxplane; |
435 |
|
if ( section.Contains("YO") ) return fYOen_maxplane; |
436 |
|
if ( section.Contains("YE") ) return fYEen_maxplane; |
437 |
|
}; |
438 |
|
return(-1); |
439 |
|
} |
440 |
|
|
441 |
|
Float_t CaloEnergy::GetEncol(Int_t i){ |
442 |
|
if ( fColumn > -1 && (((fXesel || fXosel)&&i==1) || ((fYesel || fYosel)&&i==0)) ){ |
443 |
|
Int_t t = -1; |
444 |
|
if ( i == 0 ){ |
445 |
|
if ( fColumn == 0 || fColumn == 3 || fColumn == 6 ) t = 0; |
446 |
|
if ( fColumn == 1 || fColumn == 4 || fColumn == 7 ) t = 1; |
447 |
|
if ( fColumn == 2 || fColumn == 5 || fColumn == 8 ) t = 2; |
448 |
|
} else { |
449 |
|
if ( fColumn == 0 || fColumn == 1 || fColumn == 2 ) t = 0; |
450 |
|
if ( fColumn == 3 || fColumn == 4 || fColumn == 5 ) t = 1; |
451 |
|
if ( fColumn == 6 || fColumn == 7 || fColumn == 8 ) t = 2; |
452 |
|
}; |
453 |
|
if ( debug ) printf(" encol: i %i t %i encol %f \n",i,t,encol[i][t]); |
454 |
|
return encol[i][t]; |
455 |
|
}; |
456 |
|
return(-1.); |
457 |
|
} |
458 |
|
|
459 |
|
Float_t CaloEnergy::GetMaxEnergy(){ |
460 |
|
printf(" WARNING: OBSOLETE METHOD, use GetMipEnergy() instead! \n"); |
461 |
|
return (this->GetMipEnergy()); |
462 |
|
} |
463 |
|
|
464 |
|
Float_t CaloEnergy::GetMipEnergy(){ |
465 |
|
if ( fLong ){ |
466 |
|
if ( debug ) printf(" oh! call process! with asntr %s and sntr %s \n",asntr.Data(),sntr.Data()); |
467 |
|
this->Process(asntr); |
468 |
|
}; |
469 |
|
return((fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane)); |
470 |
} |
} |
471 |
|
|
472 |
|
|
533 |
enstrip[view][plane][strip]=en; |
enstrip[view][plane][strip]=en; |
534 |
}; |
}; |
535 |
// |
// |
536 |
|
if ( debug && ((fM1+0.122-0.244*(Float_t)fRad) < 0.) ) printf("Error: (fM1+0.122-0.244*(Float_t)fRad) < 0. fM1 %f fRad %i %f \n",fM1,fRad,(fM1+0.122-0.244*(Float_t)fRad)); |
537 |
|
// |
538 |
|
// inclination factor (stolen from Daniele's code) |
539 |
|
// |
540 |
|
Float_t ytgx = 0; |
541 |
|
Float_t ytgy = 0; |
542 |
|
ytgx = 0.76 * cl2->tanx[0]; |
543 |
|
ytgy = 0.76 * cl2->tany[0]; |
544 |
|
X0pl = sqrt( pow(0.76,2.) + pow(ytgx,2.) + pow(ytgy,2.) ); |
545 |
|
// |
546 |
// sum energy plane by plane for each sections |
// sum energy plane by plane for each sections |
547 |
// |
// |
548 |
|
Float_t fen_xep[11]; |
549 |
|
Float_t fen_xop[11]; |
550 |
|
Float_t fen_yep[11]; |
551 |
|
Float_t fen_yop[11]; |
552 |
|
memset(fen_xep,0,11*sizeof(Float_t)); |
553 |
|
memset(fen_xop,0,11*sizeof(Float_t)); |
554 |
|
memset(fen_yep,0,11*sizeof(Float_t)); |
555 |
|
memset(fen_yop,0,11*sizeof(Float_t)); |
556 |
|
// |
557 |
for (Int_t i=0;i<11;i++){ |
for (Int_t i=0;i<11;i++){ |
558 |
for(strip=0; strip<96; strip++) { |
for(strip=0; strip<96; strip++) { |
559 |
|
fen_xep[i] += enstrip[1][2*i][strip]; |
560 |
|
fen_yop[i] += enstrip[0][2*i][strip]; |
561 |
|
fen_xop[i] += enstrip[1][(2*i)+1][strip]; |
562 |
|
fen_yep[i] += enstrip[0][(2*i)+1][strip]; |
563 |
if ( fRad < 0 ){ |
if ( fRad < 0 ){ |
564 |
// |
// |
565 |
// run over all the strips of the plane |
// run over all the strips of the plane |
566 |
// |
// |
567 |
en_xep[i] += enstrip[1][2*i][strip]; |
en_xep[i] += enstrip[1][2*i][strip]; |
568 |
en_yop[i] += enstrip[0][2*i][strip]; |
en_yop[i] += enstrip[0][2*i][strip]; |
572 |
// |
// |
573 |
// use only the strips inside a cylinder of given radius fRad |
// use only the strips inside a cylinder of given radius fRad |
574 |
// |
// |
575 |
if ( strip >= cl2->cibar[2*i][1]-1-fRad && strip <= cl2->cibar[2*i][1]-1+fRad ) en_xep[i] += enstrip[1][2*i][strip]; |
if ( cl2->cibar[2*i][1] >= 1 && cl2->cibar[2*i][1] <= 96 && |
576 |
if ( strip >= cl2->cibar[2*i][0]-1-fRad && strip <= cl2->cibar[2*i][0]-1+fRad ) en_yop[i] += enstrip[0][2*i][strip]; |
(strip >= (cl2->cibar[2*i][1]-1-fRad)) && (strip <= (cl2->cibar[2*i][1]-1+fRad)) ) en_xep[i] += enstrip[1][2*i][strip]; |
577 |
if ( strip >= cl2->cibar[(2*i)+1][1]-1-fRad && strip <= cl2->cibar[(2*i)+1][1]-1+fRad ) en_xop[i] += enstrip[1][(2*i)+1][strip]; |
|
578 |
if ( strip >= cl2->cibar[(2*i)+1][0]-1-fRad && strip <= cl2->cibar[(2*i)+1][0]-1+fRad ) en_yep[i] += enstrip[0][(2*i)+1][strip]; |
if ( cl2->cibar[2*i][0] >= 1 && cl2->cibar[2*i][0] <= 96 && |
579 |
|
(strip >= (cl2->cibar[2*i][0]-1-fRad)) && (strip <= (cl2->cibar[2*i][0]-1+fRad)) ) en_yop[i] += enstrip[0][2*i][strip]; |
580 |
|
|
581 |
|
if ( cl2->cibar[(2*i)+1][1] >= 1 && cl2->cibar[(2*i)+1][1] <= 96 && |
582 |
|
(strip >= (cl2->cibar[(2*i)+1][1]-1-fRad)) && (strip <= (cl2->cibar[(2*i)+1][1]-1+fRad)) ) en_xop[i] += enstrip[1][(2*i)+1][strip]; |
583 |
|
|
584 |
|
if ( cl2->cibar[(2*i)+1][0] >= 1 && cl2->cibar[(2*i)+1][0] <= 96 && |
585 |
|
(strip >= (cl2->cibar[(2*i)+1][0]-1-fRad)) && (strip <= (cl2->cibar[(2*i)+1][0]-1+fRad)) ) en_yep[i] += enstrip[0][(2*i)+1][strip]; |
586 |
}; |
}; |
587 |
}; |
}; |
588 |
|
if ( debug ) printf(" ex_xep[%i] %f cibar %i \n",i,en_xep[i],cl2->cibar[2*i][1]); |
589 |
|
if ( debug ) printf(" ex_xop[%i] %f cibar %i \n",i,en_xop[i],cl2->cibar[(2*i)+1][1]); |
590 |
|
if ( debug ) printf(" ex_yep[%i] %f cibar %i \n",i,en_yep[i],cl2->cibar[(2*i)+1][0]); |
591 |
|
if ( debug ) printf(" ex_yop[%i] %f cibar %i \n",i,en_yop[i],cl2->cibar[2*i][0]); |
592 |
energyxe += en_xep[i]; |
energyxe += en_xep[i]; |
593 |
energyyo += en_yop[i]; |
energyyo += en_yop[i]; |
594 |
energyxo += en_xop[i]; |
energyxo += en_xop[i]; |
602 |
Int_t yon = 0; |
Int_t yon = 0; |
603 |
Int_t xon = 0; |
Int_t xon = 0; |
604 |
Int_t yen = 0; |
Int_t yen = 0; |
605 |
|
Float_t en = 0.; |
606 |
// |
// |
607 |
if ( section.Contains("XE") ){ |
if ( section.Contains("XE") ){ |
608 |
yon++; |
yon++; |
609 |
xon++; |
xon++; |
610 |
yen++; |
yen++; |
611 |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
612 |
if(en_xep[ipl] > xemax_en) { |
en = fen_xep[ipl]; |
613 |
xemax_en = en_xep[ipl]; |
if ( !fAllpl ) en = en_xep[ipl]; |
614 |
|
if(en > xemax_en) { |
615 |
|
xemax_en = en; |
616 |
fMax_planexe = ipl; |
fMax_planexe = ipl; |
617 |
}; |
}; |
618 |
}; |
}; |
622 |
xon++; |
xon++; |
623 |
yen++; |
yen++; |
624 |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
625 |
if(en_yop[ipl] > yomax_en) { |
en = fen_yop[ipl]; |
626 |
yomax_en = en_yop[ipl]; |
if ( !fAllpl ) en = en_yop[ipl]; |
627 |
|
if(en > yomax_en) { |
628 |
|
yomax_en = en; |
629 |
fMax_planeyo = ipl; |
fMax_planeyo = ipl; |
630 |
}; |
}; |
631 |
}; |
}; |
634 |
if ( section.Contains("XO") ){ |
if ( section.Contains("XO") ){ |
635 |
yen++; |
yen++; |
636 |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
637 |
if(en_xop[ipl] > xomax_en) { |
en = fen_xop[ipl]; |
638 |
xomax_en = en_xop[ipl]; |
if ( !fAllpl ) en = en_xop[ipl]; |
639 |
|
if(en > xomax_en) { |
640 |
|
xomax_en = en; |
641 |
fMax_planexo = ipl; |
fMax_planexo = ipl; |
642 |
}; |
}; |
643 |
}; |
}; |
645 |
// |
// |
646 |
if ( section.Contains("YE") ){ |
if ( section.Contains("YE") ){ |
647 |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
for (Int_t ipl =0; ipl < 11; ipl ++) { |
648 |
if(en_yep[ipl] > yemax_en) { |
en = fen_yep[ipl]; |
649 |
yemax_en = en_yep[ipl]; |
if ( !fAllpl ) en = en_yep[ipl]; |
650 |
|
if(en > yemax_en) { |
651 |
|
yemax_en = en; |
652 |
fMax_planeye = ipl; |
fMax_planeye = ipl; |
653 |
}; |
}; |
654 |
}; |
}; |
658 |
// |
// |
659 |
if ( x0max > 0. ){ |
if ( x0max > 0. ){ |
660 |
if ( debug ) printf(" CALCULATE MAX PLANE GIVEN X0 ASSUMING PERPENDICULAR TRACK \n"); |
if ( debug ) printf(" CALCULATE MAX PLANE GIVEN X0 ASSUMING PERPENDICULAR TRACK \n"); |
661 |
Int_t wpl = (Int_t)roundf(x0max/0.76); |
// Int_t wpl = (Int_t)roundf(x0max/0.76); |
662 |
|
Int_t wpl = (Int_t)roundf(x0max/X0pl); |
663 |
Bool_t isY = false; |
Bool_t isY = false; |
664 |
if ( ((x0max/0.76)-(Float_t)wpl) > 0. ) isY = true; |
// if ( ((x0max/0.76)-(Float_t)wpl) > 0. ) isY = true; |
665 |
|
if ( ((x0max/X0pl)-(Float_t)wpl) > 0. ) isY = true; |
666 |
xomax_en = 0.; |
xomax_en = 0.; |
667 |
yemax_en = 0.; |
yemax_en = 0.; |
668 |
xemax_en = 0.; |
xemax_en = 0.; |
703 |
if ( section.Contains("YO") ) yomax_en = 10.; |
if ( section.Contains("YO") ) yomax_en = 10.; |
704 |
}; |
}; |
705 |
}; |
}; |
706 |
if ( debug ) printf(" x0max %f x0max/0.76 %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/0.76),wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo); |
// if ( debug ) printf(" x0max %f x0max/0.76 %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/0.76),wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo); |
707 |
|
if ( debug ) printf(" x0max %f x0max/X0pl %f X0pl %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/X0pl),X0pl,wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo); |
708 |
}; |
}; |
709 |
// |
// |
710 |
Int_t nPl = fPl; |
Int_t nPl = fPl; |
776 |
// |
// |
777 |
// for each plane of the calorimeter find the position of the track in the direction along the strip (where we do not have a measurement from the selected plane) by looking at the plane above/below of the other view and extrapolating the trajectory to the given plane |
// for each plane of the calorimeter find the position of the track in the direction along the strip (where we do not have a measurement from the selected plane) by looking at the plane above/below of the other view and extrapolating the trajectory to the given plane |
778 |
// |
// |
|
Float_t track_coordx[22][2]; |
|
|
Float_t track_coordy[22][2]; |
|
779 |
// |
// |
780 |
Float_t tgx_cl2; |
Float_t tgx_cl2; |
781 |
Float_t tgy_cl2; |
Float_t tgy_cl2; |
785 |
for (Int_t p=0; p<22; p++){ |
for (Int_t p=0; p<22; p++){ |
786 |
track_coordy[p][1] = cl2->cbar[p][1]; |
track_coordy[p][1] = cl2->cbar[p][1]; |
787 |
track_coordx[p][1] = cl2->cbar[p][0] - fabs(trk_z[p][1]-trk_z[p][0])*tgx_cl2; |
track_coordx[p][1] = cl2->cbar[p][0] - fabs(trk_z[p][1]-trk_z[p][0])*tgx_cl2; |
788 |
|
// track_coordx[p][1] = cl2->cbar[p][0] + fabs(trk_z[p][1]-trk_z[p][0])*tgx_cl2; |
789 |
track_coordx[p][0] = cl2->cbar[p][0]; |
track_coordx[p][0] = cl2->cbar[p][0]; |
790 |
track_coordy[p][0] = cl2->cbar[p][1] - fabs(trk_z[p][1]-trk_z[p][0])*tgy_cl2; |
track_coordy[p][0] = cl2->cbar[p][1] - fabs(trk_z[p][1]-trk_z[p][0])*tgy_cl2; |
791 |
|
// track_coordy[p][0] = cl2->cbar[p][1] + fabs(trk_z[p][1]-trk_z[p][0])*tgy_cl2; |
792 |
if ( debug ) printf(" p %i track_coordy[p][1] %f track_coordx[p][1] %f track_coordx[p][0] %f track_coordy[p][0] %f \n",p,track_coordy[p][1],track_coordx[p][1],track_coordx[p][0],track_coordy[p][0]); |
if ( debug ) printf(" p %i track_coordy[p][1] %f track_coordx[p][1] %f track_coordx[p][0] %f track_coordy[p][0] %f \n",p,track_coordy[p][1],track_coordx[p][1],track_coordx[p][0],track_coordy[p][0]); |
793 |
}; |
}; |
794 |
// |
// |
796 |
// |
// |
797 |
if ( section.Contains("XO") ){ |
if ( section.Contains("XO") ){ |
798 |
// |
// |
799 |
// check event is inside XO acceptance |
// find the column hit in the first plane |
800 |
|
// |
801 |
|
Int_t ix = -1; |
802 |
|
Int_t iy = -1; |
803 |
|
if ( track_coordx[(2*0)+1][1] >= (-12.054+fM) && track_coordx[(2*0)+1][1] <= (-4.246-fM) ) ix = 0; |
804 |
|
if ( track_coordx[(2*0)+1][1] >= ( -4.004+fM) && track_coordx[(2*0)+1][1] <= ( 3.804-fM) ) ix = 1; |
805 |
|
if ( track_coordx[(2*0)+1][1] >= ( 4.046+fM) && track_coordx[(2*0)+1][1] <= (11.854-fM) ) ix = 2; |
806 |
|
if ( cl2->cbar[(2*0)+1][1] >= (xo1 + fM1) && cl2->cbar[(2*0)+1][1] <= (xo2 - fM1) ) iy = 0; |
807 |
|
if ( cl2->cbar[(2*0)+1][1] >= (xo3 + fM1) && cl2->cbar[(2*0)+1][1] <= (xo4 - fM1) ) iy = 1; |
808 |
|
if ( cl2->cbar[(2*0)+1][1] >= (xo5 + fM1) && cl2->cbar[(2*0)+1][1] <= (xo6 - fM1) ) iy = 2; |
809 |
|
if ( ix > -1 && iy > -1 ) fColXO = ix + iy*3; |
810 |
|
// |
811 |
|
// check event is inside XO acceptance, if multicol is false (SingleColumn mode) then the track must be contained in a column. |
812 |
// |
// |
813 |
for (Int_t i=0; i<11; i++) { |
for (Int_t i=0; i<11; i++) { |
814 |
if |
if (( |
815 |
( |
( track_coordx[(2*i)+1][1] >= (-12.054+fM) && track_coordx[(2*i)+1][1] <= (-4.246-fM) && (ix == 0 || multicol) ) || |
816 |
(((track_coordx[(2*i)+1][1]>=(-12.054+fM))&& |
( track_coordx[(2*i)+1][1] >= ( -4.004+fM) && track_coordx[(2*i)+1][1] <= ( 3.804-fM) && (ix == 1 || multicol) ) || |
817 |
(track_coordx[(2*i)+1][1]<=(-4.246-fM)))&& |
( track_coordx[(2*i)+1][1] >= ( 4.046+fM) && track_coordx[(2*i)+1][1] <= (11.854-fM) && (ix == 2 || multicol) ) |
818 |
(((cl2->cbar[(2*i)+1][1]>=xo1 + fM1)&& |
) && ( |
819 |
(cl2->cbar[(2*i)+1][1]<=xo2 - fM1 ))|| |
( cl2->cbar[(2*i)+1][1] >= (xo1 + fM1) && cl2->cbar[(2*i)+1][1] <= (xo2 - fM1) && (iy == 0 || multicol) ) || |
820 |
((cl2->cbar[(2*i)+1][1]>=xo3 + fM1)&& |
( cl2->cbar[(2*i)+1][1] >= (xo3 + fM1) && cl2->cbar[(2*i)+1][1] <= (xo4 - fM1) && (iy == 1 || multicol) ) || |
821 |
(cl2->cbar[(2*i)+1][1]<=xo4 - fM1 ))|| |
( cl2->cbar[(2*i)+1][1] >= (xo5 + fM1) && cl2->cbar[(2*i)+1][1] <= (xo6 - fM1) && (iy == 2 || multicol) ) |
822 |
((cl2->cbar[(2*i)+1][1]>=xo5 + fM1)&& |
)){ |
|
(cl2->cbar[(2*i)+1][1]<=xo6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordx[(2*i)+1][1]>=(-4.004+fM))&& |
|
|
(track_coordx[(2*i)+1][1]<=(3.804-fM)))&& |
|
|
(((cl2->cbar[(2*i)+1][1]>=xo1 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][1]<=xo2 - fM1 ))|| |
|
|
((cl2->cbar[(2*i)+1][1]>=xo3 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][1]<=xo4 - fM1))|| |
|
|
((cl2->cbar[(2*i)+1][1]>=xo5 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][1]<=xo6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordx[(2*i)+1][1]>=(4.046+fM))&& |
|
|
(track_coordx[(2*i)+1][1]<=(11.854-fM)))&& |
|
|
(((cl2->cbar[(2*i)+1][1]>=xo1 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][1]<=xo2 - fM1))|| |
|
|
((cl2->cbar[(2*i)+1][1]>=xo3 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][1]<=xo4 - fM1 ))|| |
|
|
((cl2->cbar[(2*i)+1][1]>=xo5 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][1]<=xo6 - fM1 )))) |
|
|
){ |
|
823 |
fXosel = true; |
fXosel = true; |
824 |
fXoout = i; |
fXoout = i; |
825 |
} else { |
} else { |
839 |
// event is contained (or partially contained) hence we can integrate energy up to the maximum and calculate the energy as measured by this section |
// event is contained (or partially contained) hence we can integrate energy up to the maximum and calculate the energy as measured by this section |
840 |
// |
// |
841 |
if ( fXosel ){ |
if ( fXosel ){ |
842 |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planexo+nPl)) ;iplm++){ |
for (Int_t iplm=0; iplm<=TMath::Min(10,(Int_t)(fMax_planexo+nPl)); iplm++){ |
843 |
fXOen_maxplane += en_xop[iplm]; |
fXOen_maxplane += en_xop[iplm]; |
844 |
if ( debug ) printf(" XO iplm %i fXOen_maxplane %f en_xop[iplm] %f\n",iplm,fXOen_maxplane,en_xop[iplm]); |
if ( debug ) printf(" XO iplm %i fXOen_maxplane %f en_xop[iplm] %f\n",iplm,fXOen_maxplane,en_xop[iplm]); |
845 |
}; |
}; |
846 |
fEnergyxo = fXOen_maxplane/fConv_rxo; |
fEnergyxo = fXOen_maxplane/fConv_rxo; |
847 |
|
// |
848 |
|
for (Int_t i=0;i<11;i++){ |
849 |
|
for(strip=0; strip<96; strip++) { |
850 |
|
// |
851 |
|
// run over all the strips of the plane |
852 |
|
// |
853 |
|
if ( strip >= 0 && strip < 32 ) encol[1][0] += enstrip[1][(2*i)+1][strip]; |
854 |
|
if ( strip >= 32 && strip < 64 ) encol[1][1] += enstrip[1][(2*i)+1][strip]; |
855 |
|
if ( strip >= 64 && strip < 96 ) encol[1][2] += enstrip[1][(2*i)+1][strip]; |
856 |
|
entot[1] += enstrip[1][(2*i)+1][strip]; |
857 |
|
// |
858 |
|
}; |
859 |
|
}; |
860 |
}; |
}; |
861 |
}; |
}; |
862 |
// |
// |
863 |
if ( section.Contains("XE") ){ |
if ( section.Contains("XE") ){ |
864 |
|
// |
865 |
|
// find the column hit in the first plane |
866 |
|
// |
867 |
|
Int_t ix = -1; |
868 |
|
Int_t iy = -1; |
869 |
|
if ( track_coordx[(2*0)][1] >= (-11.854+fM) && track_coordx[(2*0)][1] <= (-4.046-fM) ) ix = 0; |
870 |
|
if ( track_coordx[(2*0)][1] >= ( -3.804+fM) && track_coordx[(2*0)][1] <= ( 4.004-fM) ) ix = 1; |
871 |
|
if ( track_coordx[(2*0)][1] >= ( 4.246+fM) && track_coordx[(2*0)][1] <= (12.054-fM) ) ix = 2; |
872 |
|
if ( cl2->cbar[(2*0)][1] >= (xe1 + fM1) && cl2->cbar[(2*0)][1] <= (xe2 - fM1) ) iy = 0; |
873 |
|
if ( cl2->cbar[(2*0)][1] >= (xe3 + fM1) && cl2->cbar[(2*0)][1] <= (xe4 - fM1) ) iy = 1; |
874 |
|
if ( cl2->cbar[(2*0)][1] >= (xe5 + fM1) && cl2->cbar[(2*0)][1] <= (xe6 - fM1) ) iy = 2; |
875 |
|
if ( ix > -1 && iy > -1 ) fColXE = ix + iy*3; |
876 |
|
// |
877 |
|
// check event is inside XO acceptance |
878 |
|
// |
879 |
for (Int_t i=0; i<11; i++) { |
for (Int_t i=0; i<11; i++) { |
880 |
if |
if (( |
881 |
( |
( track_coordx[(2*i)][1] >= (-11.854+fM) && track_coordx[(2*i)][1] <= (-4.046-fM) && (ix == 0 || multicol) ) || |
882 |
(((track_coordx[2*i][1]>=(-11.854+fM))&& |
( track_coordx[(2*i)][1] >= ( -3.804+fM) && track_coordx[(2*i)][1] <= ( 4.004-fM) && (ix == 1 || multicol) ) || |
883 |
(track_coordx[2*i][1]<=(-4.046-fM)))&& |
( track_coordx[(2*i)][1] >= ( 4.246+fM) && track_coordx[(2*i)][1] <= (12.054-fM) && (ix == 2 || multicol) ) |
884 |
(((cl2->cbar[2*i][1]>=xe1 + fM1)&& |
) && ( |
885 |
(cl2->cbar[2*i][1]<=xe2 - fM1 ))|| |
( cl2->cbar[(2*i)][1] >= (xe1 + fM1) && cl2->cbar[(2*i)][1] <= (xe2 - fM1) && (iy == 0 || multicol) ) || |
886 |
((cl2->cbar[2*i][1]>=xe3 + fM1)&& |
( cl2->cbar[(2*i)][1] >= (xe3 + fM1) && cl2->cbar[(2*i)][1] <= (xe4 - fM1) && (iy == 1 || multicol) ) || |
887 |
(cl2->cbar[2*i][1]<=xe4 - fM1 ))|| |
( cl2->cbar[(2*i)][1] >= (xe5 + fM1) && cl2->cbar[(2*i)][1] <= (xe6 - fM1) && (iy == 2 || multicol) ) |
888 |
((cl2->cbar[2*i][1]>=xe5 + fM1)&& |
)){ |
|
(cl2->cbar[2*i][1]<=xe6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordx[2*i][1]>=(-3.804+fM))&& |
|
|
(track_coordx[2*i][1]<=(4.004-fM)))&& |
|
|
(((cl2->cbar[2*i][1]>=xe1 + fM1)&& |
|
|
(cl2->cbar[2*i][1]<=xe2 - fM1 ))|| |
|
|
((cl2->cbar[2*i][1]>=xe3 + fM1)&& |
|
|
(cl2->cbar[2*i][1]<=xe4 - fM1))|| |
|
|
((cl2->cbar[2*i][1]>=xe5 + fM1)&& |
|
|
(cl2->cbar[2*i][1]<=xe6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordx[2*i][1]>=(4.246+fM))&& |
|
|
(track_coordx[2*i][1]<=(12.054-fM)))&& |
|
|
(((cl2->cbar[2*i][1]>=xe1 + fM1)&& |
|
|
(cl2->cbar[2*i][1]<=xe2 - fM1))|| |
|
|
((cl2->cbar[2*i][1]>=xe3 + fM1)&& |
|
|
(cl2->cbar[2*i][1]<=xe4 - fM1 ))|| |
|
|
((cl2->cbar[2*i][1]>=xe5 + fM1)&& |
|
|
(cl2->cbar[2*i][1]<=xe6 - fM1 )))) |
|
|
){ |
|
889 |
fXesel = true; |
fXesel = true; |
890 |
fXeout = i; |
fXeout = i; |
891 |
} else { |
} else { |
893 |
break; |
break; |
894 |
}; |
}; |
895 |
}; |
}; |
896 |
|
// |
897 |
if ( !fXesel && fXeout >= fXemin && fXeout >= (Int_t)(fMax_planexe+nPl) ){ |
if ( !fXesel && fXeout >= fXemin && fXeout >= (Int_t)(fMax_planexe+nPl) ){ |
898 |
if ( debug ) printf(" XE: this event is only partially contained: fXeout %i fXemin %i fMax_planexe + nPl %i \n",fXeout,fXemin,(Int_t)(fMax_planexe+nPl)); |
if ( debug ) printf(" XE: this event is only partially contained: fXeout %i fXemin %i fMax_planexe + nPl %i \n",fXeout,fXemin,(Int_t)(fMax_planexe+nPl)); |
899 |
fPartsel = true; |
fPartsel = true; |
905 |
if ( debug ) printf(" XE iplm %i fXOen_maxplane %f en_xop[iplm] %f\n",iplm,fXEen_maxplane,en_xep[iplm]); |
if ( debug ) printf(" XE iplm %i fXOen_maxplane %f en_xop[iplm] %f\n",iplm,fXEen_maxplane,en_xep[iplm]); |
906 |
}; |
}; |
907 |
fEnergyxe = fXEen_maxplane/fConv_rxe; |
fEnergyxe = fXEen_maxplane/fConv_rxe; |
908 |
|
// |
909 |
|
for (Int_t i=0;i<11;i++){ |
910 |
|
for(strip=0; strip<96; strip++) { |
911 |
|
// |
912 |
|
// run over all the strips of the plane |
913 |
|
// |
914 |
|
if ( strip >= 0 && strip < 32 ) encol[1][0] += enstrip[1][(2*i)][strip]; |
915 |
|
if ( strip >= 32 && strip < 64 ) encol[1][1] += enstrip[1][(2*i)][strip]; |
916 |
|
if ( strip >= 64 && strip < 96 ) encol[1][2] += enstrip[1][(2*i)][strip]; |
917 |
|
entot[1] += enstrip[1][(2*i)][strip]; |
918 |
|
// |
919 |
|
}; |
920 |
|
}; |
921 |
}; |
}; |
922 |
}; |
}; |
923 |
// |
// |
924 |
if ( section.Contains("YE") ){ |
if ( section.Contains("YE") ){ |
925 |
|
// |
926 |
|
// find the column hit in the first plane |
927 |
|
// |
928 |
|
Int_t ix = -1; |
929 |
|
Int_t iy = -1; |
930 |
|
if ( track_coordy[(2*0)+1][0] >= (-12.154+fM) && track_coordy[(2*0)+1][0] <= (-4.346-fM) ) iy = 0; |
931 |
|
if ( track_coordy[(2*0)+1][0] >= ( -4.104+fM) && track_coordy[(2*0)+1][0] <= ( 3.704-fM) ) iy = 1; |
932 |
|
if ( track_coordy[(2*0)+1][0] >= ( 3.946+fM) && track_coordy[(2*0)+1][0] <= (11.754-fM) ) iy = 2; |
933 |
|
if ( cl2->cbar[(2*0)+1][0] >= (ye1 + fM1) && cl2->cbar[(2*0)+1][0] <= (ye2 - fM1) ) ix = 0; |
934 |
|
if ( cl2->cbar[(2*0)+1][0] >= (ye3 + fM1) && cl2->cbar[(2*0)+1][0] <= (ye4 - fM1) ) ix = 1; |
935 |
|
if ( cl2->cbar[(2*0)+1][0] >= (ye5 + fM1) && cl2->cbar[(2*0)+1][0] <= (ye6 - fM1) ) ix = 2; |
936 |
|
if ( ix > -1 && iy > -1 ) fColYE = ix + iy*3; |
937 |
|
// |
938 |
|
// check event is inside XO acceptance |
939 |
|
// |
940 |
for (Int_t i=0; i<11; i++) { |
for (Int_t i=0; i<11; i++) { |
941 |
if |
if (( |
942 |
( |
( track_coordy[(2*i)+1][0] >= (-12.154+fM) && track_coordy[(2*i)+1][0] <= (-4.346-fM) && (iy == 0 || multicol) ) || |
943 |
(((track_coordy[(2*i)+1][0]>=(-12.154+fM))&& |
( track_coordy[(2*i)+1][0] >= ( -4.104+fM) && track_coordy[(2*i)+1][0] <= ( 3.704-fM) && (iy == 1 || multicol) ) || |
944 |
(track_coordy[(2*i)+1][0]<=(-4.346-fM)))&& |
( track_coordy[(2*i)+1][0] >= ( 3.946+fM) && track_coordy[(2*i)+1][0] <= (11.754-fM) && (iy == 2 || multicol) ) |
945 |
(((cl2->cbar[(2*i)+1][0]>=ye1 + fM1)&& |
) && ( |
946 |
(cl2->cbar[(2*i)+1][0]<=ye2 - fM1 ))|| |
( cl2->cbar[(2*i)+1][0] >= (ye1 + fM1) && cl2->cbar[(2*i)+1][0] <= (ye2 - fM1) && (ix == 0 || multicol) ) || |
947 |
((cl2->cbar[(2*i)+1][0]>=ye3 + fM1)&& |
( cl2->cbar[(2*i)+1][0] >= (ye3 + fM1) && cl2->cbar[(2*i)+1][0] <= (ye4 - fM1) && (ix == 1 || multicol) ) || |
948 |
(cl2->cbar[(2*i)+1][0]<=ye4 - fM1 ))|| |
( cl2->cbar[(2*i)+1][0] >= (ye5 + fM1) && cl2->cbar[(2*i)+1][0] <= (ye6 - fM1) && (ix == 2 || multicol) ) |
949 |
((cl2->cbar[(2*i)+1][0]>=ye5 + fM1)&& |
)){ |
|
(cl2->cbar[(2*i)+1][0]<=ye6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordy[(2*i)+1][0]>=(-4.104+fM))&& |
|
|
(track_coordy[(2*i)+1][0]<=(3.704-fM)))&& |
|
|
(((cl2->cbar[(2*i)+1][0]>=ye1 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][0]<=ye2 - fM1 ))|| |
|
|
((cl2->cbar[(2*i)+1][0]>=ye3 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][0]<=ye4 - fM1))|| |
|
|
((cl2->cbar[(2*i)+1][0]>=ye5 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][0]<=ye6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordy[(2*i)+1][0]>=(3.946+fM))&& |
|
|
(track_coordy[(2*i)+1][0]<=(11.754-fM)))&& |
|
|
(((cl2->cbar[(2*i)+1][0]>=ye1 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][0]<=ye2 - fM1))|| |
|
|
((cl2->cbar[(2*i)+1][0]>=ye3 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][0]<=ye4 - fM1 ))|| |
|
|
((cl2->cbar[(2*i)+1][0]>=ye5 + fM1)&& |
|
|
(cl2->cbar[(2*i)+1][0]<=ye6 - fM1 )))) |
|
|
){ |
|
950 |
fYesel = true; |
fYesel = true; |
951 |
fYeout = i; |
fYeout = i; |
952 |
} else { |
} else { |
954 |
break; |
break; |
955 |
}; |
}; |
956 |
}; |
}; |
957 |
|
// |
958 |
if ( !fYesel && fYeout >= fYemin && fYeout >= (Int_t)(fMax_planeye+nPl) ){ |
if ( !fYesel && fYeout >= fYemin && fYeout >= (Int_t)(fMax_planeye+nPl) ){ |
959 |
if ( debug ) printf(" YE: this event is only partially contained: fYeout %i fYemin %i fMax_planeye + nPl %i \n",fYeout,fYemin,(Int_t)(fMax_planeye+nPl)); |
if ( debug ) printf(" YE: this event is only partially contained: fYeout %i fYemin %i fMax_planeye + nPl %i \n",fYeout,fYemin,(Int_t)(fMax_planeye+nPl)); |
960 |
fPartsel = true; |
fPartsel = true; |
963 |
if ( fYesel ){ |
if ( fYesel ){ |
964 |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeye+nPl)) ;iplm++) fYEen_maxplane += en_yep[iplm]; |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeye+nPl)) ;iplm++) fYEen_maxplane += en_yep[iplm]; |
965 |
fEnergyye = fYEen_maxplane/fConv_rye; |
fEnergyye = fYEen_maxplane/fConv_rye; |
966 |
|
// |
967 |
|
for (Int_t i=0;i<11;i++){ |
968 |
|
for(strip=0; strip<96; strip++) { |
969 |
|
// |
970 |
|
// run over all the strips of the plane |
971 |
|
// |
972 |
|
if ( strip >= 0 && strip < 32 ) encol[0][0] += enstrip[0][(2*i)+1][strip]; |
973 |
|
if ( strip >= 32 && strip < 64 ) encol[0][1] += enstrip[0][(2*i)+1][strip]; |
974 |
|
if ( strip >= 64 && strip < 96 ) encol[0][2] += enstrip[0][(2*i)+1][strip]; |
975 |
|
entot[0] += enstrip[0][(2*i)+1][strip]; |
976 |
|
// |
977 |
|
}; |
978 |
|
}; |
979 |
|
// |
980 |
}; |
}; |
981 |
}; |
}; |
982 |
// |
// |
983 |
if ( section.Contains("YO") ){ |
if ( section.Contains("YO") ){ |
984 |
|
// |
985 |
|
// find the column hit in the first plane |
986 |
|
// |
987 |
|
Int_t ix = -1; |
988 |
|
Int_t iy = -1; |
989 |
|
if ( track_coordy[(2*0)][0] >= (-11.954+fM) && track_coordy[(2*0)][0] <= (-4.146-fM) ) iy = 0; |
990 |
|
if ( track_coordy[(2*0)][0] >= ( -3.904+fM) && track_coordy[(2*0)][0] <= ( 3.904-fM) ) iy = 1; |
991 |
|
if ( track_coordy[(2*0)][0] >= ( 4.146+fM) && track_coordy[(2*0)][0] <= (11.954-fM) ) iy = 2; |
992 |
|
if ( cl2->cbar[(2*0)][0] >= (yo1 + fM1) && cl2->cbar[(2*0)][0] <= (yo2 - fM1) ) ix = 0; |
993 |
|
if ( cl2->cbar[(2*0)][0] >= (yo3 + fM1) && cl2->cbar[(2*0)][0] <= (yo4 - fM1) ) ix = 1; |
994 |
|
if ( cl2->cbar[(2*0)][0] >= (yo5 + fM1) && cl2->cbar[(2*0)][0] <= (yo6 - fM1) ) ix = 2; |
995 |
|
if ( ix > -1 && iy > -1 ) fColYO = ix + iy*3; |
996 |
|
// |
997 |
|
// check event is inside XO acceptance |
998 |
|
// |
999 |
for (Int_t i=0; i<11; i++) { |
for (Int_t i=0; i<11; i++) { |
1000 |
if ( debug ) printf(" i %i track_coordy[2*i][0] %f cl2->cbar[2*i][0] %f \n",i,track_coordy[2*i][0],cl2->cbar[2*i][0]); |
if (( |
1001 |
if ( debug ) printf(" i %i fm %f fm1 %f yo1 %g yo2 %f yo3 %f yo4 %f yo5 %f yo6 %f \n",i,fM,fM1,yo1,yo2,yo3,yo4,yo5,yo6); |
( track_coordy[(2*i)][0] >= (-11.954+fM) && track_coordy[(2*i)][0] <= (-4.146-fM) && (iy == 0 || multicol) ) || |
1002 |
if |
( track_coordy[(2*i)][0] >= ( -3.904+fM) && track_coordy[(2*i)][0] <= ( 3.904-fM) && (iy == 1 || multicol) ) || |
1003 |
( |
( track_coordy[(2*i)][0] >= ( 4.146+fM) && track_coordy[(2*i)][0] <= (11.954-fM) && (iy == 2 || multicol) ) |
1004 |
(((track_coordy[2*i][0]>=(-11.954+fM))&& |
) && ( |
1005 |
(track_coordy[2*i][0]<=(-4.146-fM)))&& |
( cl2->cbar[(2*i)][0] >= (yo1 + fM1) && cl2->cbar[(2*i)][0] <= (yo2 - fM1) && (ix == 0 || multicol) ) || |
1006 |
(((cl2->cbar[2*i][0]>=yo1 + fM1)&& |
( cl2->cbar[(2*i)][0] >= (yo3 + fM1) && cl2->cbar[(2*i)][0] <= (yo4 - fM1) && (ix == 1 || multicol) ) || |
1007 |
(cl2->cbar[2*i][0]<=yo2 - fM1 ))|| |
( cl2->cbar[(2*i)][0] >= (yo5 + fM1) && cl2->cbar[(2*i)][0] <= (yo6 - fM1) && (ix == 2 || multicol) ) |
1008 |
((cl2->cbar[2*i][0]>=yo3 + fM1)&& |
)){ |
|
(cl2->cbar[2*i][0]<=yo4 - fM1 ))|| |
|
|
((cl2->cbar[2*i][0]>=yo5 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordy[2*i][0]>=(-3.904+fM))&& |
|
|
(track_coordy[2*i][0]<=(+3.904-fM)))&& |
|
|
(((cl2->cbar[2*i][0]>=yo1 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo2 - fM1 ))|| |
|
|
((cl2->cbar[2*i][0]>=yo3 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo4 - fM1))|| |
|
|
((cl2->cbar[2*i][0]>=yo5 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo6 - fM1 ))))|| |
|
|
|
|
|
(((track_coordy[2*i][0]>=(4.146+fM))&& |
|
|
(track_coordy[2*i][0]<=(11.954-fM)))&& |
|
|
(((cl2->cbar[2*i][0]>=yo1 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo2 - fM1))|| |
|
|
((cl2->cbar[2*i][0]>=yo3 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo4 - fM1 ))|| |
|
|
((cl2->cbar[2*i][0]>=yo5 + fM1)&& |
|
|
(cl2->cbar[2*i][0]<=yo6 - fM1 )))) |
|
|
){ |
|
1009 |
fYosel = true; |
fYosel = true; |
1010 |
fYoout = i; |
fYoout = i; |
1011 |
} else { |
} else { |
1013 |
break; |
break; |
1014 |
}; |
}; |
1015 |
}; |
}; |
1016 |
|
// |
1017 |
if ( !fYosel && fYoout >= fYomin && fYoout >= (Int_t)(fMax_planeyo+nPl) ){ |
if ( !fYosel && fYoout >= fYomin && fYoout >= (Int_t)(fMax_planeyo+nPl) ){ |
1018 |
if ( debug ) printf(" YO: this event is only partially contained: fYoout %i fYomin %i fMax_planeyo + nPl %i \n",fYoout,fYomin,(Int_t)(fMax_planeyo+nPl)); |
if ( debug ) printf(" YO: this event is only partially contained: fYoout %i fYomin %i fMax_planeyo + nPl %i \n",fYoout,fYomin,(Int_t)(fMax_planeyo+nPl)); |
1019 |
fPartsel = true; |
fPartsel = true; |
1022 |
if ( fYosel ){ |
if ( fYosel ){ |
1023 |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeyo+nPl)) ;iplm++) fYOen_maxplane += en_yop[iplm]; |
for (Int_t iplm=0;iplm<=TMath::Min(10,(Int_t)(fMax_planeyo+nPl)) ;iplm++) fYOen_maxplane += en_yop[iplm]; |
1024 |
fEnergyyo = fYOen_maxplane/fConv_ryo; |
fEnergyyo = fYOen_maxplane/fConv_ryo; |
1025 |
|
// |
1026 |
|
for (Int_t i=0;i<11;i++){ |
1027 |
|
for(strip=0; strip<96; strip++) { |
1028 |
|
// |
1029 |
|
// run over all the strips of the plane |
1030 |
|
// |
1031 |
|
if ( strip >= 0 && strip < 32 ) encol[0][0] += enstrip[0][(2*i)][strip]; |
1032 |
|
if ( strip >= 32 && strip < 64 ) encol[0][1] += enstrip[0][(2*i)][strip]; |
1033 |
|
if ( strip >= 64 && strip < 96 ) encol[0][2] += enstrip[0][(2*i)][strip]; |
1034 |
|
entot[0] += enstrip[0][(2*i)][strip]; |
1035 |
|
// |
1036 |
|
}; |
1037 |
|
}; |
1038 |
}; |
}; |
1039 |
}; |
}; |
1040 |
// |
// |
1060 |
fSel = true; |
fSel = true; |
1061 |
}; |
}; |
1062 |
fEnergy = (fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane)/fConv_rxo; |
fEnergy = (fXEen_maxplane+fYOen_maxplane+fYEen_maxplane+fXOen_maxplane)/fConv_rxo; |
1063 |
|
if ( fSel ){ |
1064 |
|
if ( fXesel ) fColumn = fColXE; |
1065 |
|
if ( fXosel ){ |
1066 |
|
if ( fColXO != fColumn && fColumn > -1 ){ |
1067 |
|
printf(" ERROR! mismatch in column number between different sections! fColumn %i fColXO %i \n",fColumn,fColXO); |
1068 |
|
} else { |
1069 |
|
fColumn = fColXO; |
1070 |
|
}; |
1071 |
|
}; |
1072 |
|
if ( fYesel ){ |
1073 |
|
if ( fColYE != fColumn && fColumn > -1 ){ |
1074 |
|
printf(" ERROR! mismatch in column number between different sections! fColumn %i fColYE %i \n",fColumn,fColYE); |
1075 |
|
} else { |
1076 |
|
fColumn = fColYE; |
1077 |
|
}; |
1078 |
|
}; |
1079 |
|
if ( fYosel ){ |
1080 |
|
if ( fColYO != fColumn && fColumn > -1 ){ |
1081 |
|
printf(" ERROR! mismatch in column number between different sections! fColumn %i fColYO %i \n",fColumn,fColYO); |
1082 |
|
} else { |
1083 |
|
fColumn = fColYO; |
1084 |
|
}; |
1085 |
|
}; |
1086 |
|
}; |
1087 |
}; |
}; |
1088 |
// |
// |
1089 |
if ( debug ) printf("sel %i indep %i fMax_plane %f conv_r %f en_maxplane %f encalo %f \n",fSel,indep,fMax_plane,fConv_rxo,fXOen_maxplane,fEnergy); |
if ( debug ) printf("sel %i indep %i fMax_plane %f conv_r %f en_maxplane %f encalo %f \n",fSel,indep,fMax_plane,fConv_rxo,fXOen_maxplane,fEnergy); |