1 /**************************************************************************
2 * Copyright(c) 2004, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 ---
17 #include <TClonesArray.h>
22 // --- AliRoot header files ---
23 #include "AliESDEvent.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 #include "AliFMDRawReader.h"
32 #include "AliRawReader.h"
33 #include "AliFMDAltroMapping.h"
34 #include "AliFMDDebug.h"
36 //_____________________________________________________________________
37 // This is the class that collects the QA data for the FMD during
40 // The following data types are picked up:
44 // Author : Hans Hjersing Dalsgaard, hans.dalsgaard@cern.ch
45 //_____________________________________________________________________
47 ClassImp(AliFMDQADataMakerRec)
49 ; // For Emacs - do not delete!
52 //_____________________________________________________________________
53 AliFMDQADataMakerRec::AliFMDQADataMakerRec() :
54 AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kFMD),
55 "FMD Quality Assurance Data Maker"),
56 fRecPointsArray("AliFMDRecPoint", 1000)
62 //_____________________________________________________________________
63 AliFMDQADataMakerRec::AliFMDQADataMakerRec(const AliFMDQADataMakerRec& qadm)
64 : AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kFMD),
65 "FMD Quality Assurance Data Maker"),
66 fRecPointsArray(qadm.fRecPointsArray)
70 // qadm Object to copy from
73 //_____________________________________________________________________
75 AliFMDQADataMakerRec::operator = (const AliFMDQADataMakerRec& qadm )
77 fRecPointsArray = qadm.fRecPointsArray;
81 //_____________________________________________________________________
82 AliFMDQADataMakerRec::~AliFMDQADataMakerRec()
88 //_____________________________________________________________________
91 AliFMDQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task,
94 // Detector specific actions at end of cycle
96 AliLog::Message(5,"FMD: end of detector cycle",
97 "AliFMDQADataMakerRec","AliFMDQADataMakerRec",
98 "AliFMDQADataMakerRec::EndOfDetectorCycle",
99 "AliFMDQADataMakerRec.cxx",95);
100 AliQAChecker::Instance()->Run(AliQAv1::kFMD, task, list);
103 //_____________________________________________________________________
104 void AliFMDQADataMakerRec::InitESDs()
106 // create Digits histograms in Digits subdir
107 const Bool_t expert = kTRUE ;
108 const Bool_t image = kTRUE ;
110 TH1F* hEnergyOfESD = new TH1F("hEnergyOfESD","Energy distribution",100,0,3);
111 hEnergyOfESD->SetXTitle("Edep/Emip");
112 hEnergyOfESD->SetYTitle("Counts");
113 Add2ESDsList(hEnergyOfESD, 0, !expert, image);
117 //_____________________________________________________________________
118 void AliFMDQADataMakerRec::InitDigits()
120 // create Digits histograms in Digits subdir
121 const Bool_t expert = kTRUE ;
122 const Bool_t image = kTRUE ;
124 TH1I* hADCCounts = new TH1I("hADCCounts",
125 "Dist of ADC counts;ADC counts;Counts",
127 Add2DigitsList(hADCCounts, 0, !expert, image);
131 //_____________________________________________________________________
132 void AliFMDQADataMakerRec::InitRecPoints()
134 // create Reconstructed Points histograms in RecPoints subdir
135 const Bool_t expert = kTRUE ;
136 const Bool_t image = kTRUE ;
138 TH1F* hEnergyOfRecpoints = new TH1F("hEnergyOfRecpoints",
139 "Energy Distribution",100,0,3);
140 hEnergyOfRecpoints->SetXTitle("Edep/Emip");
141 hEnergyOfRecpoints->SetYTitle("Counts");
142 Add2RecPointsList(hEnergyOfRecpoints,0, !expert, image);
145 //_____________________________________________________________________
146 void AliFMDQADataMakerRec::InitRaws()
148 // create Raws histograms in Raws subdir
149 const Bool_t expert = kTRUE ;
150 const Bool_t saveCorr = kTRUE ;
151 const Bool_t image = kTRUE ;
152 Int_t colors[3] = {kRed,kGreen,kBlue};
154 for(Int_t det = 1; det<=3; det++) {
155 Int_t firstring = (det==1 ? 1 : 0);
156 for(Int_t iring = firstring;iring<=1;iring++) {
157 Char_t ring = (iring == 1 ? 'I' : 'O');
159 new TH1I(Form("hADCCounts_FMD%d%c", det, ring),
160 Form("FMD%d%c ADC counts;Amplitude [ADC counts];Counts",
161 det,ring),1024,0,1023);
162 hADCCounts->SetFillStyle(3001);
163 hADCCounts->SetFillColor(colors[det-1]+2+iring);
164 hADCCounts->SetDrawOption("LOGY");
166 Int_t index1 = GetHalfringIndex(det, ring, 0,1);
167 Add2RawsList(hADCCounts, index1, !expert, image, !saveCorr);
169 for(Int_t b = 0; b<=1;b++) {
171 //Hexadecimal board numbers 0x0, 0x1, 0x10, 0x11;
172 UInt_t board = (iring == 1 ? 0 : 1);
173 board = board + b*16;
177 new TH1I(Form("hADCCounts_FMD%d%c_0x%02x", det, ring, board),
178 Form("FMD%d%c_0x%x ADC counts;Amplitude [ADC counts];Counts",
179 det,ring, board),1024,0,1023);
180 hADCCounts->SetXTitle("ADC counts");
181 hADCCounts->SetYTitle("");
182 hADCCounts->SetFillStyle(3001);
183 hADCCounts->SetFillColor(colors[det-1]+2+iring);
184 hADCCounts->SetDrawOption("LOGY");
185 Int_t index2 = GetHalfringIndex(det, ring, board/16,0);
186 Add2RawsList(hADCCounts, index2, expert, !image, !saveCorr);
194 struct FillESDHist : public AliESDFMD::ForOne
196 FillESDHist(AliFMDQADataMakerRec* m) : fM(m) {}
197 FillESDHist(const FillESDHist& o) : fM(o.fM) {}
198 FillESDHist& operator=(const FillESDHist& o) { fM = o.fM; return *this; }
199 Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, Float_t m, Float_t)
201 // Float_t mult = fmd->Multiplicity(det,ring,sec,strip);
202 if(m == AliESDFMD::kInvalidMult) return true;
204 fM->GetESDsData(0)->Fill(m);
207 AliFMDQADataMakerRec* fM;
211 //_____________________________________________________________________
212 void AliFMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
215 AliError("FMD ESD object not found!!") ;
218 AliFMDDebug(2, ("Will loop over ESD data and fill histogram"));
220 AliESDFMD* fmd = esd->GetFMDData();
230 // FIXME - we should use AliESDFMD::ForOne subclass to do this!
231 for(UShort_t det=1;det<=3;det++) {
232 UShort_t nrng = (det == 1 ? 1 : 2);
233 for (UShort_t ir = 0; ir < nrng; ir++) {
234 Char_t ring = (ir == 0 ? 'I' : 'O');
235 UShort_t nsec = (ir == 0 ? 20 : 40);
236 UShort_t nstr = (ir == 0 ? 512 : 256);
237 for(UShort_t sec =0; sec < nsec; sec++) {
238 for(UShort_t strip = 0; strip < nstr; strip++) {
239 Float_t mult = fmd->Multiplicity(det,ring,sec,strip);
240 if(mult == AliESDFMD::kInvalidMult) continue;
242 GetESDsData(0)->Fill(mult);
251 //_____________________________________________________________________
252 void AliFMDQADataMakerRec::MakeDigits()
254 // makes data from Digits
256 AliError("FMD Digit object not found!!") ;
260 for(Int_t i=0;i<fDigitsArray->GetEntriesFast();i++) {
262 AliFMDDigit* digit = static_cast<AliFMDDigit*>(fDigitsArray->At(i));
263 GetDigitsData(0)->Fill(digit->Counts());
267 //_____________________________________________________________________
268 void AliFMDQADataMakerRec::MakeDigits(TTree * digitTree)
272 fDigitsArray->Clear();
274 fDigitsArray = new TClonesArray("AliFMDDigit", 1000);
276 TBranch* branch = digitTree->GetBranch("FMD");
278 AliWarning("FMD branch in Digit Tree not found") ;
281 branch->SetAddress(&fDigitsArray);
286 //_____________________________________________________________________
287 void AliFMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
290 AliFMDRawReader fmdReader(rawReader,0);
293 fDigitsArray->Clear();
295 fDigitsArray = new TClonesArray("AliFMDDigit", 1000);
297 TClonesArray* digitsAddress = fDigitsArray;
301 digitsAddress->Clear();
302 fmdReader.ReadAdcs(digitsAddress);
303 for(Int_t i=0;i<digitsAddress->GetEntriesFast();i++) {
305 AliFMDDigit* digit = static_cast<AliFMDDigit*>(digitsAddress->At(i));
306 UShort_t det = digit->Detector();
307 Char_t ring = digit->Ring();
308 UShort_t sec = digit->Sector();
309 // UShort_t strip = digit->Strip();
310 AliFMDParameters* pars = AliFMDParameters::Instance();
311 Short_t board = pars->GetAltroMap()->Sector2Board(ring, sec);
313 Int_t index1 = GetHalfringIndex(det, ring, 0, 1);
314 GetRawsData(index1)->Fill(digit->Counts());
315 Int_t index2 = GetHalfringIndex(det, ring, board/16,0);
316 GetRawsData(index2)->Fill(digit->Counts());
321 //_____________________________________________________________________
322 void AliFMDQADataMakerRec::MakeRecPoints(TTree* clustersTree)
324 // makes data from RecPoints
326 AliFMDParameters* pars = AliFMDParameters::Instance();
327 fRecPointsArray.Clear();
328 TBranch *fmdbranch = clustersTree->GetBranch("FMD");
330 AliError("can't get the branch with the FMD recpoints !");
334 TClonesArray* RecPointsAddress = &fRecPointsArray;
336 fmdbranch->SetAddress(&RecPointsAddress);
337 fmdbranch->GetEntry(0);
338 TIter next(RecPointsAddress) ;
339 AliFMDRecPoint * rp ;
340 while ((rp = static_cast<AliFMDRecPoint*>(next()))) {
342 GetRecPointsData(0)->Fill(rp->Edep()/pars->GetEdepMip()) ;
348 //_____________________________________________________________________
349 void AliFMDQADataMakerRec::StartOfDetectorCycle()
355 //_____________________________________________________________________
356 Int_t AliFMDQADataMakerRec::GetHalfringIndex(UShort_t det,
361 UShort_t iring = (ring == 'I' ? 1 : 0);
363 Int_t index = ( ((det-1) << 3) | (iring << 2) | (board << 1) | (monitor << 0));
369 //_____________________________________________________________________