2 |
* FTrkScanQlook_EXPERT.cxx |
* FTrkScanQlook_EXPERT.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v2r00 |
* version v2r08 |
6 |
* Parameters: |
* Parameters: |
7 |
* file - the path to the root file to analyze |
* file - the path to the root file to analyze |
8 |
* outdir - total path of output file |
* outdir - total path of output file |
9 |
* event - the number of the single event to analyze |
* event - the number of the single event to analyze |
10 |
* va1 - the number of the single va1 to analyze (dsp*100+va1) |
* va1 - the number of the single va1 to analyze (dsp*100+va1) |
11 |
* value - the number of events to store in the .dat file with va1 values |
* value - the number of events to store in the .dat file with va1 values*10 + # of calibration (if #of cal = 0, select every calibration) |
12 |
* outfile - extension of output file (pdf,ps,gif,jpg) |
* outfile - extension of output file (pdf,ps,gif,jpg) |
13 |
* |
* |
14 |
*/ |
*/ |
16 |
#include <iostream> |
#include <iostream> |
17 |
#include <sstream> |
#include <sstream> |
18 |
#include <fstream> |
#include <fstream> |
19 |
|
#include <TAxis.h> |
20 |
|
#include <TH1F.h> |
21 |
|
#include <TROOT.h> |
22 |
|
#include <cstdlib> |
23 |
// |
// |
24 |
#include <TPaveText.h> |
#include <TPaveText.h> |
25 |
#include <TLatex.h> |
#include <TLatex.h> |
33 |
#include <physics/tracker/TrackerEvent.h> |
#include <physics/tracker/TrackerEvent.h> |
34 |
#include <PscuHeader.h> |
#include <PscuHeader.h> |
35 |
#include <EventHeader.h> |
#include <EventHeader.h> |
36 |
|
#include <CalibTrk1Event.h> |
37 |
|
#include <CalibTrk2Event.h> |
38 |
#include <RunHeaderEvent.h> |
#include <RunHeaderEvent.h> |
39 |
// |
// |
40 |
|
|
41 |
using namespace std; |
using namespace std; |
42 |
|
|
43 |
typedef struct trkword{ |
struct trkword{ |
44 |
int type; |
int type; |
45 |
int decode; |
int decode; |
46 |
}; |
}; |
149 |
pamela::EventHeader *eh=0,*eH=0; |
pamela::EventHeader *eh=0,*eH=0; |
150 |
pamela::PscuHeader *ph=0,*pH=0; |
pamela::PscuHeader *ph=0,*pH=0; |
151 |
pamela::RunHeaderEvent *reh=0; |
pamela::RunHeaderEvent *reh=0; |
152 |
|
pamela::CalibTrk1Event *trk1 = 0; |
153 |
|
pamela::CalibTrk2Event *trk2 = 0; |
154 |
pamela::EventCounter *cod=0; |
pamela::EventCounter *cod=0; |
155 |
|
|
156 |
pamela::PacketType *pctp=0; |
pamela::PacketType *pctp=0; |
171 |
otr->SetBranchAddress("Header",&eH); |
otr->SetBranchAddress("Header",&eH); |
172 |
otr->SetBranchAddress("RunHeader",&reh); |
otr->SetBranchAddress("RunHeader",&reh); |
173 |
|
|
174 |
|
TTree *otr1 = (TTree*)trackerFile->Get("CalibTrk1"); |
175 |
|
otr1->SetBranchAddress("CalibTrk1", &trk1); |
176 |
|
TTree *otr2 = (TTree*)trackerFile->Get("CalibTrk2"); |
177 |
|
otr2->SetBranchAddress("CalibTrk2", &trk2); |
178 |
|
|
179 |
// Define variables |
// Define variables |
180 |
Int_t nevents = tr->GetEntries(); |
Int_t nevents = tr->GetEntries(); |
181 |
Int_t neventH = otr->GetEntries(); |
Int_t neventH = otr->GetEntries(); |
182 |
|
Int_t cevents=0; |
183 |
|
|
184 |
|
if(otr1->GetEntries()==otr2->GetEntries()) |
185 |
|
cevents = otr1->GetEntries(); |
186 |
|
else printf("\nWARNING: CalibTrk1 entries is different from CalibTrk2 entries\n\n"); |
187 |
|
|
188 |
|
|
189 |
if ( nevents <= 0 ) { |
if ( nevents <= 0 ) { |
190 |
trackerFile->Close(); |
trackerFile->Close(); |
191 |
printf("The file is empty, exiting...\n"); |
printf("The file is empty, exiting...\n"); |
204 |
// |
// |
205 |
printf("\n Number of Entries: %d\n",nevents); |
printf("\n Number of Entries: %d\n",nevents); |
206 |
printf(" Number of Header Entries: %d\n",neventH); |
printf(" Number of Header Entries: %d\n",neventH); |
207 |
|
printf(" Number of Calibration: %d\n",cevents); |
208 |
|
|
209 |
Long64_t obt=0; |
ULong_t obt=0,hobt[neventH]; |
210 |
Int_t eve,cin=0; |
Int_t eve,cin=0,cins=0,br=0,olrun[neventH]; |
211 |
TString cal=""; |
TString cal=""; |
212 |
|
|
213 |
eve=3; |
eve=3; |
214 |
for(Int_t i=0;i<neventH;i++){ |
for(Int_t i=0;i<neventH;i++){ |
215 |
otr->GetEntry(i); |
otr->GetEntry(i); |
216 |
pH = eH->GetPscuHeader(); |
pH = eH->GetPscuHeader(); |
217 |
if(reh->TRK_CALIB_USED!=104){ |
hobt[i]= pH->GetOrbitalTime(); |
218 |
|
olrun[i]=reh->TRK_CALIB_USED; |
219 |
|
if(reh->TRK_CALIB_USED!=104 && br==0){ |
220 |
obt = pH->GetOrbitalTime(); |
obt = pH->GetOrbitalTime(); |
221 |
cal="Event with online calibration"; |
cal="Event with online calibration"; |
222 |
break; |
br=1; |
223 |
} |
} |
224 |
if(i==neventH-1){ |
if(i==neventH-1 && br==0){ |
225 |
cal="***** ONLINE CALIBRATION NOT FOUND IN THIS FILE *****"; |
cal="***** ONLINE CALIBRATION NOT FOUND IN THIS FILE *****"; |
226 |
eve=2; |
eve=2; |
227 |
} |
} |
228 |
} |
} |
|
if(eve==3){ |
|
|
for(Int_t i=0;i<nevents;i++){ |
|
|
tr->GetEntry(i); |
|
|
ph = eh->GetPscuHeader(); |
|
|
cod = eh->GetCounter(); |
|
|
if(i==0) cin=cod->Get(pctp->CalibTrk1); |
|
|
if(reh->TRK_CALIB_USED==104) continue; |
|
|
if(event<0 && cod->Get(pctp->CalibTrk1)==cin+1){ |
|
|
eve=i+3; |
|
|
break; |
|
|
} |
|
|
else if(event>=0 && ph->GetOrbitalTime()>obt){ |
|
|
eve=i+3; |
|
|
break; |
|
|
} |
|
|
} |
|
|
} |
|
229 |
|
|
230 |
int tot=2,totvalues=0,TOT=0; |
int tot=2,totvalues=0,TOT=0,totsum=0; |
231 |
if(event<0) tot=abs(event); |
if(event<0) tot=abs(event); |
232 |
if(value>0) totvalues=value; |
if(value>0) totvalues=value/10; |
233 |
|
|
234 |
TOT= tot > totvalues ? tot : totvalues; |
TOT= tot > totvalues ? tot : totvalues; |
235 |
|
if(value>0) totsum= (value/10) > 100 ? (value/10) : 100; |
236 |
|
else totsum=100; |
237 |
|
|
238 |
|
Int_t calev=cevents; |
239 |
|
|
240 |
TH1F *histomax[12][tot]; //histos of max signals |
TH1F *histomax[12][tot]; //histos of max signals |
241 |
TH1F *histocomp[12][tot]; //histos of compressed data |
TH1F *histocomp[12][tot]; //histos of compressed data |
242 |
TH1F *histofull[12][tot]; //histos of full data |
TH1F *histofull[12][tot]; //histos of full data |
243 |
TCanvas *c1[tot]; |
TCanvas *c1[tot],*csum[calev]; |
244 |
TH1F *histova[tot]; //histos of va1 signals |
TH1F *histova[tot]; |
245 |
|
TGraph*sum[12][calev][totsum/2]; |
246 |
TCanvas *cva[tot]; |
TCanvas *cva[tot]; |
247 |
|
|
248 |
|
TPad *trkpad[12]; //pad for histos |
249 |
|
|
250 |
|
stringstream fromfile,isfile,title,hid,message; |
251 |
|
TString figsa="",figsav="",figsava=""; |
252 |
|
Int_t canvasx=1200; |
253 |
|
Int_t canvasy=900; |
254 |
|
|
255 |
TLatex *t=new TLatex(); |
TLatex *t=new TLatex(); |
256 |
t->SetTextFont(32); |
t->SetTextFont(32); |
258 |
t->SetTextAlign(12); |
t->SetTextAlign(12); |
259 |
t->SetTextSize(0.02); |
t->SetTextSize(0.02); |
260 |
|
|
261 |
|
gStyle->SetLabelSize(0.06,"x"); |
262 |
|
gStyle->SetLabelSize(0.06,"y"); |
263 |
|
gStyle->SetTitleFontSize(0.1); |
264 |
|
gStyle->SetFillColor(10); |
265 |
|
gStyle->SetTitleFillColor(10); |
266 |
|
gStyle->SetOptStat(0); |
267 |
|
|
268 |
|
Int_t p[12],psum=0; |
269 |
|
Int_t whistostrip[3072]; |
270 |
|
for(int s=0;s<3072;s++) whistostrip[s]=s+1; |
271 |
|
|
272 |
|
if(eve==3){ |
273 |
|
for(Int_t i=0;i<nevents;i++){ |
274 |
|
tr->GetEntry(i); |
275 |
|
ph = eh->GetPscuHeader(); |
276 |
|
cod = eh->GetCounter(); |
277 |
|
if(olrun[cod->Get(pctp->RunHeader)-1]==104) continue; |
278 |
|
if(value%10!=0){ |
279 |
|
if(i==0) cin=cod->Get(pctp->CalibTrk1)+((value%10)-1); |
280 |
|
if(cod->Get(pctp->CalibTrk1)==cin+1){ |
281 |
|
eve=i+3; |
282 |
|
cin=100; |
283 |
|
} |
284 |
|
} |
285 |
|
else { |
286 |
|
if(i==0) cin=cod->Get(pctp->CalibTrk1); |
287 |
|
if(event<0 && cod->Get(pctp->CalibTrk1)==cin+1){ |
288 |
|
eve=i+3; |
289 |
|
cin=100; |
290 |
|
} |
291 |
|
else if(event>=0 && (ULong_t)ph->GetOrbitalTime()>obt && obt>0){ |
292 |
|
eve=i+3; |
293 |
|
obt=0; |
294 |
|
} |
295 |
|
} |
296 |
|
if(i==0) cins=cod->Get(pctp->CalibTrk1); |
297 |
|
if(cod->Get(pctp->CalibTrk1)>cins){ |
298 |
|
psum=1; |
299 |
|
cins=cod->Get(pctp->CalibTrk1); |
300 |
|
} |
301 |
|
if(psum==1){ |
302 |
|
figsa=out+ffile+"_FTrkScanQLook_EXPERT_sumof"; |
303 |
|
figsa+=totsum; |
304 |
|
figsa+="events_cal"; |
305 |
|
figsa+=cins; |
306 |
|
csum[cins-1] = new TCanvas(figsa.Data(),"FTrkQLookSCAN",canvasx,canvasy); |
307 |
|
csum[cins-1]->SetFillColor(10); |
308 |
|
csum[cins-1]->Range(0,0,1,1); |
309 |
|
fromfile.str(""); |
310 |
|
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Sum of "<<totsum<<" events after the "<<cins<<" calibration at OBT= "<<hobt[(cod->Get(pctp->RunHeader))-1]; |
311 |
|
t->SetTextSize(0.02); |
312 |
|
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
313 |
|
|
314 |
|
// draw pads |
315 |
|
Double_t posy = 0.95; // up y-coord - top pads |
316 |
|
Double_t hpad = 0.15; // pad height |
317 |
|
Double_t posx1=0; // left x-coord - pad column |
318 |
|
Double_t posx0=0; // x-coord - column division |
319 |
|
Double_t wrel = 0.6; // relative x size of first sub-column |
320 |
|
Double_t marg = 0.004; // margin among pads |
321 |
|
|
322 |
|
hpad = (posy-marg*11)/6; |
323 |
|
wrel = (1-marg*4)/2; |
324 |
|
for(Int_t n = 0; n<12; n++){ |
325 |
|
if ( (n+1)%2 ) { |
326 |
|
if(n>1)posy = posy-(marg*2+hpad); |
327 |
|
posx1 = marg; |
328 |
|
posx0 = posx1 + wrel; |
329 |
|
} else { |
330 |
|
posx1 = posx0 + 2*marg; |
331 |
|
posx0 = posx1 + wrel; |
332 |
|
}; |
333 |
|
|
334 |
|
/* -----------> pad for histograms */ |
335 |
|
trkpad[n] = new TPad("pad"," ",posx1,posy-hpad,posx0,posy,18,0,0); |
336 |
|
trkpad[n]->SetFillColor(19); |
337 |
|
trkpad[n]->SetFrameFillColor(10); |
338 |
|
|
339 |
|
p[n]=0; |
340 |
|
} |
341 |
|
} |
342 |
|
if(psum>0 && psum<=totsum){ |
343 |
|
// = = = = = = = = = = = = = = = = = = = = = = = = = |
344 |
|
// create header dump retrieving event info |
345 |
|
// = = = = = = = = = = = = = = = = = = = = = = = = = |
346 |
|
Int_t ndsp=0; |
347 |
|
|
348 |
|
// Int_t whistomax[3072]; |
349 |
|
Int_t whisto[3072]; |
350 |
|
// Int_t whistocomp[3072]; |
351 |
|
Int_t whistofull[3072]; |
352 |
|
|
353 |
|
// // transmitted words |
354 |
|
Int_t word = 0; |
355 |
|
Int_t iword = 0; |
356 |
|
// Int_t ii=0,ifull[12],icomp[12],imax[12],nn=0; |
357 |
|
Int_t ii=0,nn=0; |
358 |
|
trkword thisword; |
359 |
|
|
360 |
|
Int_t address,ladder; |
361 |
|
|
362 |
|
|
363 |
|
for(Int_t n = 0; n<12; n++){ |
364 |
|
ndsp = trk->DSPnumber[n]; |
365 |
|
nn = ndsp-1; |
366 |
|
// ifull[nn]=0; |
367 |
|
// icomp[nn]=0; |
368 |
|
// imax[nn]=0; |
369 |
|
|
370 |
|
for(Int_t vi = 0; vi< 3072; vi++){ |
371 |
|
// whistomax[vi] = -200; |
372 |
|
// whistocomp[vi] = -200; |
373 |
|
whistofull[vi] = -200; |
374 |
|
whisto[vi] = -200; |
375 |
|
} |
376 |
|
|
377 |
|
// // trasmitted data |
378 |
|
// |
379 |
|
address = 0; |
380 |
|
ladder = 1; |
381 |
|
for(Int_t vi = 0; vi < trk->DATAlength[n] ; vi++){ |
382 |
|
word = trk->TrackerData.At(iword); |
383 |
|
thisword = datadecode(word); |
384 |
|
iword++; |
385 |
|
switch (thisword.type){ |
386 |
|
|
387 |
|
case 0: //ADC value |
388 |
|
whisto[address] = thisword.decode; |
389 |
|
address++; |
390 |
|
// cout << " adr " << address << "\n"; |
391 |
|
break; |
392 |
|
|
393 |
|
case 1: //address |
394 |
|
address = 1024*(ladder-1) + thisword.decode; |
395 |
|
// cout << " adr " << address << "\n"; |
396 |
|
break; |
397 |
|
|
398 |
|
case 2: //end-of-ladder |
399 |
|
ladder = thisword.decode; |
400 |
|
// cout << "Ladder " << ladder << "\n"; |
401 |
|
if(ladder==3){ |
402 |
|
// end of compressed data - FILL HISTO |
403 |
|
//cout << ">>> COMPRESSED data" << "\n"; |
404 |
|
for(ii = 0; ii < 3072; ii++){ |
405 |
|
// whistocomp[ii]=whisto[ii]; |
406 |
|
whisto[ii] = -200; |
407 |
|
} |
408 |
|
address = 0; |
409 |
|
ladder = 1; |
410 |
|
}else if(ladder==6){ |
411 |
|
// end of full data - FILL HISTO |
412 |
|
//cout << ">>> FULL data" << "\n"; |
413 |
|
for(ii = 0; ii < 3072; ii++){ |
414 |
|
whistofull[ii]=whisto[ii]; |
415 |
|
whisto[ii] = -200; |
416 |
|
} |
417 |
|
address = 0; |
418 |
|
ladder = 1; |
419 |
|
}else{ |
420 |
|
if(ladder>3) ladder=ladder-3; |
421 |
|
address= ladder*1024; |
422 |
|
ladder = ladder + 1; |
423 |
|
} |
424 |
|
} |
425 |
|
} |
426 |
|
|
427 |
|
for(Int_t vi = 0; vi < 3072; vi++){ |
428 |
|
if(whistofull[vi]>-200){ |
429 |
|
if(vi==3071){ |
430 |
|
sum[nn][cins-1][p[nn]] = new TGraph(3072,whistostrip,whistofull); |
431 |
|
p[nn]++; |
432 |
|
} |
433 |
|
if(value%10==0 && value>0){ |
434 |
|
if(whistofull[vi]>0){ |
435 |
|
if(vi==0) values << (short int)ndsp << " "; |
436 |
|
values <<(short int) whistofull[vi] <<" "; |
437 |
|
if(vi==3071) values << endl; |
438 |
|
} |
439 |
|
} |
440 |
|
} |
441 |
|
} |
442 |
|
} |
443 |
|
psum++; |
444 |
|
if(value%10==0 && value>0) values << (short int)0 << endl << endl; |
445 |
|
} |
446 |
|
if(psum==totsum+1){ |
447 |
|
if(value%10==0 && value>0) values << (short int)(-cins) << endl << endl; |
448 |
|
psum=0; |
449 |
|
for(Int_t vi = 0; vi < 12; vi++){ |
450 |
|
if(p[vi]>0){ |
451 |
|
|
452 |
|
TLine li,liva1; |
453 |
|
li.SetLineColor(38); |
454 |
|
li.SetLineStyle(4); |
455 |
|
li.SetLineWidth(2); |
456 |
|
liva1.SetLineColor(42); |
457 |
|
liva1.SetLineStyle(3); |
458 |
|
liva1.SetLineWidth(1); |
459 |
|
|
460 |
|
Float_t va1x=0; |
461 |
|
csum[cins-1]->cd(); |
462 |
|
trkpad[vi]->Draw(); |
463 |
|
trkpad[vi]->cd(); |
464 |
|
trkpad[vi]->SetFillColor(10); |
465 |
|
|
466 |
|
stringstream tit,hid; |
467 |
|
tit.str(""); |
468 |
|
hid.str(""); |
469 |
|
tit<<"DSP "<<vi+1; |
470 |
|
sum[vi][cins-1][0]->SetTitle(tit.str().c_str()); |
471 |
|
sum[vi][cins-1][0]->GetXaxis()->SetTitle("channel id"); |
472 |
|
sum[vi][cins-1][0]->GetXaxis()->CenterTitle(); |
473 |
|
sum[vi][cins-1][0]->GetXaxis()->SetRangeUser(0,3073); |
474 |
|
sum[vi][cins-1][0]->GetXaxis()->SetTitleSize(0.06); |
475 |
|
sum[vi][cins-1][0]->GetXaxis()->SetTitleOffset(0.8); |
476 |
|
sum[vi][cins-1][0]->GetYaxis()->SetTitle("ADC value"); |
477 |
|
sum[vi][cins-1][0]->GetYaxis()->SetTitleSize(0.09); |
478 |
|
sum[vi][cins-1][0]->GetYaxis()->SetTitleOffset(0.4); |
479 |
|
sum[vi][cins-1][0]->GetYaxis()->CenterTitle(); |
480 |
|
sum[vi][cins-1][0]->GetYaxis()->SetRangeUser(0,4500); |
481 |
|
sum[vi][cins-1][0]->SetMarkerStyle(20); |
482 |
|
sum[vi][cins-1][0]->SetMarkerSize(0.1); |
483 |
|
sum[vi][cins-1][0]->SetMarkerColor(4); |
484 |
|
sum[vi][cins-1][0]->Draw("ap"); |
485 |
|
hid<<p[vi]<<" events"; |
486 |
|
t->SetTextSize(0.11); |
487 |
|
t->DrawLatex(2400,4250,hid.str().c_str()); |
488 |
|
for(int pi=1; pi<p[vi];pi++){ |
489 |
|
sum[vi][cins-1][pi]->SetMarkerStyle(21); |
490 |
|
sum[vi][cins-1][pi]->SetMarkerSize(0.1); |
491 |
|
sum[vi][cins-1][pi]->SetMarkerColor(4); |
492 |
|
sum[vi][cins-1][pi]->Draw("p"); |
493 |
|
} |
494 |
|
for(int va=1; va<24; va++){ |
495 |
|
va1x=128*va; |
496 |
|
liva1.DrawLine(va1x,0.,va1x,4500.); |
497 |
|
} |
498 |
|
li.DrawLine(1024.5,0.,1024.5,4500.); |
499 |
|
li.DrawLine(2048.5,0.,2048.5,4500.); |
500 |
|
csum[cins-1]->Update(); |
501 |
|
hid.str(""); |
502 |
|
} |
503 |
|
} |
504 |
|
if(strcmp(outfile.Data(),"ps") && strcmp(outfile.Data(),"pdf")){ |
505 |
|
figsa+="."+outfile; |
506 |
|
csum[cins-1]->Print(figsa.Data()); |
507 |
|
} |
508 |
|
if(cod->Get(pctp->CalibTrk1)==cevents) break; |
509 |
|
} |
510 |
|
} |
511 |
|
} |
512 |
|
|
513 |
|
if(!strcmp(outfile.Data(),"ps") || !strcmp(outfile.Data(),"pdf")){ |
514 |
|
stringstream nom1,nom2,nom3; |
515 |
|
nom1<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps("; |
516 |
|
nom2<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps"; |
517 |
|
nom3<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps)"; |
518 |
|
|
519 |
|
for(int va=0; va<cevents; va++){ |
520 |
|
if(cevents==1) csum[va]->Print(nom2.str().c_str(),"Landscape"); |
521 |
|
if(cevents>2){ |
522 |
|
if(va==0) csum[va]->Print(nom1.str().c_str(),"Landscape"); |
523 |
|
if(va>0 && va<cevents-1) csum[va]->Print(nom2.str().c_str(),"Landscape"); |
524 |
|
if(va>0 && va==cevents-1) csum[va]->Print(nom3.str().c_str(),"Landscape"); |
525 |
|
} |
526 |
|
} |
527 |
|
} |
528 |
|
// |
529 |
|
// Convert ps to pdf if required |
530 |
|
if(!strcmp(outfile.Data(),"pdf")){ |
531 |
|
stringstream com; |
532 |
|
com<<"ps2pdf13 "<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps "<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.pdf"; |
533 |
|
system(com.str().c_str()); |
534 |
|
printf("\n---> ps file converted in pdf format!\n"); |
535 |
|
com.str(""); |
536 |
|
com<<"rm -f "<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps "; |
537 |
|
system(com.str().c_str()); |
538 |
|
printf("---> ps file removed!\n\n"); |
539 |
|
com.str(""); |
540 |
|
} |
541 |
|
|
542 |
|
|
543 |
|
|
544 |
for(Int_t e=0;e<TOT;e++){ |
for(Int_t e=0;e<TOT;e++){ |
545 |
if(event<=0) |
if(event<=0 || (value%10!=0 && e==0)) |
546 |
event=eve; |
event=eve; |
547 |
else { |
else { |
548 |
event=event+1; |
event=event+1; |
556 |
tr->GetEntry(event); |
tr->GetEntry(event); |
557 |
//============================================================================ |
//============================================================================ |
558 |
|
|
559 |
gStyle->SetLabelSize(0.06,"x"); |
// draw display area |
560 |
gStyle->SetLabelSize(0.06,"y"); |
|
561 |
//gStyle->SetTitleFontSize(0.1); |
TPad *trkpad[12],*pad=0; //pad for histos |
562 |
gStyle->SetFillColor(10); |
TPaveText *trkpadtext[12]; //pad for header |
563 |
gStyle->SetTitleFillColor(10); |
Int_t ndsp=0; |
564 |
gStyle->SetTitleOffset(-1,"Y"); |
|
565 |
gStyle->SetOptStat(0); |
if((tot<TOT && e<tot) || tot==TOT){ |
|
|
|
|
// draw display area |
|
|
|
|
|
stringstream fromfile,title,hid,message; |
|
|
TString figsa="",figsav="",figsava=""; |
|
|
|
|
|
TPad *trkpad[12],*pad=0; //pad for histos |
|
|
TPaveText *trkpadtext[12]; //pad for header |
|
|
|
|
|
if((tot<TOT && e<tot) || tot==TOT){ |
|
|
Int_t canvasx=1200; |
|
|
Int_t canvasy=900; |
|
566 |
figsav=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
figsav=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
567 |
figsav+=event+1; |
figsav+=event+1; |
568 |
c1[e] = new TCanvas(figsav.Data(),"FTrkQLookSCAN",canvasx,canvasy); |
c1[e] = new TCanvas(figsav.Data(),"FTrkQLookSCAN",canvasx,canvasy); |
570 |
c1[e]->Range(0,0,1,1); |
c1[e]->Range(0,0,1,1); |
571 |
fromfile.str(""); |
fromfile.str(""); |
572 |
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event; |
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event; |
573 |
|
t->SetTextSize(0.02); |
574 |
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
575 |
t->DrawLatex(0.60,0.98,cal.Data()); |
t->DrawLatex(0.60,0.98,cal.Data()); |
576 |
|
|
577 |
|
isfile<<"Chips with white line at least once showed anomalous behaviour"; |
578 |
|
t->SetTextColor(17); |
579 |
|
t->SetTextSize(0.018); |
580 |
|
t->DrawLatex(.01,.96,isfile.str().c_str()); |
581 |
|
isfile.str(""); |
582 |
|
|
583 |
|
isfile<<"Boxes so colored point out anomalous chips (not necessarily broken)"; |
584 |
|
t->SetTextColor(107); |
585 |
|
// t->DrawLatex(.5,.96,isfile.str().c_str()); |
586 |
|
isfile.str(""); |
587 |
|
t->SetTextColor(1); |
588 |
|
|
589 |
if(va1!=0){ |
if(va1!=0){ |
590 |
figsava=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
figsava=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
630 |
/* -----------> pad for header dump */ |
/* -----------> pad for header dump */ |
631 |
trkpadtext[n] = new TPaveText((posx0+marg),(posy-hpad),posx2,posy); |
trkpadtext[n] = new TPaveText((posx0+marg),(posy-hpad),posx2,posy); |
632 |
/* -----------> HISTOGRAMS */ |
/* -----------> HISTOGRAMS */ |
633 |
|
|
634 |
|
title.str(""); |
635 |
title<<"DSP "<<n+1; |
title<<"DSP "<<n+1; |
636 |
hid<<"h"<<n+e*100; |
hid<<"h"<<n+e*100; |
637 |
histomax[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
histomax[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
653 |
pad = new TPad("padva"," ",0,0,1,0.97,18,0,0); |
pad = new TPad("padva"," ",0,0,1,0.97,18,0,0); |
654 |
pad->SetFillColor(19); |
pad->SetFillColor(19); |
655 |
pad->SetFrameFillColor(10); |
pad->SetFrameFillColor(10); |
656 |
} |
} |
657 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
658 |
// create header dump retrieving event info |
// create header dump retrieving event info |
659 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
660 |
Int_t ndsp=0; |
|
661 |
|
|
662 |
Double_t whistomax[3072]; |
Double_t whistomax[3072]; |
663 |
Double_t whisto[3072]; |
Double_t whisto[3072]; |
823 |
if(whistomax[i]>-200) imax[nn]=1; |
if(whistomax[i]>-200) imax[nn]=1; |
824 |
if(whistocomp[i]>-200) icomp[nn]=1; |
if(whistocomp[i]>-200) icomp[nn]=1; |
825 |
if(whistofull[i]>-200) ifull[nn]=1; |
if(whistofull[i]>-200) ifull[nn]=1; |
826 |
if(totvalues==TOT || (totvalues<TOT && e<totvalues)){ |
if(value%10!=0){ |
827 |
if(whistofull[i]>0){ |
if((totvalues<TOT && e<totvalues) || totvalues==TOT){ |
828 |
if(i==0) values << ndsp << " "; |
if(whistofull[i]>0){ |
829 |
values <<(short int) whistofull[i] <<" "; |
if(i==0) values << (short int)ndsp << " "; |
830 |
if(i==3071) values << endl; |
values <<(short int) whistofull[i] <<" "; |
831 |
|
if(i==3071) values << endl; |
832 |
|
} |
833 |
} |
} |
834 |
} |
} |
835 |
if((tot<TOT && e<tot) || tot==TOT){ |
if((tot<TOT && e<tot) || tot==TOT){ |
853 |
liva1.SetLineStyle(3); |
liva1.SetLineStyle(3); |
854 |
liva1.SetLineWidth(1); |
liva1.SetLineWidth(1); |
855 |
|
|
856 |
Float_t va1x=0; |
Float_t va1x=0,maxhist=4500.,minhist=-500.; |
857 |
|
|
858 |
if(va1!=0){ |
if(va1!=0){ |
859 |
cva[e]->cd(); |
cva[e]->cd(); |
898 |
if(icomp[nn]==1) histocomp[nn][e]->Draw("9bsame]["); |
if(icomp[nn]==1) histocomp[nn][e]->Draw("9bsame]["); |
899 |
if(imax[nn]==1) histomax[nn][e]->Draw("same]["); |
if(imax[nn]==1) histomax[nn][e]->Draw("same]["); |
900 |
histocomp[nn][e]->Draw("axis same"); |
histocomp[nn][e]->Draw("axis same"); |
901 |
|
|
902 |
|
b.SetFillColor(19); |
903 |
|
b.SetFillStyle(3945); |
904 |
if(nn==0){ |
if(nn==0){ |
905 |
b.SetFillColor(107); |
b.DrawBox(0.,minhist,2047.,maxhist); |
|
b.SetFillStyle(3945); |
|
|
b.DrawBox(768.,-500.,2047.,4500.); |
|
906 |
} |
} |
907 |
else if(nn==1){ |
else if(nn==1){ |
908 |
b.SetFillColor(6); |
b.DrawBox(128.,minhist,256.,maxhist); |
909 |
b.SetFillStyle(3945); |
b.DrawBox(384.,minhist,512.,maxhist); |
910 |
b.DrawBox(2944.,-500.,3060.,4500.); |
b.DrawBox(896.,minhist,1024.,maxhist); |
911 |
|
b.DrawBox(2048.,minhist,2432.,maxhist); |
912 |
b.SetFillColor(107); |
b.DrawBox(2816.,minhist,2944.,maxhist); |
913 |
b.SetFillStyle(3954); |
b.DrawBox(2944.,minhist,3070.,maxhist); |
914 |
//b.DrawBox(384.,-500.,512.,4500.); |
} |
915 |
b.DrawBox(2816.,-500.,2944.,4500.); |
else if(nn==3){ |
916 |
b.DrawBox(2048.,-500.,2176.,4500.); |
b.DrawBox(0.,minhist,256.,maxhist); |
917 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
918 |
} |
} |
919 |
else if(nn==4){ |
else if(nn==4){ |
920 |
b.SetFillColor(107); |
b.DrawBox(256.,minhist,512.,maxhist); |
921 |
b.SetFillStyle(3954); |
b.DrawBox(1792.,minhist,1920.,maxhist); |
922 |
b.DrawBox(384.,-500.,512.,4500.); |
b.DrawBox(2816.,minhist,3070.,maxhist); |
923 |
|
} |
924 |
|
else if(nn==5){ |
925 |
|
b.DrawBox(0.,minhist,256.,maxhist); |
926 |
|
b.DrawBox(896.,minhist,1024.,maxhist); |
927 |
|
b.DrawBox(1664.,minhist,1792.,maxhist); |
928 |
} |
} |
929 |
else if(nn==6){ |
else if(nn==6){ |
930 |
b.SetFillColor(6); |
b.DrawBox(512.,minhist,768.,maxhist); |
931 |
b.SetFillStyle(3945); |
b.DrawBox(1024.,minhist,1280.,maxhist); |
932 |
b.DrawBox(2560.,-500.,2816.,4500.); |
b.DrawBox(1280.,minhist,1792.,maxhist); |
933 |
b.DrawBox(1024.,-500.,1535.,4500.); |
b.DrawBox(2560.,minhist,2816.,maxhist); |
|
|
|
|
b.SetFillColor(107); |
|
|
b.SetFillStyle(3954); |
|
|
b.DrawBox(512.,-500.,768.,4500.); |
|
|
b.DrawBox(1536.,-500.,1792.,4500.); |
|
934 |
} |
} |
935 |
else if(nn==7){ |
else if(nn==7){ |
936 |
b.SetFillColor(107); |
b.DrawBox(0.,minhist,1535.,maxhist); |
937 |
b.SetFillStyle(3954); |
b.DrawBox(2024.,minhist,2280.,maxhist); |
|
b.DrawBox(512.,-500.,768.,4500.); |
|
938 |
} |
} |
939 |
else if(nn==8){ |
else if(nn==8){ |
940 |
b.SetFillColor(107); |
b.DrawBox(512.,minhist,768.,maxhist); |
|
b.SetFillStyle(3954); |
|
|
b.DrawBox(512.,-500.,768.,4500.); |
|
941 |
} |
} |
942 |
else if(nn==9){ |
else if(nn==9){ |
943 |
b.SetFillColor(107); |
b.DrawBox(0.,minhist,128.,maxhist); |
944 |
b.SetFillStyle(3954); |
b.DrawBox(256.,minhist,384.,maxhist); |
945 |
b.DrawBox(256.,-500.,384.,4500.); |
b.DrawBox(512.,minhist,640.,maxhist); |
946 |
//b.DrawBox(1280.,-500.,1408.,4500.); |
b.DrawBox(896.,minhist,1152.,maxhist); |
947 |
//b.DrawBox(1792.,-500.,1920.,4500.); |
b.DrawBox(1280.,minhist,1535.,maxhist); |
948 |
|
b.DrawBox(1664.,minhist,1920.,maxhist); |
949 |
|
b.DrawBox(2048.,minhist,2304.,maxhist); |
950 |
} |
} |
951 |
else if(nn==10){ |
else if(nn==10){ |
952 |
b.SetFillColor(107); |
b.DrawBox(0.,minhist,512.,maxhist); |
953 |
b.SetFillStyle(3954); |
b.DrawBox(1024.,minhist,3070.,maxhist); |
|
b.DrawBox(2048.,-500.,3070.,4500.); |
|
954 |
} |
} |
955 |
else if(nn==11){ |
else if(nn==11){ |
956 |
b.SetFillColor(6); |
b.DrawBox(0.,minhist,512.,maxhist); |
957 |
b.SetFillStyle(3945); |
b.DrawBox(768.,minhist,1024.,maxhist); |
958 |
b.DrawBox(768.,-500.,1024.,4500.); |
b.DrawBox(1536.,minhist,1664.,maxhist); |
959 |
|
b.DrawBox(1920.,minhist,2560.,maxhist); |
960 |
b.SetFillColor(107); |
b.DrawBox(2816.,minhist,3070.,maxhist); |
961 |
b.SetFillStyle(3954); |
} |
|
b.DrawBox(0.,-500.,512.,4500.); |
|
|
b.DrawBox(1920.,-500.,2560.,4500.); |
|
|
} |
|
962 |
for(int va=1; va<24; va++){ |
for(int va=1; va<24; va++){ |
963 |
va1x=128*va; |
va1x=128*va; |
964 |
liva1.DrawLine(va1x,-500.,va1x,4500.); |
liva1.DrawLine(va1x,minhist,va1x,maxhist); |
965 |
} |
} |
966 |
li.DrawLine(1024.5,-500.,1024.5,4500.); |
li.DrawLine(1024.5,minhist,1024.5,maxhist); |
967 |
li.DrawLine(2048.5,-500.,2048.5,4500.); |
li.DrawLine(2048.5,minhist,2048.5,maxhist); |
968 |
c1[e]->Update(); |
c1[e]->Update(); |
969 |
} |
} |
970 |
|
|
979 |
|
|
980 |
if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){ |
if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){ |
981 |
if(e==0){ |
if(e==0){ |
982 |
c1[e]->Print(nom1.str().c_str(),"Portrait"); |
c1[e]->Print(nom1.str().c_str(),"Landscape"); |
983 |
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Portrait"); |
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Landscape"); |
984 |
} |
} |
985 |
if(e>0 && tot>2){ |
if(e>0 && tot>2){ |
986 |
c1[e]->Print(nom2.str().c_str(),"Portrait"); |
c1[e]->Print(nom2.str().c_str(),"Landscape"); |
987 |
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Portrait"); |
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Landscape"); |
988 |
} |
} |
989 |
if(e==tot-1){ |
if(e==tot-1){ |
990 |
c1[e]->Print(nom2.str().c_str(),"Portrait"); |
c1[e]->Print(nom2.str().c_str(),"Landscape"); |
991 |
if(va1!=0) cva[e]->Print(nom3.str().c_str(),"Portrait"); |
if(va1!=0) cva[e]->Print(nom3.str().c_str(),"Landscape"); |
992 |
} |
} |
993 |
} |
} |
994 |
else{ |
else{ |
1000 |
} |
} |
1001 |
} |
} |
1002 |
} |
} |
1003 |
|
if(value%10!=0){ |
1004 |
if(totvalues==TOT || (totvalues<TOT && e<totvalues)) |
if(totvalues==TOT || (totvalues<TOT && e<totvalues)) |
1005 |
values << 0 << endl << endl; |
values << (short int)0 << endl << endl; |
1006 |
|
} |
1007 |
} |
} |
1008 |
stringstream com; |
stringstream com; |
1009 |
|
|