1 |
/** |
2 |
* FTrkScanQlook_EXPERT.cxx |
3 |
* |
4 |
* autor: D.Fedele |
5 |
* version v2r08 |
6 |
* Parameters: |
7 |
* file - the path to the root file to analyze |
8 |
* outdir - total path of output file |
9 |
* event - the number of the single event to analyze |
10 |
* 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*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> |
17 |
#include <sstream> |
18 |
#include <fstream> |
19 |
// |
20 |
#include <TPaveText.h> |
21 |
#include <TLatex.h> |
22 |
#include <TCanvas.h> |
23 |
#include <TGraph.h> |
24 |
#include <TFile.h> |
25 |
#include <TTree.h> |
26 |
#include <TStyle.h> |
27 |
#include <TString.h> |
28 |
// |
29 |
#include <physics/tracker/TrackerEvent.h> |
30 |
#include <PscuHeader.h> |
31 |
#include <EventHeader.h> |
32 |
#include <CalibTrk1Event.h> |
33 |
#include <CalibTrk2Event.h> |
34 |
#include <RunHeaderEvent.h> |
35 |
// |
36 |
|
37 |
using namespace std; |
38 |
|
39 |
typedef struct trkword{ |
40 |
int type; |
41 |
int decode; |
42 |
}; |
43 |
|
44 |
/*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-* |
45 |
* |
46 |
| The function "datadecode" decodes the tracker words. |
47 |
* |
48 |
| Tracker words are of three types: |
49 |
* - ADC data |
50 |
| - strip address |
51 |
* - end of ladder |
52 |
| |
53 |
* The function returns a struct variable (trkword) |
54 |
| that contains two quantities, type and decode, which are assigned |
55 |
* the following values: |
56 |
| |
57 |
* type decode |
58 |
| ---------------------------------------------------- |
59 |
* -1 error 0 |
60 |
| 0 data 0-4095 ADC values |
61 |
* 1 address 0-1023 strip address (relative to ladder) |
62 |
| 2 end-of-ladder 1,2,3 ladder number (compressed acq mode) |
63 |
* 4,5,6 ladder number + 3 (full acq mode) |
64 |
| |
65 |
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*/ |
66 |
trkword datadecode(int word){ |
67 |
int type =0; |
68 |
int data =0; |
69 |
int nodata = word&0xf000; |
70 |
int zero = 0; |
71 |
|
72 |
trkword thisword; |
73 |
|
74 |
switch (nodata>>12){ |
75 |
|
76 |
case 0: |
77 |
|
78 |
thisword.decode = word; |
79 |
thisword.type = 0; |
80 |
// cout << thisword.decode << "\n"; |
81 |
return (thisword); //>>>>> ADC data (0) |
82 |
|
83 |
case 1: |
84 |
|
85 |
type = word&0xC00; |
86 |
data = word&0x3ff; |
87 |
|
88 |
switch(type>>10){ |
89 |
|
90 |
case 0: |
91 |
thisword.decode = data; |
92 |
thisword.type = 1; //>>> address (1) |
93 |
return (thisword); |
94 |
|
95 |
case 2: |
96 |
// if(data>=4)data = data-3; |
97 |
if(data>6){ |
98 |
printf("Error on data \n"); |
99 |
thisword.decode = zero; |
100 |
thisword.type = -1; |
101 |
return (thisword); //>>>>> error (-1) |
102 |
} |
103 |
thisword.decode = data; |
104 |
thisword.type = 2; |
105 |
return (thisword); //>>> end-of-ladder |
106 |
|
107 |
default: |
108 |
printf("Error on data \n"); |
109 |
thisword.decode = zero; |
110 |
thisword.type = -1; |
111 |
return (thisword); //>>>>> error (-1) |
112 |
|
113 |
} |
114 |
|
115 |
default: |
116 |
|
117 |
printf("Error on data \n"); |
118 |
thisword.decode = zero; |
119 |
thisword.type = -1; |
120 |
return (thisword); //>>>>> error (-1) |
121 |
} |
122 |
} |
123 |
|
124 |
|
125 |
void FTrkScanQLook_EXPERT(TString file, TString outdir,Int_t event, Int_t va1, Int_t value, TString outfile) |
126 |
{ |
127 |
|
128 |
// |
129 |
// obtain information about the data file and select the output file |
130 |
Int_t dwpos = file.Last('/'); |
131 |
Int_t dwpos1 = file.Last('.'); |
132 |
TString base,ffile ; |
133 |
ffile=file(dwpos+1,dwpos1-(dwpos+1)); |
134 |
if(dwpos>0) base=file(0,dwpos); |
135 |
|
136 |
TString out; |
137 |
if(outdir.Length()==0){ |
138 |
out = base; |
139 |
}else{ |
140 |
out = outdir; |
141 |
} |
142 |
if(out.Last('/')+1<out.Length()) out+="/"; |
143 |
|
144 |
pamela::tracker::TrackerEvent *trk=0; |
145 |
pamela::EventHeader *eh=0,*eH=0; |
146 |
pamela::PscuHeader *ph=0,*pH=0; |
147 |
pamela::RunHeaderEvent *reh=0; |
148 |
pamela::CalibTrk1Event *trk1 = 0; |
149 |
pamela::CalibTrk2Event *trk2 = 0; |
150 |
pamela::EventCounter *cod=0; |
151 |
|
152 |
pamela::PacketType *pctp=0; |
153 |
// open files |
154 |
TFile *trackerFile = new TFile(file); |
155 |
if ( !trackerFile ){ |
156 |
trackerFile->Close(); |
157 |
printf("No tracker file! \n"); |
158 |
return; |
159 |
} |
160 |
|
161 |
//Takes the tree and branches |
162 |
TTree *tr = (TTree*)trackerFile->Get("Physics"); |
163 |
tr->SetBranchAddress("Tracker",&trk); |
164 |
tr->SetBranchAddress("Header",&eh); |
165 |
|
166 |
TTree *otr = (TTree*)trackerFile->Get("RunHeader"); |
167 |
otr->SetBranchAddress("Header",&eH); |
168 |
otr->SetBranchAddress("RunHeader",&reh); |
169 |
|
170 |
TTree *otr1 = (TTree*)trackerFile->Get("CalibTrk1"); |
171 |
otr1->SetBranchAddress("CalibTrk1", &trk1); |
172 |
TTree *otr2 = (TTree*)trackerFile->Get("CalibTrk2"); |
173 |
otr2->SetBranchAddress("CalibTrk2", &trk2); |
174 |
|
175 |
// Define variables |
176 |
Int_t nevents = tr->GetEntries(); |
177 |
Int_t neventH = otr->GetEntries(); |
178 |
Int_t cevents=0; |
179 |
|
180 |
if(otr1->GetEntries()==otr2->GetEntries()) |
181 |
cevents = otr1->GetEntries(); |
182 |
else printf("\nWARNING: CalibTrk1 entries is different from CalibTrk2 entries\n\n"); |
183 |
|
184 |
|
185 |
if ( nevents <= 0 ) { |
186 |
trackerFile->Close(); |
187 |
printf("The file is empty, exiting...\n"); |
188 |
return; |
189 |
} |
190 |
|
191 |
// |
192 |
// open the output text files for the alarms |
193 |
TString fname = out+ffile; |
194 |
int slen=fname.Length(); |
195 |
int last=fname.Last('_'); |
196 |
if(last<(slen-1))fname.Append('_'); |
197 |
fname.Append("FTrk-DSP-values.dat"); |
198 |
ofstream values(fname,ios::out); |
199 |
|
200 |
// |
201 |
printf("\n Number of Entries: %d\n",nevents); |
202 |
printf(" Number of Header Entries: %d\n",neventH); |
203 |
printf(" Number of Calibration: %d\n",cevents); |
204 |
|
205 |
ULong_t obt=0,hobt[neventH]; |
206 |
Int_t eve,cin=0,cins=0,br=0,olrun[neventH]; |
207 |
TString cal=""; |
208 |
|
209 |
eve=3; |
210 |
for(Int_t i=0;i<neventH;i++){ |
211 |
otr->GetEntry(i); |
212 |
pH = eH->GetPscuHeader(); |
213 |
hobt[i]= pH->GetOrbitalTime(); |
214 |
olrun[i]=reh->TRK_CALIB_USED; |
215 |
if(reh->TRK_CALIB_USED!=104 && br==0){ |
216 |
obt = pH->GetOrbitalTime(); |
217 |
cal="Event with online calibration"; |
218 |
br=1; |
219 |
} |
220 |
if(i==neventH-1 && br==0){ |
221 |
cal="***** ONLINE CALIBRATION NOT FOUND IN THIS FILE *****"; |
222 |
eve=2; |
223 |
} |
224 |
} |
225 |
|
226 |
int tot=2,totvalues=0,TOT=0,totsum=0; |
227 |
if(event<0) tot=abs(event); |
228 |
if(value>0) totvalues=value/10; |
229 |
|
230 |
TOT= tot > totvalues ? tot : totvalues; |
231 |
if(value>0) totsum= (value/10) > 100 ? (value/10) : 100; |
232 |
else totsum=100; |
233 |
|
234 |
Int_t calev=cevents; |
235 |
|
236 |
TH1F *histomax[12][tot]; //histos of max signals |
237 |
TH1F *histocomp[12][tot]; //histos of compressed data |
238 |
TH1F *histofull[12][tot]; //histos of full data |
239 |
TCanvas *c1[tot],*csum[calev]; |
240 |
TH1F *histova[tot]; |
241 |
TGraph*sum[12][calev][totsum/2]; |
242 |
TCanvas *cva[tot]; |
243 |
|
244 |
TPad *trkpad[12]; //pad for histos |
245 |
|
246 |
stringstream fromfile,isfile,title,hid,message; |
247 |
TString figsa="",figsav="",figsava=""; |
248 |
Int_t canvasx=1200; |
249 |
Int_t canvasy=900; |
250 |
|
251 |
TLatex *t=new TLatex(); |
252 |
t->SetTextFont(32); |
253 |
t->SetTextColor(1); |
254 |
t->SetTextAlign(12); |
255 |
t->SetTextSize(0.02); |
256 |
|
257 |
gStyle->SetLabelSize(0.06,"x"); |
258 |
gStyle->SetLabelSize(0.06,"y"); |
259 |
gStyle->SetTitleFontSize(0.1); |
260 |
gStyle->SetFillColor(10); |
261 |
gStyle->SetTitleFillColor(10); |
262 |
gStyle->SetOptStat(0); |
263 |
|
264 |
Int_t p[12],psum=0; |
265 |
Int_t whistostrip[3072]; |
266 |
for(int s=0;s<3072;s++) whistostrip[s]=s+1; |
267 |
|
268 |
if(eve==3){ |
269 |
for(Int_t i=0;i<nevents;i++){ |
270 |
tr->GetEntry(i); |
271 |
ph = eh->GetPscuHeader(); |
272 |
cod = eh->GetCounter(); |
273 |
if(olrun[cod->Get(pctp->RunHeader)-1]==104) continue; |
274 |
if(value%10!=0){ |
275 |
if(i==0) cin=cod->Get(pctp->CalibTrk1)+((value%10)-1); |
276 |
if(cod->Get(pctp->CalibTrk1)==cin+1){ |
277 |
eve=i+3; |
278 |
cin=100; |
279 |
} |
280 |
} |
281 |
else { |
282 |
if(i==0) cin=cod->Get(pctp->CalibTrk1); |
283 |
if(event<0 && cod->Get(pctp->CalibTrk1)==cin+1){ |
284 |
eve=i+3; |
285 |
cin=100; |
286 |
} |
287 |
else if(event>=0 && (ULong_t)ph->GetOrbitalTime()>obt && obt>0){ |
288 |
eve=i+3; |
289 |
obt=0; |
290 |
} |
291 |
} |
292 |
if(i==0) cins=cod->Get(pctp->CalibTrk1); |
293 |
if(cod->Get(pctp->CalibTrk1)>cins){ |
294 |
psum=1; |
295 |
cins=cod->Get(pctp->CalibTrk1); |
296 |
} |
297 |
if(psum==1){ |
298 |
figsa=out+ffile+"_FTrkScanQLook_EXPERT_sumof"; |
299 |
figsa+=totsum; |
300 |
figsa+="events_cal"; |
301 |
figsa+=cins; |
302 |
csum[cins-1] = new TCanvas(figsa.Data(),"FTrkQLookSCAN",canvasx,canvasy); |
303 |
csum[cins-1]->SetFillColor(10); |
304 |
csum[cins-1]->Range(0,0,1,1); |
305 |
fromfile.str(""); |
306 |
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Sum of "<<totsum<<" events after the "<<cins<<" calibration at OBT= "<<hobt[(cod->Get(pctp->RunHeader))-1]; |
307 |
t->SetTextSize(0.02); |
308 |
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
309 |
|
310 |
// draw pads |
311 |
Double_t posy = 0.95; // up y-coord - top pads |
312 |
Double_t hpad = 0.15; // pad height |
313 |
Double_t posx1=0; // left x-coord - pad column |
314 |
Double_t posx0=0; // x-coord - column division |
315 |
Double_t wrel = 0.6; // relative x size of first sub-column |
316 |
Double_t marg = 0.004; // margin among pads |
317 |
|
318 |
hpad = (posy-marg*11)/6; |
319 |
wrel = (1-marg*4)/2; |
320 |
for(Int_t n = 0; n<12; n++){ |
321 |
if ( (n+1)%2 ) { |
322 |
if(n>1)posy = posy-(marg*2+hpad); |
323 |
posx1 = marg; |
324 |
posx0 = posx1 + wrel; |
325 |
} else { |
326 |
posx1 = posx0 + 2*marg; |
327 |
posx0 = posx1 + wrel; |
328 |
}; |
329 |
|
330 |
/* -----------> pad for histograms */ |
331 |
trkpad[n] = new TPad("pad"," ",posx1,posy-hpad,posx0,posy,18,0,0); |
332 |
trkpad[n]->SetFillColor(19); |
333 |
trkpad[n]->SetFrameFillColor(10); |
334 |
|
335 |
p[n]=0; |
336 |
} |
337 |
} |
338 |
if(psum>0 && psum<=totsum){ |
339 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
340 |
// create header dump retrieving event info |
341 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
342 |
Int_t ndsp=0; |
343 |
|
344 |
Int_t whistomax[3072]; |
345 |
Int_t whisto[3072]; |
346 |
Int_t whistocomp[3072]; |
347 |
Int_t whistofull[3072]; |
348 |
|
349 |
// // transmitted words |
350 |
Int_t word = 0; |
351 |
Int_t iword = 0; |
352 |
Int_t ii=0,ifull[12],icomp[12],imax[12],nn=0; |
353 |
trkword thisword; |
354 |
|
355 |
Int_t address,ladder; |
356 |
|
357 |
|
358 |
for(Int_t n = 0; n<12; n++){ |
359 |
ndsp = trk->DSPnumber[n]; |
360 |
nn = ndsp-1; |
361 |
ifull[nn]=0; |
362 |
icomp[nn]=0; |
363 |
imax[nn]=0; |
364 |
|
365 |
for(Int_t vi = 0; vi< 3072; vi++){ |
366 |
whistomax[vi] = -200; |
367 |
whistocomp[vi] = -200; |
368 |
whistofull[vi] = -200; |
369 |
whisto[vi] = -200; |
370 |
} |
371 |
|
372 |
// // trasmitted data |
373 |
// |
374 |
address = 0; |
375 |
ladder = 1; |
376 |
for(Int_t vi = 0; vi < trk->DATAlength[n] ; vi++){ |
377 |
word = trk->TrackerData.At(iword); |
378 |
thisword = datadecode(word); |
379 |
iword++; |
380 |
switch (thisword.type){ |
381 |
|
382 |
case 0: //ADC value |
383 |
whisto[address] = thisword.decode; |
384 |
address++; |
385 |
// cout << " adr " << address << "\n"; |
386 |
break; |
387 |
|
388 |
case 1: //address |
389 |
address = 1024*(ladder-1) + thisword.decode; |
390 |
// cout << " adr " << address << "\n"; |
391 |
break; |
392 |
|
393 |
case 2: //end-of-ladder |
394 |
ladder = thisword.decode; |
395 |
// cout << "Ladder " << ladder << "\n"; |
396 |
if(ladder==3){ |
397 |
// end of compressed data - FILL HISTO |
398 |
//cout << ">>> COMPRESSED data" << "\n"; |
399 |
for(ii = 0; ii < 3072; ii++){ |
400 |
whistocomp[ii]=whisto[ii]; |
401 |
whisto[ii] = -200; |
402 |
} |
403 |
address = 0; |
404 |
ladder = 1; |
405 |
}else if(ladder==6){ |
406 |
// end of full data - FILL HISTO |
407 |
//cout << ">>> FULL data" << "\n"; |
408 |
for(ii = 0; ii < 3072; ii++){ |
409 |
whistofull[ii]=whisto[ii]; |
410 |
whisto[ii] = -200; |
411 |
} |
412 |
address = 0; |
413 |
ladder = 1; |
414 |
}else{ |
415 |
if(ladder>3) ladder=ladder-3; |
416 |
address= ladder*1024; |
417 |
ladder = ladder + 1; |
418 |
} |
419 |
} |
420 |
} |
421 |
|
422 |
for(Int_t vi = 0; vi < 3072; vi++){ |
423 |
if(whistofull[vi]>-200){ |
424 |
if(vi==3071){ |
425 |
sum[nn][cins-1][p[nn]] = new TGraph(3072,whistostrip,whistofull); |
426 |
p[nn]++; |
427 |
} |
428 |
if(value%10==0 && value>0){ |
429 |
if(whistofull[vi]>0){ |
430 |
if(vi==0) values << (short int)ndsp << " "; |
431 |
values <<(short int) whistofull[vi] <<" "; |
432 |
if(vi==3071) values << endl; |
433 |
} |
434 |
} |
435 |
} |
436 |
} |
437 |
} |
438 |
psum++; |
439 |
if(value%10==0 && value>0) values << (short int)0 << endl << endl; |
440 |
} |
441 |
if(psum==totsum+1){ |
442 |
if(value%10==0 && value>0) values << (short int)(-cins) << endl << endl; |
443 |
psum=0; |
444 |
for(Int_t vi = 0; vi < 12; vi++){ |
445 |
if(p[vi]>0){ |
446 |
|
447 |
TLine li,liva1; |
448 |
li.SetLineColor(38); |
449 |
li.SetLineStyle(4); |
450 |
li.SetLineWidth(2); |
451 |
liva1.SetLineColor(42); |
452 |
liva1.SetLineStyle(3); |
453 |
liva1.SetLineWidth(1); |
454 |
|
455 |
Float_t va1x=0; |
456 |
csum[cins-1]->cd(); |
457 |
trkpad[vi]->Draw(); |
458 |
trkpad[vi]->cd(); |
459 |
trkpad[vi]->SetFillColor(10); |
460 |
|
461 |
stringstream tit,hid; |
462 |
tit.str(""); |
463 |
hid.str(""); |
464 |
tit<<"DSP "<<vi+1; |
465 |
sum[vi][cins-1][0]->SetTitle(tit.str().c_str()); |
466 |
sum[vi][cins-1][0]->GetXaxis()->SetTitle("channel id"); |
467 |
sum[vi][cins-1][0]->GetXaxis()->CenterTitle(); |
468 |
sum[vi][cins-1][0]->GetXaxis()->SetRangeUser(0,3073); |
469 |
sum[vi][cins-1][0]->GetXaxis()->SetTitleSize(0.06); |
470 |
sum[vi][cins-1][0]->GetXaxis()->SetTitleOffset(0.8); |
471 |
sum[vi][cins-1][0]->GetYaxis()->SetTitle("ADC value"); |
472 |
sum[vi][cins-1][0]->GetYaxis()->SetTitleSize(0.09); |
473 |
sum[vi][cins-1][0]->GetYaxis()->SetTitleOffset(0.4); |
474 |
sum[vi][cins-1][0]->GetYaxis()->CenterTitle(); |
475 |
sum[vi][cins-1][0]->GetYaxis()->SetRangeUser(0,4500); |
476 |
sum[vi][cins-1][0]->SetMarkerStyle(20); |
477 |
sum[vi][cins-1][0]->SetMarkerSize(0.1); |
478 |
sum[vi][cins-1][0]->SetMarkerColor(4); |
479 |
sum[vi][cins-1][0]->Draw("ap"); |
480 |
hid<<p[vi]<<" events"; |
481 |
t->SetTextSize(0.11); |
482 |
t->DrawLatex(2400,4250,hid.str().c_str()); |
483 |
for(int pi=1; pi<p[vi];pi++){ |
484 |
sum[vi][cins-1][pi]->SetMarkerStyle(21); |
485 |
sum[vi][cins-1][pi]->SetMarkerSize(0.1); |
486 |
sum[vi][cins-1][pi]->SetMarkerColor(4); |
487 |
sum[vi][cins-1][pi]->Draw("p"); |
488 |
} |
489 |
for(int va=1; va<24; va++){ |
490 |
va1x=128*va; |
491 |
liva1.DrawLine(va1x,0.,va1x,4500.); |
492 |
} |
493 |
li.DrawLine(1024.5,0.,1024.5,4500.); |
494 |
li.DrawLine(2048.5,0.,2048.5,4500.); |
495 |
csum[cins-1]->Update(); |
496 |
hid.str(""); |
497 |
} |
498 |
} |
499 |
if(strcmp(outfile.Data(),"ps") && strcmp(outfile.Data(),"pdf")){ |
500 |
figsa+="."+outfile; |
501 |
csum[cins-1]->Print(figsa.Data()); |
502 |
} |
503 |
if(cod->Get(pctp->CalibTrk1)==cevents) break; |
504 |
} |
505 |
} |
506 |
} |
507 |
|
508 |
if(!strcmp(outfile.Data(),"ps") || !strcmp(outfile.Data(),"pdf")){ |
509 |
stringstream nom1,nom2,nom3; |
510 |
nom1<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps("; |
511 |
nom2<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps"; |
512 |
nom3<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps)"; |
513 |
|
514 |
for(int va=0; va<cevents; va++){ |
515 |
if(cevents==1) csum[va]->Print(nom2.str().c_str(),"Landscape"); |
516 |
if(cevents>2){ |
517 |
if(va==0) csum[va]->Print(nom1.str().c_str(),"Landscape"); |
518 |
if(va>0 && va<cevents-1) csum[va]->Print(nom2.str().c_str(),"Landscape"); |
519 |
if(va>0 && va==cevents-1) csum[va]->Print(nom3.str().c_str(),"Landscape"); |
520 |
} |
521 |
} |
522 |
} |
523 |
// |
524 |
// Convert ps to pdf if required |
525 |
if(!strcmp(outfile.Data(),"pdf")){ |
526 |
stringstream com; |
527 |
com<<"ps2pdf13 "<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps "<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.pdf"; |
528 |
system(com.str().c_str()); |
529 |
printf("\n---> ps file converted in pdf format!\n"); |
530 |
com.str(""); |
531 |
com<<"rm -f "<<out<<ffile<<"_FTrkScanQLook_EXPERT_sumof"<<totsum<<"event.ps "; |
532 |
system(com.str().c_str()); |
533 |
printf("---> ps file removed!\n\n"); |
534 |
com.str(""); |
535 |
} |
536 |
|
537 |
|
538 |
|
539 |
for(Int_t e=0;e<TOT;e++){ |
540 |
if(event<=0 || (value%10!=0 && e==0)) |
541 |
event=eve; |
542 |
else { |
543 |
event=event+1; |
544 |
if(event>eve-3 && eve>2) |
545 |
cal="Event with online calibration"; |
546 |
else |
547 |
cal="***** ONLINE CALIBRATION NOT FOUND IN THIS FILE *****"; |
548 |
} |
549 |
printf("Scan of Entry %d\n",event); |
550 |
|
551 |
tr->GetEntry(event); |
552 |
//============================================================================ |
553 |
|
554 |
// draw display area |
555 |
|
556 |
TPad *trkpad[12],*pad=0; //pad for histos |
557 |
TPaveText *trkpadtext[12]; //pad for header |
558 |
Int_t ndsp=0; |
559 |
|
560 |
if((tot<TOT && e<tot) || tot==TOT){ |
561 |
figsav=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
562 |
figsav+=event+1; |
563 |
c1[e] = new TCanvas(figsav.Data(),"FTrkQLookSCAN",canvasx,canvasy); |
564 |
c1[e]->SetFillColor(10); |
565 |
c1[e]->Range(0,0,1,1); |
566 |
fromfile.str(""); |
567 |
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event; |
568 |
t->SetTextSize(0.02); |
569 |
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
570 |
t->DrawLatex(0.60,0.98,cal.Data()); |
571 |
|
572 |
isfile<<"Boxes so colored point out chips broken before the launch and never resumed"; |
573 |
t->SetTextColor(6); |
574 |
t->SetTextSize(0.018); |
575 |
t->DrawLatex(.05,.96,isfile.str().c_str()); |
576 |
isfile.str(""); |
577 |
|
578 |
isfile<<"Boxes so colored point out anomalous chips (not necessarily broken)"; |
579 |
t->SetTextColor(107); |
580 |
t->DrawLatex(.5,.96,isfile.str().c_str()); |
581 |
isfile.str(""); |
582 |
t->SetTextColor(1); |
583 |
|
584 |
if(va1!=0){ |
585 |
figsava=out+ffile+"_FTrkScanQLook_EXPERT_ev"; |
586 |
figsava+=event+1; |
587 |
figsava+="_DSP"; |
588 |
figsava+=(int)(va1/100); |
589 |
figsava+="_VA1-"; |
590 |
figsava+=va1%100; |
591 |
cva[e] = new TCanvas(figsava.Data(),"TrkQLookSCAN VA1",canvasx,canvasy); |
592 |
cva[e]->SetFillColor(10); |
593 |
cva[e]->Range(0,0,1,1); |
594 |
fromfile.str(""); |
595 |
fromfile<<"FTrkScanQLook_EXPERT File: "<<ffile<<" ----> Entry "<<event<<" --> DSP "<<(int)(va1/100)<<" --> va1 "<<va1%100; |
596 |
t->DrawLatex(0.02,0.98,fromfile.str().c_str()); |
597 |
// t->DrawLatex(0.65,0.98,cal.Data()); |
598 |
} |
599 |
|
600 |
// draw pads |
601 |
Double_t posy = 0.95; // up y-coord - top pads |
602 |
Double_t hpad = 0.15; // pad height |
603 |
Double_t posx1=0; // left x-coord - pad column |
604 |
Double_t posx2=0; // right x-coord - pad olumn |
605 |
Double_t posx0=0; // x-coord - column division |
606 |
Double_t wrel = 0.6; // relative x size of first sub-column |
607 |
Double_t marg = 0.004; // margin among pads |
608 |
|
609 |
for(Int_t n = 0; n<12; n++){ |
610 |
if ( (n+1)%2 ) { |
611 |
if(n>1)posy = posy-(marg*2+hpad); |
612 |
posx1 = marg; |
613 |
posx2 = 0.5 - marg; |
614 |
posx0 = 0.5*wrel; |
615 |
} else { |
616 |
posx1 = posx1 + 0.5; |
617 |
posx2 = posx2 + 0.5; |
618 |
posx0 = posx0 + 0.5; |
619 |
}; |
620 |
|
621 |
/* -----------> pad for histograms */ |
622 |
trkpad[n] = new TPad("pad"," ",posx1,posy-hpad,posx0-marg,posy,18,0,0); |
623 |
trkpad[n]->SetFillColor(19); |
624 |
trkpad[n]->SetFrameFillColor(10); |
625 |
/* -----------> pad for header dump */ |
626 |
trkpadtext[n] = new TPaveText((posx0+marg),(posy-hpad),posx2,posy); |
627 |
/* -----------> HISTOGRAMS */ |
628 |
|
629 |
title.str(""); |
630 |
title<<"DSP "<<n+1; |
631 |
hid<<"h"<<n+e*100; |
632 |
histomax[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
633 |
hid<<"hh"<<n+e*100; |
634 |
histocomp[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
635 |
hid<<"hhh"<<n+e*100; |
636 |
histofull[n][e] = new TH1F(hid.str().c_str(),title.str().c_str(),3073,-0.5,3072.5); |
637 |
title.str(""); |
638 |
hid.str(""); |
639 |
} |
640 |
if(va1!=0){ |
641 |
title.str(""); |
642 |
title<<"DSP "<<(int)(va1/100)<<" -- va1 "<<va1%100; |
643 |
hid<<"va"<<e; |
644 |
histova[e] = new TH1F(hid.str().c_str(),title.str().c_str(),130,0,129); |
645 |
} //end loop on views |
646 |
|
647 |
/* -----------> pad for histograms */ |
648 |
pad = new TPad("padva"," ",0,0,1,0.97,18,0,0); |
649 |
pad->SetFillColor(19); |
650 |
pad->SetFrameFillColor(10); |
651 |
} |
652 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
653 |
// create header dump retrieving event info |
654 |
// = = = = = = = = = = = = = = = = = = = = = = = = = |
655 |
|
656 |
|
657 |
Double_t whistomax[3072]; |
658 |
Double_t whisto[3072]; |
659 |
Double_t whistocomp[3072]; |
660 |
Double_t whistofull[3072]; |
661 |
|
662 |
//============================================= |
663 |
// transmitted words |
664 |
Int_t word = 0; |
665 |
Int_t iword = 0; |
666 |
Int_t TOTDATAlength_check = 0; |
667 |
Int_t ii=0,ifull[12],icomp[12],imax[12],nn=0; |
668 |
trkword thisword; |
669 |
|
670 |
Int_t address,ladder; |
671 |
|
672 |
|
673 |
for(Int_t n = 0; n<12; n++){ |
674 |
|
675 |
ndsp = trk->DSPnumber[n]; |
676 |
nn = ndsp-1; |
677 |
ifull[nn]=0; |
678 |
icomp[nn]=0; |
679 |
imax[nn]=0; |
680 |
if(ndsp>0){ |
681 |
if(ndsp<13){ |
682 |
if((tot<TOT && e<tot) || tot==TOT){ |
683 |
|
684 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
685 |
* |
686 |
* Write event LEVEL0 report |
687 |
* |
688 |
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*/ |
689 |
|
690 |
trkpadtext[nn]->SetTextFont(40); |
691 |
trkpadtext[nn]->SetFillColor(33); |
692 |
trkpadtext[nn]->SetTextSize(0.012); |
693 |
trkpadtext[nn]->SetTextAlign(13); |
694 |
|
695 |
trkpadtext[nn]->AddText(" "); |
696 |
message<<"DAQ mode --------> "<<trk->DAQmode[n]; |
697 |
trkpadtext[nn]->AddText(message.str().c_str()); |
698 |
message.str(""); |
699 |
message<<"Event number --------> "<<trk->eventn[n]; |
700 |
trkpadtext[nn]->AddText(message.str().c_str()); |
701 |
message.str(""); |
702 |
message<<"13-bit words --------> "<<trk->DATAlength[n]; |
703 |
trkpadtext[nn]->AddText(message.str().c_str()); |
704 |
message.str(""); |
705 |
if (!(nn%2)&&trk->signcluster[n][0]!=0) message<<"L1 add: "<<trk->addrcluster[n][0]<<" - sign: "<<1024-(trk->signcluster[n][0]); |
706 |
else message<<"L1 add: "<<trk->addrcluster[n][0]<<" - sign: "<<(trk->signcluster[n][0]); |
707 |
trkpadtext[nn]->AddText(message.str().c_str()); |
708 |
message.str(""); |
709 |
if (!(nn%2)&&trk->signcluster[n][1]!=0) message<<"L2 add: "<<trk->addrcluster[n][1]<<" - sign: "<<1024-(trk->signcluster[n][1]); |
710 |
else message<<"L2 add: "<<trk->addrcluster[n][1]<<" - sign: "<<trk->signcluster[n][1]; |
711 |
trkpadtext[nn]->AddText(message.str().c_str()); |
712 |
message.str(""); |
713 |
if (!(nn%2)&&trk->signcluster[n][2]!=0) message<<"L3 add: "<<trk->addrcluster[n][2]<<" - sign: "<<1024-(trk->signcluster[n][2]); |
714 |
else message<<"L3 add: "<<trk->addrcluster[n][2]<<" - sign: "<<trk->signcluster[n][2]; |
715 |
trkpadtext[nn]->AddText(message.str().c_str()); |
716 |
message.str(""); |
717 |
message<<"NCLUST "<<trk->nclust[n]<<" CUTC "<<trk->cutc[n]<<" CUTCL "<<trk->cutcl[n]; |
718 |
trkpadtext[nn]->AddText(message.str().c_str()); |
719 |
message.str(""); |
720 |
message<<"Comp. time "<<trk->compressiontime[n]<<" x 0.051ms = "<<0.051*trk->compressiontime[n]<<" ms"; |
721 |
trkpadtext[nn]->AddText(message.str().c_str()); |
722 |
message.str(""); |
723 |
trkpadtext[nn]->AddText(" "); |
724 |
message<<"CRC --> "<<trk->crc[n]; |
725 |
trkpadtext[nn]->AddText(message.str().c_str()); |
726 |
message.str(""); |
727 |
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]; |
728 |
trkpadtext[nn]->AddText(message.str().c_str()); |
729 |
message.str(""); |
730 |
trkpadtext[nn]->AddText(" "); |
731 |
trkpadtext[nn]->AddLine(0,0,0,0); |
732 |
message<<"PNum "<<trk->pnum[n]<<" - BId "<<trk->bid[n]<<" ALARM "<<trk->alarm[n]; |
733 |
trkpadtext[nn]->AddText(message.str().c_str()); |
734 |
message.str(""); |
735 |
message<<"Cmd "<<trk->cmdnum[n]<<" --- Answer length "<<trk->aswr[n]<<" byte"; |
736 |
trkpadtext[nn]->AddText(message.str().c_str()); |
737 |
message.str(""); |
738 |
trkpadtext[nn]->AddText(" "); |
739 |
|
740 |
TOTDATAlength_check = TOTDATAlength_check + trk->DATAlength[n]; |
741 |
} |
742 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
743 |
* |
744 |
* Plot event LEVEL0 histo |
745 |
* |
746 |
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*/ |
747 |
|
748 |
//============================================= |
749 |
} |
750 |
for(Int_t i = 0; i< 3072; i++){ |
751 |
whistomax[i] = -200; |
752 |
whistocomp[i] = -200; |
753 |
whistofull[i] = -200; |
754 |
whisto[i] = -200; |
755 |
} |
756 |
|
757 |
// =============== |
758 |
// trasmitted data |
759 |
// =============== |
760 |
|
761 |
address = 0; |
762 |
ladder = 1; |
763 |
for(Int_t i = 0; i < trk->DATAlength[n] ; i++){ |
764 |
word = trk->TrackerData.At(iword); |
765 |
thisword = datadecode(word); |
766 |
iword++; |
767 |
switch (thisword.type){ |
768 |
|
769 |
case 0: //ADC value |
770 |
whisto[address] = thisword.decode; |
771 |
address++; |
772 |
// cout << " adr " << address << "\n"; |
773 |
break; |
774 |
|
775 |
case 1: //address |
776 |
address = 1024*(ladder-1) + thisword.decode; |
777 |
// cout << " adr " << address << "\n"; |
778 |
break; |
779 |
|
780 |
case 2: //end-of-ladder |
781 |
ladder = thisword.decode; |
782 |
// cout << "Ladder " << ladder << "\n"; |
783 |
if(ladder==3){ |
784 |
// end of compressed data - FILL HISTO |
785 |
//cout << ">>> COMPRESSED data" << "\n"; |
786 |
for(ii = 0; ii < 3072; ii++){ |
787 |
whistocomp[ii]=whisto[ii]; |
788 |
whisto[ii] = -200; |
789 |
} |
790 |
address = 0; |
791 |
ladder = 1; |
792 |
}else if(ladder==6){ |
793 |
// end of full data - FILL HISTO |
794 |
//cout << ">>> FULL data" << "\n"; |
795 |
for(ii = 0; ii < 3072; ii++){ |
796 |
whistofull[ii]=whisto[ii]; |
797 |
whisto[ii] = -200; |
798 |
} |
799 |
address = 0; |
800 |
ladder = 1; |
801 |
}else{ |
802 |
if(ladder>3) ladder=ladder-3; |
803 |
address= ladder*1024; |
804 |
ladder = ladder + 1; |
805 |
} |
806 |
} |
807 |
} |
808 |
|
809 |
|
810 |
// =============== |
811 |
// maximum signals |
812 |
// =============== |
813 |
if(trk->signcluster[nn][0]!=0) whistomax[ 0+(int)trk->addrcluster[nn][0]] = whistocomp[ 0+(int)trk->addrcluster[nn][0]]; |
814 |
if(trk->signcluster[nn][1]!=0) whistomax[1024+(int)trk->addrcluster[nn][1]] = whistocomp[1024+(int)trk->addrcluster[nn][1]]; |
815 |
if(trk->signcluster[nn][2]!=0) whistomax[2048+(int)trk->addrcluster[nn][2]] = whistocomp[2048+(int)trk->addrcluster[nn][2]]; |
816 |
|
817 |
for(Int_t i = 0; i < 3072; i++){ |
818 |
if(whistomax[i]>-200) imax[nn]=1; |
819 |
if(whistocomp[i]>-200) icomp[nn]=1; |
820 |
if(whistofull[i]>-200) ifull[nn]=1; |
821 |
if(value%10!=0){ |
822 |
if((totvalues<TOT && e<totvalues) || totvalues==TOT){ |
823 |
if(whistofull[i]>0){ |
824 |
if(i==0) values << (short int)ndsp << " "; |
825 |
values <<(short int) whistofull[i] <<" "; |
826 |
if(i==3071) values << endl; |
827 |
} |
828 |
} |
829 |
} |
830 |
if((tot<TOT && e<tot) || tot==TOT){ |
831 |
histomax[nn][e]->Fill((Float_t)i,whistomax[i]); |
832 |
histocomp[nn][e]->Fill((Float_t)i,whistocomp[i]); |
833 |
histofull[nn][e]->Fill((Float_t)i,whistofull[i]); |
834 |
if(va1!=0 && ndsp==(int)(va1/100) && i>=128*((va1%100)-1)){ |
835 |
if(i<128*(va1%100)) |
836 |
histova[e]->Fill((Float_t)(i-128*((va1%100)-1)+1),whistofull[i]); |
837 |
} |
838 |
} |
839 |
} |
840 |
|
841 |
if((tot<TOT && e<tot) || tot==TOT){ |
842 |
TBox b; |
843 |
TLine li,liva1; |
844 |
li.SetLineColor(38); |
845 |
li.SetLineStyle(4); |
846 |
li.SetLineWidth(2); |
847 |
liva1.SetLineColor(42); |
848 |
liva1.SetLineStyle(3); |
849 |
liva1.SetLineWidth(1); |
850 |
|
851 |
Float_t va1x=0; |
852 |
|
853 |
if(va1!=0){ |
854 |
cva[e]->cd(); |
855 |
pad->Draw(); |
856 |
pad->cd(); |
857 |
pad->SetFillColor(10); |
858 |
histova[e]->SetTitleSize(0.01); |
859 |
histova[e]->GetYaxis()->SetRangeUser(1500,4500); |
860 |
histova[e]->SetLineColor(40); |
861 |
histova[e]->SetFillColor(40); |
862 |
histova[e]->SetLineWidth(1); |
863 |
histova[e]->SetLineStyle(2); |
864 |
// histova[e]->GetYaxis()->SetLabelSize(0.03); |
865 |
histova[e]->Draw(""); |
866 |
cva[e]->Update(); |
867 |
} |
868 |
|
869 |
c1[e]->cd(); |
870 |
trkpadtext[nn]->Draw(); |
871 |
trkpad[nn]->Draw(); |
872 |
trkpad[nn]->cd(); |
873 |
trkpad[nn]->SetFillColor(10); |
874 |
|
875 |
histocomp[nn][e]->SetTitleSize(0.1); |
876 |
histocomp[nn][e]->GetYaxis()->SetRangeUser(-500,4500); |
877 |
histocomp[nn][e]->SetLineStyle(1); |
878 |
histocomp[nn][e]->SetLineColor(38); |
879 |
histocomp[nn][e]->SetFillColor(38); |
880 |
histocomp[nn][e]->SetLineWidth(1); |
881 |
histocomp[nn][e]->Draw(""); |
882 |
|
883 |
histofull[nn][e]->SetLineColor(40); |
884 |
histofull[nn][e]->SetFillColor(40); |
885 |
histofull[nn][e]->SetLineWidth(1); |
886 |
histofull[nn][e]->SetLineStyle(2); |
887 |
|
888 |
histomax[nn][e]->SetLineColor(2); |
889 |
histomax[nn][e]->SetLineWidth(1); |
890 |
histomax[nn][e]->SetLineStyle(3); |
891 |
|
892 |
if(ifull[nn]==1) histofull[nn][e]->Draw("9bsame]["); |
893 |
if(icomp[nn]==1) histocomp[nn][e]->Draw("9bsame]["); |
894 |
if(imax[nn]==1) histomax[nn][e]->Draw("same]["); |
895 |
histocomp[nn][e]->Draw("axis same"); |
896 |
if(nn==0){ |
897 |
b.SetFillColor(107); |
898 |
b.SetFillStyle(3945); |
899 |
b.DrawBox(256.,-500.,384.,4500); |
900 |
b.DrawBox(768.,-500.,2047.,4500.); |
901 |
} |
902 |
else if(nn==1){ |
903 |
b.SetFillColor(6); |
904 |
b.SetFillStyle(3945); |
905 |
b.DrawBox(2944.,-500.,3060.,4500.); |
906 |
|
907 |
b.SetFillColor(107); |
908 |
b.SetFillStyle(3954); |
909 |
b.DrawBox(384.,-500.,512.,4500.); |
910 |
b.DrawBox(896.,-500.,1024.,4500.); |
911 |
b.DrawBox(2048.,-500.,2432.,4500.); |
912 |
b.DrawBox(2816.,-500.,2944.,4500.); |
913 |
} |
914 |
else if(nn==3){ |
915 |
b.SetFillColor(107); |
916 |
b.SetFillStyle(3954); |
917 |
b.DrawBox(2816.,-500.,3070.,4500.); |
918 |
} |
919 |
else if(nn==4){ |
920 |
b.SetFillColor(107); |
921 |
b.SetFillStyle(3954); |
922 |
b.DrawBox(256.,-500.,512.,4500.); |
923 |
b.DrawBox(2816.,-500.,3070.,4500.); |
924 |
} |
925 |
else if(nn==5){ |
926 |
b.SetFillColor(107); |
927 |
b.SetFillStyle(3954); |
928 |
b.DrawBox(896.,-500.,1024.,4500.); |
929 |
// b.DrawBox(1664.,-500.,1792.,4500.); |
930 |
} |
931 |
else if(nn==6){ |
932 |
b.SetFillColor(107); |
933 |
b.SetFillStyle(3954); |
934 |
b.DrawBox(512.,-500.,768.,4500.); |
935 |
b.DrawBox(1024.,-500.,1280.,4500.); |
936 |
b.DrawBox(1280.,-500.,1792.,4500.); |
937 |
b.DrawBox(2560.,-500.,2816.,4500.); |
938 |
} |
939 |
else if(nn==7){ |
940 |
b.SetFillColor(107); |
941 |
b.SetFillStyle(3954); |
942 |
b.DrawBox(512.,-500.,768.,4500.); |
943 |
b.DrawBox(896.,-500.,1535.,4500.); |
944 |
b.DrawBox(2024.,-500.,2280.,4500.); |
945 |
} |
946 |
else if(nn==8){ |
947 |
b.SetFillColor(107); |
948 |
b.SetFillStyle(3954); |
949 |
b.DrawBox(512.,-500.,768.,4500.); |
950 |
} |
951 |
else if(nn==9){ |
952 |
b.SetFillColor(107); |
953 |
b.SetFillStyle(3954); |
954 |
b.DrawBox(0.,-500.,128.,4500.); |
955 |
b.DrawBox(256.,-500.,384.,4500.); |
956 |
b.DrawBox(512.,-500.,640.,4500.); |
957 |
b.DrawBox(896.,-500.,1024.,4500.); |
958 |
b.DrawBox(1280.,-500.,1535.,4500.); |
959 |
b.DrawBox(1664.,-500.,1920.,4500.); |
960 |
} |
961 |
else if(nn==10){ |
962 |
b.SetFillColor(107); |
963 |
b.SetFillStyle(3954); |
964 |
b.DrawBox(0.,-500.,512.,4500.); |
965 |
b.DrawBox(1024.,-500.,3070.,4500.); |
966 |
} |
967 |
else if(nn==11){ |
968 |
b.SetFillColor(6); |
969 |
b.SetFillStyle(3945); |
970 |
b.DrawBox(768.,-500.,1024.,4500.); |
971 |
|
972 |
b.SetFillColor(107); |
973 |
b.SetFillStyle(3954); |
974 |
b.DrawBox(0.,-500.,512.,4500.); |
975 |
b.DrawBox(1920.,-500.,2560.,4500.); |
976 |
b.DrawBox(2816.,-500.,3070.,4500.); |
977 |
} |
978 |
for(int va=1; va<24; va++){ |
979 |
va1x=128*va; |
980 |
liva1.DrawLine(va1x,-500.,va1x,4500.); |
981 |
} |
982 |
li.DrawLine(1024.5,-500.,1024.5,4500.); |
983 |
li.DrawLine(2048.5,-500.,2048.5,4500.); |
984 |
c1[e]->Update(); |
985 |
} |
986 |
|
987 |
} |
988 |
}//end loop on views |
989 |
|
990 |
if((tot<TOT && e<tot) || tot==TOT){ |
991 |
stringstream nom1,nom2,nom3; |
992 |
nom1<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps("; |
993 |
nom2<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps"; |
994 |
nom3<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps)"; |
995 |
|
996 |
if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){ |
997 |
if(e==0){ |
998 |
c1[e]->Print(nom1.str().c_str(),"Landscape"); |
999 |
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Landscape"); |
1000 |
} |
1001 |
if(e>0 && tot>2){ |
1002 |
c1[e]->Print(nom2.str().c_str(),"Landscape"); |
1003 |
if(va1!=0) cva[e]->Print(nom2.str().c_str(),"Landscape"); |
1004 |
} |
1005 |
if(e==tot-1){ |
1006 |
c1[e]->Print(nom2.str().c_str(),"Landscape"); |
1007 |
if(va1!=0) cva[e]->Print(nom3.str().c_str(),"Landscape"); |
1008 |
} |
1009 |
} |
1010 |
else{ |
1011 |
figsav+="."+outfile; |
1012 |
c1[e]->Print(figsav.Data()); |
1013 |
if(va1!=0){ |
1014 |
figsava+="."+outfile; |
1015 |
cva[e]->Print(figsava.Data()); |
1016 |
} |
1017 |
} |
1018 |
} |
1019 |
if(value%10!=0){ |
1020 |
if(totvalues==TOT || (totvalues<TOT && e<totvalues)) |
1021 |
values << (short int)0 << endl << endl; |
1022 |
} |
1023 |
} |
1024 |
stringstream com; |
1025 |
|
1026 |
values.close(); |
1027 |
if(value==0){ |
1028 |
com<<"rm -f "<<fname; |
1029 |
system(com.str().c_str()); |
1030 |
} |
1031 |
// |
1032 |
// Convert ps to pdf if required |
1033 |
if(!strcmp(outfile.Data(),"pdf")){ |
1034 |
com<<"ps2pdf13 "<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps "<<out<<ffile<<"_FTrkScanQLook_EXPERT.pdf"; |
1035 |
system(com.str().c_str()); |
1036 |
printf("\n---> ps file converted in pdf format!\n"); |
1037 |
com.str(""); |
1038 |
com<<"rm -f "<<out<<ffile<<"_FTrkScanQLook_EXPERT.ps "; |
1039 |
system(com.str().c_str()); |
1040 |
printf("---> ps file removed!\n\n"); |
1041 |
com.str(""); |
1042 |
} |
1043 |
|
1044 |
return; |
1045 |
} |