1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
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
21 Revision 1.8 2007/05/03 09:45:09 decaro
22 Coding convention: RN11 violation -> suppression
24 Revision 1.7 2007/05/02 14:09:39 arcelli
25 Retrieval of Env. Temperature removed (will get it from the GRP)
27 Revision 1.6 2007/04/04 17:19:19 arcelli
28 Moved some printout to debug level
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
33 Revision 1.4 2007/02/19 15:41:55 decaro
34 Coding convention: few corrections
36 Revision 1.3 2007/01/24 11:19:58 arcelli
37 Modify ProcessData to return a logical (CZ)
39 Revision 1.2 2006/12/18 18:17:38 arcelli
40 Updated Aliases for DCS TOF datapoints (C.Zampolli)
42 Revision 1.1 2006/10/26 09:10:52 arcelli
43 Class for handling the TOF DCS data in the Shuttle (C.Zampolli)
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
54 #include "TTimeStamp.h"
58 #include "AliDCSValue.h"
61 #include "AliTOFDataDCS.h"
62 #include "AliTOFFormatDCS.h"
68 ClassImp(AliTOFDataDCS)
70 //---------------------------------------------------------------
71 AliTOFDataDCS::AliTOFDataDCS():
76 fStartTimeDCSQuery(0),
84 for(int i=0;i<kNHV;i++) {
94 //---------------------------------------------------------------
95 AliTOFDataDCS::AliTOFDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime, UInt_t startTimeDCSQuery, UInt_t endTimeDCSQuery):
98 fStartTime(startTime),
100 fStartTimeDCSQuery(startTimeDCSQuery),
101 fEndTimeDCSQuery(endTimeDCSQuery),
102 fIsProcessed(kFALSE),
106 // constructor with arguments
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()));
118 //---------------------------------------------------------------
120 AliTOFDataDCS::AliTOFDataDCS(const AliTOFDataDCS & data):
123 fStartTime(data.fStartTime),
124 fEndTime(data.fEndTime),
125 fStartTimeDCSQuery(data.fStartTimeDCSQuery),
126 fEndTimeDCSQuery(data.fEndTimeDCSQuery),
127 fIsProcessed(data.fIsProcessed),
134 for(int i=0;i<kNAliases;i++) {
135 fAliasNames[i]=data.fAliasNames[i];
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];
147 //---------------------------------------------------------------
149 AliTOFDataDCS& AliTOFDataDCS:: operator=(const AliTOFDataDCS & data) {
151 // assignment operator
156 TObject::operator=(data);
158 fStartTime=data.GetStartTime();
159 fEndTime=data.GetEndTime();
160 fStartTimeDCSQuery=data.GetStartTimeDCSQuery();
161 fEndTimeDCSQuery=data.GetEndTimeDCSQuery();
163 for(int i=0;i<kNAliases;i++) {
164 fAliasNames[i]=data.GetAliasName(i);
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);
177 //---------------------------------------------------------------
178 AliTOFDataDCS::~AliTOFDataDCS() {
182 for(int i=0;i<kNHV;i++) {
195 //---------------------------------------------------------------
196 Bool_t AliTOFDataDCS::ProcessData(TMap& aliasMap){
198 // method to process the data
200 if(!(fAliasNames[0])) Init();
209 Float_t timedelta[2];
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));
216 if (fEndTime==fStartTime){
217 AliError(Form(" Run with null time length: start time = %i = end time = %i",fStartTime,fEndTime));
223 AliDCSValue* aValue0;
224 AliDCSValue* aValue1;
226 // starting loop on aliases
227 for(int j=0; j<kNAliases; j++){
228 for (Int_t k=0;k<2;k++) {
232 //AliInfo(Form("j = %i, with alias = %s",j,fAliasNames[j].Data()));
233 aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
235 AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
237 return kFALSE; // returning only in case we are not in a FDR run
244 Introduce(j, aliasArr);
246 if(aliasArr->GetEntries()==0){
247 AliError(Form("Alias %s has no entries! Nothing will be stored",
248 fAliasNames[j].Data()));
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());
258 fHVvp[j]->SetFloat(0,val0);
259 fHVvp[j]->SetTimeStampFloat(0,time0);
262 fHVvn[j-kNHV]->SetFloat(0,val0);
263 fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
266 fHVip[j-2*kNHV]->SetFloat(0,val0);
267 fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
270 fHVin[j-3*kNHV]->SetFloat(0,val0);
271 fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
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());
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);
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);
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);
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);
319 TIter iterarray(aliasArr);
321 Int_t nentries = aliasArr->GetEntries();
322 Int_t deltaTimeStamp = (Int_t) nentries/3;
323 Int_t deltaTimeStamp1 = (Int_t) nentries/2;
325 // filling aliases with 10 floats+1 Usign
327 for (Int_t k=0;k<3;k++){
328 index = deltaTimeStamp*k;
333 index=deltaTimeStamp1;
338 aValue = (AliDCSValue*) aliasArr->At(index);
339 val = aValue->GetFloat();
340 time = (Float_t) (aValue->GetTimeStamp());
342 fHVvp[j]->SetFloat(k,val);
343 fHVvp[j]->SetTimeStampFloat(k,time);
346 fHVvn[j-kNHV]->SetFloat(k,val);
347 fHVvn[j-kNHV]->SetTimeStampFloat(k,time);
350 fHVip[j-2*kNHV]->SetFloat(k,val);
351 fHVip[j-2*kNHV]->SetTimeStampFloat(k,time);
354 fHVin[j-3*kNHV]->SetFloat(k,val);
355 fHVin[j-3*kNHV]->SetTimeStampFloat(k,time);
359 //computing the most significant variations
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;
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;
375 if (delta[1]<=delta[0]) {
376 Float_t temp = delta[1];
377 Float_t timetemp = timedelta[1];
380 timedelta[1]=timedelta[0];
381 timedelta[0]=timetemp;
385 for (Int_t kk=0;kk<2;kk++){
387 fHVvp[j]->SetDelta(kk,delta[kk]);
388 fHVvp[j]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
391 fHVvn[j-kNHV]->SetDelta(kk,delta[kk]);
392 fHVvn[j-kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
395 fHVip[j-2*kNHV]->SetDelta(kk,delta[kk]);
396 fHVip[j-2*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
399 fHVin[j-3*kNHV]->SetDelta(kk,delta[kk]);
400 fHVin[j-3*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
410 //---------------------------------------------------------------
411 void AliTOFDataDCS::Init(){
413 // initialization of aliases and DCS data
416 for(int i=0;i<kNAliases;i++){
419 fAliasNames[i] = "tof_hv_vp_";
420 sindex.Form("%02i",i);
421 fAliasNames[i] += sindex;
422 fHVvp[i] = new AliTOFFormatDCS();
425 fAliasNames[i] = "tof_hv_vn_";
426 sindex.Form("%02i",i-kNHV);
427 fAliasNames[i] += sindex;
428 fHVvn[i-kNHV] = new AliTOFFormatDCS();
432 fAliasNames[i] = "tof_hv_ip_";
433 sindex.Form("%02i",i-2*kNHV);
434 fAliasNames[i] += sindex;
435 fHVip[i-2*kNHV] = new AliTOFFormatDCS();
438 fAliasNames[i] = "tof_hv_in_";
439 sindex.Form("%02i",i-3*kNHV);
440 fAliasNames[i] += sindex;
441 fHVin[i-3*kNHV] = new AliTOFFormatDCS();
447 //---------------------------------------------------------------
448 void AliTOFDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const
451 // method to introduce new aliases
454 entries = aliasArr->GetEntries();
457 AliDebug(2,Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
458 AliDebug(2,Form(" %d DP values collected",entries));
462 //---------------------------------------------------------------
463 void AliTOFDataDCS::Draw(const Option_t* /*option*/)
465 // Draw all histos and graphs
467 if(!fIsProcessed) return;
470 TString canvasHistoName="Histos";
471 ch=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);