]>
Commit | Line | Data |
---|---|---|
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 | Produces the data needed to calculate the quality assurance. | |
19 | All data must be mergeable objects. | |
20 | */ | |
21 | ||
22 | // --- ROOT system --- | |
23 | #include <TClonesArray.h> | |
24 | #include <TFile.h> | |
25 | #include <TH1F.h> | |
26 | #include <TH1I.h> | |
27 | #include <TH2I.h> | |
28 | #include <TParameter.h> | |
29 | ||
30 | // --- Standard library --- | |
31 | ||
32 | // --- AliRoot header files --- | |
33 | #include "AliESDEvent.h" | |
34 | #include "AliLog.h" | |
35 | #include "AliCDBManager.h" | |
36 | #include "AliCDBStorage.h" | |
37 | #include "AliCDBEntry.h" | |
38 | #include "AliVZEROQADataMakerRec.h" | |
39 | #include "AliQAChecker.h" | |
40 | #include "AliRawReader.h" | |
41 | #include "AliVZERORawStream.h" | |
42 | #include "AliVZEROReconstructor.h" | |
43 | ||
44 | ||
45 | ClassImp(AliVZEROQADataMakerRec) | |
46 | ||
47 | //____________________________________________________________________________ | |
48 | AliVZEROQADataMakerRec::AliVZEROQADataMakerRec() : | |
49 | AliQADataMakerRec(AliQA::GetDetName(AliQA::kVZERO), "VZERO Quality Assurance Data Maker"), | |
50 | fCalibData(GetCalibData()), | |
51 | fEvent(0) | |
52 | ||
53 | { | |
54 | // constructor | |
55 | ||
56 | for(Int_t i=0; i<64; i++){ | |
57 | fEven[i] = 0; | |
58 | fOdd[i] = 0; } | |
59 | ||
60 | for(Int_t i=0; i<128; i++){ | |
61 | fADC_Mean[i] = 0.0; } | |
62 | } | |
63 | ||
64 | //____________________________________________________________________________ | |
65 | AliVZEROQADataMakerRec::AliVZEROQADataMakerRec(const AliVZEROQADataMakerRec& qadm) : | |
66 | AliQADataMakerRec(), | |
67 | fCalibData(GetCalibData()), | |
68 | fEvent(0) | |
69 | ||
70 | { | |
71 | //copy constructor | |
72 | SetName((const char*)qadm.GetName()) ; | |
73 | SetTitle((const char*)qadm.GetTitle()); | |
74 | } | |
75 | ||
76 | //__________________________________________________________________ | |
77 | AliVZEROQADataMakerRec& AliVZEROQADataMakerRec::operator = (const AliVZEROQADataMakerRec& qadm ) | |
78 | { | |
79 | // Equal operator | |
80 | this->~AliVZEROQADataMakerRec(); | |
81 | new(this) AliVZEROQADataMakerRec(qadm); | |
82 | return *this; | |
83 | } | |
84 | ||
85 | //____________________________________________________________________________ | |
86 | AliVZEROCalibData* AliVZEROQADataMakerRec::GetCalibData() const | |
87 | ||
88 | { | |
89 | AliCDBManager *man = AliCDBManager::Instance(); | |
90 | ||
91 | AliCDBEntry *entry=0; | |
92 | ||
93 | entry = man->Get("VZERO/Calib/Data"); | |
94 | ||
95 | // Retrieval of data in directory VZERO/Calib/Data: | |
96 | ||
97 | AliVZEROCalibData *calibdata = 0; | |
98 | ||
99 | if (entry) calibdata = (AliVZEROCalibData*) entry->GetObject(); | |
100 | if (!calibdata) AliFatal("No calibration data from calibration database !"); | |
101 | ||
102 | return calibdata; | |
103 | } | |
104 | ||
105 | //____________________________________________________________________________ | |
106 | void AliVZEROQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list) | |
107 | { | |
108 | //Detector specific actions at end of cycle | |
109 | // do the QA checking | |
110 | AliQAChecker::Instance()->Run(AliQA::kVZERO, task, list) ; | |
111 | } | |
112 | ||
113 | //____________________________________________________________________________ | |
114 | void AliVZEROQADataMakerRec::InitESDs() | |
115 | { | |
116 | //Create histograms to control ESD | |
117 | ||
118 | TH1I * h1 = new TH1I("hVZERONbPMA", "Number of PMs fired in V0A", 80, 0, 80); | |
119 | h1->Sumw2() ; | |
120 | Add2ESDsList(h1, 0) ; | |
121 | ||
122 | TH1I * h2 = new TH1I("hVZERONbPMC", "Number of PMs fired in V0C", 80, 0, 80); | |
123 | h2->Sumw2() ; | |
124 | Add2ESDsList(h2, 1) ; | |
125 | ||
126 | TH1I * h3 = new TH1I("hVZEROMultA", "Multiplicity in V0A", 50, 0, 50) ; | |
127 | h3->Sumw2() ; | |
128 | Add2ESDsList(h3, 2) ; | |
129 | ||
130 | TH1I * h4 = new TH1I("hVZEROMultC", "Multiplicity in V0C", 50, 0, 50) ; | |
131 | h4->Sumw2() ; | |
132 | Add2ESDsList(h4, 3) ; | |
133 | ||
134 | TH2F * h5 = new TH2F("fVzeroMult", "Vzero multiplicity", | |
135 | 64, -0.5, 63.5,1000, -0.5, 99.5); | |
136 | h5->GetXaxis()->SetTitle("Vzero PMT"); | |
137 | h5->GetYaxis()->SetTitle("Multiplicity"); | |
138 | h5->Sumw2() ; | |
139 | Add2ESDsList(h5, 4) ; | |
140 | TH1F * h6 = new TH1F("fBBA","BB Vzero A", 32, -0.5,31.5); | |
141 | h6->Sumw2(); | |
142 | Add2ESDsList(h6, 5) ; | |
143 | TH1F * h7 = new TH1F("fBGA","BG Vzero A", 32, -0.5,31.5); | |
144 | h7->Sumw2(); | |
145 | Add2ESDsList(h7, 6) ; | |
146 | TH1F * h8 = new TH1F("fBBC","BB Vzero C", 32, -0.5,31.5); | |
147 | h8->Sumw2(); | |
148 | Add2ESDsList(h8, 7) ; | |
149 | TH1F * h9 = new TH1F("fBGC","BG Vzero C", 32, -0.5,31.5); | |
150 | h9->Sumw2(); | |
151 | Add2ESDsList(h9, 8) ; | |
152 | ||
153 | TH2F *h10 = new TH2F("fVzeroAdc", "Vzero Adc", | |
154 | 64, -0.5, 63.5,1024, -0.5, 1023.5); | |
155 | h10->GetXaxis()->SetTitle("Vzero PMT"); | |
156 | h10->GetYaxis()->SetTitle("Adc counts"); | |
157 | h10->Sumw2() ; | |
158 | Add2ESDsList(h10, 9); | |
159 | ||
160 | TH2F *h11 = new TH2F("fVzeroTime", "Vzero Time", | |
161 | 64, -0.5, 63.5,300, -0.5, 149.5); | |
162 | h11->GetXaxis()->SetTitle("Vzero PMT"); | |
163 | h11->GetYaxis()->SetTitle("Time [100 ps]"); | |
164 | h11->Sumw2() ; | |
165 | Add2ESDsList(h11,10); | |
166 | ||
167 | } | |
168 | ||
169 | //____________________________________________________________________________ | |
170 | void AliVZEROQADataMakerRec::InitRaws() | |
171 | { | |
172 | // create Raws histograms in Raws subdir | |
173 | ||
174 | char ADCname[12]; | |
175 | char texte[40]; | |
176 | TH1I *hRawADC0[64]; | |
177 | TH1I *hRawADC1[64]; | |
178 | ||
179 | Bool_t expert = kTRUE ; | |
180 | Bool_t saveCorr = kTRUE ; | |
181 | ||
182 | TH2I *h0 = new TH2I("hCellADCMap0","ADC vs Cell for EVEN Integrator", 70, 0, 70, 512, 0, 1024); | |
183 | h0->Sumw2(); | |
184 | Add2RawsList(h0,0, !expert, !saveCorr) ; | |
185 | TH2I *h1 = new TH2I("hCellADCMap1","ADC vs Cell for ODD Integrator", 70, 0, 70, 512, 0, 1024); | |
186 | h1->Sumw2(); | |
187 | Add2RawsList(h1,1, !expert, !saveCorr) ; | |
188 | ||
189 | TH2F *hMeanADC0 = new TH2F("hCellMeanADCMap0","Mean ADC vs cell for EVEN integrator",70,-0.5,69.5,512,-0.5,511.5); | |
190 | Add2RawsList(hMeanADC0,2, !expert, !saveCorr); | |
191 | TH2F *hMeanADC1 = new TH2F("hCellMeanADCMap1","Mean ADC vs cell for ODD integrator",70,-0.5,69.5,512,-0.5,511.5); | |
192 | Add2RawsList(hMeanADC1,3, !expert, !saveCorr); | |
193 | ||
194 | TH1I *hMulV0A = new TH1I("hMulV0A", "Multiplicity in V0A", 40, 0, 40) ; | |
195 | Add2RawsList(hMulV0A,4, !expert, saveCorr); | |
196 | TH1I *hMulV0C = new TH1I("hMulV0C", "Multiplicity in V0C", 40, 0, 40) ; | |
197 | Add2RawsList(hMulV0C,5, !expert, saveCorr); | |
198 | ||
199 | for (Int_t i=0; i<64; i++) | |
200 | { | |
201 | sprintf(ADCname,"hRaw0ADC%d",i); | |
202 | sprintf(texte,"Raw ADC in cell %d for even integrator",i); | |
203 | hRawADC0[i]= new TH1I(ADCname,texte,1024,0,1024); | |
204 | Add2RawsList(hRawADC0[i],i+6, expert, !saveCorr); | |
205 | ||
206 | sprintf(ADCname,"hRaw1ADC%d",i); | |
207 | sprintf(texte,"Raw ADC in cell %d for odd integrator",i); | |
208 | hRawADC1[i]= new TH1I(ADCname,texte,1024,0,1024); | |
209 | Add2RawsList(hRawADC1[i],i+6+64, expert, !saveCorr); | |
210 | } | |
211 | } | |
212 | ||
213 | //____________________________________________________________________________ | |
214 | void AliVZEROQADataMakerRec::MakeESDs(AliESDEvent * esd) | |
215 | { | |
216 | // make QA data from ESDs | |
217 | ||
218 | AliESDVZERO *esdVZERO=esd->GetVZEROData(); | |
219 | ||
220 | if (esdVZERO) { | |
221 | GetESDsData(0)->Fill(esdVZERO->GetNbPMV0A()); | |
222 | GetESDsData(1)->Fill(esdVZERO->GetNbPMV0C()); | |
223 | GetESDsData(2)->Fill(esdVZERO->GetMTotV0A()); | |
224 | GetESDsData(3)->Fill(esdVZERO->GetMTotV0C()); | |
225 | for(Int_t i=0;i<64;i++) { | |
226 | GetESDsData(4)->Fill((Float_t) i,(Float_t) esdVZERO->GetMultiplicity(i)); | |
227 | GetESDsData(9)->Fill((Float_t) i,(Float_t) esdVZERO->GetAdc(i)); | |
228 | GetESDsData(10)->Fill((Float_t) i,(Float_t) esdVZERO->GetTime(i)); | |
229 | } | |
230 | for(Int_t i=0;i<32;i++) { | |
231 | if(esdVZERO->BBTriggerV0A(i)) | |
232 | GetESDsData(5)->Fill((Float_t) i); | |
233 | if(esdVZERO->BGTriggerV0A(i)) | |
234 | GetESDsData(6)->Fill((Float_t) i); | |
235 | if(esdVZERO->BBTriggerV0C(i)) | |
236 | GetESDsData(7)->Fill((Float_t) i); | |
237 | if(esdVZERO->BGTriggerV0C(i)) | |
238 | GetESDsData(8)->Fill((Float_t) i); | |
239 | } | |
240 | } | |
241 | ||
242 | } | |
243 | ||
244 | //____________________________________________________________________________ | |
245 | void AliVZEROQADataMakerRec::MakeRaws(AliRawReader* rawReader) | |
246 | { | |
247 | //Fill histograms with Raws, computes average ADC values dynamically (pedestal subtracted) | |
248 | ||
249 | rawReader->Reset() ; | |
250 | AliVZERORawStream* rawStream = new AliVZERORawStream(rawReader); | |
251 | rawStream->Next(); | |
252 | ||
253 | Float_t ChargeEoI, Threshold; // for pedestal subtraction | |
254 | ||
255 | // for(Int_t i=0; i<128; i++) { printf(" i = %d pedestal = %f sigma = %f \n\n", | |
256 | // i, fCalibData->GetPedestal(i),fCalibData->GetSigma(i) );} | |
257 | ||
258 | GetRawsData(2)->Reset(); // to keep only the last value of the ADC average | |
259 | GetRawsData(3)->Reset(); | |
260 | ||
261 | Int_t mulV0A = 0 ; | |
262 | Int_t mulV0C = 0 ; | |
263 | ||
264 | Int_t thresholV0A = 50 ; | |
265 | Int_t thresholV0C = 50 ; | |
266 | ||
267 | for(Int_t j=0; j<64; j++) { | |
268 | Int_t i = rawStream->GetOfflineChannel(j); // Convert Online to Offline channel number | |
269 | if(!rawStream->GetIntegratorFlag(i,10)) | |
270 | { | |
271 | // even integrator - fills index6 to 69 | |
272 | GetRawsData(0)->Fill(i,rawStream->GetADC(i)) ; | |
273 | GetRawsData(i+6)->Fill(rawStream->GetADC(i)) ; | |
274 | ChargeEoI = (float)rawStream->GetADC(i) - fCalibData->GetPedestal(i); | |
275 | Threshold = 3.0 * fCalibData->GetSigma(i); | |
276 | if(rawStream->GetADC(i)<1023 && ChargeEoI > Threshold) { | |
277 | fADC_Mean[i] = ((fADC_Mean[i]*fEven[i]) + rawStream->GetADC(i)) / (fEven[i]+1); | |
278 | GetRawsData(2)->Fill((Float_t) i, fADC_Mean[i]); | |
279 | fEven[i]+=1; } | |
280 | } | |
281 | else if(rawStream->GetIntegratorFlag(i,10)) | |
282 | { | |
283 | // odd integrator - fills index 70 to 133 | |
284 | GetRawsData(1)->Fill(i,rawStream->GetADC(i)) ; | |
285 | GetRawsData(i+6+64)->Fill(rawStream->GetADC(i)) ; | |
286 | ChargeEoI = (float)rawStream->GetADC(i) - fCalibData->GetPedestal(i+64) ; | |
287 | Threshold = 3.0 * fCalibData->GetSigma(i+64); | |
288 | if(rawStream->GetADC(i)<1023 && ChargeEoI > Threshold) { | |
289 | fADC_Mean[i+64] = ((fADC_Mean[i+64]*fOdd[i]) + rawStream->GetADC(i)) / (fOdd[i]+1); | |
290 | GetRawsData(3)->Fill((Float_t) i, fADC_Mean[i+64]); | |
291 | fOdd[i]+=1; } | |
292 | } | |
293 | if( i < 32 ){ | |
294 | if (rawStream->GetADC(i) > thresholV0C) | |
295 | mulV0C++ ; | |
296 | }else{ | |
297 | if (rawStream->GetADC(i) > thresholV0A) | |
298 | mulV0A++ ; | |
299 | } | |
300 | } | |
301 | ||
302 | fEvent++; | |
303 | GetRawsData(4)->Fill(mulV0A) ; | |
304 | TParameter<double> * p = dynamic_cast<TParameter<double>*>(GetParameterList()->FindObject(Form("%s_%s_%s", GetName(), AliQA::GetTaskName(AliQA::kRAWS).Data(), GetRawsData(4)->GetName()))) ; | |
305 | p->SetVal(mulV0A) ; | |
306 | GetRawsData(5)->Fill(mulV0C) ; | |
307 | p = dynamic_cast<TParameter<double>*>(GetParameterList()->FindObject(Form("%s_%s_%s", GetName(), AliQA::GetTaskName(AliQA::kRAWS).Data(), GetRawsData(5)->GetName()))) ; | |
308 | p->SetVal(mulV0C) ; | |
309 | } | |
310 | ||
311 | //____________________________________________________________________________ | |
312 | void AliVZEROQADataMakerRec::StartOfDetectorCycle() | |
313 | { | |
314 | // Detector specific actions at start of cycle | |
315 | ||
316 | } |