Updated alignment code (Mikolaj)
[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"
86c6cec8 31#include "AliFMDRawReader.h"
3ceaa9ad 32#include "AliRawReader.h"
c9dd1c4d 33
34//_____________________________________________________________________
35// This is the class that collects the QA data for the FMD during
36// reconstruction.
37//
38// The following data types are picked up:
c9dd1c4d 39// - rec points
40// - esd data
c9dd1c4d 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),
a7e41e8d 53 "FMD Quality Assurance Data Maker"),
198942c7 54 fDigitsArray("AliFMDDigit", 0),
55 fRecPointsArray("AliFMDRecPoint", 1000)
c9dd1c4d 56{
57 // ctor
56236ce9 58
c9dd1c4d 59}
60
61//_____________________________________________________________________
a7e41e8d 62AliFMDQADataMakerRec::AliFMDQADataMakerRec(const AliFMDQADataMakerRec& qadm)
63 : AliQADataMakerRec(AliQA::GetDetName(AliQA::kFMD),
64 "FMD Quality Assurance Data Maker"),
65 fDigitsArray(qadm.fDigitsArray),
66 fRecPointsArray(qadm.fRecPointsArray)
c9dd1c4d 67{
ffa78f64 68 // copy ctor
c9dd1c4d 69 // Parameters:
70 // qadm Object to copy from
71
72}
9bd2ccc2 73//_____________________________________________________________________
a7e41e8d 74AliFMDQADataMakerRec& AliFMDQADataMakerRec::operator = (const AliFMDQADataMakerRec& qadm )
75{
76 fDigitsArray = qadm.fDigitsArray;
77 fRecPointsArray = qadm.fRecPointsArray;
78
79 return *this;
80}
81//_____________________________________________________________________
9bd2ccc2 82AliFMDQADataMakerRec::~AliFMDQADataMakerRec()
83{
56236ce9 84
9bd2ccc2 85}
c9dd1c4d 86
87
88//_____________________________________________________________________
89
90void
92a357bf 91AliFMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task,
c9dd1c4d 92 TObjArray * list)
93{
94 // Detector specific actions at end of cycle
95 // do the QA checking
ffa78f64 96 AliLog::Message(5,"FMD: end of detector cycle",
97 "AliFMDQADataMakerRec","AliFMDQADataMakerRec",
98 "AliFMDQADataMakerRec::EndOfDetectorCycle",
99 "AliFMDQADataMakerRec.cxx",95);
c9dd1c4d 100 AliQAChecker::Instance()->Run(AliQA::kFMD, task, list);
101}
102
103//_____________________________________________________________________
104void AliFMDQADataMakerRec::InitESDs()
105{
106 // create Digits histograms in Digits subdir
107 TH1F* hEnergyOfESD = new TH1F("hEnergyOfESD","Energy distribution",100,0,3);
108 hEnergyOfESD->SetXTitle("Edep/Emip");
109 hEnergyOfESD->SetYTitle("Counts");
110 Add2ESDsList(hEnergyOfESD, 0);
111
112}
113
114//_____________________________________________________________________
86c6cec8 115/*void AliFMDQADataMakerRec::InitDigits()
c9dd1c4d 116{
117 // create Digits histograms in Digits subdir
118 TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",
119 1024,0,1024);
c9dd1c4d 120 hADCCounts->SetXTitle("ADC counts");
121 hADCCounts->SetYTitle("");
c9dd1c4d 122 Add2DigitsList(hADCCounts, 0);
c9dd1c4d 123
124}
86c6cec8 125*/
c9dd1c4d 126//_____________________________________________________________________
127void AliFMDQADataMakerRec::InitRecPoints()
128{
129 TH1F* hEnergyOfRecpoints = new TH1F("hEnergyOfRecpoints",
130 "Energy Distribution",100,0,3);
131 hEnergyOfRecpoints->SetXTitle("Edep/Emip");
132 hEnergyOfRecpoints->SetYTitle("");
133 Add2RecPointsList(hEnergyOfRecpoints,0);
134}
135
136//_____________________________________________________________________
137void AliFMDQADataMakerRec::InitRaws()
138{
86c6cec8 139 TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",
198942c7 140 1024,0,1023);
86c6cec8 141 hADCCounts->SetXTitle("ADC counts");
142 hADCCounts->SetYTitle("");
143 Add2RawsList(hADCCounts, 0);
144
c9dd1c4d 145}
146
147//_____________________________________________________________________
148void AliFMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
149{
150 if(!esd) {
151 AliError("FMD ESD object not found!!") ;
152 return;
153 }
154 AliESDFMD* fmd = esd->GetFMDData();
155 if (!fmd) return;
156
157 for(UShort_t det=1;det<=3;det++) {
158 for (UShort_t ir = 0; ir < 2; ir++) {
159 Char_t ring = (ir == 0 ? 'I' : 'O');
160 UShort_t nsec = (ir == 0 ? 20 : 40);
161 UShort_t nstr = (ir == 0 ? 512 : 256);
162 for(UShort_t sec =0; sec < nsec; sec++) {
163 for(UShort_t strip = 0; strip < nstr; strip++) {
164 Float_t mult = fmd->Multiplicity(det,ring,sec,strip);
165 if(mult == AliESDFMD::kInvalidMult) continue;
166
167 GetESDsData(0)->Fill(mult);
168 }
169 }
170 }
171 }
172}
173
86c6cec8 174/*
c9dd1c4d 175//_____________________________________________________________________
176void AliFMDQADataMakerRec::MakeDigits(TClonesArray * digits)
177{
178 // makes data from Digits
179 if(!digits) {
ffa78f64 180 AliError("FMD Digit object not found!!") ;
181 return;
c9dd1c4d 182 }
56236ce9 183 for(Int_t i=0;i<digits->GetEntriesFast();i++) {
c9dd1c4d 184 //Raw ADC counts
56236ce9 185 AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
c9dd1c4d 186 GetDigitsData(0)->Fill(digit->Counts());
187 }
188}
189
190//_____________________________________________________________________
191void AliFMDQADataMakerRec::MakeDigits(TTree * digitTree)
192{
193
56236ce9 194 fDigitsArray.Clear();
ffa78f64 195 TBranch* branch = digitTree->GetBranch("FMD");
196 if (!branch) {
c9dd1c4d 197 AliWarning("FMD branch in Digit Tree not found") ;
198 return;
199 }
56236ce9 200 TClonesArray* digitsAddress = &fDigitsArray;
201 branch->SetAddress(&digitsAddress);
ffa78f64 202 branch->GetEntry(0);
56236ce9 203 MakeDigits(digitsAddress);
c9dd1c4d 204}
86c6cec8 205*/
c9dd1c4d 206//_____________________________________________________________________
86c6cec8 207void AliFMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
c9dd1c4d 208{
198942c7 209
86c6cec8 210 AliFMDRawReader fmdReader(rawReader,0);
211 TClonesArray* digitsAddress = &fDigitsArray;
3ceaa9ad 212
213 rawReader->Reset();
214
215 while(rawReader->NextEvent()) {
216
198942c7 217 digitsAddress->Clear();
218 fmdReader.ReadAdcs(digitsAddress);
219 for(Int_t i=0;i<digitsAddress->GetEntriesFast();i++) {
3ceaa9ad 220 //Raw ADC counts
198942c7 221 AliFMDDigit* digit = static_cast<AliFMDDigit*>(digitsAddress->At(i));
3ceaa9ad 222 GetRawsData(0)->Fill(digit->Counts());
223 }
86c6cec8 224 }
c9dd1c4d 225
226}
227
228//_____________________________________________________________________
229void AliFMDQADataMakerRec::MakeRecPoints(TTree* clustersTree)
230{
231 // makes data from RecPoints
198942c7 232 AliFMDParameters* pars = AliFMDParameters::Instance();
56236ce9 233 fRecPointsArray.Clear();
c9dd1c4d 234 TBranch *fmdbranch = clustersTree->GetBranch("FMD");
235 if (!fmdbranch) {
236 AliError("can't get the branch with the FMD recpoints !");
237 return;
238 }
239
56236ce9 240 TClonesArray* RecPointsAddress = &fRecPointsArray;
9bd2ccc2 241
56236ce9 242 fmdbranch->SetAddress(&RecPointsAddress);
c9dd1c4d 243 fmdbranch->GetEntry(0);
56236ce9 244 TIter next(RecPointsAddress) ;
c9dd1c4d 245 AliFMDRecPoint * rp ;
246 while ((rp = static_cast<AliFMDRecPoint*>(next()))) {
198942c7 247
248 GetRecPointsData(0)->Fill(rp->Edep()/pars->GetEdepMip()) ;
249
c9dd1c4d 250 }
c9dd1c4d 251
252}
253
254//_____________________________________________________________________
255void AliFMDQADataMakerRec::StartOfDetectorCycle()
256{
257 // What
258 // to
259 // do?
260}
261//_____________________________________________________________________
262//
263// EOF
264//