1. Re-organization of function calls in TestPreprocessor to make it more meaningful.
[u/mrichter/AliRoot.git] / SHUTTLE / test / AliTPCDataDCS.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 //
18 // This is an example of a container class
19 // of the data retreieved from the DCS archive DB.
20 // It is called by the detector's Preprocessor and
21 // it is stored in the CDB.
22 //
23
24 #include "AliTPCDataDCS.h"
25
26 #include "AliCDBMetaData.h"
27 #include "AliDCSValue.h"
28 #include "AliLog.h"
29
30 #include <TTimeStamp.h>
31 #include <TObjString.h>
32 #include <TH2F.h>
33 #include <TProfile.h>
34 #include <TGraph.h>
35 #include <TDatime.h>
36 #include <TStyle.h>
37 #include <TCanvas.h>
38
39 ClassImp(AliTPCDataDCS)
40
41 //---------------------------------------------------------------
42 AliTPCDataDCS::AliTPCDataDCS():
43         TObject(),
44         fRun(0),
45         fStartTime(0),
46         fEndTime(0),
47         fGraphs("TGraph",kNGraphs),
48         fIsProcessed(kFALSE)
49 {
50 // default constructor
51
52         for(int i=0;i<kNHistos;i++) fHv[i]=0;
53         fFunc=0;
54 }
55
56 //---------------------------------------------------------------
57 AliTPCDataDCS::AliTPCDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime):
58         TObject(),
59         fRun(nRun),
60         fStartTime(startTime),
61         fEndTime(endTime),
62         fGraphs("TGraph",kNGraphs),
63         fIsProcessed(kFALSE)
64 {
65 // constructor
66
67         AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", nRun,
68         TTimeStamp(startTime).AsString(),
69         TTimeStamp(endTime).AsString()));
70
71         for(int i=0;i<kNHistos;i++) fHv[i]=0;
72         fFunc=0;
73
74 }
75
76 //---------------------------------------------------------------
77 AliTPCDataDCS::~AliTPCDataDCS()
78 {
79 // destructor
80
81         Clean();
82 }
83
84 //---------------------------------------------------------------
85 void AliTPCDataDCS::Clean()
86 {
87 // destructor
88
89         for(int i=0;i<kNHistos;i++) {
90                 if(fHv[i]) delete fHv[i];
91                 fHv[i]=0;
92         }
93         fGraphs.Clear("C");
94         if(fFunc) delete fFunc;
95         fFunc=0;
96 }
97
98 //---------------------------------------------------------------
99 void AliTPCDataDCS::Init()
100 {
101 // Init alias names and histos
102
103         TH1::AddDirectory(kFALSE);
104
105         fGraphs.SetOwner(1);
106
107         for(int i=0;i<kNAliases;i++){
108                 fAliasNames[i] = "TpcHvSect0";
109                 fAliasNames[i] += i;
110                 fAliasNames[i] += ".FloatValue";
111         }
112
113         for(int i=0;i<kNHistos;i++){
114                 fHv[i] = new TH1F(fAliasNames[i].Data(),fAliasNames[i].Data(), 20, kHvMin, kHvMax);
115                 fHv[i]->GetXaxis()->SetTitle("Hv");
116         }
117 }
118
119 //---------------------------------------------------------------
120 void AliTPCDataDCS::ProcessData(TMap& aliasMap)
121 {
122 // process TMap of alias values retrieved from the DCS archive DB
123
124         Clean();
125         Init();
126
127         TObjArray *aliasArr;
128         AliDCSValue* aValue;
129         for(int j=0; j<kNAliases; j++){
130                 aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
131                 if(!aliasArr){
132                         AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
133                         continue;
134                 }
135                 Introduce(j, aliasArr);
136
137                 if(aliasArr->GetEntries()<2){
138                         AliError(Form("Alias %s has just %d entries!",
139                                         fAliasNames[j].Data(),aliasArr->GetEntries()));
140                         continue;
141                 }
142
143                 TIter iterarray(aliasArr);
144
145                 Double_t *time = new Double_t[aliasArr->GetEntries()];
146                 Double_t *val = new Double_t[aliasArr->GetEntries()];
147
148                 UInt_t ne=0;
149                 while ((aValue = (AliDCSValue*) iterarray.Next())) {
150                 val[ne] = aValue->GetFloat();
151                 time[ne] = (Double_t) (aValue->GetTimeStamp());
152                 // fill histos (alias 0-2)
153                 if(j < 3) fHv[j]->Fill(val[ne]);
154                 ne++;
155                 }
156                 // fill graphs (alias 3-5)
157                 if(j >= 3) CreateGraph(j, aliasArr->GetEntries(), time, val);
158                 delete[] val;
159                 delete[] time;
160         }
161
162         // calculate mean and rms of the first two histos
163         for(int i=0;i<kNHistos;i++){
164                 fMean[i] = fHv[i]->GetMean();
165                 fWidth[i] = fHv[i]->GetRMS();
166         }
167
168         // pol1 fit of the first graph
169 //      if(fGraphs.GetEntries() > 0){
170 //              ((TGraph*) fGraphs.UncheckedAt(0))->Fit("pol1");
171 //              fFunc = ((TGraph*) fGraphs.UncheckedAt(0))->GetFunction("pol1");
172 //      }
173
174         fIsProcessed=kTRUE;
175
176
177 }
178
179 //---------------------------------------------------------------
180 void AliTPCDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)
181 {
182 // Say something about alias value array being processed
183
184         int entries=aliasArr->GetEntries();
185         AliInfo(Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
186         AliInfo(Form("          %d DP values collected",entries));
187
188 }
189
190 //---------------------------------------------------------------
191 void AliTPCDataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y)
192 {
193 // create TGraph to be filled with DP values
194
195         TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);
196
197         gr->GetXaxis()->SetTimeDisplay(1);
198         gr->SetTitle(fAliasNames[i].Data());
199
200         AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));
201
202
203 }
204
205 //---------------------------------------------------------------
206 void AliTPCDataDCS::Draw(const Option_t* /*option*/)
207 {
208 // Draw all histos and graphs
209
210   if(!fIsProcessed) return;
211
212   TCanvas *ch;
213   TString canvasHistoName="Histos";
214   ch=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
215   ch->Divide(2,2);
216   ch->cd(1);
217   fHv[0]->Draw();
218   ch->cd(2);
219   fHv[1]->Draw();
220   ch->cd(3);
221   fHv[2]->Draw();
222
223
224   if(fGraphs.GetEntries() == 0) return;
225
226   TCanvas *cg;
227   TString canvasGraphName="Graphs";
228   cg=new TCanvas(canvasGraphName,canvasGraphName,40,40,600,600);
229   cg->Divide(2,2);
230   cg->cd(1);
231   ((TGraph*) fGraphs.UncheckedAt(0))->Draw("alp");
232
233   cg->cd(2);
234   ((TGraph*) fGraphs.UncheckedAt(1))->Draw("alp");
235   cg->cd(3);
236   ((TGraph*) fGraphs.UncheckedAt(2))->Draw("alp");
237
238   if(fFunc){
239         cg->cd(4);
240         fFunc->Draw("l");
241   }
242
243 }
244