4646#include " TClonesArray.h"
4747#include " TMath.h"
4848
49- double fpos_range_min_FrsSci[2 ] = { -4 , 40 };
50- double fpos_range_max_FrsSci[2 ] = { 4 , 50 };
51- double ftof_range_min_FrsSci[2 ] = {1375 , 1240 };
52- double ftof_range_max_FrsSci[2 ] = {1400 , 1275 };
53-
49+ double fpos_range_min_FrsSci[2 ] = { -4 , 40 };
50+ double fpos_range_max_FrsSci[2 ] = { 4 , 50 };
51+ double ftof_range_min_FrsSci[2 ] = { 1375 , 1240 };
52+ double ftof_range_max_FrsSci[2 ] = { 1400 , 1275 };
5453
5554R3BAnalysisIncomingID::R3BAnalysisIncomingID ()
5655 : R3BAnalysisIncomingID(" AnalysisIncomingID" , 1 )
@@ -283,14 +282,17 @@ void R3BAnalysisIncomingID::Exec(Option_t*)
283282 0.5 * (FrsSciTraw[i * fFrsSciNbPmts][hitr] + FrsSciTraw[i * fFrsSciNbPmts + 1][hitl]) -
284283 (FrsSciTraw[i * fFrsSciNbPmts + 2][0] + 8192. * 5.);
285284 }*/
286-
287- StartTraw_atTcal = fTimeStitch ->GetTime ( 0.5 * (FrsSciTraw[i * fFrsSciNbPmts ][hitr] + FrsSciTraw[i * fFrsSciNbPmts + 1 ][hitl]) -
288- FrsSciTraw[i * fFrsSciNbPmts + 2 ][0 ], " vftx" , " vftx" );
289-
290-
285+
286+ StartTraw_atTcal = fTimeStitch ->GetTime (
287+ 0.5 * (FrsSciTraw[i * fFrsSciNbPmts ][hitr] + FrsSciTraw[i * fFrsSciNbPmts + 1 ][hitl]) -
288+ FrsSciTraw[i * fFrsSciNbPmts + 2 ][0 ],
289+ " vftx" ,
290+ " vftx" );
291+
291292 auto TofRaw = fTimeStitch ->GetTime (fHeader ->GetTStartMaster () - StartTraw_atTcal, " vftx" , " vftx" );
292293 auto PosRaw = FrsSciTraw[i * fFrsSciNbPmts ][hitr] - FrsSciTraw[i * fFrsSciNbPmts + 1 ][hitl];
293- if (ftof_range_min_FrsSci[i] < TofRaw && TofRaw < ftof_range_max_FrsSci[i] && fpos_range_min_FrsSci[i] < PosRaw && PosRaw < fpos_range_max_FrsSci[i])
294+ if (ftof_range_min_FrsSci[i] < TofRaw && TofRaw < ftof_range_max_FrsSci[i] &&
295+ fpos_range_min_FrsSci[i] < PosRaw && PosRaw < fpos_range_max_FrsSci[i])
294296 {
295297 indexr[i] = hitr;
296298 indexl[i] = hitl;
@@ -299,28 +301,29 @@ void R3BAnalysisIncomingID::Exec(Option_t*)
299301 }
300302 }
301303
302- if (indexr[i] >= 0 && indexl[i] >= 0 && multTofRaw[i] == 1 && i== 1 )
304+ if (indexr[i] >= 0 && indexl[i] >= 0 && multTofRaw[i] == 1 && i == 1 )
303305 {
304- /* if ((FrsSciTC[i * fFrsSciNbPmts][indexr[i]] < FrsSciTC[i * fFrsSciNbPmts + 2][0]) &&
305- (FrsSciTC[i * fFrsSciNbPmts + 1][indexl[i]] < FrsSciTC[i * fFrsSciNbPmts + 2][0]))
306- {
307- StartTraw_atTcal = 0.5 * (FrsSciTraw[i * fFrsSciNbPmts][indexr[i]] +
308- FrsSciTraw[i * fFrsSciNbPmts + 1][indexl[i]]) -
309- FrsSciTraw[i * fFrsSciNbPmts + 2][0] - (i == 1 ? 130. : 0.);
310- }
311- else
312- {
313- StartTraw_atTcal = 0.5 * (FrsSciTraw[i * fFrsSciNbPmts][indexr[i]] +
314- FrsSciTraw[i * fFrsSciNbPmts + 1][indexl[i]]) -
315- (FrsSciTraw[i * fFrsSciNbPmts + 2][0] + 8192 * 5.) - (i == 1 ? 130. : 0.);
316- }*/
317-
318-
319- StartTraw_atTcal = fTimeStitch ->GetTime ( 0.5 * (FrsSciTraw[i * fFrsSciNbPmts ][indexr[i]] +
320- FrsSciTraw[i * fFrsSciNbPmts + 1 ][indexl[i]]) -
321- FrsSciTraw[i * fFrsSciNbPmts + 2 ][0 ] , " vftx" , " vftx" ) - (i == 1 ? 130 . : 0 .);
322-
323-
306+ /* if ((FrsSciTC[i * fFrsSciNbPmts][indexr[i]] < FrsSciTC[i * fFrsSciNbPmts + 2][0]) &&
307+ (FrsSciTC[i * fFrsSciNbPmts + 1][indexl[i]] < FrsSciTC[i * fFrsSciNbPmts + 2][0]))
308+ {
309+ StartTraw_atTcal = 0.5 * (FrsSciTraw[i * fFrsSciNbPmts][indexr[i]] +
310+ FrsSciTraw[i * fFrsSciNbPmts + 1][indexl[i]]) -
311+ FrsSciTraw[i * fFrsSciNbPmts + 2][0] - (i == 1 ? 130. : 0.);
312+ }
313+ else
314+ {
315+ StartTraw_atTcal = 0.5 * (FrsSciTraw[i * fFrsSciNbPmts][indexr[i]] +
316+ FrsSciTraw[i * fFrsSciNbPmts + 1][indexl[i]]) -
317+ (FrsSciTraw[i * fFrsSciNbPmts + 2][0] + 8192 * 5.) - (i == 1 ? 130. : 0.);
318+ }*/
319+
320+ StartTraw_atTcal = fTimeStitch ->GetTime (0.5 * (FrsSciTraw[i * fFrsSciNbPmts ][indexr[i]] +
321+ FrsSciTraw[i * fFrsSciNbPmts + 1 ][indexl[i]]) -
322+ FrsSciTraw[i * fFrsSciNbPmts + 2 ][0 ],
323+ " vftx" ,
324+ " vftx" ) -
325+ (i == 1 ? 130 . : 0 .);
326+
324327 auto TofRaw = fTimeStitch ->GetTime (fHeader ->GetTStartMaster () - StartTraw_atTcal, " vftx" , " vftx" );
325328 // fh2_Tstop_vs_Tstart_Zgt5[i]->Fill(StartTraw_atTcal, fHeader->GetTStartMaster());
326329 auto PosRaw = FrsSciTraw[i * fFrsSciNbPmts ][indexr[i]] - FrsSciTraw[i * fFrsSciNbPmts + 1 ][indexl[i]] -
0 commit comments