]>
Commit | Line | Data |
---|---|---|
76b6018a | 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 | ||
80663cb8 | 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 <TH1F.h> | |
76b6018a | 29 | |
30 | class TH2; | |
31 | class AliCDBMetaData; | |
32 | class TDatime; | |
33 | ||
34 | // AliVZERODataDCS class | |
35 | // main aim to introduce the aliases for the VZERO DCS | |
36 | // data points to be then | |
37 | // stored in the OCDB, and to process them. | |
38 | // ProcessData() method called by VZEROPreprocessor | |
39 | ||
40 | ClassImp(AliVZERODataDCS) | |
41 | ||
42 | //_____________________________________________________________________________ | |
43 | AliVZERODataDCS::AliVZERODataDCS(): | |
44 | TObject(), | |
45 | fRun(0), | |
46 | fStartTime(0), | |
47 | fEndTime(0), | |
48 | fGraphs("TGraph",kNGraphs), | |
49 | fIsProcessed(kFALSE) | |
50 | { | |
51 | // Default constructor | |
7495d2be | 52 | for(int i=0;i<kNHvChannel;i++) fDeadChannel[i] = kFALSE; |
76b6018a | 53 | } |
54 | ||
55 | //_____________________________________________________________________________ | |
56 | AliVZERODataDCS::AliVZERODataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime): | |
57 | TObject(), | |
58 | fRun(nRun), | |
59 | fStartTime(startTime), | |
60 | fEndTime(endTime), | |
61 | fGraphs("TGraph",kNGraphs), | |
62 | fIsProcessed(kFALSE) | |
63 | { | |
64 | ||
65 | // constructor with arguments | |
7495d2be | 66 | for(int i=0;i<kNHvChannel;i++) fDeadChannel[i] = kFALSE; |
76b6018a | 67 | |
68 | AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", nRun, | |
69 | TTimeStamp(startTime).AsString(), | |
70 | TTimeStamp(endTime).AsString())); | |
71 | ||
72 | Init(); | |
73 | ||
74 | } | |
75 | ||
76 | //_____________________________________________________________________________ | |
77 | AliVZERODataDCS::~AliVZERODataDCS() { | |
78 | ||
79 | // destructor | |
80 | fGraphs.Clear("C"); | |
81 | ||
82 | } | |
83 | ||
84 | //_____________________________________________________________________________ | |
85 | void AliVZERODataDCS::ProcessData(TMap& aliasMap){ | |
86 | ||
87 | // method to process the data | |
88 | ||
89 | if(!(fAliasNames[0])) Init(); | |
90 | ||
91 | TObjArray *aliasArr; | |
92 | AliDCSValue* aValue; | |
93 | ||
94 | // starting loop on aliases | |
95 | for(int iAlias=0; iAlias<kNAliases; iAlias++){ | |
96 | ||
97 | aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[iAlias].Data()); | |
98 | if(!aliasArr){ | |
99 | AliError(Form("Alias %s not found!", fAliasNames[iAlias].Data())); | |
100 | return; | |
101 | } | |
102 | ||
fad64858 | 103 | //Introduce(iAlias, aliasArr); |
76b6018a | 104 | |
105 | if(aliasArr->GetEntries()<2){ | |
106 | AliError(Form("Alias %s has just %d entries!", | |
107 | fAliasNames[iAlias].Data(),aliasArr->GetEntries())); | |
108 | continue; | |
109 | } | |
110 | ||
111 | TIter iterarray(aliasArr); | |
112 | ||
113 | Int_t nentries = aliasArr->GetEntries(); | |
114 | ||
80663cb8 | 115 | Double_t *times = new Double_t[nentries]; |
116 | Double_t *values = new Double_t[nentries]; | |
76b6018a | 117 | |
118 | UInt_t iValue=0; | |
d2f86b51 | 119 | Float_t variation = 0.0; |
120 | ||
76b6018a | 121 | while((aValue = (AliDCSValue*) iterarray.Next())) { |
80663cb8 | 122 | values[iValue] = aValue->GetFloat(); |
7495d2be | 123 | if(iValue>0) { |
7495d2be | 124 | if(values[iValue-1]>0.) variation = TMath::Abs(values[iValue]-values[iValue-1])/values[iValue-1]; |
125 | if(variation > 0.10) fDeadChannel[GetOfflineChannel(iAlias)] = kTRUE; | |
126 | } | |
80663cb8 | 127 | times[iValue] = (Double_t) (aValue->GetTimeStamp()); |
128 | fHv[iAlias]->Fill(values[iValue]); | |
7495d2be | 129 | printf("%s %f Dead=%d\n",fAliasNames[iAlias].Data(),values[iValue],fDeadChannel[GetOfflineChannel(iAlias)]); |
76b6018a | 130 | iValue++; |
131 | } | |
80663cb8 | 132 | CreateGraph(iAlias, aliasArr->GetEntries(), times, values); // fill graphs |
76b6018a | 133 | |
80663cb8 | 134 | delete[] values; |
135 | delete[] times; | |
76b6018a | 136 | } |
137 | ||
138 | // calculate mean and rms of the first two histos | |
ebee0db9 | 139 | // and convert index to aliroot channel |
76b6018a | 140 | for(int i=0;i<kNAliases;i++){ |
ebee0db9 | 141 | Int_t iChannel = GetOfflineChannel(i); |
142 | fMeanHV[iChannel] = fHv[i]->GetMean(); | |
143 | fWidthHV[iChannel] = fHv[i]->GetRMS(); | |
76b6018a | 144 | } |
76b6018a | 145 | |
146 | fIsProcessed=kTRUE; | |
147 | } | |
148 | ||
149 | //_____________________________________________________________________________ | |
150 | void AliVZERODataDCS::Init(){ | |
151 | ||
152 | // initialization of aliases and DCS data | |
153 | ||
154 | TString sindex; | |
155 | int iAlias = 0; | |
156 | ||
157 | for(int iSide = 0; iSide<2 ; iSide++){ | |
158 | for(int iRing = 0; iRing<4 ; iRing++){ | |
159 | for(int iSector = 0; iSector<8 ; iSector++){ | |
160 | if(iSide == 0) fAliasNames[iAlias] = "V00/HV/V0A/SECTOR"; | |
161 | else fAliasNames[iAlias] = "V00/HV/V0C/SECTOR"; | |
162 | sindex.Form("%d/RING%d",iSector,iRing); | |
163 | fAliasNames[iAlias] += sindex; | |
164 | ||
165 | fHv[iAlias] = new TH1F(fAliasNames[iAlias].Data(),fAliasNames[iAlias].Data(), 2000, kHvMin, kHvMax); | |
166 | fHv[iAlias]->GetXaxis()->SetTitle("Hv"); | |
167 | iAlias++; | |
168 | } | |
169 | } | |
170 | } | |
171 | if(iAlias!=kNAliases) | |
172 | AliError(Form("Number of DCS Aliases defined not correct")); | |
173 | ||
174 | } | |
175 | ||
176 | //_____________________________________________________________________________ | |
177 | void AliVZERODataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const | |
178 | { | |
179 | ||
180 | // method to introduce new aliases | |
181 | ||
c61a7285 | 182 | int entries=aliasArr->GetEntries(); |
76b6018a | 183 | AliInfo(Form("************ Alias: %s **********",fAliasNames[numAlias].Data())); |
184 | AliInfo(Form(" %d DP values collected",entries)); | |
185 | ||
186 | } | |
187 | ||
188 | //_____________________________________________________________________________ | |
189 | void AliVZERODataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y) | |
190 | { | |
191 | ||
192 | // Create graphics | |
193 | ||
194 | TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y); | |
195 | ||
196 | gr->GetXaxis()->SetTimeDisplay(1); | |
197 | gr->SetTitle(fAliasNames[i].Data()); | |
198 | ||
199 | AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast())); | |
200 | ||
201 | } | |
202 | ||
203 | ||
204 | //_____________________________________________________________________________ | |
205 | void AliVZERODataDCS::Draw(const Option_t* /*option*/) | |
206 | { | |
207 | // Draw all histos and graphs | |
208 | ||
209 | if(!fIsProcessed) return; | |
210 | ||
211 | if(fGraphs.GetEntries()==0) return; | |
212 | ||
80663cb8 | 213 | TString canvasName; |
76b6018a | 214 | TCanvas *cHV[8]; |
215 | ||
216 | for(int iSide = 0 ;iSide<2;iSide++){ | |
217 | for(int iRing=0;iRing<4;iRing++){ | |
80663cb8 | 218 | if(iSide == 0) canvasName = "V0A_Ring"; |
219 | else canvasName = "V0C_Ring"; | |
220 | canvasName += iRing; | |
76b6018a | 221 | int iCanvas = iSide*4 + iRing; |
80663cb8 | 222 | cHV[iCanvas] = new TCanvas(canvasName,canvasName); |
76b6018a | 223 | cHV[iCanvas]->Divide(3,3); |
224 | for(int iSector=0;iSector<8;iSector++){ | |
225 | cHV[iCanvas]->cd(iSector+1); | |
226 | int iChannel = iSide*32 + iRing*8 + iSector; | |
227 | ((TGraph*) fGraphs.UncheckedAt(iChannel))->SetMarkerStyle(20); | |
228 | ((TGraph*) fGraphs.UncheckedAt(iChannel))->Draw("ALP"); | |
229 | ||
230 | } | |
231 | ||
232 | } | |
233 | } | |
234 | ||
235 | } | |
236 |