Added macro to lauch TOF QA task with extended functionality via plugin
[u/mrichter/AliRoot.git] / TOF / AliTOFDataDCS.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /*
17 $Log: AliTOFDataDCS.cxx,v $
18 Revision 1.9  2007/05/04 14:02:45  decaro
19 AliTOFDataDCS::Draw(Option_t *) method declared const: compiling warning suppression
20
21 Revision 1.8  2007/05/03 09:45:09  decaro
22 Coding convention: RN11 violation -> suppression
23
24 Revision 1.7  2007/05/02 14:09:39  arcelli
25 Retrieval of Env. Temperature removed (will get it from the GRP)
26
27 Revision 1.6  2007/04/04 17:19:19  arcelli
28 Moved some printout to debug level
29
30 Revision 1.5  2007/02/20 15:57:00  decaro
31 Raw data update: to read the TOF raw data defined in UNPACKED mode
32
33 Revision 1.4  2007/02/19 15:41:55  decaro
34 Coding convention: few corrections
35
36 Revision 1.3  2007/01/24 11:19:58  arcelli
37 Modify ProcessData to return a logical (CZ)
38
39 Revision 1.2  2006/12/18 18:17:38  arcelli
40 Updated Aliases for DCS TOF datapoints (C.Zampolli)
41
42 Revision 1.1  2006/10/26 09:10:52  arcelli
43 Class for handling the TOF DCS data in the Shuttle (C.Zampolli)
44
45 */  
46
47 // AliTOFDataDCS class
48 // main aim to introduce the aliases for the TOF DCS
49 // data points to be then
50 // stored in the OCDB, and to process them. 
51 // Process() method called by TOFPrepr
52
53 #include "TString.h"
54 #include "TTimeStamp.h"
55 #include "TMap.h"
56 #include "TCanvas.h"
57
58 #include "AliDCSValue.h"
59 #include "AliLog.h"
60
61 #include "AliTOFDataDCS.h"
62 #include "AliTOFFormatDCS.h"
63
64 class TH2;
65 class AliCDBMetaData;
66 class TDatime;
67
68 ClassImp(AliTOFDataDCS)
69
70 //---------------------------------------------------------------
71 AliTOFDataDCS::AliTOFDataDCS():
72         TObject(),
73         fRun(0),
74         fStartTime(0),
75         fEndTime(0),
76         fStartTimeDCSQuery(0),
77         fEndTimeDCSQuery(0),
78         fIsProcessed(kFALSE),
79         fFDR(kFALSE)
80 {
81
82   // main constructor 
83
84   for(int i=0;i<kNHV;i++) {
85     fHVvp[i]=0x0;
86     fHVvn[i]=0x0;
87     fHVip[i]=0x0;
88     fHVin[i]=0x0;
89   }
90     
91   
92 }
93
94 //---------------------------------------------------------------
95 AliTOFDataDCS::AliTOFDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime, UInt_t startTimeDCSQuery, UInt_t endTimeDCSQuery):
96         TObject(),
97         fRun(nRun),
98         fStartTime(startTime),
99         fEndTime(endTime),
100         fStartTimeDCSQuery(startTimeDCSQuery),
101         fEndTimeDCSQuery(endTimeDCSQuery),
102         fIsProcessed(kFALSE),
103         fFDR(kFALSE)
104 {
105
106   // constructor with arguments
107
108         AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s \n\tStartTime DCS Query %s \n\tEndTime DCS Query %s", nRun,
109         TTimeStamp(startTime).AsString(),
110         TTimeStamp(endTime).AsString(), 
111         TTimeStamp(startTimeDCSQuery).AsString(), 
112         TTimeStamp(endTimeDCSQuery).AsString()));
113
114         Init();
115
116 }
117
118 //---------------------------------------------------------------
119
120 AliTOFDataDCS::AliTOFDataDCS(const AliTOFDataDCS & data):
121   TObject(data), 
122   fRun(data.fRun),
123   fStartTime(data.fStartTime),
124   fEndTime(data.fEndTime),
125   fStartTimeDCSQuery(data.fStartTimeDCSQuery),
126   fEndTimeDCSQuery(data.fEndTimeDCSQuery),
127   fIsProcessed(data.fIsProcessed),
128   fFDR(data.fFDR)
129
130 {
131
132   // copy constructor
133
134   for(int i=0;i<kNAliases;i++) {
135     fAliasNames[i]=data.fAliasNames[i];
136   }
137  
138   for(int i=0;i<kNHV;i++) {
139     fHVvp[i]=data.fHVvp[i];
140     fHVvn[i]=data.fHVvn[i];
141     fHVip[i]=data.fHVip[i];
142     fHVin[i]=data.fHVin[i];
143   }
144   
145     
146 }
147 //---------------------------------------------------------------
148
149 AliTOFDataDCS& AliTOFDataDCS:: operator=(const AliTOFDataDCS & data) { 
150
151   // assignment operator
152
153   if (this == &data)
154     return *this;
155
156   TObject::operator=(data);
157   fRun=data.GetRun();
158   fStartTime=data.GetStartTime();
159   fEndTime=data.GetEndTime();
160   fStartTimeDCSQuery=data.GetStartTimeDCSQuery();
161   fEndTimeDCSQuery=data.GetEndTimeDCSQuery();
162
163   for(int i=0;i<kNAliases;i++) {
164     fAliasNames[i]=data.GetAliasName(i);
165   }
166
167   for(int i=0;i<kNHV;i++) {
168     fHVvp[i]=data.GetHVvp(i);
169     fHVvn[i]=data.GetHVvn(i);
170     fHVip[i]=data.GetHVip(i);
171     fHVin[i]=data.GetHVin(i);
172   }
173
174
175   return *this;
176 }
177 //---------------------------------------------------------------
178 AliTOFDataDCS::~AliTOFDataDCS() {
179
180   // destructor
181
182   for(int i=0;i<kNHV;i++) {
183     delete fHVvp[i];
184     fHVvp[i]=0;
185     delete fHVvn[i];
186     fHVvn[i]=0;
187     delete fHVip[i];
188     fHVip[i]=0;
189     delete fHVin[i];
190     fHVin[i]=0;
191   }
192   
193 }
194
195 //---------------------------------------------------------------
196 Bool_t AliTOFDataDCS::ProcessData(TMap& aliasMap){
197
198   // method to process the data
199
200   if(!(fAliasNames[0])) Init();
201
202   Float_t val=0;
203   Float_t val0=0;
204   Float_t val1=0;
205   Float_t time=0; 
206   Float_t time0=0; 
207   Float_t time1=0; 
208   Float_t delta[2];
209   Float_t timedelta[2];
210
211   AliInfo(Form(" Start Time = %i",fStartTime));
212   AliInfo(Form(" End Time = %i",fEndTime));
213   AliInfo(Form(" Start Time DCS Query= %i",fStartTimeDCSQuery));
214   AliInfo(Form(" End Time DCS Query= %i",fEndTimeDCSQuery));
215
216   if (fEndTime==fStartTime){
217     AliError(Form(" Run with null time length: start time = %i = end time = %i",fStartTime,fEndTime));
218     return kFALSE;
219   }
220
221   TObjArray *aliasArr;
222   AliDCSValue* aValue;
223   AliDCSValue* aValue0;
224   AliDCSValue* aValue1;
225
226   // starting loop on aliases
227   for(int j=0; j<kNAliases; j++){
228     for (Int_t k=0;k<2;k++) {
229       delta[k]=0;
230       timedelta[k]=0;
231     }
232     //AliInfo(Form("j = %i, with alias = %s",j,fAliasNames[j].Data()));
233     aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
234     if(!aliasArr){
235       AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
236       if (!fFDR){
237         return kFALSE;    // returning only in case we are not in a FDR run
238       }
239       else {
240         continue;
241       }
242     }
243
244     Introduce(j, aliasArr);
245     
246     if(aliasArr->GetEntries()==0){
247       AliError(Form("Alias %s has no entries! Nothing will be stored",
248                     fAliasNames[j].Data()));
249       continue;
250     }
251     if(aliasArr->GetEntries() == 1){
252       AliWarning(Form("Alias %s has just %d entries! Only first value will be stored",
253                     fAliasNames[j].Data(),aliasArr->GetEntries()));
254       aValue0 = (AliDCSValue*) aliasArr->At(0);
255       val0 = aValue0->GetFloat();
256       time0 = (Float_t) (aValue0->GetTimeStamp());
257       if (j<kNHV){
258         fHVvp[j]->SetFloat(0,val0);
259         fHVvp[j]->SetTimeStampFloat(0,time0);
260       }
261       else if (j<kNHV*2){
262         fHVvn[j-kNHV]->SetFloat(0,val0);
263         fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
264       }
265       else if (j<kNHV*3){
266         fHVip[j-2*kNHV]->SetFloat(0,val0);
267         fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
268       }
269       else if (j<kNHV*4){
270         fHVin[j-3*kNHV]->SetFloat(0,val0);
271         fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
272       }
273       continue;
274     }
275     if(aliasArr->GetEntries()==2){
276       AliWarning(Form("Alias %s has just %d entries! Storing first and second entry only and the delta between these two",
277                     fAliasNames[j].Data(),aliasArr->GetEntries()));
278       aValue0 = (AliDCSValue*) aliasArr->At(0);
279       val0 = aValue0->GetFloat();
280       time0 = (Float_t) (aValue0->GetTimeStamp());
281       aValue1 = (AliDCSValue*) aliasArr->At(1);
282       val1 = aValue1->GetFloat();
283       time1 = (Float_t) (aValue1->GetTimeStamp());
284       if (j<kNHV){
285         fHVvp[j]->SetFloat(0,val0);
286         fHVvp[j]->SetTimeStampFloat(0,time0);
287         fHVvp[j]->SetFloat(1,val1);
288         fHVvp[j]->SetTimeStampFloat(1,time1);
289         fHVvp[j]->SetDelta(0,TMath::Abs(val1-val0));
290         fHVvp[j]->SetTimeStampDelta(0,time1-time0);
291       }
292       else if (j<kNHV*2){
293         fHVvn[j-kNHV]->SetFloat(0,val0);
294         fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
295         fHVvn[j-kNHV]->SetFloat(1,val1);
296         fHVvn[j-kNHV]->SetTimeStampFloat(1,time1);
297         fHVvn[j-kNHV]->SetDelta(0,TMath::Abs(val1-val0));
298         fHVvn[j-kNHV]->SetTimeStampDelta(0,time1-time0);
299       }
300       else if (j<kNHV*3){
301         fHVip[j-2*kNHV]->SetFloat(0,val0);
302         fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
303         fHVip[j-2*kNHV]->SetFloat(1,val1);
304         fHVip[j-2*kNHV]->SetTimeStampFloat(1,time1);
305         fHVip[j-2*kNHV]->SetDelta(0,TMath::Abs(val1-val0));
306         fHVip[j-2*kNHV]->SetTimeStampDelta(0,time1-time0);
307       }
308       else if (j<kNHV*4){
309         fHVin[j-3*kNHV]->SetFloat(0,val0);
310         fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
311         fHVin[j-3*kNHV]->SetFloat(1,val1);
312         fHVin[j-3*kNHV]->SetTimeStampFloat(1,time1);
313         fHVin[j-3*kNHV]->SetDelta(0,TMath::Abs(val1-val0));
314         fHVin[j-3*kNHV]->SetTimeStampDelta(0,time1-time0);
315       }
316       continue;
317     }
318     
319     TIter iterarray(aliasArr);
320     
321     Int_t nentries = aliasArr->GetEntries();
322     Int_t deltaTimeStamp = (Int_t) nentries/3;
323     Int_t deltaTimeStamp1 = (Int_t) nentries/2;
324
325     // filling aliases with 10 floats+1 Usign
326     Int_t index = 0;
327     for (Int_t k=0;k<3;k++){
328       index = deltaTimeStamp*k;
329       if (k==0) {
330         index=0;
331       }
332       else if (k==1) {
333         index=deltaTimeStamp1;
334       } 
335       else if (k==2) {
336         index=nentries-1; 
337       }
338       aValue = (AliDCSValue*) aliasArr->At(index);
339       val = aValue->GetFloat();
340       time = (Float_t) (aValue->GetTimeStamp());
341       if (j<kNHV){
342         fHVvp[j]->SetFloat(k,val);
343         fHVvp[j]->SetTimeStampFloat(k,time);
344       }
345       else if (j<kNHV*2){
346         fHVvn[j-kNHV]->SetFloat(k,val);
347         fHVvn[j-kNHV]->SetTimeStampFloat(k,time);
348       }
349       else if (j<kNHV*3){
350         fHVip[j-2*kNHV]->SetFloat(k,val);
351         fHVip[j-2*kNHV]->SetTimeStampFloat(k,time);
352       }
353       else if (j<kNHV*4){
354         fHVin[j-3*kNHV]->SetFloat(k,val);
355         fHVin[j-3*kNHV]->SetTimeStampFloat(k,time);
356       }
357     }
358   
359     //computing the most significant variations
360
361     //Float_t timeDiff = (Float_t)(fEndTime-fStartTime);
362     Float_t timeDiff = (Float_t)(fEndTimeDCSQuery-fStartTimeDCSQuery);
363     Int_t deltamin = (Int_t)(60/timeDiff*nentries); //sampling every minute
364     Int_t klast = nentries-deltamin;
365       
366     for (Int_t k=0;k<klast;k++){
367       aValue = (AliDCSValue*) aliasArr->At(k);
368       aValue1 = (AliDCSValue*) aliasArr->At(k+deltamin);
369       val = aValue->GetFloat();
370       val1 = aValue1->GetFloat();
371       if (delta[0]<=TMath::Abs(val1-val)) {
372         delta[0]=TMath::Abs(val1-val);
373         timedelta[0] = (Float_t)k;
374       }
375       if (delta[1]<=delta[0]) {
376         Float_t temp = delta[1];
377         Float_t timetemp = timedelta[1];
378         delta[1]=delta[0];
379         delta[0]=temp;
380         timedelta[1]=timedelta[0];
381         timedelta[0]=timetemp;
382       }
383     }
384     
385     for (Int_t kk=0;kk<2;kk++){
386       if (j<kNHV){
387         fHVvp[j]->SetDelta(kk,delta[kk]);
388         fHVvp[j]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
389       }
390       else if (j<kNHV*2){
391         fHVvn[j-kNHV]->SetDelta(kk,delta[kk]);
392         fHVvn[j-kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
393       }
394       else if (j<kNHV*3){
395         fHVip[j-2*kNHV]->SetDelta(kk,delta[kk]);
396         fHVip[j-2*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
397       }
398       else if (j<kNHV*4){
399         fHVin[j-3*kNHV]->SetDelta(kk,delta[kk]);
400         fHVin[j-3*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
401       }
402     }
403   }
404   
405   fIsProcessed=kTRUE;
406
407   return kTRUE;
408 }
409
410 //---------------------------------------------------------------
411 void AliTOFDataDCS::Init(){
412
413   // initialization of aliases and DCS data
414
415   TString sindex;
416   for(int i=0;i<kNAliases;i++){
417     //HV, v
418     if (i<kNHV){
419         fAliasNames[i] = "tof_hv_vp_";
420         sindex.Form("%02i",i);
421         fAliasNames[i] += sindex;
422         fHVvp[i] = new AliTOFFormatDCS();
423     }
424     else if (i<kNHV*2){
425         fAliasNames[i] = "tof_hv_vn_";
426         sindex.Form("%02i",i-kNHV);
427         fAliasNames[i] += sindex;
428         fHVvn[i-kNHV] = new AliTOFFormatDCS();
429     }
430     //HV, i
431     else if (i<kNHV*3){
432         fAliasNames[i] = "tof_hv_ip_";
433         sindex.Form("%02i",i-2*kNHV);
434         fAliasNames[i] += sindex;
435         fHVip[i-2*kNHV] = new AliTOFFormatDCS();
436     }
437     else if (i<kNHV*4){
438         fAliasNames[i] = "tof_hv_in_";
439         sindex.Form("%02i",i-3*kNHV);
440         fAliasNames[i] += sindex;
441         fHVin[i-3*kNHV] = new AliTOFFormatDCS();
442     }
443   }
444 }
445
446
447 //---------------------------------------------------------------
448 void AliTOFDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const
449 {
450
451   // method to introduce new aliases
452
453   int entries=0;
454   entries = aliasArr->GetEntries();
455   int nal=0;
456   nal=numAlias;
457   AliDebug(2,Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
458   AliDebug(2,Form("     %d DP values collected",entries));
459
460 }
461
462 //---------------------------------------------------------------
463 void AliTOFDataDCS::Draw(const Option_t* /*option*/)
464 {
465   // Draw all histos and graphs
466
467   if(!fIsProcessed) return;
468
469   TCanvas *ch;
470   TString canvasHistoName="Histos";
471   ch=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
472   ch->cd();
473
474   // to be implemented
475
476 }
477