| 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 |
// |
// |