Introducing event specie in QA (Yves)
[u/mrichter/AliRoot.git] / FMD / AliFMDQADataMakerSim.cxx
CommitLineData
c9dd1c4d 1/**************************************************************************
ffa78f64 2 * Copyright(c) 2004, ALICE Experiment at CERN, All rights reserved. *
c9dd1c4d 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// --- ROOT system ---
16#include <iostream>
17#include <TClonesArray.h>
18#include <TFile.h>
19#include <TH1F.h>
20#include <TH1I.h>
21
22// --- AliRoot header files ---
23#include "AliESDEvent.h"
24#include "AliLog.h"
25#include "AliFMDQADataMakerSim.h"
26#include "AliFMDDigit.h"
27#include "AliFMDHit.h"
28#include "AliQAChecker.h"
29#include "AliFMDParameters.h"
028cb80b 30#include "AliFMDSDigit.h"
c9dd1c4d 31
32//_____________________________________________________________________
33// This is the class that collects the QA data for the FMD during simulation.
34// The following data types are picked up:
35// - hits
36// - digits
37// The following data types are not supported (yet):
38// - raws
39// - sdigits
40// Author : Hans Hjersing Dalsgaard, Niels Bohr Institute, hans.dalsgaard@cern.ch
41//_____________________________________________________________________
42
43ClassImp(AliFMDQADataMakerSim)
ffa78f64 44#if 0
45; // This line is for Emacs - do not delete!
46#endif
c9dd1c4d 47//_____________________________________________________________________
ffa78f64 48AliFMDQADataMakerSim::AliFMDQADataMakerSim()
49 : AliQADataMakerSim(AliQA::GetDetName(AliQA::kFMD),
9bd2ccc2 50 "FMD Quality Assurance Data Maker"),
028cb80b 51 fSDigitsArray("AliFMDSDigit", 1000),
56236ce9 52 fDigitsArray("AliFMDDigit", 1000),
53 fHitsArray("AliFMDHit", 10)
c9dd1c4d 54{
55 // ctor
56236ce9 56
c9dd1c4d 57}
58
59//_____________________________________________________________________
a7e41e8d 60AliFMDQADataMakerSim::AliFMDQADataMakerSim(const AliFMDQADataMakerSim& qadm)
61 : AliQADataMakerSim(),
028cb80b 62 fSDigitsArray(qadm.fSDigitsArray),
a7e41e8d 63 fDigitsArray(qadm.fDigitsArray),
64 fHitsArray(qadm.fHitsArray)
c9dd1c4d 65{
66 //copy ctor
a7e41e8d 67
c9dd1c4d 68 // Parameters:
69 // qadm Object to copy from
70
71}
9bd2ccc2 72//_____________________________________________________________________
a7e41e8d 73AliFMDQADataMakerSim& AliFMDQADataMakerSim::operator = (const AliFMDQADataMakerSim& qadm )
74{
028cb80b 75 fSDigitsArray = qadm.fSDigitsArray;
a7e41e8d 76 fDigitsArray = qadm.fDigitsArray;
77 fHitsArray = qadm.fHitsArray;
78
79 return *this;
80}
81//_____________________________________________________________________
9bd2ccc2 82AliFMDQADataMakerSim::~AliFMDQADataMakerSim()
83{
56236ce9 84
9bd2ccc2 85}
c9dd1c4d 86
87//_____________________________________________________________________
92a357bf 88void AliFMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task,
57acd2d2 89 TObjArray ** list)
c9dd1c4d 90{
91 //Detector specific actions at end of cycle
92 // do the QA checking
ffa78f64 93 AliLog::Message(5,"FMD: end of detector cycle",
94 "AliFMDQADataMakerSim","AliFMDQADataMakerSim",
95 "AliFMDQADataMakerSim::EndOfDetectorCycle",
96 "AliFMDQADataMakerSim.cxx",83);
c9dd1c4d 97 AliQAChecker::Instance()->Run(AliQA::kFMD, task, list) ;
98
99}
028cb80b 100//_____________________________________________________________________
101void AliFMDQADataMakerSim::InitSDigits()
102{
103 // create SDigits histograms in SDigits subdir
104 TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",1024,0,1024);
105 hADCCounts->SetXTitle("ADC counts");
106 Add2SDigitsList(hADCCounts, 0);
107}
c9dd1c4d 108
109//____________________________________________________________________
110void AliFMDQADataMakerSim::InitHits()
111{
112 // create Digits histograms in Digits subdir
ffa78f64 113 TH1F* hEnergyOfHits = new TH1F("hEnergyOfHits","Energy distribution",100,0,3);
c9dd1c4d 114 hEnergyOfHits->SetXTitle("Edep");
115 hEnergyOfHits->SetYTitle("Counts");
116 Add2HitsList(hEnergyOfHits, 0);
117}
118
119//_____________________________________________________________________
120void AliFMDQADataMakerSim::InitDigits()
121{
122 // create Digits histograms in Digits subdir
123 TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",1024,0,1024);
c9dd1c4d 124 hADCCounts->SetXTitle("ADC counts");
c9dd1c4d 125 Add2DigitsList(hADCCounts, 0);
c9dd1c4d 126}
127
128//_____________________________________________________________________
129void AliFMDQADataMakerSim::MakeHits(TClonesArray * hits)
130{
a7e41e8d 131 TIter next(hits);
c9dd1c4d 132 AliFMDHit * hit;
133 while ((hit = static_cast<AliFMDHit *>(next())))
134 GetHitsData(0)->Fill(hit->Edep()/hit->Length()*0.032);
135}
136
137//_____________________________________________________________________
138void AliFMDQADataMakerSim::MakeHits(TTree * hitTree)
139{
140 // make QA data from Hit Tree
141
56236ce9 142 fHitsArray.Clear();
143
c9dd1c4d 144 TBranch * branch = hitTree->GetBranch("FMD") ;
145 if (!branch) {
146 AliWarning("FMD branch in Hit Tree not found") ;
147 return;
148 }
56236ce9 149
150 TClonesArray* hitsAddress = &fHitsArray;
151
152 branch->SetAddress(&hitsAddress) ;
c9dd1c4d 153
154 for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
155 branch->GetEntry(ientry);
56236ce9 156 MakeHits(hitsAddress); //tmp);
c9dd1c4d 157 }
9bd2ccc2 158
c9dd1c4d 159}
160
161//_____________________________________________________________________
162void AliFMDQADataMakerSim::MakeDigits(TClonesArray * digits)
163{
164 // makes data from Digits
165 if(!digits) return;
166
56236ce9 167 for(Int_t i = 0 ; i < digits->GetEntriesFast() ; i++) {
c9dd1c4d 168 //Raw ADC counts
169 AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
170 GetDigitsData(0)->Fill(digit->Counts());
171 }
172}
173
174//_____________________________________________________________________
175void AliFMDQADataMakerSim::MakeDigits(TTree * digitTree)
176{
177
56236ce9 178 fDigitsArray.Clear();
c9dd1c4d 179 TBranch * branch = digitTree->GetBranch("FMD") ;
180 if (!branch) {
181 AliWarning("FMD branch in Digit Tree not found") ;
182 return;
183 }
56236ce9 184 TClonesArray* digitAddress = &fDigitsArray;
185 branch->SetAddress(&digitAddress) ;
c9dd1c4d 186 branch->GetEntry(0) ;
56236ce9 187 MakeDigits(digitAddress) ;
c9dd1c4d 188}
189
028cb80b 190//_____________________________________________________________________
191void AliFMDQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
192{
193 // makes data from Digits
194 if(!sdigits) return;
195
196 for(Int_t i = 0 ; i < sdigits->GetEntriesFast() ; i++) {
197 //Raw ADC counts
198 AliFMDSDigit* sdigit = static_cast<AliFMDSDigit*>(sdigits->At(i));
199 GetSDigitsData(0)->Fill(sdigit->Counts());
200 }
201}
202
203//_____________________________________________________________________
204void AliFMDQADataMakerSim::MakeSDigits(TTree * sdigitTree)
205{
206
207 fSDigitsArray.Clear();
208 TBranch * branch = sdigitTree->GetBranch("FMD") ;
209 if (!branch) {
210 AliWarning("FMD branch in SDigit Tree not found") ;
211 return;
212 }
213 TClonesArray* sdigitAddress = &fSDigitsArray;
214 branch->SetAddress(&sdigitAddress) ;
215 branch->GetEntry(0) ;
216 MakeSDigits(sdigitAddress) ;
217}
218
c9dd1c4d 219//_____________________________________________________________________
220void AliFMDQADataMakerSim::StartOfDetectorCycle()
221{
222
223}
224//_____________________________________________________________________
225//
226// EOF
227//