/[PAMELA software]/quicklook/tracker/flight/macros/FTrkQLook_BASIC.cxx
ViewVC logotype

Contents of /quicklook/tracker/flight/macros/FTrkQLook_BASIC.cxx

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.11 - (show annotations) (download)
Thu Jul 20 11:09:33 2006 UTC (18 years, 4 months ago) by pam-fi
Branch: MAIN
CVS Tags: v1r07, v1r08, v1r09
Changes since 1.10: +15 -11 lines
v1r07

1 /**
2 * FTrkQLook_BASIC.cxx
3 *
4 * autor: D.Fedele
5 * version v1r06
6 * Parameters:
7 * file - the data file to analyze
8 * fromevent - first event to analyze
9 * toevent - last event to analyze
10 * outdir - total path of output file
11 * outfile - extension of output file (pdf,ps,gif,jpg)
12 *
13 */
14 //
15 #include <TPaveText.h>
16 #include <TLatex.h>
17 #include <TCanvas.h>
18 #include <TGraph.h>
19 #include <TStyle.h>
20 #include <TFile.h>
21 #include <TTree.h>
22 #include <TArrow.h>
23 //
24 #include <physics/tracker/TrackerEvent.h>
25 #include <PscuHeader.h>
26 #include <EventHeader.h>
27 #include <RunHeaderEvent.h>
28 #include <EventCounter.h>
29 #include <PacketType.h>
30 //
31 #define MAXSTORAGE 50000
32
33 void FTrkQLook_BASIC(TString file,Int_t fromevent,Int_t toevent, TString outdir,TString outfile)
34 {
35 //
36 // obtain information about the data file and select the output dir
37 Int_t dwpos = file.Last('/');
38 Int_t dwpos1 = file.Last('.');
39 TString base,ffile ;
40 ffile=file(dwpos+1,dwpos1-(dwpos+1));
41 if(dwpos>0) base=file(0,dwpos);
42
43 TString out;
44 if(outdir.Length()==0){
45 out = base;
46 }else{
47 out = outdir;
48 }
49 if(out.Last('/')+1<out.Length()) out+="/";
50
51 //
52 // inizialise the variables and open the file
53 pamela::tracker::TrackerEvent *te=0;
54 pamela::EventHeader *eh=0,*eH=0,*ceh=0;
55 pamela::RunHeaderEvent *reh=0;
56 pamela::PscuHeader *ph=0,*pH=0;
57 pamela::EventCounter *cod=0;
58
59 pamela::PacketType *pctp=0;
60
61 TFile *datafile = new TFile(file);
62 TTree *otr = (TTree*)datafile->Get("RunHeader");
63 otr->SetBranchAddress("Header",&eH);
64 otr->SetBranchAddress("RunHeader",&reh);
65 TTree *tr = (TTree*)datafile->Get("Physics");
66 tr->SetBranchAddress("Tracker",&te);
67 tr->SetBranchAddress("Header",&eh);
68 TTree *ctr = (TTree*)datafile->Get("CalibTrk1");
69 ctr->SetBranchAddress("Header",&ceh);
70
71 Long64_t neventC = ctr->GetEntries();
72 Long64_t nevent = tr->GetEntries();
73 Long64_t neventH = otr->GetEntries();
74 Int_t minevent=0;
75 Int_t maxevent=0;
76
77 printf("Number of total events: %lld\nNumber of total header events: %lld\n",nevent,neventH);
78 printf("Number of calibration events: %lld\n",neventC);
79
80 if (nevent<=0){
81 datafile->Close();
82 return;
83 }
84 if ( fromevent > toevent && toevent>0 ){
85 printf("It must be fromevent < toevent \n");
86 return;
87 }
88 if ( fromevent > nevent || fromevent < 0 ) {
89 printf("You can choose fromevent between 0 (all) and %lld \n",nevent);
90 return;
91 }
92 if ( toevent > nevent || toevent < 0 ) {
93 printf("You can choose toevent between 0 (all) and %lld \n",nevent);
94 return;
95 }
96 if ( fromevent == 0 ) {
97 minevent = 0;
98 maxevent = nevent;
99 } else {
100 minevent = fromevent;
101 if ( toevent > 0 ){
102 maxevent = toevent+1;
103 } else if (toevent > nevent) {
104 maxevent = nevent;
105 } else {
106 maxevent = toevent+1;
107 }
108 nevent=maxevent-minevent ;
109 }
110
111 //
112 // other variables definitions
113 stringstream oss,fromfile,isfile;
114 //
115 // information about the RunHeader
116 ULong64_t HOBT[neventH];
117 Int_t trk_cal_us[neventH];
118 for (Int_t vi=0; vi<neventH;vi++){
119 HOBT[vi]=0;
120 trk_cal_us[vi]=0;
121 }
122
123 //
124 // information about RunHeader
125 Int_t countnboot=1;
126 for (Int_t ev=0; ev<neventH; ev++){
127 otr->GetEntry(ev);
128 pH = eH->GetPscuHeader();
129 HOBT[ev]= pH->GetOrbitalTime();
130 trk_cal_us[ev]=reh->TRK_CALIB_USED;
131 if((HOBT[ev]<HOBT[ev-1]) && ev>0)
132 countnboot+=1;
133 }
134 countnboot+=2*(Int_t)nevent/MAXSTORAGE;
135 // printf("\ncountnboot=%d\n",countnboot);
136
137 //
138 // information about calibration OBT
139 ULong64_t COBT[neventC];
140 for (Int_t vi=0; vi<neventC;vi++){
141 COBT[vi]=0;
142 }
143 for (Int_t ev=0; ev<neventC; ev++){
144 ctr->GetEntry(ev);
145 pH = ceh->GetPscuHeader();
146 COBT[ev]= pH->GetOrbitalTime();
147 }
148
149 //
150 // Style options
151 gStyle->SetLabelSize(0.06,"x");
152 gStyle->SetLabelSize(0.06,"y");
153 gStyle->SetStatFontSize(0.075);
154 gStyle->SetOptStat(10);
155 gStyle->SetFillColor(10);
156 gStyle->SetTitleFontSize(0.1);
157 gStyle->SetTitleFillColor(10);
158 gStyle->SetTitleOffset(0.8,"y");
159 gStyle->SetTitleOffset(0.9,"x");
160 gStyle->SetTitleSize(0.06,"y");
161 gStyle->SetTitleSize(0.055,"x");
162
163 //***************************************************************************************
164 // LOOP on each event
165 //***************************************************************************************
166
167 if (fromevent!=0)
168 printf("\n Scan of events from %i to %i ... \n",minevent,maxevent-1);
169 else
170 printf("\n Scan of events from %i to %i ... \n",minevent+1,maxevent);
171
172
173 Int_t minev=minevent,maxev=maxevent,hin=0,hfin=0,cin=0,cfin=0;
174 TPad *pad[12][countnboot];
175 TGraph *dataletime[12][countnboot],*dataletime1[12][countnboot];
176 TCanvas *DataTimeCanv[countnboot];
177 for(Int_t ii=0; ii<countnboot;ii++){
178 fromfile<<"FTrkQLook_BASIC File: "<<ffile;
179 isfile<<"DATALENGTH vs. OBT pag"<<ii+1;
180 DataTimeCanv[ii]=new TCanvas(isfile.str().c_str(),isfile.str().c_str(),900,1200);
181 DataTimeCanv[ii]->SetFillColor(10);
182 DataTimeCanv[ii]->Range(0,0,100,100);
183
184 TLatex *t=new TLatex();
185 t->SetTextFont(32);
186 t->SetTextColor(1);
187 t->SetTextAlign(12);
188 t->SetTextSize(0.02);
189 t->DrawLatex(2.,98.7,fromfile.str().c_str());
190 TLatex *t1=new TLatex();
191 t1->SetTextFont(32);
192 t1->SetTextColor(1);
193 t1->SetTextAlign(12);
194 t1->SetTextSize(0.02);
195 t1->DrawLatex(70.,98.7,isfile.str().c_str());
196 fromfile.str("");
197 isfile.str("");
198
199 fromfile<<"D = Default Calibration";
200 isfile<<"O = OnLine Calibration";
201 t->SetTextColor(6);
202 t->SetTextSize(0.018);
203 t->DrawLatex(70.,97.,fromfile.str().c_str());
204 t->SetTextColor(3);
205 t->DrawLatex(70.,96.,isfile.str().c_str());
206 fromfile.str("");
207 isfile.str("");
208
209 fromfile<<"The green arrow (if present) points out the time of the online calibration";
210 t->DrawLatex(7.,96.,fromfile.str().c_str());
211 fromfile.str("");
212
213 //*************************************************************************************
214 //book pads and histos
215 //***************************************************************************************
216
217
218 Float_t posy = 0.95; // up y-coord - top pads
219 Float_t hpad = 0; // pad height
220 Float_t posx1=0; // left x-coord - pad column
221 Float_t posx0=0; // x-coord - column division
222 Float_t wrel = 0; // relative x size of first sub-column
223 Float_t marg = 0.004; // margin among pads
224
225 hpad = (posy-marg*11)/6;
226 wrel = (1-marg*4)/2;
227 stringstream title;
228 stringstream hid;
229
230 for(Int_t n = 0; n<12; n++) {
231 if ( (n+1)%2==1 ) {
232 if(n>1) posy = posy-(marg*2+hpad);
233 posx1 = marg;
234 posx0 = posx1 + wrel;
235 }
236 else {
237 posx1 = posx0 + 2*marg;
238 posx0 = posx1 + wrel;
239 }
240
241 /* -----------> pad for histograms */
242 oss<<"pad"<<n*100+ii;
243 pad[n][ii]=new TPad(oss.str().c_str()," ",posx1,posy-hpad,posx0,posy,18,0,0);
244 oss.str("");
245 };
246
247 TLine li;
248 li.SetLineColor(1);
249 li.SetLineStyle(1);
250 li.SetLineWidth(1);
251
252 TArrow ar;
253 ar.SetLineColor(3);
254 stringstream calus;
255
256 TLatex *t2=new TLatex();
257 t2->SetTextFont(32);
258 t2->SetTextColor(1);
259 t2->SetTextAlign(13);
260 t2->SetTextSize(0.08);
261
262 Int_t i=0;
263 Long64_t x[MAXSTORAGE], xb[MAXSTORAGE];
264 Float_t yyd[MAXSTORAGE][12],yyb[MAXSTORAGE][12];
265 Int_t countbad[12];
266 Float_t perc=0,xMIN=0.,xMAX=0.;
267 for (Int_t n=0; n<12 ; n++)
268 countbad[n]=0;
269
270 //
271 // obtain values of the datalenght
272 for (Int_t ev=minev; ev<maxevent; ev++){
273 tr->GetEntry(ev);
274 ph = eh->GetPscuHeader();
275 cod = eh->GetCounter();
276
277 if(ev==minev){
278 if(cod->Get(pctp->CalibTrk1)>0) cin=cod->Get(pctp->CalibTrk1)-1;
279 else cin=cod->Get(pctp->CalibTrk1);
280 if(cin==cfin-1) cin+=1;
281
282 if(cod->Get(pctp->RunHeader)>0) hin=cod->Get(pctp->RunHeader)-1;
283 else hin=cod->Get(pctp->RunHeader);
284 if(hin==hfin-1) hin+=1;
285 }
286
287 if(ev==maxevent-1) maxev=maxevent-1;
288
289 if((ph->GetOrbitalTime()<x[ev-minev-1] && ev-minev!=0) || ev-minev==MAXSTORAGE){
290 maxev=ev;
291 break;
292 }
293 else{
294 cfin=cod->Get(pctp->CalibTrk1);
295 hfin=cod->Get(pctp->RunHeader);
296 x[(ev-minev)]= ph->GetOrbitalTime();
297 i=0;
298
299 for (Int_t n=0; n<12 ; n++){
300 yyb[countbad[n]][i]=0;
301 xb[countbad[n]]= 0;
302
303 i=te->DSPnumber[n]-1;
304
305 yyd[(ev-minev)][i]=te->DATAlength[n];
306 if(i==6){
307 if(yyd[(ev-minev)][i]>1500){
308 if(yyd[(ev-minev)][i]<3075){
309 yyb[countbad[i]][i]= yyd[(ev-minev)][i];
310 xb[countbad[i]]= x[(ev-minev)];
311 countbad[i]+=1;
312 }
313 }
314 }
315 else{
316 if(yyd[(ev-minev)][i]>750){
317 if(yyd[(ev-minev)][i]<3075){
318 yyb[countbad[i]][i]= yyd[(ev-minev)][i];
319 xb[countbad[i]]= x[(ev-minev)];
320 countbad[i]+=1;
321 }
322 }
323 }
324 }
325 }
326 }
327
328 //
329 // define limit for the Xaxis of the graphs
330
331 xMAX=x[maxev-minev-1]+(x[maxev-minev-1]-x[0])/10;
332 if(xMAX>1000000) xMIN=x[0]-(x[maxev-minev-1]-x[0])/10;
333 if(xMAX<1000000 || xMIN<0) xMIN=0.;
334
335
336 //
337 // Draw Histos
338 for (Int_t i=0; i<12 ; i++){
339
340 Double_t yd[maxev-minev],ybd[maxev-minev];
341 Double_t xd[maxev-minev],xbd[maxev-minev];
342
343 for(Int_t v=0;v<maxev-minev;v++){
344 yd[v]=(Double_t)yyd[v][i];
345 ybd[v]=(Double_t)yyb[v][i];
346 xd[v]=(Double_t)x[v];
347 xbd[v]=(Double_t)xb[v];
348 }
349
350 if(xMAX<1000000){
351 xd[maxev-minev-1]=0;
352 yd[maxev-minev-1]=0.;
353 }
354
355 if((maxev-minev)>1000){
356 perc=(countbad[i]*100)/(maxev-minev);
357 if(perc>10) pad[i][ii]->SetFillColor(2);
358 else pad[i][ii]->SetFillColor(10);
359 }
360 else{
361 if(countbad[i]>=100) pad[i][ii]->SetFillColor(2);
362 else pad[i][ii]->SetFillColor(10);
363 }
364
365 oss<<"DSP "<<i+1;
366 DataTimeCanv[ii]->cd();
367 pad[i][ii]->SetFrameFillColor(10);
368 pad[i][ii]->Draw();
369 pad[i][ii]->cd();
370 dataletime[i][ii]= new TGraph((maxev-minev),xd,yd);
371 dataletime[i][ii]->SetTitle(oss.str().c_str());
372 dataletime[i][ii]->GetXaxis()->SetTitle("OBT (ms)");
373 dataletime[i][ii]->GetXaxis()->CenterTitle();
374 dataletime[i][ii]->GetXaxis()->SetRangeUser(xMIN,xMAX);
375 dataletime[i][ii]->GetYaxis()->SetTitle("datalength (Word 13 bit)");
376 dataletime[i][ii]->GetYaxis()->CenterTitle();
377 if(i==6) dataletime[i][ii]->GetYaxis()->SetRangeUser(0,4500);
378 else dataletime[i][ii]->GetYaxis()->SetRangeUser(0,4000);
379 dataletime[i][ii]->SetMarkerStyle(21);
380 if((maxev-minev)<50) dataletime[i][ii]->SetMarkerSize(0.5);
381 else dataletime[i][ii]->SetMarkerSize(0.1);
382 dataletime[i][ii]->SetMarkerColor(4);
383 dataletime[i][ii]->Draw("ap");
384
385
386 if((maxev-minev)>1000 && perc>10){
387 dataletime1[i][ii]= new TGraph(countbad[i],xbd,ybd);
388 dataletime1[i][ii]->SetMarkerStyle(21);
389 if((maxev-minev)<50) dataletime1[i][ii]->SetMarkerSize(0.5);
390 else dataletime1[i][ii]->SetMarkerSize(0.1);
391 dataletime1[i][ii]->SetMarkerColor(2);
392 dataletime1[i][ii]->Draw("psame");
393 }
394 else if((maxev-minev)<1000 && countbad[i]>=100){
395 dataletime1[i][ii]= new TGraph(countbad[i],xbd,ybd);
396 dataletime1[i][ii]->SetMarkerStyle(21);
397 if((maxev-minev)<50) dataletime1[i][ii]->SetMarkerSize(0.5);
398 else dataletime1[i][ii]->SetMarkerSize(0.1);
399 dataletime1[i][ii]->SetMarkerColor(2);
400 dataletime1[i][ii]->Draw("psame");
401 }
402 li.SetLineColor(1);
403 li.SetLineStyle(1);
404 li.SetLineWidth(1);
405 if(i!=6) li.DrawLine(xMIN,750,xMAX,750);
406 else li.DrawLine(xMIN,1500,xMAX,1500);
407 li.DrawLine(xMIN,3075,xMAX,3075);
408
409 li.SetLineColor(12);
410 li.SetLineStyle(4);
411 li.SetLineWidth(1);
412 for(Int_t j=hin;j<hfin;j++){
413 if(i==6) li.DrawLine(HOBT[j],0.,HOBT[j],4500.);
414 else li.DrawLine(HOBT[j],0.,HOBT[j],4000.);
415 if(trk_cal_us[j]==104){
416 calus<<"D";
417 t2->SetTextColor(6);
418 if(i==6) t2->DrawLatex(HOBT[j],4350.,calus.str().c_str());
419 else t2->DrawLatex(HOBT[j],3850.,calus.str().c_str());
420 calus.str("");
421 }
422 else{
423 calus<<"O";
424 t2->SetTextColor(3);
425 if(i==6) t2->DrawLatex(HOBT[j],4350.,calus.str().c_str());
426 else t2->DrawLatex(HOBT[j],3850.,calus.str().c_str());
427 calus.str("");
428 }
429 }
430 for(Int_t j=cin;j<cfin;j++){
431 if(i==6) ar.DrawArrow(COBT[j],1700.,COBT[j],2700.,0.01,"<");
432 else ar.DrawArrow(COBT[j],1000.,COBT[j],2000.,0.01,"<");
433 }
434
435 oss.str("");
436 DataTimeCanv[ii]->Update();
437 }
438
439 minev=maxev;
440 if(maxev==maxevent-1) {
441 countnboot=ii+1;
442 break;
443 }
444 }
445 printf("... end of packets. \n");
446
447 //*************************************************************************
448 // Save output Files
449 //*************************************************************************
450 stringstream nom1,nom2,nom3;
451
452 for(Int_t fl=0;fl<countnboot;fl++){
453 if(countnboot==1){
454 nom1<<ffile<<"_FTrkQLook_BASIC."<<outfile.Data();
455 DataTimeCanv[fl]->Print(out+nom1.str().c_str());
456 nom1.str("");
457 }
458
459 if(countnboot>=2){
460 if(!strcmp(outfile.Data(),"ps") || !strcmp(outfile.Data(),"pdf")){
461 nom1.str("");
462 nom2.str("");
463 nom3.str("");
464 nom1<<ffile<<"_FTrkQLook_BASIC.ps(";
465 nom2<<ffile<<"_FTrkQLook_BASIC.ps";
466 nom3<<ffile<<"_FTrkQLook_BASIC.ps)";
467 if(fl==0) DataTimeCanv[fl]->Print(out+nom1.str().c_str(),"portrait");
468 else if(fl==countnboot-1) DataTimeCanv[fl]->Print(out+nom3.str().c_str(),"portrait");
469 else DataTimeCanv[fl]->Print(out+nom2.str().c_str(),"portrait");
470
471 }
472 else{
473 nom1.str("");
474 nom1<<ffile<<"_FTrkQLook_BASIC-pag"<<fl+1<<"."<<outfile.Data();
475 DataTimeCanv[fl]->Print(out+nom1.str().c_str());
476 }
477 }
478 }
479
480 //
481 // Convert ps to pdf if required
482 if(!strcmp(outfile.Data(),"pdf") && countnboot>=2){
483 stringstream com;
484 com<<"ps2pdf13 "<<out<<ffile<<"_FTrkQLook_BASIC.ps "<<out<<ffile<<"_FTrkQLook_BASIC.pdf";
485 system(com.str().c_str());
486 printf("\n---> ps file converted in pdf format!\n");
487 com.str("");
488 com<<"rm -f "<<out<<ffile<<"_FTrkQLook_BASIC.ps";
489 system(com.str().c_str());
490 printf("---> ps file removed!\n\n");
491 com.str("");
492 }
493
494 datafile->Close();
495 gROOT->Reset();
496 return;
497 }

  ViewVC Help
Powered by ViewVC 1.1.23