2 |
* FTrkQLook_BASIC.cxx |
* FTrkQLook_BASIC.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r08 |
* version v1r12 |
6 |
* Parameters: |
* Parameters: |
7 |
* file - the data file to analyze |
* file - the data file to analyze |
8 |
* fromevent - first event to analyze |
* fromevent - first event to analyze |
113 |
stringstream oss,fromfile,isfile; |
stringstream oss,fromfile,isfile; |
114 |
// |
// |
115 |
// information about the RunHeader |
// information about the RunHeader |
116 |
ULong64_t HOBT[neventH]; |
ULong_t HOBT[neventH]; |
117 |
Int_t trk_cal_us[neventH]; |
Int_t trk_cal_us[neventH]; |
118 |
for (Int_t vi=0; vi<neventH;vi++){ |
for (Int_t vi=0; vi<neventH;vi++){ |
119 |
HOBT[vi]=0; |
HOBT[vi]=0; |
136 |
|
|
137 |
// |
// |
138 |
// information about calibration OBT |
// information about calibration OBT |
139 |
ULong64_t COBT[neventC]; |
ULong_t COBT[neventC]; |
140 |
for (Int_t vi=0; vi<neventC;vi++){ |
for (Int_t vi=0; vi<neventC;vi++){ |
141 |
COBT[vi]=0; |
COBT[vi]=0; |
142 |
} |
} |
260 |
t2->SetTextSize(0.08); |
t2->SetTextSize(0.08); |
261 |
|
|
262 |
Int_t i=0; |
Int_t i=0; |
263 |
Long64_t x[MAXSTORAGE], xb[MAXSTORAGE]; |
ULong_t x[MAXSTORAGE], xb[MAXSTORAGE]; |
264 |
Float_t yyd[MAXSTORAGE][12],yyb[MAXSTORAGE][12]; |
Float_t yyd[MAXSTORAGE][12],yyb[MAXSTORAGE][12]; |
265 |
Int_t countbad[12]; |
Int_t countbad[12]; |
266 |
Float_t perc=0,xMIN=0.,xMAX=0.; |
Float_t perc=0; |
267 |
|
Double_t xMIN=0.,xMAX=0.; |
268 |
for (Int_t n=0; n<12 ; n++) |
for (Int_t n=0; n<12 ; n++) |
269 |
countbad[n]=0; |
countbad[n]=0; |
270 |
|
|
287 |
|
|
288 |
if(ev==maxevent-1) maxev=maxevent-1; |
if(ev==maxevent-1) maxev=maxevent-1; |
289 |
|
|
290 |
if((ph->GetOrbitalTime()<x[ev-minev-1] && ev-minev!=0) || ev-minev==MAXSTORAGE){ |
if(((ULong_t)ph->GetOrbitalTime()<x[ev-minev-1] && ev-minev!=0) || ev-minev==MAXSTORAGE){ |
291 |
maxev=ev; |
maxev=ev; |
292 |
break; |
break; |
293 |
} |
} |
334 |
|
|
335 |
xMAX=x[maxev-minev-1]+(x[maxev-minev-1]-x[0])/10; |
xMAX=x[maxev-minev-1]+(x[maxev-minev-1]-x[0])/10; |
336 |
if(xMAX>1000000) xMIN=x[0]-(x[maxev-minev-1]-x[0])/10; |
if(xMAX>1000000) xMIN=x[0]-(x[maxev-minev-1]-x[0])/10; |
337 |
if(xMAX<1000000 || xMIN<0) xMIN=0.; |
if(xMAX<1000000 || xMIN<0) xMIN=0.; |
338 |
|
if(xMIN==xMAX) xMIN=x[0] * x[0]/xMAX; |
339 |
|
|
340 |
|
|
341 |
// |
// |
368 |
else pad[n][ii]->SetFillColor(10); |
else pad[n][ii]->SetFillColor(10); |
369 |
} |
} |
370 |
|
|
371 |
oss<<"DSP "<<i+1; |
oss<<"DSP "<<n+1; |
372 |
DataTimeCanv[ii]->cd(); |
DataTimeCanv[ii]->cd(); |
373 |
pad[n][ii]->SetFrameFillColor(10); |
pad[n][ii]->SetFrameFillColor(10); |
374 |
pad[n][ii]->Draw(); |
pad[n][ii]->Draw(); |
375 |
pad[n][ii]->cd(); |
pad[n][ii]->cd(); |
376 |
dataletime[n][ii]= new TGraph((maxev-minev),xd,yd); |
dataletime[n][ii]= new TGraph((maxev-minev),xd,yd); |
377 |
|
if(maxev-minev==1) dataletime[n][ii]->GetXaxis()->SetNdivisions(-501); |
378 |
dataletime[n][ii]->SetTitle(oss.str().c_str()); |
dataletime[n][ii]->SetTitle(oss.str().c_str()); |
379 |
dataletime[n][ii]->GetXaxis()->SetTitle("OBT (ms)"); |
dataletime[n][ii]->GetXaxis()->SetTitle("OBT (ms)"); |
380 |
dataletime[n][ii]->GetXaxis()->CenterTitle(); |
dataletime[n][ii]->GetXaxis()->CenterTitle(); |