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