Fix to HV datapoints retrieval in case the HV is switched off before run-completed...
[u/mrichter/AliRoot.git] / VZERO / AliVZERODataDCS.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 #include "AliVZERODataDCS.h"
18
19 #include "AliDCSValue.h"
20 #include "AliLog.h"
21
22 #include <TGraph.h>
23 #include <TAxis.h>
24 #include <TCanvas.h>
25 #include <TTimeStamp.h>
26 #include <TMap.h>
27 #include <TString.h>
28 #include <TObjString.h>
29 #include <TH1F.h>
30
31 class TH2;
32 class AliCDBMetaData;
33 class TDatime;
34
35 // AliVZERODataDCS class
36 // main aim to introduce the aliases for the VZERO DCS
37 // data points to be then
38 // stored in the OCDB, and to process them. 
39 // ProcessData() method called by VZEROPreprocessor
40
41 ClassImp(AliVZERODataDCS)
42
43 //_____________________________________________________________________________
44 AliVZERODataDCS::AliVZERODataDCS():
45         TObject(),
46         fRun(0),
47         fStartTime(0),
48         fEndTime(0),
49         fDaqStartTime(0),
50         fDaqEndTime(0),
51     fGraphs("TGraph",kNGraphs),
52         fFEEParameters(NULL),
53         fIsProcessed(kFALSE)
54
55 {
56   // Default constructor
57         for(int i=0;i<kNHvChannel;i++) {
58                 fDeadChannel[i] = kFALSE;
59                 fMeanHV[i]      = 100.0;
60                 fWidthHV[i]     = 0.0; 
61         }
62 }
63
64 //_____________________________________________________________________________
65 AliVZERODataDCS::AliVZERODataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime, UInt_t daqStartTime, UInt_t daqEndTime):
66         TObject(),
67         fRun(nRun),
68         fStartTime(startTime),
69         fEndTime(endTime),
70         fDaqStartTime(daqStartTime),
71         fDaqEndTime(daqEndTime),
72         fGraphs("TGraph",kNGraphs),
73         fFEEParameters(new TMap()),
74         fIsProcessed(kFALSE)
75
76 {
77
78   // constructor with arguments
79         for(int i=0;i<kNHvChannel;i++) {
80          fDeadChannel[i] = kFALSE;        
81          fMeanHV[i]      = 100.0;
82      fWidthHV[i]     = 0.0; 
83         }
84         AliInfo(Form("\n\tRun %d \n\tTime Created %s \n\tTime Completed %s \n\tDAQ start %s \n\tDAQ end %s  ", nRun,
85                 TTimeStamp(startTime).AsString(),
86                 TTimeStamp(endTime).AsString(),
87                 TTimeStamp(daqStartTime).AsString(),
88                 TTimeStamp(daqEndTime).AsString()));
89         
90         fFEEParameters->SetOwnerValue();
91         Init();
92
93 }
94
95 //_____________________________________________________________________________
96 AliVZERODataDCS::~AliVZERODataDCS() {
97
98   // destructor
99   fGraphs.Clear("C");
100   delete fFEEParameters;
101
102 }
103
104 //_____________________________________________________________________________
105 void AliVZERODataDCS::ProcessData(TMap& aliasMap){
106
107   // method to process the data
108
109   if(!(fAliasNames[0])) Init();
110
111   TObjArray *aliasArr;
112   AliDCSValue* aValue;
113
114   // starting loop on aliases
115   for(int iAlias=0; iAlias<kNAliases; iAlias++){
116
117     aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[iAlias].Data());
118     if(!aliasArr){
119       AliError(Form("Alias %s not found!", fAliasNames[iAlias].Data()));
120       continue;
121     }
122
123     //Introduce(iAlias, aliasArr);
124     
125     if(aliasArr->GetEntries()<2){
126       AliWarning(Form("Alias %s has just %d entries!",
127                     fAliasNames[iAlias].Data(),aliasArr->GetEntries()));
128     }
129     
130     TIter iterarray(aliasArr);
131         
132     if(iAlias<kNHvChannel){ // Treating HV values
133         Int_t nentries = aliasArr->GetEntries();
134
135         Double_t *times = new Double_t[nentries];
136         Double_t *values = new Double_t[nentries];
137
138         UInt_t iValue=0;
139         Float_t variation = 0.0;
140
141         while((aValue = (AliDCSValue*) iterarray.Next())) {
142                         values[iValue] = aValue->GetFloat();
143                         UInt_t currentTime = aValue->GetTimeStamp();
144                         times[iValue] = (Double_t) (currentTime);
145                         
146                         if(iValue>0) {
147                                 if(values[iValue-1]>0.) variation = TMath::Abs(values[iValue]-values[iValue-1])/values[iValue-1];
148                                 if(currentTime>fDaqEndTime && variation>0.10) continue;
149                                 if(variation > 0.10) fDeadChannel[GetOfflineChannel(iAlias)] = kTRUE;
150                         }
151                         fHv[iAlias]->Fill(values[iValue]);
152                         printf("%s %f Dead=%d\n",fAliasNames[iAlias].Data(),values[iValue],fDeadChannel[GetOfflineChannel(iAlias)]);
153                         iValue++;
154         }      
155         CreateGraph(iAlias, aliasArr->GetEntries(), times, values); // fill graphs 
156
157         // calculate mean and rms of the first two histos
158         // and convert index to aliroot channel
159         Int_t iChannel     = GetOfflineChannel(iAlias); 
160         fMeanHV[iChannel]  = fHv[iAlias]->GetMean();
161         fWidthHV[iChannel] = fHv[iAlias]->GetRMS();
162
163         delete[] values;
164         delete[] times; 
165         } else { // Treating FEE Parameters
166                 AliDCSValue * lastVal = NULL;
167                 while((aValue = (AliDCSValue*) iterarray.Next())) lastVal = aValue; // Take only the last value
168                 fFEEParameters->Add(new TObjString(fAliasNames[iAlias].Data()),lastVal);
169         }      
170   }
171   
172   fIsProcessed=kTRUE;
173 }
174
175 //_____________________________________________________________________________
176 void AliVZERODataDCS::Init(){
177
178   // initialization of aliases and DCS data
179
180   TString sindex;
181   int iAlias = 0;
182   
183   for(int iSide = 0; iSide<2 ; iSide++){
184         for(int iRing = 0; iRing<4 ; iRing++){
185                 for(int iSector = 0; iSector<8 ; iSector++){
186                         if(iSide == 0) fAliasNames[iAlias] = "V00/HV/V0A/SECTOR";
187                         else fAliasNames[iAlias] = "V00/HV/V0C/SECTOR";
188                         sindex.Form("%d/RING%d",iSector,iRing);
189                         fAliasNames[iAlias] += sindex;
190                         
191                         fHv[iAlias] = new TH1F(fAliasNames[iAlias].Data(),fAliasNames[iAlias].Data(), 3000, kHvMin, kHvMax);
192                         fHv[iAlias]->GetXaxis()->SetTitle("Hv");
193                         iAlias++;
194                 }
195         }
196   }
197
198  // Time Resolution Parameters
199         
200         for(int iCIU = 0; iCIU<8 ; iCIU++){
201                 fAliasNames[iAlias++] = Form("V00/FEE/CIU%d/TimeResolution",iCIU);
202                 fAliasNames[iAlias++] = Form("V00/FEE/CIU%d/WidthResolution",iCIU);
203         }
204
205   if(iAlias!=kNAliases) 
206               AliError(Form("Number of DCS Aliases defined not correct"));
207
208 }
209
210 //_____________________________________________________________________________
211 void AliVZERODataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const
212 {
213
214   // method to introduce new aliases
215
216   int entries=aliasArr->GetEntries();
217   AliInfo(Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
218   AliInfo(Form("        %d DP values collected",entries));
219
220 }
221
222 //_____________________________________________________________________________
223 void AliVZERODataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y)
224 {
225
226    // Create graphics
227    
228    TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);
229
230    gr->GetXaxis()->SetTimeDisplay(1);
231    gr->SetTitle(fAliasNames[i].Data());
232
233    AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));
234
235 }
236
237
238 //_____________________________________________________________________________
239 void AliVZERODataDCS::Draw(const Option_t* /*option*/)
240 {
241 // Draw all histos and graphs
242
243   if(!fIsProcessed) return;
244
245   if(fGraphs.GetEntries()==0)  return;
246   
247   TString canvasName;
248   TCanvas *cHV[8];
249   
250   for(int iSide = 0 ;iSide<2;iSide++){
251         for(int iRing=0;iRing<4;iRing++){
252                 if(iSide == 0)  canvasName = "V0A_Ring";
253                 else  canvasName = "V0C_Ring";
254                 canvasName += iRing;
255                 int iCanvas = iSide*4 + iRing;
256                 cHV[iCanvas] = new TCanvas(canvasName,canvasName);
257                 cHV[iCanvas]->Divide(3,3);
258                 for(int iSector=0;iSector<8;iSector++){
259                         cHV[iCanvas]->cd(iSector+1);
260                         int iChannel = iSide*32 + iRing*8 + iSector; 
261                         ((TGraph*) fGraphs.UncheckedAt(iChannel))->SetMarkerStyle(20);
262                         ((TGraph*) fGraphs.UncheckedAt(iChannel))->Draw("ALP");
263
264                 }
265                                 
266         }
267   }
268
269 }
270