1 |
/** |
2 |
* TmtcTemperature |
3 |
* author Marcelli |
4 |
* version 1.1 - 10 December 2005 |
5 |
* |
6 |
* Version 1.2 - 28 December 2004 |
7 |
* If TFile does not exist the function exit to prompt |
8 |
* If nevents = 0 the function exit to promt |
9 |
* Initialize to 0 alla array dedicated to graphs |
10 |
* |
11 |
* version 1.3 - 22 February 2005 - Nagni |
12 |
* For compatibility with batch mode excution: |
13 |
* 1) Added "include <iostream>" and "using namespace std" |
14 |
* 2) Removed gROOT->Reset() |
15 |
* |
16 |
* version 1.4 - 25 February 2005 - Marcelli |
17 |
* Canvas of thermistors divided in 12 pads |
18 |
* |
19 |
* Version 1.5 |
20 |
* Date 08 March 2005 - Nagni |
21 |
* Added "format" parameter for saving the produced image in various formats |
22 |
* (for a complete list of types refer to TPad::SaveAs method) |
23 |
* |
24 |
* version 1.6 |
25 |
* Date 24 May 2005 - Marcelli |
26 |
* Added four thermistor (files after DW_050516_007 has 16 thermistors instead of 12) |
27 |
* |
28 |
* version 1.7 |
29 |
* Date 02 February 2006 - Marcelli |
30 |
* Modified to work with new version of YODA |
31 |
* |
32 |
* version 1.8 |
33 |
* Date March 2006 - Marcelli |
34 |
* For compilation: |
35 |
* Added function "int main(int argc, char* argv[])" |
36 |
* |
37 |
* Description : To monitor information from TMTC Packet (voltage, CC,....) |
38 |
* |
39 |
* Parameters: |
40 |
* TString base - the path to the root directory for the specific Pamela unpack session |
41 |
* TString outDir - the path where to save the output image (Default = base) |
42 |
* TString format - the format which will be used fo rsave the produced images (Default = "jpg") |
43 |
* |
44 |
*/ |
45 |
|
46 |
|
47 |
|
48 |
|
49 |
#include <iostream> |
50 |
#include <sstream> |
51 |
#include "TStyle.h" |
52 |
#include "TFile.h" |
53 |
#include "TAxis.h" |
54 |
#include "TLegend.h" |
55 |
#include "TF1.h" |
56 |
#include "TList.h" |
57 |
#include "TTree.h" |
58 |
#include "TObjString.h" |
59 |
#include "TCanvas.h" |
60 |
#include "TGraph.h" |
61 |
#include "TMultiGraph.h" |
62 |
#include "TH1F.h" |
63 |
#include "TGaxis.h" |
64 |
#include "TString.h" |
65 |
#include "TPaveText.h" |
66 |
#include "tmtc/TmtcEvent.h" |
67 |
#include "tmtc/TmtcRecord.h" |
68 |
|
69 |
|
70 |
using namespace std; |
71 |
|
72 |
/////////////////////////////////////////////////////////////////////////////////////////////// |
73 |
/////////////////////////////////////FUNZIONI DI CONVERSIONE/////////////////////////////////// |
74 |
/////////////////////////////////////////////////////////////////////////////////////////////// |
75 |
|
76 |
|
77 |
//reference: Laben Document, TL 18348, pag. 79 |
78 |
float convert_th(int TH) { |
79 |
|
80 |
float a,q,deltax,deltay; |
81 |
float gradi[111],grado=-1; |
82 |
int adc[111]={4095, 3868, 3654, 3454, 3265, 3088, 2922, 2765, 2618, 2479, 2348, 2225, 2110, 2000, 1897, 1800, 1709, 1622, 1541, 1464, 1391, 1322, 1257, 1196, 1138, 1083, 1031, 981, 935, 890, 849, 809, 771, 736, 702, 670, 639, 610, 583, 557, 532, 508, 486, 464, 444, 425, 406, 389, 372, 356, 341, 327, 313, 300, 288, 276, 264, 254, 243, 233, 224, 215, 206, 198, 190, 183, 175, 169, 162, 156, 150, 144, 138, 133, 128, 123, 118, 114, 110, 106, 102, 98, 94, 91, 87, 84, 81, 78, 75, 73, 70, 67, 65, 63, 60, 58, 56, 54, 52, 50, 49, 47, 45, 44, 42, 41, 39, 38, 37, 35, 34}; |
83 |
int maxpos=-1,minpos=-1,i; |
84 |
|
85 |
for (i=0;i<111;i++) |
86 |
{ |
87 |
gradi[i]=-35+i; |
88 |
} |
89 |
|
90 |
if (TH > 4095) |
91 |
{ |
92 |
grado=100.; |
93 |
return grado; |
94 |
} |
95 |
|
96 |
if (TH <= 34) |
97 |
{ |
98 |
grado=100.; |
99 |
return grado; |
100 |
} |
101 |
|
102 |
|
103 |
if (TH==4095) |
104 |
{ |
105 |
grado=-35.; |
106 |
return grado; |
107 |
} |
108 |
else |
109 |
{ |
110 |
|
111 |
for (i=1;i<111;i++) |
112 |
{ |
113 |
if (TH>=adc[i]) |
114 |
{ |
115 |
minpos=i; |
116 |
maxpos=i-1; |
117 |
break; |
118 |
} |
119 |
} |
120 |
deltax=adc[maxpos]-adc[minpos]; |
121 |
deltay=gradi[maxpos]-gradi[minpos]; |
122 |
a=deltay/deltax; |
123 |
q=gradi[maxpos]-a*adc[maxpos]; |
124 |
grado=a*TH+q; |
125 |
} |
126 |
return (grado); |
127 |
} |
128 |
|
129 |
|
130 |
|
131 |
//reference: Laben Document, TL 18348, pag. 76 |
132 |
float convert_volt(int in) { |
133 |
|
134 |
float a,q,deltax,deltay; |
135 |
float inputvoltage[52],voltage=-1; |
136 |
float adc[52]; |
137 |
int maxpos=-1,minpos=-1,i; |
138 |
|
139 |
for (i=0;i<52;i++) |
140 |
{ |
141 |
inputvoltage[i] = (i*0.10); |
142 |
} |
143 |
for (i=0;i<52;i++) |
144 |
{ |
145 |
adc[i] = (i*5); |
146 |
} |
147 |
|
148 |
if (in==0) |
149 |
{ |
150 |
voltage= 0.00; |
151 |
} |
152 |
else |
153 |
{ |
154 |
|
155 |
for (i=0;i<52;i++) |
156 |
{ |
157 |
if (in<=adc[i]) |
158 |
{ |
159 |
minpos=i; |
160 |
maxpos=i-1; |
161 |
break; |
162 |
} |
163 |
} |
164 |
deltax=adc[maxpos]-adc[minpos]; |
165 |
deltay=inputvoltage[maxpos]-inputvoltage[minpos]; |
166 |
a=deltay/deltax; |
167 |
q=inputvoltage[maxpos]-a*adc[maxpos]; |
168 |
voltage=a*in+q; |
169 |
} |
170 |
return (voltage); |
171 |
|
172 |
} |
173 |
|
174 |
|
175 |
//////////////////////////////////////////////////////////////////////////////// |
176 |
////////////////////////// TMTC TEMPERATURE ////////////////////////////////// |
177 |
////////////////////////////////////////////////////////////////////////////////// |
178 |
|
179 |
void TmtcTemperature(TString base, TString outDir, TString format){ |
180 |
|
181 |
//------- load root file -------------- |
182 |
TFile *file = new TFile(base.Data()); |
183 |
|
184 |
if ( outDir == "" ) outDir = "."; |
185 |
|
186 |
if (!file){ |
187 |
printf("No such file in the directory has been found"); |
188 |
return; |
189 |
} |
190 |
|
191 |
//--- Takes the tree of Tmtc ----------- |
192 |
TPaveText *pt=0; |
193 |
TPaveText *pt1=0; |
194 |
|
195 |
TTree *tr = (TTree*)file->Get("Tmtc"); |
196 |
TBranch *tmtcBr = tr->GetBranch("Tmtc"); |
197 |
|
198 |
pamela::TmtcEvent *tme=0; |
199 |
pamela::TmtcRecord *tmr=0; |
200 |
|
201 |
tr->SetBranchAddress("Tmtc", &tme); |
202 |
|
203 |
//////--------to analize data-----------------////// |
204 |
Long64_t nevents = tmtcBr->GetEntries(); |
205 |
if (nevents<=0) { |
206 |
file->Close(); |
207 |
printf("nevents = %i", nevents); |
208 |
return; |
209 |
} |
210 |
|
211 |
|
212 |
string titolocc[] ={"IPM1","IPM2","IPM3","IPM4","IPM5","IPM6","KHB_HOT","KHB_COLD","IDAQ_HOT","IDAQ_COLD","VCB_STANDBY","VRL_HOT","VRL_COLD","PSB","TOFHV_HOT","TOFHV_COLD"}; |
213 |
|
214 |
string titoloth[] ={"IPM_TH","S1_TH","S4_TH","S4_ND_PLATE_TH","TRK_TH1","TRK_TH2","FLUID_IN_TH","FLUID_OUT_TH","VME_TH","DCDC_TH","CPU_TH1","CPU_TH2", "IPM_CPU_TH", "VRL_TH1", "VRL_TH2", "VME_COOL_TH"}; |
215 |
|
216 |
|
217 |
const Int_t size = nevents; |
218 |
std::stringstream oss; |
219 |
Int_t recordstot=0; |
220 |
Int_t recordstotmax=0; |
221 |
Int_t records[size]; |
222 |
Int_t ev[size]; |
223 |
Double_t OBT[size]; |
224 |
Int_t tot=0; |
225 |
Int_t m=0; |
226 |
Int_t l=0; |
227 |
|
228 |
|
229 |
for (Int_t i = 0; i < size; i++){ ///ciclo per decidere le pagine |
230 |
tmtcBr->GetEntry(i); |
231 |
Long64_t tmpSize = tme->Records->GetEntries(); |
232 |
recordstot=recordstot+tmpSize; |
233 |
if (tmpSize>0){ |
234 |
tmr = (pamela::TmtcRecord*)tme->Records->At(0); |
235 |
if(i==0){ |
236 |
ev[0]=0; |
237 |
records[0]=0; |
238 |
OBT[0]=tmr->TM_RECORD_OBT; |
239 |
}else if (fmod(i, 50)==0 || ((tmr->TM_RECORD_OBT)<OBT[i-1]) || (i==(size-1)) ){ |
240 |
//}else if ((fmod(i, 20)==0) || (i==(size-1)) ){ |
241 |
tot=tot+1; |
242 |
if (recordstotmax<recordstot) recordstotmax=recordstot; |
243 |
records[tot]=recordstot; |
244 |
OBT[tot]=tmr->TM_RECORD_OBT; |
245 |
ev[tot]=i; |
246 |
recordstot=0; |
247 |
} |
248 |
} |
249 |
} |
250 |
|
251 |
|
252 |
|
253 |
const Int_t lungmax=16*recordstotmax; |
254 |
const Int_t lungmin=6*recordstotmax; |
255 |
|
256 |
|
257 |
Double_t xrecordobtcc[lungmax], yccdiagacq[lungmax], ythana[lungmax], xrecordobtth[lungmax], xrecordobtdea[lungmin], ydea[lungmin], xrecordobtcc_1[lungmax], xrecordobtcc_1d[lungmax], xrecordobtcc_1u[lungmax], yccdiagacq_1[lungmax], yccdiagacq_1d[lungmax], yccdiagacq_1u[lungmax], ythana_1[lungmax], xrecordobtth_1[lungmax], xrecordobtdea_1[lungmin], ydea_1[lungmin], limth[lungmax], limvolt1[lungmin], limvolt2[lungmin], limvolt3[lungmin]; |
258 |
|
259 |
Double_t yccdiagacq_11u[lungmax], yccdiagacq_11d[lungmax], yccipm12u[lungmax], yccipm12d[lungmax], yccipmkk1u[lungmax], yccipmkk1d[lungmax], yccipmerror[lungmax], yccipmerror2[lungmax]; |
260 |
|
261 |
|
262 |
|
263 |
/////ciclo sulle pagine//////////////////// |
264 |
for (Int_t interval=0; interval<tot; interval++){ |
265 |
|
266 |
///-------to create canvas--------------------// |
267 |
TCanvas *Canvascc = new TCanvas("Tmtc_1", base, 2500, 2000); |
268 |
Canvascc->Divide(1,15); |
269 |
Canvascc->SetFillColor(10); |
270 |
|
271 |
TCanvas *Canvasthdea = new TCanvas("Tmtc_4", base, 1280, 1024); |
272 |
Canvasthdea->SetFillColor(10); |
273 |
Canvasthdea->Divide(4,4); |
274 |
|
275 |
|
276 |
TCanvas *Canvasvoltdea = new TCanvas("Tmtc_5", base, 1280, 1024); |
277 |
Canvasvoltdea->SetFillColor(10); |
278 |
Canvasvoltdea->Divide(2,3); |
279 |
|
280 |
|
281 |
recordstot=0; |
282 |
|
283 |
for (Int_t i = ev[interval]; i < ev[interval+1]; i++){ //ciclo su un sottogruppo |
284 |
// cout<<"ciclo da "<<ev[interval]<<" a "<< ev[interval+1]<<"\n"; |
285 |
tmtcBr->GetEntry(i); |
286 |
l=0; |
287 |
m=0; |
288 |
Long64_t tmpSize = tme->Records->GetEntries(); |
289 |
Int_t size_b = tmpSize; |
290 |
for (Int_t j = 0; j < size_b; j++){ |
291 |
tmr = (pamela::TmtcRecord*)tme->Records->At(j); |
292 |
|
293 |
for (Int_t k =0; k <16; k++){ |
294 |
yccdiagacq[16*recordstot+16*j+k] = (((tmr->TM_DIAG_AND_BILEVEL_ACQ)>>(15-k))&0x0001); |
295 |
xrecordobtcc[16*recordstot+16*j+k] = tmr->TM_RECORD_OBT; |
296 |
} |
297 |
|
298 |
for (Int_t k =0; k <16; k++){ |
299 |
ythana[16*recordstot+16*j+k] = convert_th(tmr->TM_TH_ANA[k]); |
300 |
xrecordobtth[16*recordstot+16*j+k] = tmr->TM_RECORD_OBT; |
301 |
} |
302 |
|
303 |
for (Int_t k =0; k <6; k++){ |
304 |
ydea[6*recordstot+6*j+k] = convert_volt(tmr->TM_DEA_ANA[k]); |
305 |
xrecordobtdea[6*recordstot+6*j+k] = tmr->TM_RECORD_OBT; |
306 |
} |
307 |
} |
308 |
recordstot=recordstot+tmpSize; |
309 |
} |
310 |
|
311 |
TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString(); |
312 |
filename = ((TObjString*)filename.Tokenize('.')->First())->GetString(); |
313 |
|
314 |
TLegend *leg1 = new TLegend(0.87,0.80,0.98,0.94, ""); |
315 |
leg1->SetFillColor(10); |
316 |
leg1->SetBorderSize(0); |
317 |
|
318 |
Canvascc->cd(); |
319 |
pt = new TPaveText (.87,.95,.98,.98); |
320 |
oss.str(""); |
321 |
oss<<filename.Data(); |
322 |
pt->AddText(oss.str().c_str()); |
323 |
pt->SetTextColor(1); |
324 |
//pt->SetTextSize(1); |
325 |
pt->SetFillColor(10); |
326 |
pt->SetBorderSize(0); |
327 |
pt->Draw(); |
328 |
|
329 |
/* |
330 |
TLine *line = new TLine; |
331 |
line->SetLineColor(4); |
332 |
line->SetLineWidth(2); |
333 |
line->DrawLine(0, 0.60, 1, 0.60); |
334 |
line->Draw(); |
335 |
line->DrawLine(0, 0.33, 1, 0.33); |
336 |
line->Draw(); |
337 |
line->DrawLine(0, 0.20, 1, 0.20); |
338 |
line->Draw(); |
339 |
line->DrawLine(0, 0.13, 1, 0.13); |
340 |
line->Draw(); |
341 |
*/ |
342 |
|
343 |
|
344 |
|
345 |
//CC Graph |
346 |
for (Int_t k =0; k<16; k++){ |
347 |
for (Int_t i = 0; i < recordstot; i++){ |
348 |
if (k>0) { |
349 |
yccdiagacq_11u[i]=yccdiagacq_1u[i]; |
350 |
yccdiagacq_11d[i]=yccdiagacq_1d[i]; |
351 |
} |
352 |
yccdiagacq_1d[i]=0; |
353 |
yccdiagacq_1u[i]=0; |
354 |
if (yccdiagacq[16*i+k]==0) { |
355 |
yccdiagacq_1d[i]=-1; |
356 |
xrecordobtcc_1[i]= xrecordobtcc[16*i+k]; |
357 |
} |
358 |
else if (yccdiagacq[16*i+k]==1){ |
359 |
yccdiagacq_1u[i]=1; |
360 |
xrecordobtcc_1[i]= xrecordobtcc[16*i+k]; |
361 |
} |
362 |
if (k==1) { |
363 |
if (yccdiagacq_11u[i]==1 && yccdiagacq_1u[i]==1) yccipm12u[i]=1; |
364 |
else yccipm12u[i]=0; |
365 |
if (yccdiagacq_11d[i]==-1 && yccdiagacq_1d[i]==-1) yccipm12d[i]=-1; |
366 |
else yccipm12d[i]=0; |
367 |
} |
368 |
if (k==1 || k==3 || k==5 || k==7 || k==9) { |
369 |
if (yccdiagacq_11u[i]==1 && yccdiagacq_1u[i]==1) yccipmkk1u[i]=1; |
370 |
else yccipmkk1u[i]=0; |
371 |
if (yccdiagacq_11d[i]==-1 && yccdiagacq_1d[i]==-1) yccipmkk1d[i]=-1; |
372 |
else yccipmkk1d[i]=0; |
373 |
if (k==7 || k==9) { |
374 |
if (yccipmkk1u[i]==1 && yccipm12d[i]==0 && yccipm12u[i]==0) yccipmerror[i]=1; |
375 |
else if (yccipmkk1d[i]==-1 && yccipm12d[i]==0 && yccipm12u[i]==0 || (yccipmkk1d[i]==-1 && yccipm12u[i]==1)) yccipmerror[i]=-1; |
376 |
else yccipmerror[i]=0; |
377 |
} |
378 |
} |
379 |
if (k==6 || k==7 || k==8 || k==9) { |
380 |
if (yccdiagacq_1d[i]==-1 && yccipm12u[i]==1) yccipmerror2[i]=-1; |
381 |
else yccipmerror2[i]=0; |
382 |
} |
383 |
if (k==13) { |
384 |
if (yccdiagacq_1d[i]==-1 && yccipm12d[i]==0 && yccipm12u[i]==0) yccdiagacq_1d[i]==-1; //psb on |
385 |
if (yccdiagacq_1u[i]==1 && yccipm12u[i]==1) yccdiagacq_1u[i]==1; //psb off |
386 |
if (yccdiagacq_1d[i]==-1 && yccipm12d[i]==-1) yccipmerror[i]=-1; |
387 |
else if (yccdiagacq_1u[i]==1 && yccipm12d[i]==0 && yccipm12u[i]==0) yccipmerror[i]=1; |
388 |
else yccipmerror[i]=0; //psb alarm |
389 |
} |
390 |
if (k==15) { |
391 |
yccdiagacq_11u[i]=0; |
392 |
yccdiagacq_11d[i]=0; |
393 |
} |
394 |
} |
395 |
|
396 |
//CC IPM |
397 |
if (k==0 || k==1 ||k==2 || k==3 || k==4 || k==5 ){ |
398 |
Canvascc->cd(m+1); |
399 |
Canvascc->cd(m+1)->SetGridx(); |
400 |
gPad->SetFillColor(10); |
401 |
gPad->SetFrameBorderSize(0); |
402 |
gPad->SetLineColor(1); |
403 |
gPad->SetFrameBorderMode(0); |
404 |
|
405 |
|
406 |
pt = new TPaveText (.01,.30,.10,.70); |
407 |
oss.str(""); |
408 |
oss<<titolocc[k]; |
409 |
pt->AddText(oss.str().c_str()); |
410 |
pt->SetFillColor(10); |
411 |
pt->SetBorderSize(0); |
412 |
pt->Draw(); |
413 |
|
414 |
|
415 |
TGraph *ccrecordd = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1d); |
416 |
ccrecordd->SetTitle(""); |
417 |
ccrecordd->SetFillColor(8); |
418 |
ccrecordd->GetXaxis()->SetLabelSize(0); |
419 |
ccrecordd->GetXaxis()->SetTickLength(0); |
420 |
ccrecordd->GetYaxis()->SetLabelSize(0); |
421 |
ccrecordd->GetYaxis()->SetTickLength(0); |
422 |
ccrecordd->SetMaximum(1.3); |
423 |
ccrecordd->SetMinimum(-1.3); |
424 |
ccrecordd->GetXaxis()->SetAxisColor(0); |
425 |
ccrecordd->Draw("AB"); |
426 |
if (k==0) leg1->AddEntry(ccrecordd,"ON","f"); |
427 |
|
428 |
TGraph *ccrecordu = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1u); |
429 |
ccrecordu->SetTitle(""); |
430 |
ccrecordu->SetFillColor(17); |
431 |
ccrecordu->GetXaxis()->SetLabelSize(0); |
432 |
ccrecordu->GetXaxis()->SetTickLength(0); |
433 |
ccrecordu->GetYaxis()->SetLabelSize(0); |
434 |
ccrecordu->GetYaxis()->SetTickLength(0); |
435 |
ccrecordu->GetXaxis()->SetAxisColor(0); |
436 |
ccrecordu->SetMaximum(1.3); |
437 |
ccrecordu->SetMinimum(-1.3); |
438 |
ccrecordu->Draw("Bsame"); |
439 |
if (k==0) leg1->AddEntry(ccrecordu,"OFF","f"); |
440 |
|
441 |
|
442 |
if (k==1 || k==3 || k==5){ |
443 |
Canvascc->cd(m); |
444 |
TGraph *ccrecorderrord = new TGraph(recordstot, xrecordobtcc_1, yccipmkk1d); |
445 |
ccrecorderrord->SetTitle(""); |
446 |
ccrecorderrord->SetFillColor(2); |
447 |
ccrecorderrord->GetXaxis()->SetLabelSize(0); |
448 |
ccrecorderrord->GetXaxis()->SetTickLength(0); |
449 |
ccrecorderrord->GetYaxis()->SetLabelSize(0); |
450 |
ccrecorderrord->GetYaxis()->SetTickLength(0); |
451 |
ccrecorderrord->SetMaximum(1.3); |
452 |
ccrecorderrord->SetMinimum(-1.3); |
453 |
ccrecorderrord->GetXaxis()->SetAxisColor(0); |
454 |
ccrecorderrord->Draw("Bsame"); |
455 |
|
456 |
TGraph *ccrecorderroru = new TGraph(recordstot, xrecordobtcc_1, yccipmkk1u); |
457 |
ccrecorderroru->SetTitle(""); |
458 |
ccrecorderroru->SetFillColor(2); |
459 |
ccrecorderroru->GetXaxis()->SetLabelSize(0); |
460 |
ccrecorderroru->GetXaxis()->SetTickLength(0); |
461 |
ccrecorderroru->GetYaxis()->SetLabelSize(0); |
462 |
ccrecorderroru->GetYaxis()->SetTickLength(0); |
463 |
ccrecorderroru->SetMaximum(1.3); |
464 |
ccrecorderroru->SetMinimum(-1.3); |
465 |
ccrecorderroru->GetXaxis()->SetAxisColor(0); |
466 |
ccrecorderroru->Draw("Bsame"); |
467 |
|
468 |
Canvascc->cd(m+1); |
469 |
ccrecorderrord->Draw("Bsame"); |
470 |
ccrecorderroru->Draw("Bsame"); |
471 |
} |
472 |
|
473 |
|
474 |
Double_t xmin; |
475 |
Double_t xmax; |
476 |
xmin=ccrecordu->GetXaxis()->GetXmin(); |
477 |
xmax=ccrecordu->GetXaxis()->GetXmax(); |
478 |
TGaxis *axis = new TGaxis (xmin, 0, xmax, 0, xmin, xmax, 510, "+-", 0); |
479 |
axis->Draw("same"); |
480 |
/*TF1 *func1 = new TF1("func1", "0"); |
481 |
func1->SetRange(xmin,xrecordobtcc_1[recordstot-1]); |
482 |
func1->SetLineColor(1); |
483 |
func1->SetLineStyle(1); |
484 |
func1->SetLineWidth(2); |
485 |
func1->Draw("same");*/ |
486 |
|
487 |
m=m+1; |
488 |
} |
489 |
|
490 |
//CC idaq-khb-tof |
491 |
if (k==6 || k==7 || k==8 || k==9 || k==14 || k==15){ |
492 |
Canvascc->cd(l+7); |
493 |
gPad->SetFillColor(10); |
494 |
gPad->SetFrameBorderSize(0); |
495 |
|
496 |
pt = new TPaveText (.00,.30,.09,.70); |
497 |
oss.str(""); |
498 |
oss<<titolocc[k]; |
499 |
pt->AddText(oss.str().c_str()); |
500 |
pt->SetFillColor(10); |
501 |
pt->SetBorderSize(0); |
502 |
pt->Draw(); |
503 |
|
504 |
|
505 |
TGraph *ccrecordd = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1d); |
506 |
ccrecordd->SetTitle(""); |
507 |
if (k==14 || k==15) ccrecordd->SetFillColor(2); |
508 |
else ccrecordd->SetFillColor(8); |
509 |
ccrecordd->GetXaxis()->SetLabelSize(0); |
510 |
ccrecordd->GetXaxis()->SetTickLength(0); |
511 |
ccrecordd->GetYaxis()->SetLabelSize(0); |
512 |
ccrecordd->GetYaxis()->SetTickLength(0); |
513 |
ccrecordd->SetMaximum(1.3); |
514 |
ccrecordd->SetMinimum(-1.3); |
515 |
ccrecordd->GetXaxis()->SetAxisColor(0); |
516 |
ccrecordd->Draw("AB"); |
517 |
|
518 |
TGraph *ccrecordu = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1u); |
519 |
ccrecordu->SetTitle(""); |
520 |
if (k==14 || k==15) { |
521 |
ccrecordu->SetFillColor(7); |
522 |
if (k==14) leg1->AddEntry(ccrecordu,"ON-OFF-LATCHUP","f"); |
523 |
} |
524 |
else ccrecordu->SetFillColor(17); |
525 |
ccrecordu->GetXaxis()->SetAxisColor(0); |
526 |
ccrecordu->GetXaxis()->SetLabelSize(0); |
527 |
ccrecordu->GetXaxis()->SetTickLength(0); |
528 |
ccrecordu->GetYaxis()->SetLabelSize(0); |
529 |
ccrecordu->GetYaxis()->SetTickLength(0); |
530 |
ccrecordu->SetMaximum(1.3); |
531 |
ccrecordu->SetMinimum(-1.3); |
532 |
ccrecordu->Draw("Bsame"); |
533 |
|
534 |
if (k==6 || k==7 || k==8 || k==9){ |
535 |
Canvascc->cd(l+7); |
536 |
TGraph *ccrecorderroru = new TGraph(recordstot, xrecordobtcc_1, yccipmerror2); |
537 |
ccrecorderroru->SetTitle(""); |
538 |
ccrecorderroru->SetFillColor(2); |
539 |
ccrecorderroru->GetXaxis()->SetLabelSize(0); |
540 |
ccrecorderroru->GetXaxis()->SetTickLength(0); |
541 |
ccrecorderroru->GetYaxis()->SetLabelSize(0); |
542 |
ccrecorderroru->GetYaxis()->SetTickLength(0); |
543 |
ccrecorderroru->SetMaximum(1.3); |
544 |
ccrecorderroru->SetMinimum(-1.3); |
545 |
ccrecorderroru->GetXaxis()->SetAxisColor(0); |
546 |
ccrecorderroru->Draw("Bsame"); |
547 |
} |
548 |
|
549 |
|
550 |
if (k==7 || k==9){ |
551 |
Canvascc->cd(l+7-1); |
552 |
TGraph *ccrecorderroru = new TGraph(recordstot, xrecordobtcc_1, yccipmerror); |
553 |
ccrecorderroru->SetTitle(""); |
554 |
ccrecorderroru->SetFillColor(2); |
555 |
ccrecorderroru->GetXaxis()->SetLabelSize(0); |
556 |
ccrecorderroru->GetXaxis()->SetTickLength(0); |
557 |
ccrecorderroru->GetYaxis()->SetLabelSize(0); |
558 |
ccrecorderroru->GetYaxis()->SetTickLength(0); |
559 |
ccrecorderroru->SetMaximum(1.3); |
560 |
ccrecorderroru->SetMinimum(-1.3); |
561 |
ccrecorderroru->GetXaxis()->SetAxisColor(0); |
562 |
if (k==7) leg1->AddEntry(ccrecorderroru,"ALARM","f"); |
563 |
ccrecorderroru->Draw("Bsame"); |
564 |
|
565 |
Canvascc->cd(l+7); |
566 |
ccrecorderroru->Draw("Bsame"); |
567 |
} |
568 |
|
569 |
|
570 |
Double_t xmin; |
571 |
Double_t xmax; |
572 |
xmin=ccrecordu->GetXaxis()->GetXmin(); |
573 |
xmax=ccrecordu->GetXaxis()->GetXmax(); |
574 |
TGaxis *axis = new TGaxis (xmin, 0, xmax, 0, xmin, xmax, 510, "+-", 1.5); |
575 |
axis->Draw("same"); |
576 |
/*TF1 *func1 = new TF1("func1", "0"); |
577 |
func1->SetRange(xmin,xrecordobtcc_1[recordstot-1]); |
578 |
func1->SetLineColor(1); |
579 |
func1->SetLineStyle(1); |
580 |
func1->SetLineWidth(2); |
581 |
func1->Draw("same");*/ |
582 |
l=l+1; |
583 |
|
584 |
} |
585 |
|
586 |
//CC vrl |
587 |
if (k==11 || k==12 || k==13 ){ |
588 |
if (k==13) Canvascc->cd(13); |
589 |
if (k==11) Canvascc->cd(14); |
590 |
if (k==12) Canvascc->cd(15); |
591 |
gPad->SetFillColor(10); |
592 |
gPad->SetFrameBorderSize(0); |
593 |
|
594 |
pt = new TPaveText (.01,.30,.09,.70); |
595 |
oss.str(""); |
596 |
oss<<titolocc[k]; |
597 |
pt->AddText(oss.str().c_str()); |
598 |
pt->SetFillColor(10); |
599 |
pt->SetBorderSize(0); |
600 |
pt->Draw(); |
601 |
|
602 |
if (k==11 || k==12){ |
603 |
TGraph *ccrecordd = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1d); |
604 |
ccrecordd->SetTitle(""); |
605 |
ccrecordd->SetFillColor(8); |
606 |
ccrecordd->GetXaxis()->SetAxisColor(0); |
607 |
ccrecordd->GetXaxis()->SetLabelSize(0); |
608 |
ccrecordd->GetXaxis()->SetTickLength(0); |
609 |
ccrecordd->GetYaxis()->SetLabelSize(0); |
610 |
ccrecordd->GetYaxis()->SetTickLength(0); |
611 |
ccrecordd->SetMaximum(1.3); |
612 |
ccrecordd->SetMinimum(-1.3); |
613 |
ccrecordd->Draw("AB"); |
614 |
|
615 |
|
616 |
TGraph *ccrecordu = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1u); |
617 |
ccrecordu->SetTitle(""); |
618 |
ccrecordu->SetFillColor(17); |
619 |
ccrecordu->GetXaxis()->SetAxisColor(0); |
620 |
ccrecordu->GetXaxis()->SetLabelSize(0); |
621 |
ccrecordu->GetXaxis()->SetTickLength(0); |
622 |
ccrecordu->GetYaxis()->SetLabelSize(0); |
623 |
ccrecordu->GetYaxis()->SetTickLength(0); |
624 |
ccrecordu->SetMaximum(1.3); |
625 |
ccrecordu->SetMinimum(-1.3); |
626 |
ccrecordu->Draw("Bsame"); |
627 |
|
628 |
Double_t xmin; |
629 |
Double_t xmax; |
630 |
xmin=ccrecordu->GetXaxis()->GetXmin(); |
631 |
xmax=ccrecordu->GetXaxis()->GetXmax(); |
632 |
if (k==12){ |
633 |
TGaxis *axis = new TGaxis (xmin, -1, xmax, -1, xmin, xmax, 510, "", 0); |
634 |
axis->SetLabelSize(0.24); |
635 |
axis->SetLabelColor(1); |
636 |
axis->SetTitle("OBT (ms)"); |
637 |
//axis->CenterTitle(); |
638 |
axis->SetTitleOffset(0.4); |
639 |
axis->SetTitleSize(0.21); |
640 |
axis->Draw(); |
641 |
TF1 *func1 = new TF1("func1", "0"); |
642 |
func1->SetRange(xmin,xmax); |
643 |
func1->SetLineColor(1); |
644 |
func1->SetLineStyle(1); |
645 |
func1->SetLineWidth(1); |
646 |
func1->Draw("same"); |
647 |
} |
648 |
if (k==11){ |
649 |
TGaxis *axis = new TGaxis (xmin, 0, xmax, 0, xmin, xmax, 510, "+-", 0); |
650 |
axis->SetLineWidth(1); |
651 |
axis->SetTickSize(2); |
652 |
axis->Draw("same"); |
653 |
} |
654 |
} |
655 |
|
656 |
if (k==13){ |
657 |
Canvascc->cd(13); |
658 |
|
659 |
TGraph *ccrecordd = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1d); |
660 |
ccrecordd->SetTitle(""); |
661 |
ccrecordd->SetFillColor(8); |
662 |
ccrecordd->GetXaxis()->SetAxisColor(0); |
663 |
ccrecordd->GetXaxis()->SetLabelSize(0); |
664 |
ccrecordd->GetXaxis()->SetTickLength(0); |
665 |
ccrecordd->GetYaxis()->SetLabelSize(0); |
666 |
ccrecordd->GetYaxis()->SetTickLength(0); |
667 |
ccrecordd->SetMaximum(1.3); |
668 |
ccrecordd->SetMinimum(-1.3); |
669 |
ccrecordd->Draw("AB"); |
670 |
|
671 |
TGraph *ccrecordu = new TGraph(recordstot, xrecordobtcc_1, yccdiagacq_1u); |
672 |
ccrecordu->SetTitle(""); |
673 |
ccrecordu->SetFillColor(17); |
674 |
ccrecordu->GetXaxis()->SetAxisColor(0); |
675 |
ccrecordu->GetXaxis()->SetLabelSize(0); |
676 |
ccrecordu->GetXaxis()->SetTickLength(0); |
677 |
ccrecordu->GetYaxis()->SetLabelSize(0); |
678 |
ccrecordu->GetYaxis()->SetTickLength(0); |
679 |
ccrecordu->SetMaximum(1.3); |
680 |
ccrecordu->SetMinimum(-1.3); |
681 |
ccrecordu->Draw("Bsame"); |
682 |
|
683 |
TGraph *ccrecorderroru = new TGraph(recordstot, xrecordobtcc_1, yccipmerror); |
684 |
ccrecorderroru->SetTitle(""); |
685 |
ccrecorderroru->SetFillColor(2); |
686 |
ccrecorderroru->GetXaxis()->SetLabelSize(0); |
687 |
ccrecorderroru->GetXaxis()->SetTickLength(0); |
688 |
ccrecorderroru->GetYaxis()->SetLabelSize(0); |
689 |
ccrecorderroru->GetYaxis()->SetTickLength(0); |
690 |
ccrecorderroru->SetMaximum(1.3); |
691 |
ccrecorderroru->SetMinimum(-1.3); |
692 |
ccrecorderroru->GetXaxis()->SetAxisColor(0); |
693 |
ccrecorderroru->Draw("Bsame"); |
694 |
|
695 |
|
696 |
Double_t xmin=ccrecordu->GetXaxis()->GetXmin(); |
697 |
Double_t xmax=ccrecordu->GetXaxis()->GetXmax(); |
698 |
TGaxis *axis = new TGaxis (xmin, 0, xmax, 0, xmin, xmax, 510, "+-", 0); |
699 |
axis->SetLineWidth(1); |
700 |
axis->SetTickSize(2); |
701 |
axis->Draw("same"); |
702 |
} |
703 |
|
704 |
} |
705 |
|
706 |
} |
707 |
|
708 |
Canvascc->cd(); |
709 |
leg1->Draw(); |
710 |
Canvascc->Update(); |
711 |
|
712 |
// TH Graph |
713 |
Canvasthdea->cd(); |
714 |
//thermistors number 12 |
715 |
if (tr->GetBranch("Records.TM_TH_ANA[12]")){ |
716 |
for (Int_t k =0; k<12; k++){ |
717 |
for (Int_t i = 0; i < recordstot; i++){ |
718 |
ythana_1[i]= ythana[12*i+k]; |
719 |
xrecordobtth_1[i]= xrecordobtth[12*i+k]; |
720 |
limth[i] = 45; |
721 |
} |
722 |
Canvasthdea->cd(k+1); |
723 |
TGraph *graph = new TGraph(recordstot, xrecordobtth_1, limth); |
724 |
graph->SetMarkerColor(50); |
725 |
graph->SetLineColor(50); |
726 |
graph->SetMarkerStyle(21); |
727 |
graph->SetMarkerSize(0.3); |
728 |
graph->SetMinimum(-10); |
729 |
graph->SetMaximum(60); |
730 |
graph->GetXaxis()->SetTitle("OBT"); |
731 |
graph->GetXaxis()->CenterTitle(); |
732 |
graph->GetYaxis()->SetTitle("Temperature value (°C)"); |
733 |
graph->GetYaxis()->CenterTitle(); |
734 |
oss.str(""); |
735 |
oss << filename.Data() <<": " <<titoloth[k] ; |
736 |
graph->SetTitle(oss.str().c_str()); |
737 |
|
738 |
TPaveText *pt = new TPaveText (.45,.7,.55,.77); |
739 |
pt->AddText("OFF"); |
740 |
pt->SetBorderSize(1); |
741 |
pt->Draw(); |
742 |
graph->Draw("ACP"); |
743 |
pt->Draw(); |
744 |
|
745 |
TGraph *threcord = new TGraph(recordstot, xrecordobtth_1, ythana_1); |
746 |
threcord->SetLineColor(kBlue); |
747 |
threcord->SetMarkerColor(kBlue); |
748 |
threcord->SetMarkerStyle(21); |
749 |
threcord->Draw("C"); |
750 |
} |
751 |
} |
752 |
//thermistors number=16 |
753 |
if (tr->GetBranch("Records.TM_TH_ANA[16]")){ |
754 |
for (Int_t k =0; k<16; k++){ |
755 |
for (Int_t i = 0; i < recordstot; i++){ |
756 |
ythana_1[i]= ythana[16*i+k]; |
757 |
xrecordobtth_1[i]= xrecordobtth[16*i+k]; |
758 |
limth[i] = 45; |
759 |
} |
760 |
Canvasthdea->cd(k+1); |
761 |
TGraph *graph = new TGraph(recordstot, xrecordobtth_1, limth); |
762 |
graph->SetMarkerColor(50); |
763 |
graph->SetLineColor(50); |
764 |
graph->SetMarkerStyle(21); |
765 |
graph->SetMarkerSize(0.3); |
766 |
graph->SetMinimum(-10); |
767 |
graph->SetMaximum(60); |
768 |
graph->GetXaxis()->SetTitle("OBT"); |
769 |
graph->GetXaxis()->CenterTitle(); |
770 |
graph->GetYaxis()->SetTitle("Temperature value (°C)"); |
771 |
graph->GetYaxis()->CenterTitle(); |
772 |
oss.str(""); |
773 |
oss << filename.Data() <<": " <<titoloth[k] ; |
774 |
graph->SetTitle(oss.str().c_str()); |
775 |
|
776 |
TPaveText *pt = new TPaveText (.45,.7,.55,.77); |
777 |
pt->AddText("OFF"); |
778 |
pt->SetBorderSize(1); |
779 |
pt->Draw(); |
780 |
graph->Draw("ACP"); |
781 |
pt->Draw(); |
782 |
|
783 |
TGraph *threcord = new TGraph(recordstot, xrecordobtth_1, ythana_1); |
784 |
threcord->SetLineColor(kBlue); |
785 |
threcord->SetMarkerColor(kBlue); |
786 |
threcord->SetMarkerStyle(21); |
787 |
threcord->Draw("C"); |
788 |
} |
789 |
} |
790 |
|
791 |
// DEA Graph |
792 |
Canvasvoltdea->cd(); |
793 |
for (Int_t k =0; k<6; k++){ |
794 |
for (Int_t i = 0; i < recordstot; i++){ |
795 |
ydea_1[i]= ydea[6*i + k]; |
796 |
xrecordobtdea_1[i]= xrecordobtdea[6*i+k]; |
797 |
limvolt1[i]= 0.65; |
798 |
limvolt2[i]= 3.15; |
799 |
limvolt3[i]= 4.25; |
800 |
} |
801 |
Canvasvoltdea->cd(k+1); |
802 |
TGraph *graph1 = new TGraph(recordstot, xrecordobtth_1, limvolt1); |
803 |
TPaveText *pt = new TPaveText (.84,.13,.92,.19); |
804 |
pt->AddText("OFF"); |
805 |
pt->SetBorderSize(1); |
806 |
pt->Draw(); |
807 |
TPaveText *pt1 = new TPaveText (.84,.67,.92,.73); |
808 |
pt1->AddText("ON"); |
809 |
pt1->SetBorderSize(1); |
810 |
pt1->Draw(); |
811 |
graph1->SetMarkerColor(50); |
812 |
graph1->SetLineColor(50); |
813 |
graph1->SetMarkerStyle(21); |
814 |
graph1->SetMarkerSize(0.3); |
815 |
graph1->SetMinimum(0.0); |
816 |
graph1->SetMaximum(5.0); |
817 |
graph1->GetXaxis()->SetTitle("OBT"); |
818 |
graph1->GetXaxis()->CenterTitle(); |
819 |
graph1->GetYaxis()->SetTitle("Voltage value (V)"); |
820 |
graph1->GetYaxis()->CenterTitle(); |
821 |
oss.str(""); |
822 |
oss << filename.Data()<<": IPM "<<k+1<<" Voltage"; |
823 |
graph1->SetTitle(oss.str().c_str()); |
824 |
graph1->Draw("ACP"); |
825 |
pt1->Draw(); |
826 |
pt->Draw(); |
827 |
|
828 |
TGraph *graph2 = new TGraph(recordstot, xrecordobtth_1, limvolt2); |
829 |
graph2->SetMarkerColor(50); |
830 |
graph2->SetLineColor(50); |
831 |
graph2->SetMarkerStyle(21); |
832 |
graph2->SetMarkerSize(0.2); |
833 |
graph2->SetMinimum(0.0); |
834 |
graph2->SetMaximum(5.0); |
835 |
graph2->Draw("CP"); |
836 |
|
837 |
TGraph *graph3 = new TGraph(recordstot, xrecordobtth_1, limvolt3); |
838 |
graph3->SetMarkerColor(50); |
839 |
graph3->SetLineColor(50); |
840 |
graph3->SetMarkerStyle(21); |
841 |
graph3->SetMarkerSize(0.2); |
842 |
graph3->SetMinimum(0.0); |
843 |
graph3->SetMaximum(5.0); |
844 |
graph3->Draw("CP"); |
845 |
|
846 |
TGraph *dearecord = new TGraph(recordstot, xrecordobtdea_1, ydea_1); |
847 |
dearecord->SetLineColor(kBlue); |
848 |
dearecord->SetLineWidth(1); |
849 |
dearecord->SetMarkerStyle(21); |
850 |
dearecord->SetMarkerSize(0.5); |
851 |
dearecord->SetMarkerColor(kBlue); |
852 |
dearecord->Draw("CPL"); |
853 |
} |
854 |
|
855 |
|
856 |
stringstream oss1, oss2, oss3, oss4; |
857 |
oss.str(""); |
858 |
oss1.str(""); |
859 |
oss2.str(""); |
860 |
oss << outDir.Data() << filename.Data() << "_Tmtc_cc_"<<interval+1<<"." << format; |
861 |
oss1 << outDir.Data() << filename.Data() << "_Tmtc_th_dea_"<<interval+1<<"." << format; |
862 |
oss2 << outDir.Data() << filename.Data() << "_Tmtc_volt_dea_"<<interval+1<<"." << format; |
863 |
Canvascc->SaveAs(oss.str().c_str()); |
864 |
Canvasthdea->SaveAs(oss1.str().c_str()); |
865 |
Canvasvoltdea->SaveAs(oss2.str().c_str()); //// fine del lavoro sul sottogruppo |
866 |
|
867 |
Canvascc->Clear(); |
868 |
Canvasthdea->Clear(); |
869 |
Canvasvoltdea->Clear(); |
870 |
//gROOT->GetListOfCanvases()->Delete(); |
871 |
//gROOT->Reset(); |
872 |
|
873 |
} |
874 |
|
875 |
|
876 |
file->Close(); |
877 |
} |
878 |
|
879 |
|
880 |
////////////////////////////////////////////////////////////////////////////////// |
881 |
///////////////////////// MAIN /////////////////////////////////////////////// |
882 |
////////////////////////////////////////////////////////////////////////////////// |
883 |
|
884 |
|
885 |
int main(int argc, char* argv[]){ |
886 |
TString path; |
887 |
TString outDir = "./"; |
888 |
TString format = "jpg"; |
889 |
|
890 |
if (argc < 2){ |
891 |
printf("You have to insert at least the file to analyze \n"); |
892 |
printf("Try '--help' for more information. \n"); |
893 |
exit(1); |
894 |
} |
895 |
|
896 |
if (!strcmp(argv[1], "--help")){ |
897 |
printf( "Usage: TmtcTemperature FILE [OPTION] \n"); |
898 |
printf( "\t --help Print this help and exit \n"); |
899 |
printf( "\t -outDir[path] Path where to put the output [default ./] \n"); |
900 |
printf( "\t -format[jpg|ps|gif] Format for output files [default 'jpg'] \n"); |
901 |
exit(1); |
902 |
} |
903 |
|
904 |
path=argv[1]; |
905 |
|
906 |
for (int i = 2; i < argc; i++){ |
907 |
|
908 |
|
909 |
if (!strcmp(argv[i], "-outDir")){ |
910 |
if (++i >= argc){ |
911 |
printf( "-outDir needs arguments. \n"); |
912 |
printf( "Try '--help' for more information. \n"); |
913 |
exit(1); |
914 |
} |
915 |
else{ |
916 |
outDir = argv[i]; |
917 |
continue; |
918 |
} |
919 |
} |
920 |
|
921 |
if (!strcmp(argv[i], "-format")){ |
922 |
if (++i >= argc){ |
923 |
printf( "-format needs arguments. \n"); |
924 |
printf( "Try '--help' for more information. \n"); |
925 |
exit(1); |
926 |
} |
927 |
else{ |
928 |
format = argv[i]; |
929 |
continue; |
930 |
} |
931 |
} |
932 |
|
933 |
|
934 |
} |
935 |
|
936 |
TmtcTemperature(argv[1], outDir, format); |
937 |
|
938 |
} |