Removal of the simulated CDH which is not present in the CTP raw-data payload. This...
[u/mrichter/AliRoot.git] / FMD / AliFMDQADataMakerRec.cxx
CommitLineData
c9dd1c4d 1/**************************************************************************
2 * Copyright(c) 2004, 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// --- 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 "AliFMDQADataMakerRec.h"
26#include "AliFMDDigit.h"
27#include "AliFMDRecPoint.h"
28#include "AliQAChecker.h"
29#include "AliESDFMD.h"
30#include "AliFMDParameters.h"
31
32//_____________________________________________________________________
33// This is the class that collects the QA data for the FMD during
34// reconstruction.
35//
36// The following data types are picked up:
37// - digits
38// - rec points
39// - esd data
40// The following data types are not supported (yet):
41// - raws
42// Author : Hans Hjersing Dalsgaard, hans.dalsgaard@cern.ch
43//_____________________________________________________________________
44
45ClassImp(AliFMDQADataMakerRec)
46#if 0
47; // For Emacs - do not delete!
48#endif
49
50//_____________________________________________________________________
51AliFMDQADataMakerRec::AliFMDQADataMakerRec() :
52 AliQADataMakerRec(AliQA::GetDetName(AliQA::kFMD),
53 "FMD Quality Assurance Data Maker")
54{
55 // ctor
56}
57
58//_____________________________________________________________________
c61a7285 59AliFMDQADataMakerRec::AliFMDQADataMakerRec(const AliFMDQADataMakerRec& /*qadm*/)
ffa78f64 60 : AliQADataMakerRec()
c9dd1c4d 61{
ffa78f64 62 // copy ctor
c9dd1c4d 63 // Parameters:
64 // qadm Object to copy from
65
66}
67
68
69//_____________________________________________________________________
70
71void
92a357bf 72AliFMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task,
c9dd1c4d 73 TObjArray * list)
74{
75 // Detector specific actions at end of cycle
76 // do the QA checking
ffa78f64 77 AliLog::Message(5,"FMD: end of detector cycle",
78 "AliFMDQADataMakerRec","AliFMDQADataMakerRec",
79 "AliFMDQADataMakerRec::EndOfDetectorCycle",
80 "AliFMDQADataMakerRec.cxx",95);
c9dd1c4d 81 AliQAChecker::Instance()->Run(AliQA::kFMD, task, list);
82}
83
84//_____________________________________________________________________
85void AliFMDQADataMakerRec::InitESDs()
86{
87 // create Digits histograms in Digits subdir
88 TH1F* hEnergyOfESD = new TH1F("hEnergyOfESD","Energy distribution",100,0,3);
89 hEnergyOfESD->SetXTitle("Edep/Emip");
90 hEnergyOfESD->SetYTitle("Counts");
91 Add2ESDsList(hEnergyOfESD, 0);
92
93}
94
95//_____________________________________________________________________
96void AliFMDQADataMakerRec::InitDigits()
97{
98 // create Digits histograms in Digits subdir
99 TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",
100 1024,0,1024);
c9dd1c4d 101 hADCCounts->SetXTitle("ADC counts");
102 hADCCounts->SetYTitle("");
c9dd1c4d 103 Add2DigitsList(hADCCounts, 0);
c9dd1c4d 104
105}
106
107//_____________________________________________________________________
108void AliFMDQADataMakerRec::InitRecPoints()
109{
110 TH1F* hEnergyOfRecpoints = new TH1F("hEnergyOfRecpoints",
111 "Energy Distribution",100,0,3);
112 hEnergyOfRecpoints->SetXTitle("Edep/Emip");
113 hEnergyOfRecpoints->SetYTitle("");
114 Add2RecPointsList(hEnergyOfRecpoints,0);
115}
116
117//_____________________________________________________________________
118void AliFMDQADataMakerRec::InitRaws()
119{
120
121}
122
123//_____________________________________________________________________
124void AliFMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
125{
126 if(!esd) {
127 AliError("FMD ESD object not found!!") ;
128 return;
129 }
130 AliESDFMD* fmd = esd->GetFMDData();
131 if (!fmd) return;
132
133 for(UShort_t det=1;det<=3;det++) {
134 for (UShort_t ir = 0; ir < 2; ir++) {
135 Char_t ring = (ir == 0 ? 'I' : 'O');
136 UShort_t nsec = (ir == 0 ? 20 : 40);
137 UShort_t nstr = (ir == 0 ? 512 : 256);
138 for(UShort_t sec =0; sec < nsec; sec++) {
139 for(UShort_t strip = 0; strip < nstr; strip++) {
140 Float_t mult = fmd->Multiplicity(det,ring,sec,strip);
141 if(mult == AliESDFMD::kInvalidMult) continue;
142
143 GetESDsData(0)->Fill(mult);
144 }
145 }
146 }
147 }
148}
149
150//_____________________________________________________________________
151void AliFMDQADataMakerRec::MakeDigits(TClonesArray * digits)
152{
153 // makes data from Digits
154 if(!digits) {
ffa78f64 155 AliError("FMD Digit object not found!!") ;
156 return;
c9dd1c4d 157 }
158 for(Int_t i=0;i<digits->GetEntries();i++) {
159 //Raw ADC counts
160 AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
161 GetDigitsData(0)->Fill(digit->Counts());
162 }
163}
164
165//_____________________________________________________________________
166void AliFMDQADataMakerRec::MakeDigits(TTree * digitTree)
167{
168
ffa78f64 169 TClonesArray* digits = new TClonesArray("AliFMDDigit", 1000);
170 TBranch* branch = digitTree->GetBranch("FMD");
171 if (!branch) {
c9dd1c4d 172 AliWarning("FMD branch in Digit Tree not found") ;
173 return;
174 }
175
ffa78f64 176 branch->SetAddress(&digits);
177 branch->GetEntry(0);
178 MakeDigits(digits);
c9dd1c4d 179}
180
181//_____________________________________________________________________
c61a7285 182void AliFMDQADataMakerRec::MakeRaws(AliRawReader* /*rawReader*/)
c9dd1c4d 183{
184
185}
186
187//_____________________________________________________________________
188void AliFMDQADataMakerRec::MakeRecPoints(TTree* clustersTree)
189{
190 // makes data from RecPoints
191 TBranch *fmdbranch = clustersTree->GetBranch("FMD");
192 if (!fmdbranch) {
193 AliError("can't get the branch with the FMD recpoints !");
194 return;
195 }
196
ffa78f64 197 TClonesArray * fmdrecpoints = new TClonesArray("AliFMDRecPoint", 1000);
c9dd1c4d 198 fmdbranch->SetAddress(&fmdrecpoints);
199 fmdbranch->GetEntry(0);
200
201 TIter next(fmdrecpoints) ;
202 AliFMDRecPoint * rp ;
203 while ((rp = static_cast<AliFMDRecPoint*>(next()))) {
204 GetRecPointsData(0)->Fill(rp->Particles()) ;
205 }
206 fmdrecpoints->Delete();
207 delete fmdrecpoints;
208
209}
210
211//_____________________________________________________________________
212void AliFMDQADataMakerRec::StartOfDetectorCycle()
213{
214 // What
215 // to
216 // do?
217}
218//_____________________________________________________________________
219//
220// EOF
221//