2 |
* FTrkQLook_BASIC.cxx |
* FTrkQLook_BASIC.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r10 |
* 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 |
20 |
#include <TFile.h> |
#include <TFile.h> |
21 |
#include <TTree.h> |
#include <TTree.h> |
22 |
#include <TArrow.h> |
#include <TArrow.h> |
23 |
|
#include <TAxis.h> |
24 |
|
#include <TROOT.h> |
25 |
|
#include <cstdlib> |
26 |
|
#include <sys/stat.h> |
27 |
// |
// |
28 |
#include <physics/tracker/TrackerEvent.h> |
#include <physics/tracker/TrackerEvent.h> |
29 |
#include <PscuHeader.h> |
#include <PscuHeader.h> |
117 |
stringstream oss,fromfile,isfile; |
stringstream oss,fromfile,isfile; |
118 |
// |
// |
119 |
// information about the RunHeader |
// information about the RunHeader |
120 |
ULong64_t HOBT[neventH]; |
ULong_t HOBT[neventH]; |
121 |
Int_t trk_cal_us[neventH]; |
Int_t trk_cal_us[neventH]; |
122 |
for (Int_t vi=0; vi<neventH;vi++){ |
for (Int_t vi=0; vi<neventH;vi++){ |
123 |
HOBT[vi]=0; |
HOBT[vi]=0; |
140 |
|
|
141 |
// |
// |
142 |
// information about calibration OBT |
// information about calibration OBT |
143 |
ULong64_t COBT[neventC]; |
ULong_t COBT[neventC]; |
144 |
for (Int_t vi=0; vi<neventC;vi++){ |
for (Int_t vi=0; vi<neventC;vi++){ |
145 |
COBT[vi]=0; |
COBT[vi]=0; |
146 |
} |
} |
264 |
t2->SetTextSize(0.08); |
t2->SetTextSize(0.08); |
265 |
|
|
266 |
Int_t i=0; |
Int_t i=0; |
267 |
Long64_t x[MAXSTORAGE], xb[MAXSTORAGE]; |
ULong_t x[MAXSTORAGE], xb[MAXSTORAGE]; |
268 |
Float_t yyd[MAXSTORAGE][12],yyb[MAXSTORAGE][12]; |
Float_t yyd[MAXSTORAGE][12],yyb[MAXSTORAGE][12]; |
269 |
Int_t countbad[12]; |
Int_t countbad[12]; |
270 |
Float_t perc=0,xMIN=0.,xMAX=0.; |
Float_t perc=0; |
271 |
|
Double_t xMIN=0.,xMAX=0.; |
272 |
for (Int_t n=0; n<12 ; n++) |
for (Int_t n=0; n<12 ; n++) |
273 |
countbad[n]=0; |
countbad[n]=0; |
274 |
|
|
291 |
|
|
292 |
if(ev==maxevent-1) maxev=maxevent-1; |
if(ev==maxevent-1) maxev=maxevent-1; |
293 |
|
|
294 |
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){ |
295 |
maxev=ev; |
maxev=ev; |
296 |
break; |
break; |
297 |
} |
} |
338 |
|
|
339 |
xMAX=x[maxev-minev-1]+(x[maxev-minev-1]-x[0])/10; |
xMAX=x[maxev-minev-1]+(x[maxev-minev-1]-x[0])/10; |
340 |
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; |
341 |
if(xMAX<1000000 || xMIN<0) xMIN=0.; |
if(xMAX<1000000 || xMIN<0) xMIN=0.; |
342 |
|
if(xMIN==xMAX) xMIN=x[0] * x[0]/xMAX; |
343 |
|
|
344 |
|
|
345 |
// |
// |