Fixes for bug #49914: Compilation breaks in trunk, and bug #48629: Trunk cannot read...
[u/mrichter/AliRoot.git] / PMD / AliPMDQADataMakerSim.cxx
CommitLineData
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
42ClassImp(AliPMDQADataMakerSim)
43
44//____________________________________________________________________________
45AliPMDQADataMakerSim::AliPMDQADataMakerSim() :
4e25ac79 46 AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kPMD), "PMD Quality Assurance Data Maker")
a9523c3e 47{
48 // ctor
49}
50
51//____________________________________________________________________________
52AliPMDQADataMakerSim::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//__________________________________________________________________
61AliPMDQADataMakerSim& AliPMDQADataMakerSim::operator = (const AliPMDQADataMakerSim& qadm )
62{
63 // Assign operator.
64 this->~AliPMDQADataMakerSim();
65 new(this) AliPMDQADataMakerSim(qadm);
66 return *this;
67}
68
69//____________________________________________________________________________
70void 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//____________________________________________________________________________
92void 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//____________________________________________________________________________
115void 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//____________________________________________________________________________
138void 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//____________________________________________________________________________
184void 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//____________________________________________________________________________
218void 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//____________________________________________________________________________
249void 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//____________________________________________________________________________
273void 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//____________________________________________________________________________
302void 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
329void AliPMDQADataMakerSim::StartOfDetectorCycle()
330{
331 //Detector specific actions at start of cycle
332
333}
334//____________________________________________________________________________
335
4e25ac79 336void AliPMDQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
a9523c3e 337{
338 //Detector specific actions at end of cycle
339 // do the QA checking
4e25ac79 340 AliQAChecker::Instance()->Run(AliQAv1::kPMD, task, list) ;
a9523c3e 341}