Corrected code to get rid of compilation warnings
[u/mrichter/AliRoot.git] / FMD / AliFMDQADataMakerSim.cxx
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 "AliFMDQADataMakerSim.h"
26 #include "AliFMDDigit.h"
27 #include "AliFMDHit.h"
28 #include "AliQAChecker.h"
29 #include "AliFMDParameters.h"
30
31 //_____________________________________________________________________
32 // This is the class that collects the QA data for the FMD during simulation.
33 // The following data types are picked up:
34 // - hits
35 // - digits
36 // The following data types are not supported (yet):
37 // - raws
38 // - sdigits
39 // Author : Hans Hjersing Dalsgaard, Niels Bohr Institute, hans.dalsgaard@cern.ch
40 //_____________________________________________________________________
41
42 ClassImp(AliFMDQADataMakerSim)
43 #if 0
44 ; // This line is for Emacs - do not delete!
45 #endif
46 //_____________________________________________________________________
47 AliFMDQADataMakerSim::AliFMDQADataMakerSim() 
48   :  AliQADataMakerSim(AliQA::GetDetName(AliQA::kFMD),
49                        "FMD Quality Assurance Data Maker"),
50      fDigitsArray(0),
51      fHitsArray(0)
52 {
53   // ctor
54   fDigitsArray = new TClonesArray("AliFMDDigit", 1000) ; 
55   fHitsArray   = new TClonesArray("AliFMDHit", 10);
56 }
57
58 //_____________________________________________________________________
59 AliFMDQADataMakerSim::AliFMDQADataMakerSim(const AliFMDQADataMakerSim& qadm) 
60   : AliQADataMakerSim(),
61     fDigitsArray(qadm.fDigitsArray),
62     fHitsArray(qadm.fHitsArray)
63 {
64   //copy ctor 
65   
66   // Parameters: 
67   //    qadm    Object to copy from
68   
69 }
70 //_____________________________________________________________________
71 AliFMDQADataMakerSim& AliFMDQADataMakerSim::operator = (const AliFMDQADataMakerSim& qadm ) 
72 {
73   fDigitsArray = qadm.fDigitsArray;
74   fHitsArray = qadm.fHitsArray;
75   
76   return *this;
77 }
78 //_____________________________________________________________________
79 AliFMDQADataMakerSim::~AliFMDQADataMakerSim()
80 {
81   delete fDigitsArray;
82   delete fHitsArray;
83 }
84
85 //_____________________________________________________________________
86 void AliFMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, 
87                                               TObjArray * list)
88 {
89   //Detector specific actions at end of cycle
90   // do the QA checking
91   AliLog::Message(5,"FMD: end of detector cycle",
92                   "AliFMDQADataMakerSim","AliFMDQADataMakerSim",
93                   "AliFMDQADataMakerSim::EndOfDetectorCycle",
94                   "AliFMDQADataMakerSim.cxx",83);
95   AliQAChecker::Instance()->Run(AliQA::kFMD, task, list) ;  
96   
97 }
98
99 //____________________________________________________________________ 
100 void AliFMDQADataMakerSim::InitHits()
101 {
102   // create Digits histograms in Digits subdir
103   TH1F* hEnergyOfHits = new TH1F("hEnergyOfHits","Energy distribution",100,0,3);
104   hEnergyOfHits->SetXTitle("Edep");
105   hEnergyOfHits->SetYTitle("Counts");
106   Add2HitsList(hEnergyOfHits, 0);
107 }
108
109 //_____________________________________________________________________
110 void AliFMDQADataMakerSim::InitDigits()
111 {
112   // create Digits histograms in Digits subdir
113   TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",1024,0,1024);
114   hADCCounts->SetXTitle("ADC counts");
115   Add2DigitsList(hADCCounts, 0);
116 }
117
118 //_____________________________________________________________________
119 void AliFMDQADataMakerSim::MakeHits(TClonesArray * hits)
120 {
121   TIter next(hits);
122   AliFMDHit * hit;
123   while ((hit = static_cast<AliFMDHit *>(next()))) 
124     GetHitsData(0)->Fill(hit->Edep()/hit->Length()*0.032);
125 }
126
127 //_____________________________________________________________________
128 void AliFMDQADataMakerSim::MakeHits(TTree * hitTree)
129 {
130   // make QA data from Hit Tree
131   
132   TBranch * branch = hitTree->GetBranch("FMD") ;
133   if (!branch) {
134     AliWarning("FMD branch in Hit Tree not found") ; 
135     return;
136   }
137   fHitsArray->Clear();
138   branch->SetAddress(&fHitsArray) ;
139   
140   for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
141     branch->GetEntry(ientry);
142     MakeHits(fHitsArray);   //tmp); 
143   }     
144   
145 }
146
147 //_____________________________________________________________________
148 void AliFMDQADataMakerSim::MakeDigits(TClonesArray * digits)
149 {
150   // makes data from Digits
151   if(!digits) return;
152
153   for(Int_t i = 0 ; i < fDigitsArray->GetEntriesFast() ; i++) {
154     //Raw ADC counts
155     AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
156     GetDigitsData(0)->Fill(digit->Counts());
157   }
158 }
159
160 //_____________________________________________________________________
161 void AliFMDQADataMakerSim::MakeDigits(TTree * digitTree)
162 {
163   
164   fDigitsArray->Clear();
165   TBranch * branch = digitTree->GetBranch("FMD") ;
166   if (!branch)    {
167       AliWarning("FMD branch in Digit Tree not found") ; 
168       return;
169   } 
170   branch->SetAddress(&fDigitsArray) ;
171   branch->GetEntry(0) ; 
172   MakeDigits(fDigitsArray) ; 
173 }
174
175 //_____________________________________________________________________ 
176 void AliFMDQADataMakerSim::StartOfDetectorCycle()
177 {
178    
179 }
180 //_____________________________________________________________________ 
181 //
182 // EOF
183 //