2 |
* FTrkScanQlook_EXPERT.cxx |
* FTrkScanQlook_EXPERT.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version v1r10 |
* 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 |
* DSPprint - the number of a particular DSP to draw (0 don't draw) |
* va1 - the number of the single va1 to analyze (dsp*100+va1) |
11 |
* outfile - extension of output file (pdf,ps,gif,jpg) |
* 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) |
13 |
* |
* |
14 |
*/ |
*/ |
15 |
// |
// |
16 |
#include <iostream> |
#include <iostream> |
17 |
#include <sstream> |
#include <sstream> |
18 |
|
#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 |
}; |
}; |
126 |
} |
} |
127 |
|
|
128 |
|
|
129 |
void FTrkScanQLook_EXPERT(TString file, TString outdir,Int_t event, Int_t DSPprint, TString outfile) |
void FTrkScanQLook_EXPERT(TString file, TString outdir,Int_t event, Int_t va1, Int_t value, TString outfile) |
130 |
{ |
{ |
131 |
|
|
132 |
// |
// |
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; |
155 |
|
|
156 |
|
pamela::PacketType *pctp=0; |
157 |
// open files |
// open files |
158 |
TFile *trackerFile = new TFile(file); |
TFile *trackerFile = new TFile(file); |
159 |
if ( !trackerFile ){ |
if ( !trackerFile ){ |
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"); |
192 |
return; |
return; |
193 |
} |
} |
194 |
|
|
195 |
|
// |
196 |
|
// open the output text files for the alarms |
197 |
|
TString fname = out+ffile; |
198 |
|
int slen=fname.Length(); |
199 |
|
int last=fname.Last('_'); |
200 |
|
if(last<(slen-1))fname.Append('_'); |
201 |
|
fname.Append("FTrk-DSP-values.dat"); |
202 |
|
ofstream values(fname,ios::out); |
203 |
|
|
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 ev[2]; |
Int_t eve,cin=0,cins=0,br=0,olrun[neventH]; |
211 |
TString cal=""; |
TString cal=""; |
212 |
|
|
213 |
ev[0]=0; |
eve=3; |
|
ev[1]=0; |
|
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 |
ev[0]=2; |
eve=2; |
|
ev[1]=3; |
|
227 |
} |
} |
228 |
} |
} |
229 |
if(ev[0]==0){ |
|
230 |
|
int tot=2,totvalues=0,TOT=0,totsum=0; |
231 |
|
if(event<0) tot=abs(event); |
232 |
|
if(value>0) totvalues=value/10; |
233 |
|
|
234 |
|
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 |
241 |
|
TH1F *histocomp[12][tot]; //histos of compressed data |
242 |
|
TH1F *histofull[12][tot]; //histos of full data |
243 |
|
TCanvas *c1[tot],*csum[calev]; |
244 |
|
TH1F *histova[tot]; |
245 |
|
TGraph*sum[12][calev][totsum/2]; |
246 |
|
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(); |
256 |
|
t->SetTextFont(32); |
257 |
|
t->SetTextColor(1); |
258 |
|
t->SetTextAlign(12); |
259 |
|
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++){ |
for(Int_t i=0;i<nevents;i++){ |
274 |
tr->GetEntry(i); |
tr->GetEntry(i); |
275 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
276 |
if(ph->GetOrbitalTime()>obt){ |
cod = eh->GetCounter(); |
277 |
ev[0]=i+2; |
if(olrun[cod->Get(pctp->RunHeader)-1]==104) continue; |
278 |
ev[1]=i+3; |
if(value%10!=0){ |
279 |
break; |
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 |
|
|
|
TH1F *histomax[12][2]; //histos of max signals |
|
|
TH1F *histocomp[12][2]; //histos of compressed data |
|
|
TH1F *histofull[12][2]; //histos of full data |
|
|
TCanvas *c1[2]; |
|
542 |
|
|
543 |
for(Int_t e=0;e<2;e++){ |
|
544 |
event=ev[e]; |
for(Int_t e=0;e<TOT;e++){ |
545 |
|
if(event<=0 || (value%10!=0 && e==0)) |
546 |
|
event=eve; |
547 |
|
else { |
548 |
|
event=event+1; |
549 |
|
if(event>eve-3 && eve>2) |
550 |
|
cal="Event with online calibration"; |
551 |
|
else |
552 |
|
cal="***** ONLINE CALIBRATION NOT FOUND IN THIS FILE *****"; |
553 |
|
} |
554 |
printf("Scan of Entry %d\n",event); |
printf("Scan of Entry %d\n",event); |
555 |
|
|
556 |
tr->GetEntry(event); |
tr->GetEntry(event); |
557 |
//============================================================================ |
//============================================================================ |
558 |
|
|
|
gStyle->SetLabelSize(0.06,"x"); |
|
|
gStyle->SetLabelSize(0.06,"y"); |
|
|
gStyle->SetTitleFontSize(0.1); |
|
|
gStyle->SetFillColor(10); |
|
|
gStyle->SetTitleFillColor(10); |
|
|
gStyle->SetTitleOffset(-1,"Y"); |
|
|
gStyle->SetOptStat(0); |
|
|
|
|
559 |
// draw display area |
// draw display area |
560 |
|
|
561 |
Int_t canvasx=1200; |
TPad *trkpad[12],*pad=0; //pad for histos |
|
Int_t canvasy=900; |
|
|
stringstream figsav; |
|
|
figsav.str(""); |
|
|
figsav<<out<<ffile<<"_FTrkScanQLook_EXPERT_ev"<<event+1<<"."<<outfile.Data(); |
|
|
c1[e] = new TCanvas(figsav.str().c_str(),"FTrkQLookSCAN",canvasx,canvasy); |
|
|
c1[e]->SetFillColor(10); |
|
|
c1[e]->Range(0,0,1,1); |
|
|
stringstream fromfile; |
|
|
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event; |
|
|
TLatex *t=new TLatex(); |
|
|
t->SetTextFont(32); |
|
|
t->SetTextColor(1); |
|
|
t->SetTextAlign(12); |
|
|
t->SetTextSize(0.02); |
|
|
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
|
|
t->DrawLatex(0.60,0.98,cal.Data()); |
|
|
|
|
|
// draw pads |
|
|
TPad *trkpad[12]; //pad for histos |
|
562 |
TPaveText *trkpadtext[12]; //pad for header |
TPaveText *trkpadtext[12]; //pad for header |
563 |
Double_t posy = 0.95; // up y-coord - top pads |
Int_t ndsp=0; |
|
Double_t hpad = 0.15; // pad height |
|
|
Double_t posx1=0; // left x-coord - pad column |
|
|
Double_t posx2=0; // right x-coord - pad olumn |
|
|
Double_t posx0=0; // x-coord - column division |
|
|
Double_t wrel = 0.6; // relative x size of first sub-column |
|
|
Double_t marg = 0.004; // margin among pads |
|
564 |
|
|
565 |
|
if((tot<TOT && e<tot) || tot==TOT){ |
566 |
|
figsav=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
567 |
|
figsav+=event+1; |
568 |
|
c1[e] = new TCanvas(figsav.Data(),"FTrkQLookSCAN",canvasx,canvasy); |
569 |
|
c1[e]->SetFillColor(10); |
570 |
|
c1[e]->Range(0,0,1,1); |
571 |
|
fromfile.str(""); |
572 |
|
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event; |
573 |
|
t->SetTextSize(0.02); |
574 |
|
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
575 |
|
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){ |
590 |
|
figsava=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
591 |
|
figsava+=event+1; |
592 |
|
figsava+="_DSP"; |
593 |
|
figsava+=(int)(va1/100); |
594 |
|
figsava+="_VA1-"; |
595 |
|
figsava+=va1%100; |
596 |
|
cva[e] = new TCanvas(figsava.Data(),"TrkQLookSCAN VA1",canvasx,canvasy); |
597 |
|
cva[e]->SetFillColor(10); |
598 |
|
cva[e]->Range(0,0,1,1); |
599 |
|
fromfile.str(""); |
600 |
|
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event<<" --> DSP "<<(int)(va1/100)<<" --> va1 "<<va1%100; |
601 |
|
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
602 |
|
// t->DrawLatex(0.65,0.98,cal.Data()); |
603 |
|
} |
604 |
|
|
605 |
stringstream title; |
// draw pads |
606 |
stringstream hid; |
Double_t posy = 0.95; // up y-coord - top pads |
607 |
for(Int_t n = 0; n<12; n++){ |
Double_t hpad = 0.15; // pad height |
608 |
if ( (n+1)%2 ) { |
Double_t posx1=0; // left x-coord - pad column |
609 |
if(n>1)posy = posy-(marg*2+hpad); |
Double_t posx2=0; // right x-coord - pad olumn |
610 |
posx1 = marg; |
Double_t posx0=0; // x-coord - column division |
611 |
posx2 = 0.5 - marg; |
Double_t wrel = 0.6; // relative x size of first sub-column |
612 |
posx0 = 0.5*wrel; |
Double_t marg = 0.004; // margin among pads |
613 |
} else { |
|
614 |
posx1 = posx1 + 0.5; |
for(Int_t n = 0; n<12; n++){ |
615 |
posx2 = posx2 + 0.5; |
if ( (n+1)%2 ) { |
616 |
posx0 = posx0 + 0.5; |
if(n>1)posy = posy-(marg*2+hpad); |
617 |
}; |
posx1 = marg; |
618 |
|
posx2 = 0.5 - marg; |
619 |
|
posx0 = 0.5*wrel; |
620 |
|
} else { |
621 |
|
posx1 = posx1 + 0.5; |
622 |
|
posx2 = posx2 + 0.5; |
623 |
|
posx0 = posx0 + 0.5; |
624 |
|
}; |
625 |
|
|
626 |
|
/* -----------> pad for histograms */ |
627 |
|
trkpad[n] = new TPad("pad"," ",posx1,posy-hpad,posx0-marg,posy,18,0,0); |
628 |
|
trkpad[n]->SetFillColor(19); |
629 |
|
trkpad[n]->SetFrameFillColor(10); |
630 |
|
/* -----------> pad for header dump */ |
631 |
|
trkpadtext[n] = new TPaveText((posx0+marg),(posy-hpad),posx2,posy); |
632 |
|
/* -----------> HISTOGRAMS */ |
633 |
|
|
634 |
|
title.str(""); |
635 |
|
title<<"DSP "<<n+1; |
636 |
|
hid<<"h"<<n+e*100; |
637 |
|
histomax[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
638 |
|
hid<<"hh"<<n+e*100; |
639 |
|
histocomp[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
640 |
|
hid<<"hhh"<<n+e*100; |
641 |
|
histofull[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
642 |
|
title.str(""); |
643 |
|
hid.str(""); |
644 |
|
} |
645 |
|
if(va1!=0){ |
646 |
|
title.str(""); |
647 |
|
title<<"DSP "<<(int)(va1/100)<<" -- va1 "<<va1%100; |
648 |
|
hid<<"va"<<e; |
649 |
|
histova[e] = new TH1F(hid.str().c_str(),title.str().c_str(),130,0,129); |
650 |
|
} //end loop on views |
651 |
|
|
652 |
/* -----------> pad for histograms */ |
/* -----------> pad for histograms */ |
653 |
trkpad[n] = new TPad("pad"," ",posx1,posy-hpad,posx0-marg,posy,18,0,0); |
pad = new TPad("padva"," ",0,0,1,0.97,18,0,0); |
654 |
trkpad[n]->SetFillColor(19); |
pad->SetFillColor(19); |
655 |
trkpad[n]->SetFrameFillColor(10); |
pad->SetFrameFillColor(10); |
656 |
/* -----------> pad for header dump */ |
} |
|
trkpadtext[n] = new TPaveText((posx0+marg),(posy-hpad),posx2,posy); |
|
|
/* -----------> HISTOGRAMS */ |
|
|
|
|
|
title<<"DSP "<<n+1; |
|
|
hid<<"h"<<n+e*100; |
|
|
histomax[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
|
|
hid<<"hh"<<n+e*100; |
|
|
histocomp[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
|
|
hid<<"hhh"<<n+e*100; |
|
|
histofull[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
|
|
title.str(""); |
|
|
hid.str(""); |
|
|
} //end loop on views |
|
|
|
|
657 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
658 |
// create header dump retrieving event info |
// create header dump retrieving event info |
659 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
660 |
Int_t ndsp=0; |
|
|
stringstream message; |
|
661 |
|
|
662 |
Double_t whistomax[3072]; |
Double_t whistomax[3072]; |
663 |
Double_t whisto[3072]; |
Double_t whisto[3072]; |
673 |
trkword thisword; |
trkword thisword; |
674 |
|
|
675 |
Int_t address,ladder; |
Int_t address,ladder; |
676 |
|
|
677 |
|
|
678 |
for(Int_t n = 0; n<12; n++){ |
for(Int_t n = 0; n<12; n++){ |
679 |
|
|
680 |
ndsp = trk->DSPnumber[n]; |
ndsp = trk->DSPnumber[n]; |
682 |
ifull[nn]=0; |
ifull[nn]=0; |
683 |
icomp[nn]=0; |
icomp[nn]=0; |
684 |
imax[nn]=0; |
imax[nn]=0; |
685 |
|
if(ndsp>0){ |
686 |
|
if(ndsp<13){ |
687 |
|
if((tot<TOT && e<tot) || tot==TOT){ |
688 |
|
|
689 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
690 |
* |
* |
691 |
* Write event LEVEL0 report |
* Write event LEVEL0 report |
692 |
* |
* |
693 |
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*/ |
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*/ |
694 |
|
|
695 |
trkpadtext[nn]->SetTextFont(40); |
trkpadtext[nn]->SetTextFont(40); |
696 |
trkpadtext[nn]->SetFillColor(33); |
trkpadtext[nn]->SetFillColor(33); |
697 |
trkpadtext[nn]->SetTextSize(0.012); |
trkpadtext[nn]->SetTextSize(0.012); |
698 |
trkpadtext[nn]->SetTextAlign(13); |
trkpadtext[nn]->SetTextAlign(13); |
699 |
|
|
700 |
trkpadtext[nn]->AddText(" "); |
trkpadtext[nn]->AddText(" "); |
701 |
message<<"DAQ mode --------> "<<trk->DAQmode[n]; |
message<<"DAQ mode --------> "<<trk->DAQmode[n]; |
702 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
703 |
message.str(""); |
message.str(""); |
704 |
message<<"Event number --------> "<<trk->eventn[n]; |
message<<"Event number --------> "<<trk->eventn[n]; |
705 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
706 |
message.str(""); |
message.str(""); |
707 |
message<<"13-bit words --------> "<<trk->DATAlength[n]; |
message<<"13-bit words --------> "<<trk->DATAlength[n]; |
708 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
709 |
message.str(""); |
message.str(""); |
710 |
if (!(nn%2)&&trk->signcluster[n][0]!=0) message<<"L1 add: "<<trk->addrcluster[n][0]<<" - sign: "<<1024-(trk->signcluster[n][0]); |
if (!(nn%2)&&trk->signcluster[n][0]!=0) message<<"L1 add: "<<trk->addrcluster[n][0]<<" - sign: "<<1024-(trk->signcluster[n][0]); |
711 |
else message<<"L1 add: "<<trk->addrcluster[n][0]<<" - sign: "<<(trk->signcluster[n][0]); |
else message<<"L1 add: "<<trk->addrcluster[n][0]<<" - sign: "<<(trk->signcluster[n][0]); |
712 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
713 |
message.str(""); |
message.str(""); |
714 |
if (!(nn%2)&&trk->signcluster[n][1]!=0) message<<"L2 add: "<<trk->addrcluster[n][1]<<" - sign: "<<1024-(trk->signcluster[n][1]); |
if (!(nn%2)&&trk->signcluster[n][1]!=0) message<<"L2 add: "<<trk->addrcluster[n][1]<<" - sign: "<<1024-(trk->signcluster[n][1]); |
715 |
else message<<"L2 add: "<<trk->addrcluster[n][1]<<" - sign: "<<trk->signcluster[n][1]; |
else message<<"L2 add: "<<trk->addrcluster[n][1]<<" - sign: "<<trk->signcluster[n][1]; |
716 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
717 |
message.str(""); |
message.str(""); |
718 |
if (!(nn%2)&&trk->signcluster[n][2]!=0) message<<"L3 add: "<<trk->addrcluster[n][2]<<" - sign: "<<1024-(trk->signcluster[n][2]); |
if (!(nn%2)&&trk->signcluster[n][2]!=0) message<<"L3 add: "<<trk->addrcluster[n][2]<<" - sign: "<<1024-(trk->signcluster[n][2]); |
719 |
else message<<"L3 add: "<<trk->addrcluster[n][2]<<" - sign: "<<trk->signcluster[n][2]; |
else message<<"L3 add: "<<trk->addrcluster[n][2]<<" - sign: "<<trk->signcluster[n][2]; |
720 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
721 |
message.str(""); |
message.str(""); |
722 |
message<<"NCLUST "<<trk->nclust[n]<<" CUTC "<<trk->cutc[n]<<" CUTCL "<<trk->cutcl[n]; |
message<<"NCLUST "<<trk->nclust[n]<<" CUTC "<<trk->cutc[n]<<" CUTCL "<<trk->cutcl[n]; |
723 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
724 |
message.str(""); |
message.str(""); |
725 |
message<<"Comp. time "<<trk->compressiontime[n]<<" x 0.051ms = "<<0.051*trk->compressiontime[n]<<" ms"; |
message<<"Comp. time "<<trk->compressiontime[n]<<" x 0.051ms = "<<0.051*trk->compressiontime[n]<<" ms"; |
726 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
727 |
message.str(""); |
message.str(""); |
728 |
trkpadtext[nn]->AddText(" "); |
trkpadtext[nn]->AddText(" "); |
729 |
message<<"CRC --> "<<trk->crc[n]; |
message<<"CRC --> "<<trk->crc[n]; |
730 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
731 |
message.str(""); |
message.str(""); |
732 |
message<<"FL1-6 --> "<<trk->fl1[n]<<" "<<trk->fl2[n]<<" "<<trk->fl3[n]<<" "<<trk->fl4[n]<<" "<<trk->fl5[n]<<" "<<trk->fl6[n]<<" FC "<<trk->fc[n]; |
message<<"FL1-6 --> "<<trk->fl1[n]<<" "<<trk->fl2[n]<<" "<<trk->fl3[n]<<" "<<trk->fl4[n]<<" "<<trk->fl5[n]<<" "<<trk->fl6[n]<<" FC "<<trk->fc[n]; |
733 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
734 |
message.str(""); |
message.str(""); |
735 |
trkpadtext[nn]->AddText(" "); |
trkpadtext[nn]->AddText(" "); |
736 |
trkpadtext[nn]->AddLine(0,0,0,0); |
trkpadtext[nn]->AddLine(0,0,0,0); |
737 |
message<<"PNum "<<trk->pnum[n]<<" - BId "<<trk->bid[n]<<" ALARM "<<trk->alarm[n]; |
message<<"PNum "<<trk->pnum[n]<<" - BId "<<trk->bid[n]<<" ALARM "<<trk->alarm[n]; |
738 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
739 |
message.str(""); |
message.str(""); |
740 |
message<<"Cmd "<<trk->cmdnum[n]<<" --- Answer length "<<trk->aswr[n]<<" byte"; |
message<<"Cmd "<<trk->cmdnum[n]<<" --- Answer length "<<trk->aswr[n]<<" byte"; |
741 |
trkpadtext[nn]->AddText(message.str().c_str()); |
trkpadtext[nn]->AddText(message.str().c_str()); |
742 |
message.str(""); |
message.str(""); |
743 |
trkpadtext[nn]->AddText(" "); |
trkpadtext[nn]->AddText(" "); |
744 |
|
|
745 |
TOTDATAlength_check = TOTDATAlength_check + trk->DATAlength[n]; |
TOTDATAlength_check = TOTDATAlength_check + trk->DATAlength[n]; |
746 |
|
} |
747 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
748 |
* |
* |
749 |
* Plot event LEVEL0 histo |
* Plot event LEVEL0 histo |
750 |
* |
* |
751 |
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*/ |
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*/ |
752 |
|
|
753 |
//============================================= |
//============================================= |
754 |
|
} |
755 |
for(Int_t i = 0; i< 3072; i++){ |
for(Int_t i = 0; i< 3072; i++){ |
756 |
whistomax[i] = -200; |
whistomax[i] = -200; |
757 |
whistocomp[i] = -200; |
whistocomp[i] = -200; |
758 |
whistofull[i] = -200; |
whistofull[i] = -200; |
759 |
whisto[i] = -200; |
whisto[i] = -200; |
760 |
} |
} |
761 |
|
|
762 |
// =============== |
// =============== |
763 |
// trasmitted data |
// trasmitted data |
764 |
// =============== |
// =============== |
765 |
|
|
766 |
address = 0; |
address = 0; |
767 |
ladder = 1; |
ladder = 1; |
768 |
for(Int_t i = 0; i < trk->DATAlength[n] ; i++){ |
for(Int_t i = 0; i < trk->DATAlength[n] ; i++){ |
769 |
word = trk->TrackerData.At(iword); |
word = trk->TrackerData.At(iword); |
770 |
thisword = datadecode(word); |
thisword = datadecode(word); |
771 |
iword++; |
iword++; |
772 |
switch (thisword.type){ |
switch (thisword.type){ |
773 |
|
|
774 |
case 0: //ADC value |
case 0: //ADC value |
775 |
whisto[address] = thisword.decode; |
whisto[address] = thisword.decode; |
776 |
address++; |
address++; |
777 |
// cout << " adr " << address << "\n"; |
// cout << " adr " << address << "\n"; |
778 |
break; |
break; |
779 |
|
|
780 |
case 1: //address |
case 1: //address |
781 |
address = 1024*(ladder-1) + thisword.decode; |
address = 1024*(ladder-1) + thisword.decode; |
782 |
// cout << " adr " << address << "\n"; |
// cout << " adr " << address << "\n"; |
783 |
break; |
break; |
784 |
|
|
785 |
case 2: //end-of-ladder |
case 2: //end-of-ladder |
786 |
ladder = thisword.decode; |
ladder = thisword.decode; |
787 |
// cout << "Ladder " << ladder << "\n"; |
// cout << "Ladder " << ladder << "\n"; |
788 |
if(ladder==3){ |
if(ladder==3){ |
789 |
// end of compressed data - FILL HISTO |
// end of compressed data - FILL HISTO |
790 |
//cout << ">>> COMPRESSED data" << "\n"; |
//cout << ">>> COMPRESSED data" << "\n"; |
791 |
for(ii = 0; ii < 3072; ii++){ |
for(ii = 0; ii < 3072; ii++){ |
792 |
whistocomp[ii]=whisto[ii]; |
whistocomp[ii]=whisto[ii]; |
793 |
whisto[ii] = -200; |
whisto[ii] = -200; |
794 |
|
} |
795 |
|
address = 0; |
796 |
|
ladder = 1; |
797 |
|
}else if(ladder==6){ |
798 |
|
// end of full data - FILL HISTO |
799 |
|
//cout << ">>> FULL data" << "\n"; |
800 |
|
for(ii = 0; ii < 3072; ii++){ |
801 |
|
whistofull[ii]=whisto[ii]; |
802 |
|
whisto[ii] = -200; |
803 |
|
} |
804 |
|
address = 0; |
805 |
|
ladder = 1; |
806 |
|
}else{ |
807 |
|
if(ladder>3) ladder=ladder-3; |
808 |
|
address= ladder*1024; |
809 |
|
ladder = ladder + 1; |
810 |
} |
} |
811 |
address = 0; |
} |
|
ladder = 1; |
|
|
}else if(ladder==6){ |
|
|
// end of full data - FILL HISTO |
|
|
//cout << ">>> FULL data" << "\n"; |
|
|
for(ii = 0; ii < 3072; ii++){ |
|
|
whistofull[ii]=whisto[ii]; |
|
|
whisto[ii] = -200; |
|
|
} |
|
|
address = 0; |
|
|
ladder = 1; |
|
|
}else{ |
|
|
if(ladder>3) ladder=ladder-3; |
|
|
address= ladder*1024; |
|
|
ladder = ladder + 1; |
|
|
} |
|
812 |
} |
} |
|
} |
|
813 |
|
|
814 |
|
|
815 |
// =============== |
// =============== |
816 |
// maximum signals |
// maximum signals |
817 |
// =============== |
// =============== |
818 |
if(trk->signcluster[nn][0]!=0) whistomax[ 0+(int)trk->addrcluster[nn][0]] = whistocomp[ 0+(int)trk->addrcluster[nn][0]]; |
if(trk->signcluster[nn][0]!=0) whistomax[ 0+(int)trk->addrcluster[nn][0]] = whistocomp[ 0+(int)trk->addrcluster[nn][0]]; |
819 |
if(trk->signcluster[nn][1]!=0) whistomax[1024+(int)trk->addrcluster[nn][1]] = whistocomp[1024+(int)trk->addrcluster[nn][1]]; |
if(trk->signcluster[nn][1]!=0) whistomax[1024+(int)trk->addrcluster[nn][1]] = whistocomp[1024+(int)trk->addrcluster[nn][1]]; |
820 |
if(trk->signcluster[nn][2]!=0) whistomax[2048+(int)trk->addrcluster[nn][2]] = whistocomp[2048+(int)trk->addrcluster[nn][2]]; |
if(trk->signcluster[nn][2]!=0) whistomax[2048+(int)trk->addrcluster[nn][2]] = whistocomp[2048+(int)trk->addrcluster[nn][2]]; |
821 |
|
|
822 |
for(Int_t i = 0; i < 3072; i++){ |
for(Int_t i = 0; i < 3072; i++){ |
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 |
histomax[nn][e]->Fill((Float_t)i,whistomax[i]); |
if(value%10!=0){ |
827 |
histocomp[nn][e]->Fill((Float_t)i,whistocomp[i]); |
if((totvalues<TOT && e<totvalues) || totvalues==TOT){ |
828 |
histofull[nn][e]->Fill((Float_t)i,whistofull[i]); |
if(whistofull[i]>0){ |
829 |
} |
if(i==0) values << (short int)ndsp << " "; |
830 |
|
values <<(short int) whistofull[i] <<" "; |
831 |
TBox b; |
if(i==3071) values << endl; |
832 |
|
} |
833 |
c1[e]->cd(); |
} |
834 |
trkpadtext[nn]->Draw(); |
} |
835 |
trkpad[nn]->Draw(); |
if((tot<TOT && e<tot) || tot==TOT){ |
836 |
trkpad[nn]->cd(); |
histomax[nn][e]->Fill((Float_t)i,whistomax[i]); |
837 |
trkpad[nn]->SetFillColor(10); |
histocomp[nn][e]->Fill((Float_t)i,whistocomp[i]); |
838 |
|
histofull[nn][e]->Fill((Float_t)i,whistofull[i]); |
839 |
histocomp[nn][e]->GetYaxis()->SetRangeUser(-500,4500); |
if(va1!=0 && ndsp==(int)(va1/100) && i>=128*((va1%100)-1)){ |
840 |
histocomp[nn][e]->SetLineStyle(1); |
if(i<128*(va1%100)) |
841 |
histocomp[nn][e]->SetLineColor(38); |
histova[e]->Fill((Float_t)(i-128*((va1%100)-1)+1),whistofull[i]); |
842 |
histocomp[nn][e]->SetFillColor(38); |
} |
843 |
histocomp[nn][e]->SetLineWidth(1); |
} |
844 |
histocomp[nn][e]->Draw(""); |
} |
845 |
|
|
846 |
histofull[nn][e]->SetLineColor(40); |
if((tot<TOT && e<tot) || tot==TOT){ |
847 |
histofull[nn][e]->SetFillColor(40); |
TBox b; |
848 |
histofull[nn][e]->SetLineWidth(1); |
TLine li,liva1; |
849 |
histofull[nn][e]->SetLineStyle(2); |
li.SetLineColor(38); |
850 |
|
li.SetLineStyle(4); |
851 |
histomax[nn][e]->SetLineColor(2); |
li.SetLineWidth(2); |
852 |
histomax[nn][e]->SetLineWidth(1); |
liva1.SetLineColor(42); |
853 |
histomax[nn][e]->SetLineStyle(3); |
liva1.SetLineStyle(3); |
854 |
|
liva1.SetLineWidth(1); |
855 |
if(ifull[nn]==1) histofull[nn][e]->Draw("9bsame]["); |
|
856 |
if(icomp[nn]==1) histocomp[nn][e]->Draw("9bsame]["); |
Float_t va1x=0,maxhist=4500.,minhist=-500.; |
857 |
if(imax[nn]==1) histomax[nn][e]->Draw("same]["); |
|
858 |
histocomp[nn][e]->Draw("axis same"); |
if(va1!=0){ |
859 |
if(nn==1){ |
cva[e]->cd(); |
860 |
b.SetFillColor(6); |
pad->Draw(); |
861 |
b.SetFillStyle(3945); |
pad->cd(); |
862 |
b.DrawBox(2944.,-500.,3060.,4500.); |
pad->SetFillColor(10); |
863 |
|
histova[e]->SetTitleSize(0.01); |
864 |
b.SetFillColor(107); |
histova[e]->GetYaxis()->SetRangeUser(1500,4500); |
865 |
b.SetFillStyle(3954); |
histova[e]->SetLineColor(40); |
866 |
b.DrawBox(2816.,-500.,2944.,4500.); |
histova[e]->SetFillColor(40); |
867 |
b.DrawBox(2048.,-500.,2176.,4500.); |
histova[e]->SetLineWidth(1); |
868 |
} |
histova[e]->SetLineStyle(2); |
869 |
else if(nn==4){ |
// histova[e]->GetYaxis()->SetLabelSize(0.03); |
870 |
b.SetFillColor(107); |
histova[e]->Draw(""); |
871 |
b.SetFillStyle(3954); |
cva[e]->Update(); |
872 |
b.DrawBox(384.,-500.,512.,4500.); |
} |
873 |
} |
|
874 |
else if(nn==6){ |
c1[e]->cd(); |
875 |
b.SetFillColor(6); |
trkpadtext[nn]->Draw(); |
876 |
b.SetFillStyle(3945); |
trkpad[nn]->Draw(); |
877 |
b.DrawBox(2560.,-500.,2816.,4500.); |
trkpad[nn]->cd(); |
878 |
b.DrawBox(1024.,-500.,1535.,4500.); |
trkpad[nn]->SetFillColor(10); |
879 |
|
|
880 |
b.SetFillColor(107); |
histocomp[nn][e]->SetTitleSize(0.1); |
881 |
b.SetFillStyle(3954); |
histocomp[nn][e]->GetYaxis()->SetRangeUser(-500,4500); |
882 |
b.DrawBox(512.,-500.,768.,4500.); |
histocomp[nn][e]->SetLineStyle(1); |
883 |
b.DrawBox(1536.,-500.,1792.,4500.); |
histocomp[nn][e]->SetLineColor(38); |
884 |
} |
histocomp[nn][e]->SetFillColor(38); |
885 |
else if(nn==7){ |
histocomp[nn][e]->SetLineWidth(1); |
886 |
b.SetFillColor(107); |
histocomp[nn][e]->Draw(""); |
887 |
b.SetFillStyle(3954); |
|
888 |
b.DrawBox(512.,-500.,768.,4500.); |
histofull[nn][e]->SetLineColor(40); |
889 |
} |
histofull[nn][e]->SetFillColor(40); |
890 |
else if(nn==11){ |
histofull[nn][e]->SetLineWidth(1); |
891 |
b.SetFillColor(6); |
histofull[nn][e]->SetLineStyle(2); |
892 |
b.SetFillStyle(3945); |
|
893 |
b.DrawBox(768.,-500.,1024.,4500.); |
histomax[nn][e]->SetLineColor(2); |
894 |
|
histomax[nn][e]->SetLineWidth(1); |
895 |
b.SetFillColor(107); |
histomax[nn][e]->SetLineStyle(3); |
896 |
b.SetFillStyle(3954); |
|
897 |
b.DrawBox(0.,-500.,512.,4500); |
if(ifull[nn]==1) histofull[nn][e]->Draw("9bsame]["); |
898 |
b.DrawBox(1920.,-500.,2048.,4500); |
if(icomp[nn]==1) histocomp[nn][e]->Draw("9bsame]["); |
899 |
|
if(imax[nn]==1) histomax[nn][e]->Draw("same]["); |
900 |
|
histocomp[nn][e]->Draw("axis same"); |
901 |
|
|
902 |
|
b.SetFillColor(19); |
903 |
|
b.SetFillStyle(3945); |
904 |
|
if(nn==0){ |
905 |
|
b.DrawBox(0.,minhist,2047.,maxhist); |
906 |
|
} |
907 |
|
else if(nn==1){ |
908 |
|
b.DrawBox(128.,minhist,256.,maxhist); |
909 |
|
b.DrawBox(384.,minhist,512.,maxhist); |
910 |
|
b.DrawBox(896.,minhist,1024.,maxhist); |
911 |
|
b.DrawBox(2048.,minhist,2432.,maxhist); |
912 |
|
b.DrawBox(2816.,minhist,2944.,maxhist); |
913 |
|
b.DrawBox(2944.,minhist,3070.,maxhist); |
914 |
|
} |
915 |
|
else if(nn==3){ |
916 |
|
b.DrawBox(0.,minhist,256.,maxhist); |
917 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
918 |
|
} |
919 |
|
else if(nn==4){ |
920 |
|
b.DrawBox(256.,minhist,512.,maxhist); |
921 |
|
b.DrawBox(1792.,minhist,1920.,maxhist); |
922 |
|
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){ |
930 |
|
b.DrawBox(512.,minhist,768.,maxhist); |
931 |
|
b.DrawBox(1024.,minhist,1280.,maxhist); |
932 |
|
b.DrawBox(1280.,minhist,1792.,maxhist); |
933 |
|
b.DrawBox(2560.,minhist,2816.,maxhist); |
934 |
|
} |
935 |
|
else if(nn==7){ |
936 |
|
b.DrawBox(0.,minhist,1535.,maxhist); |
937 |
|
b.DrawBox(2024.,minhist,2280.,maxhist); |
938 |
|
} |
939 |
|
else if(nn==8){ |
940 |
|
b.DrawBox(512.,minhist,768.,maxhist); |
941 |
|
} |
942 |
|
else if(nn==9){ |
943 |
|
b.DrawBox(0.,minhist,128.,maxhist); |
944 |
|
b.DrawBox(256.,minhist,384.,maxhist); |
945 |
|
b.DrawBox(512.,minhist,640.,maxhist); |
946 |
|
b.DrawBox(896.,minhist,1152.,maxhist); |
947 |
|
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){ |
952 |
|
b.DrawBox(0.,minhist,512.,maxhist); |
953 |
|
b.DrawBox(1024.,minhist,3070.,maxhist); |
954 |
|
} |
955 |
|
else if(nn==11){ |
956 |
|
b.DrawBox(0.,minhist,512.,maxhist); |
957 |
|
b.DrawBox(768.,minhist,1024.,maxhist); |
958 |
|
b.DrawBox(1536.,minhist,1664.,maxhist); |
959 |
|
b.DrawBox(1920.,minhist,2560.,maxhist); |
960 |
|
b.DrawBox(2816.,minhist,3070.,maxhist); |
961 |
|
} |
962 |
|
for(int va=1; va<24; va++){ |
963 |
|
va1x=128*va; |
964 |
|
liva1.DrawLine(va1x,minhist,va1x,maxhist); |
965 |
|
} |
966 |
|
li.DrawLine(1024.5,minhist,1024.5,maxhist); |
967 |
|
li.DrawLine(2048.5,minhist,2048.5,maxhist); |
968 |
|
c1[e]->Update(); |
969 |
|
} |
970 |
|
|
971 |
} |
} |
|
c1[e]->Update(); |
|
|
|
|
972 |
}//end loop on views |
}//end loop on views |
973 |
|
|
974 |
|
if((tot<TOT && e<tot) || tot==TOT){ |
975 |
|
stringstream nom1,nom2,nom3; |
976 |
|
nom1<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps("; |
977 |
|
nom2<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps"; |
978 |
|
nom3<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps)"; |
979 |
|
|
980 |
|
if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){ |
981 |
|
if(e==0){ |
982 |
|
c1[e]->Print(nom1.str().c_str(),"Landscape"); |
983 |
|
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Landscape"); |
984 |
|
} |
985 |
|
if(e>0 && tot>2){ |
986 |
|
c1[e]->Print(nom2.str().c_str(),"Landscape"); |
987 |
|
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Landscape"); |
988 |
|
} |
989 |
|
if(e==tot-1){ |
990 |
|
c1[e]->Print(nom2.str().c_str(),"Landscape"); |
991 |
|
if(va1!=0) cva[e]->Print(nom3.str().c_str(),"Landscape"); |
992 |
|
} |
993 |
|
} |
994 |
|
else{ |
995 |
|
figsav+="."+outfile; |
996 |
|
c1[e]->Print(figsav.Data()); |
997 |
|
if(va1!=0){ |
998 |
|
figsava+="."+outfile; |
999 |
|
cva[e]->Print(figsava.Data()); |
1000 |
|
} |
1001 |
|
} |
1002 |
|
} |
1003 |
|
if(value%10!=0){ |
1004 |
|
if(totvalues==TOT || (totvalues<TOT && e<totvalues)) |
1005 |
|
values << (short int)0 << endl << endl; |
1006 |
|
} |
1007 |
|
} |
1008 |
|
stringstream com; |
1009 |
|
|
1010 |
c1[e]->Print(figsav.str().c_str()); |
values.close(); |
1011 |
|
if(value==0){ |
1012 |
|
com<<"rm -f "<<fname; |
1013 |
|
system(com.str().c_str()); |
1014 |
} |
} |
1015 |
|
// |
1016 |
|
// Convert ps to pdf if required |
1017 |
|
if(!strcmp(outfile.Data(),"pdf")){ |
1018 |
|
com<<"ps2pdf13 "<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps "<<out<<ffile<<"_FTrkScanQLook_EXPERT.pdf"; |
1019 |
|
system(com.str().c_str()); |
1020 |
|
printf("\n---> ps file converted in pdf format!\n"); |
1021 |
|
com.str(""); |
1022 |
|
com<<"rm -f "<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps "; |
1023 |
|
system(com.str().c_str()); |
1024 |
|
printf("---> ps file removed!\n\n"); |
1025 |
|
com.str(""); |
1026 |
|
} |
1027 |
|
|
1028 |
return; |
return; |
1029 |
} |
} |