]>
Commit | Line | Data |
---|---|---|
a9523c3e | 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 | B.K. Nandi | |
21 | */ | |
22 | ||
23 | // --- ROOT system --- | |
24 | #include <TClonesArray.h> | |
25 | #include <TFile.h> | |
26 | #include <TH1F.h> | |
27 | #include <TH1I.h> | |
28 | #include <TH2F.h> | |
29 | #include <TTree.h> | |
30 | ||
31 | // --- Standard library --- | |
32 | ||
33 | // --- AliRoot header files --- | |
34 | ||
35 | #include "AliLog.h" | |
36 | #include "AliPMDhit.h" | |
37 | #include "AliPMDsdigit.h" | |
38 | #include "AliPMDdigit.h" | |
39 | #include "AliPMDQADataMakerSim.h" | |
40 | #include "AliQAChecker.h" | |
41 | ||
42 | ClassImp(AliPMDQADataMakerSim) | |
43 | ||
44 | //____________________________________________________________________________ | |
45 | AliPMDQADataMakerSim::AliPMDQADataMakerSim() : | |
46 | AliQADataMakerSim(AliQA::GetDetName(AliQA::kPMD), "PMD Quality Assurance Data Maker") | |
47 | { | |
48 | // ctor | |
49 | } | |
50 | ||
51 | //____________________________________________________________________________ | |
52 | AliPMDQADataMakerSim::AliPMDQADataMakerSim(const AliPMDQADataMakerSim& qadm) : | |
53 | AliQADataMakerSim() | |
54 | { | |
55 | //copy ctor | |
56 | SetName((const char*)qadm.GetName()) ; | |
57 | SetTitle((const char*)qadm.GetTitle()); | |
58 | } | |
59 | ||
60 | //__________________________________________________________________ | |
61 | AliPMDQADataMakerSim& AliPMDQADataMakerSim::operator = (const AliPMDQADataMakerSim& qadm ) | |
62 | { | |
63 | // Assign operator. | |
64 | this->~AliPMDQADataMakerSim(); | |
65 | new(this) AliPMDQADataMakerSim(qadm); | |
66 | return *this; | |
67 | } | |
68 | ||
69 | //____________________________________________________________________________ | |
70 | void AliPMDQADataMakerSim::InitHits() | |
71 | { | |
72 | // create Hits histograms in Hits subdir | |
73 | ||
2c1131dd | 74 | TH1F *h0 = new TH1F("hPreHitsEdep","Hits energy distribution in (keV)PRE(PMD)", 500, 0., 500.); |
a9523c3e | 75 | h0->Sumw2() ; |
76 | Add2HitsList(h0, 0) ; | |
77 | ||
2c1131dd | 78 | TH1F *h1 = new TH1F("hCpvHitsEdep","Hits energy distribution in (keV)CPV(PMD)", 500, 0., 500.); |
a9523c3e | 79 | h1->Sumw2() ; |
80 | Add2HitsList(h1, 1) ; | |
81 | ||
2c1131dd | 82 | TH1I *h2 = new TH1I("hPreHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 3000) ; |
a9523c3e | 83 | h2->Sumw2() ; |
84 | Add2HitsList(h2, 2) ; | |
85 | ||
2c1131dd | 86 | TH1I *h3 = new TH1I("hCpvHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 3000) ; |
a9523c3e | 87 | h2->Sumw2() ; |
88 | Add2HitsList(h3, 3) ; | |
89 | } | |
90 | ||
91 | //____________________________________________________________________________ | |
92 | void AliPMDQADataMakerSim::InitSDigits() | |
93 | { | |
94 | // create SDigits histograms in SDigits subdir | |
95 | ||
2c1131dd | 96 | TH1F *h0 = new TH1F("hPreSDigitsEdep","SDigits energy distribution in(keV) PRE(PMD)", 500, 0., 500.); |
a9523c3e | 97 | h0->Sumw2(); |
98 | Add2SDigitsList(h0, 0); | |
99 | ||
2c1131dd | 100 | TH1F *h1 = new TH1F("hCpvSDigitsEdep","SDigits energy distribution in (keV)CPV(PMD)", 500, 0., 500.); |
a9523c3e | 101 | h1->Sumw2(); |
102 | Add2SDigitsList(h1, 1); | |
103 | ||
2c1131dd | 104 | TH1I *h2 = new TH1I("hPreSDigitsMult","SDigits multiplicity distribution in PRE(PMD)", 500, 0., 1000.); |
a9523c3e | 105 | h2->Sumw2(); |
106 | Add2SDigitsList(h2, 2); | |
107 | ||
2c1131dd | 108 | TH1I *h3 = new TH1I("hCpvSDigitsMult","SDigits multiplicity distribution in CPV(PMD)", 500, 0., 1000.); |
a9523c3e | 109 | h3->Sumw2(); |
110 | Add2SDigitsList(h3, 3); | |
111 | ||
112 | } | |
113 | ||
114 | //____________________________________________________________________________ | |
115 | void AliPMDQADataMakerSim::InitDigits() | |
116 | { | |
117 | // create Digits histograms in Digits subdir | |
118 | ||
119 | TH1F *h0 = new TH1F("hPreDigitsEdep","Digits energy distribution in PRE(PMD)", 100, 0., 2000.); | |
120 | h0->Sumw2(); | |
121 | Add2DigitsList(h0, 0); | |
122 | ||
123 | TH1F *h1 = new TH1F("hCpvDigitsEdep","Digits energy distribution in CPV(PMD)", 100, 0., 2000.); | |
124 | h1->Sumw2(); | |
125 | Add2DigitsList(h1, 1); | |
126 | ||
2c1131dd | 127 | TH1I *h2 = new TH1I("hPreDigitsMult","Digits multiplicity distribution in PRE(PMD)", 500, 0, 1000) ; |
a9523c3e | 128 | h2->Sumw2(); |
129 | Add2DigitsList(h2, 2); | |
130 | ||
2c1131dd | 131 | TH1I *h3 = new TH1I("hCpvDigitsMult","Digits multiplicity distribution in CPV(PMD)", 500, 0, 1000); |
a9523c3e | 132 | h3->Sumw2(); |
133 | Add2DigitsList(h3, 3); | |
134 | ||
135 | } | |
136 | ||
137 | //____________________________________________________________________________ | |
138 | void AliPMDQADataMakerSim::MakeHits(TClonesArray *hits) | |
139 | { | |
140 | //make QA data from Hits | |
141 | ||
142 | Int_t premul = 0, cpvmul = 0; | |
2c1131dd | 143 | Float_t edepkev = 0.; |
a9523c3e | 144 | TIter next(hits); |
145 | AliPMDhit * hit; | |
2c1131dd | 146 | |
a9523c3e | 147 | while ( (hit = dynamic_cast<AliPMDhit *>(next())) ) |
2c1131dd | 148 | { |
a9523c3e | 149 | if (hit->Z() > 361.5) |
2c1131dd | 150 | { |
151 | edepkev = (hit->GetEnergy())/1000.; | |
152 | GetHitsData(0)->Fill(edepkev); | |
a9523c3e | 153 | premul++; |
2c1131dd | 154 | } |
a9523c3e | 155 | else if (hit->Z() < 361.5) |
2c1131dd | 156 | { |
157 | edepkev = (hit->GetEnergy())/1000.; | |
158 | GetHitsData(1)->Fill(edepkev); | |
a9523c3e | 159 | cpvmul++; |
2c1131dd | 160 | } |
a9523c3e | 161 | } |
162 | ||
163 | if(premul <= 0) | |
164 | { | |
165 | GetHitsData(2)->Fill(-1.); | |
166 | } | |
167 | else | |
168 | { | |
169 | GetHitsData(2)->Fill(premul); | |
170 | } | |
171 | ||
172 | if(cpvmul <= 0) | |
173 | { | |
174 | GetHitsData(3)->Fill(-1.); | |
175 | } | |
176 | else | |
177 | { | |
178 | GetHitsData(3)->Fill(cpvmul); | |
179 | } | |
180 | ||
181 | } | |
182 | ||
183 | //____________________________________________________________________________ | |
184 | void AliPMDQADataMakerSim::MakeHits(TTree * hitTree) | |
185 | { | |
186 | // make QA data from Hit Tree | |
187 | ||
188 | TBranch * branch = hitTree->GetBranch("PMD") ; | |
189 | if ( ! branch ) | |
190 | { | |
191 | AliWarning("PMD branch in Hit Tree not found") ; | |
192 | return; | |
193 | } | |
194 | ||
195 | static TClonesArray statichits("AliPMDhit", 1000); | |
196 | statichits.Clear(); | |
197 | TClonesArray *hits = &statichits; | |
198 | static TClonesArray staticdummy("AliPMDhit", 1000); | |
199 | staticdummy.Clear(); | |
200 | TClonesArray *dummy = &staticdummy; | |
201 | branch->SetAddress(&dummy); | |
202 | Int_t index = 0 ; | |
203 | ||
204 | for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) { | |
205 | branch->GetEntry(ientry) ; | |
206 | for (Int_t ihit = 0 ; ihit < dummy->GetEntries() ; ihit++) { | |
207 | AliPMDhit * hit = dynamic_cast<AliPMDhit *> (dummy->At(ihit)) ; | |
208 | new((*hits)[index]) AliPMDhit(*hit) ; | |
209 | ||
210 | index++ ; | |
211 | } | |
212 | } | |
213 | ||
214 | MakeHits(hits); | |
215 | ||
216 | } | |
217 | //____________________________________________________________________________ | |
218 | void AliPMDQADataMakerSim::MakeSDigits(TClonesArray * sdigits) | |
219 | { | |
220 | // makes data from SDigits | |
221 | ||
222 | Int_t cpvmul = 0, premul = 0; | |
2c1131dd | 223 | Float_t edepkev = 0.; |
a9523c3e | 224 | |
225 | TIter next(sdigits) ; | |
226 | AliPMDsdigit * sdigit ; | |
227 | while ( (sdigit = dynamic_cast<AliPMDsdigit *>(next())) ) | |
228 | { | |
229 | if(sdigit->GetDetector() == 0) | |
230 | { | |
2c1131dd | 231 | edepkev = (sdigit->GetCellEdep())/1000.; |
232 | GetSDigitsData(0)->Fill(edepkev); | |
233 | premul++; | |
a9523c3e | 234 | } |
235 | if(sdigit->GetDetector() == 1) | |
236 | { | |
2c1131dd | 237 | edepkev = (sdigit->GetCellEdep())/1000.; |
238 | GetSDigitsData(1)->Fill(edepkev); | |
239 | cpvmul++; | |
a9523c3e | 240 | } |
241 | ||
242 | } | |
243 | if (premul > 0) GetSDigitsData(2)->Fill(premul); | |
244 | if (cpvmul > 0) GetSDigitsData(3)->Fill(cpvmul); | |
245 | ||
246 | } | |
247 | ||
248 | //____________________________________________________________________________ | |
249 | void AliPMDQADataMakerSim::MakeSDigits(TTree * sdigitTree) | |
250 | { | |
251 | // makes data from SDigit Tree | |
252 | ||
253 | TClonesArray * sdigits = new TClonesArray("AliPMDsdigit", 1000) ; | |
254 | ||
255 | TBranch * branch = sdigitTree->GetBranch("PMDSDigit") ; | |
256 | branch->SetAddress(&sdigits) ; | |
257 | ||
258 | if ( ! branch ) | |
259 | { | |
260 | AliWarning("PMD branch in SDigit Tree not found") ; | |
261 | } | |
262 | else | |
263 | { | |
264 | for (Int_t ient = 0; ient < branch->GetEntries(); ient++) | |
265 | { | |
266 | branch->GetEntry(ient) ; | |
267 | MakeSDigits(sdigits) ; | |
268 | } | |
269 | } | |
270 | } | |
271 | ||
272 | //____________________________________________________________________________ | |
273 | void AliPMDQADataMakerSim::MakeDigits(TClonesArray * digits) | |
274 | { | |
275 | // makes data from Digits | |
276 | ||
277 | Int_t cpvmul = 0, premul = 0; | |
278 | ||
279 | TIter next(digits) ; | |
280 | AliPMDdigit * digit ; | |
281 | while ( (digit = dynamic_cast<AliPMDdigit *>(next())) ) | |
282 | { | |
283 | if(digit->GetDetector() == 0) | |
284 | { | |
285 | GetDigitsData(0)->Fill( digit->GetADC()) ; | |
286 | premul++; | |
287 | } | |
288 | if(digit->GetDetector() == 1) | |
289 | { | |
290 | GetDigitsData(1)->Fill( digit->GetADC()); | |
291 | cpvmul++; | |
292 | } | |
293 | } | |
294 | ||
295 | if (premul > 0) GetDigitsData(2)->Fill(premul); | |
296 | if (cpvmul > 0) GetDigitsData(3)->Fill(cpvmul); | |
297 | ||
298 | ||
299 | } | |
300 | ||
301 | //____________________________________________________________________________ | |
302 | void AliPMDQADataMakerSim::MakeDigits(TTree * digitTree) | |
303 | { | |
304 | // makes data from Digit Tree | |
305 | ||
306 | TClonesArray * digits = new TClonesArray("AliPMDdigit", 1000) ; | |
307 | ||
308 | TBranch * branch = digitTree->GetBranch("PMDDigit") ; | |
309 | branch->SetAddress(&digits) ; | |
310 | ||
311 | if ( ! branch ) | |
312 | { | |
313 | AliWarning("PMD branch in Digit Tree not found") ; | |
314 | } | |
315 | else | |
316 | { | |
317 | for (Int_t ient = 0; ient < branch->GetEntries(); ient++) | |
318 | { | |
a9523c3e | 319 | branch->GetEntry(ient) ; |
320 | MakeDigits(digits) ; | |
321 | } | |
322 | ||
323 | } | |
324 | } | |
325 | ||
326 | ||
327 | //____________________________________________________________________________ | |
328 | ||
329 | void AliPMDQADataMakerSim::StartOfDetectorCycle() | |
330 | { | |
331 | //Detector specific actions at start of cycle | |
332 | ||
333 | } | |
334 | //____________________________________________________________________________ | |
335 | ||
336 | void AliPMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list) | |
337 | { | |
338 | //Detector specific actions at end of cycle | |
339 | // do the QA checking | |
340 | AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ; | |
341 | } |