Removed generated files@
[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 #include "AliFMDSDigit.h"
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
43 ClassImp(AliFMDQADataMakerSim)
44 #if 0
45 ; // This line is for Emacs - do not delete!
46 #endif
47 //_____________________________________________________________________
48 AliFMDQADataMakerSim::AliFMDQADataMakerSim() 
49   :  AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kFMD),
50                        "FMD Quality Assurance Data Maker")
51 {
52   // ctor
53
54 }
55
56 //_____________________________________________________________________
57 AliFMDQADataMakerSim::AliFMDQADataMakerSim(const AliFMDQADataMakerSim& /*qadm*/) 
58   : AliQADataMakerSim()
59 {
60   // copy ctor 
61   // 
62   // Parameters: 
63   //    qadm    Object to copy from
64   
65 }
66 //_____________________________________________________________________
67 AliFMDQADataMakerSim& 
68 AliFMDQADataMakerSim::operator = (const AliFMDQADataMakerSim& ) 
69 {
70   
71   return *this;
72 }
73 //_____________________________________________________________________
74 AliFMDQADataMakerSim::~AliFMDQADataMakerSim()
75 {
76
77 }
78
79 //_____________________________________________________________________
80 void AliFMDQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, 
81                                               TObjArray ** list)
82 {
83   //Detector specific actions at end of cycle
84   // do the QA checking
85   AliLog::Message(5,"FMD: end of detector cycle",
86                   "AliFMDQADataMakerSim","AliFMDQADataMakerSim",
87                   "AliFMDQADataMakerSim::EndOfDetectorCycle",
88                   "AliFMDQADataMakerSim.cxx",83);
89   AliQAChecker::Instance()->Run(AliQAv1::kFMD, task, list) ;  
90   
91 }
92 //_____________________________________________________________________
93 void AliFMDQADataMakerSim::InitSDigits()
94 {
95   // create SDigits histograms in SDigits subdir
96   const Bool_t expert   = kTRUE ; 
97   const Bool_t image    = kTRUE ; 
98   
99   TH1I* hADCCounts = new TH1I("hADCCounts",
100                               "Dist of ADC counts;ADC counts;Entries",
101                               1024,0,1024);
102   hADCCounts->SetXTitle("ADC counts");
103   Add2SDigitsList(hADCCounts, 0, !expert, image);
104 }
105
106 //____________________________________________________________________ 
107 void AliFMDQADataMakerSim::InitHits()
108 {
109   // create Digits histograms in Digits subdir
110   const Bool_t expert   = kTRUE ; 
111   const Bool_t image    = kTRUE ; 
112   
113   TH1F* hEnergyOfHits = new TH1F("hEnergyOfHits",
114                                  "Energy distribution;Energy [MeV];Counts",
115                                  100,0,3);
116   hEnergyOfHits->SetXTitle("Edep");
117   hEnergyOfHits->SetYTitle("Counts");
118   Add2HitsList(hEnergyOfHits, 0, !expert, image);
119 }
120
121 //_____________________________________________________________________
122 void AliFMDQADataMakerSim::InitDigits()
123 {
124   // create Digits histograms in Digits subdir
125   const Bool_t expert   = kTRUE ; 
126   const Bool_t image    = kTRUE ; 
127   
128   TH1I* hADCCounts = new TH1I("hADCCounts",
129                               "Dist of ADC counts; ADC counts;Entries",
130                               1024,0,1024);
131   hADCCounts->SetXTitle("ADC counts");
132   Add2DigitsList(hADCCounts, 0, !expert, image);
133 }
134
135 //_____________________________________________________________________
136 void AliFMDQADataMakerSim::MakeHits()
137 {
138   // Check id histograms already created for this Event Specie
139   if ( ! GetHitsData(0) )
140     InitHits() ;
141
142   TIter next(fHitsArray);
143   AliFMDHit * hit;
144   while ((hit = static_cast<AliFMDHit *>(next()))) 
145     GetHitsData(0)->Fill(hit->Edep()/hit->Length()*0.032);
146 }
147
148 //_____________________________________________________________________
149 void AliFMDQADataMakerSim::MakeHits(TTree * hitTree)
150 {
151   // make QA data from Hit Tree
152   // 
153   // Parameters: 
154   //   hitTree    Hits container 
155   //
156   if (!fHitsArray) 
157     fHitsArray = new TClonesArray("AliFMDHit", 1000) ; 
158   fHitsArray->Clear() ; 
159   
160   TBranch * branch = hitTree->GetBranch("FMD") ;
161   if (!branch) {
162     AliWarning("FMD branch in Hit Tree not found") ; 
163     return;
164   }
165     
166   branch->SetAddress(&fHitsArray) ;
167   
168   for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
169     branch->GetEntry(ientry);
170     MakeHits();   //tmp); 
171     fHitsArray->Clear() ; 
172   }     
173 }
174
175 //_____________________________________________________________________
176 void AliFMDQADataMakerSim::MakeDigits()
177 {
178   // makes data from Digits
179   // 
180   // Parameters: 
181   //    none
182   if(!fDigitsArray) return;
183   
184   for(Int_t i = 0 ; i < fDigitsArray->GetEntriesFast() ; i++) {
185     //Raw ADC counts
186     AliFMDDigit* digit = static_cast<AliFMDDigit*>(fDigitsArray->At(i));
187     GetDigitsData(0)->Fill(digit->Counts());
188   }
189 }
190
191 //_____________________________________________________________________
192 void AliFMDQADataMakerSim::MakeDigits(TTree * digitTree)
193 {
194   // Make data from digits. 
195   // 
196   // Parameters: 
197   //    digitTree    Tree holding digits. 
198   
199   if (!fDigitsArray) 
200     fDigitsArray = new TClonesArray("AliFMDDigit", 1000) ; 
201   fDigitsArray->Clear();
202   
203   TBranch * branch = digitTree->GetBranch("FMD") ;
204   if (!branch)    {
205       AliWarning("FMD branch in Digit Tree not found") ; 
206       return;
207   } 
208   branch->SetAddress(&fDigitsArray) ;
209
210   if (fDigitsArray) fDigitsArray->Clear();
211
212   branch->GetEntry(0) ; 
213   MakeDigits() ; 
214 }
215
216 //_____________________________________________________________________
217 void AliFMDQADataMakerSim::MakeSDigits()
218 {
219   // makes data from Digits
220   // 
221   // Parameters: 
222   //   none 
223   if(!fSDigitsArray) return;
224   
225    for(Int_t i = 0 ; i < fSDigitsArray->GetEntriesFast() ; i++) {
226     //Raw ADC counts
227     AliFMDSDigit* sdigit = static_cast<AliFMDSDigit*>(fSDigitsArray->At(i));
228     GetSDigitsData(0)->Fill(sdigit->Counts());
229   }
230 }
231
232 //_____________________________________________________________________
233 void AliFMDQADataMakerSim::MakeSDigits(TTree * sdigitTree)
234 {
235   // Make data from digits. 
236   // 
237   // Parameters: 
238   //    digitTree    Tree holding digits. 
239   
240   if (!fSDigitsArray) 
241     fSDigitsArray = new TClonesArray("AliFMDSDigit", 1000) ; 
242   fSDigitsArray->Clear() ;
243
244   TBranch * branch = sdigitTree->GetBranch("FMD") ;
245   if (!branch)    {
246     AliWarning("FMD branch in SDigit Tree not found") ; 
247     return;
248   } 
249   branch->SetAddress(&fSDigitsArray) ;
250   branch->GetEntry(0) ; 
251   MakeSDigits() ; 
252 }
253
254 //_____________________________________________________________________ 
255 void AliFMDQADataMakerSim::StartOfDetectorCycle()
256 {
257   // Does 
258   // not 
259   // do 
260   // anything 
261 }
262 //_____________________________________________________________________ 
263 //
264 // EOF
265 //