Make and print an image of QA user flagged histograms (Yves)
[u/mrichter/AliRoot.git] / MUON / AliMUONQADataMakerSim.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, 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
16 // $Id$
17
18 #include "AliMUONQADataMakerSim.h"
19 #include "AliMUONHit.h"  
20 #include "AliMUONDigit.h"  
21 #include "AliMUONVHitStore.h"
22 #include "AliMUONVDigitStore.h"
23
24 // --- AliRoot header files ---
25 #include "AliLog.h"
26 #include "AliQAChecker.h"
27
28 // --- ROOT system ---
29 #include <TClonesArray.h>
30 #include <TFile.h> 
31 #include <TH1F.h> 
32 #include <TH1I.h> 
33 #include <TH2F.h> 
34 #include <TTree.h>
35
36 //-----------------------------------------------------------------------------
37 /// \class AliMUONQADataMakerSim
38 ///
39 /// MUON base class for quality assurance data (histo) maker
40 ///
41 /// \author C. Finck
42
43 /// \cond CLASSIMP
44 ClassImp(AliMUONQADataMakerSim)
45 /// \endcond
46            
47 //____________________________________________________________________________ 
48 AliMUONQADataMakerSim::AliMUONQADataMakerSim() : 
49     AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kMUON), "MUON Quality Assurance Data Maker"),
50     fHitStore(0x0),
51     fDigitStore(0x0)   
52 {
53   /// Default constructor
54
55   AliDebug(1,"");
56 }
57
58 //____________________________________________________________________________ 
59 AliMUONQADataMakerSim::AliMUONQADataMakerSim(const AliMUONQADataMakerSim& qadm) :
60     AliQADataMakerSim(),
61   fHitStore(0x0),
62   fDigitStore(0x0)
63 {
64   /// Copy constructor
65
66   AliDebug(1,"");
67
68     if ( qadm.fHitStore ) 
69     {
70       fHitStore = static_cast<AliMUONVHitStore*>(qadm.fHitStore->Clone());
71     }
72     if ( qadm.fDigitStore ) 
73     {
74       fDigitStore = static_cast<AliMUONVDigitStore*>(qadm.fDigitStore->Clone());
75     }
76     SetName((const char*)qadm.GetName()) ; 
77     SetTitle((const char*)qadm.GetTitle()); 
78 }
79
80 //__________________________________________________________________
81 AliMUONQADataMakerSim& AliMUONQADataMakerSim::operator = (const AliMUONQADataMakerSim& qadm )
82 {
83   /// Assignment operator
84
85     AliDebug(1,"");
86
87     this->~AliMUONQADataMakerSim();
88     new(this) AliMUONQADataMakerSim(qadm);
89     return *this;
90 }
91
92 //__________________________________________________________________
93 AliMUONQADataMakerSim::~AliMUONQADataMakerSim()
94 {
95   /// Destructor
96
97   AliDebug(1,"");
98
99   delete fHitStore;
100   delete fDigitStore;
101 }
102
103 //__________________________________________________________________
104 void AliMUONQADataMakerSim::InitHits() 
105 {
106   /// Initialized hit spectra
107   const Bool_t expert   = kTRUE ; 
108   const Bool_t image    = kTRUE ; 
109   
110   TH1F* h0 = new TH1F("hHitDetElem", "DetElemId distribution in Hits", 1400, 100., 1500.); 
111   Add2HitsList(h0, 0, !expert, image);
112
113   TH1F* h1 = new TH1F("hHitPtot", "P distribution in Hits ", 300, 0., 300.); 
114   Add2HitsList(h1, 1, !expert, image);
115   return;
116
117
118 //__________________________________________________________________
119 void AliMUONQADataMakerSim::InitSDigits() 
120 {
121   /// Initialized SDigits spectra
122   const Bool_t expert   = kTRUE ; 
123   const Bool_t image    = kTRUE ; 
124   
125   TH1I* h0 = new TH1I("hSDigitsDetElem", "Detection element distribution in SDigits",  1400, 100, 1500); 
126   Add2SDigitsList(h0, 0, !expert, image);
127
128   TH1F* h1 = new TH1F("hSDigitsCharge", "Charge distribution in SDigits", 4096, 0, 4095); 
129   Add2SDigitsList(h1, 1, !expert, image);
130
131 }  
132
133 //__________________________________________________________________
134 void AliMUONQADataMakerSim::InitDigits() 
135 {
136   /// Initialized Digits spectra 
137   const Bool_t expert   = kTRUE ; 
138   const Bool_t image    = kTRUE ; 
139   
140   TH1I* h0 = new TH1I("hDigitsDetElem", "Detection element distribution in Digits",  1400, 100, 1500); 
141   Add2DigitsList(h0, 0, !expert, image);
142
143   TH1I* h1 = new TH1I("hDigitsADC", "ADC distribution in Digits", 4096, 0, 4095); 
144   Add2DigitsList(h1, 1, !expert, image);  
145
146
147
148 //__________________________________________________________________
149 void AliMUONQADataMakerSim::MakeHits(TTree* hitsTree)        
150 {
151   /// makes data from Hits
152   if (!fHitStore)
153     fHitStore = AliMUONVHitStore::Create(*hitsTree);
154   fHitStore->Connect(*hitsTree, false);
155   hitsTree->GetEvent(0);
156     
157   TIter next(fHitStore->CreateIterator());
158
159   AliMUONHit* hit = 0x0;
160
161   while ( ( hit = static_cast<AliMUONHit*>(next()) ) )
162   {
163     GetHitsData(0)->Fill(hit->DetElemId());
164     GetHitsData(1)->Fill(hit->Momentum());
165   }
166
167   
168 }
169
170 //__________________________________________________________________
171 void AliMUONQADataMakerSim::MakeSDigits(TTree* sdigitsTree)        
172 {
173   /// makes data from SDigits
174   if (!fDigitStore)
175     fDigitStore = AliMUONVDigitStore::Create(*sdigitsTree);
176   fDigitStore->Connect(*sdigitsTree, false);
177   sdigitsTree->GetEvent(0);
178     
179   TIter next(fDigitStore->CreateIterator());
180
181   AliMUONVDigit* dig = 0x0;
182
183   while ( ( dig = static_cast<AliMUONVDigit*>(next()) ) )
184   {
185     GetSDigitsData(0)->Fill(dig->DetElemId());
186     GetSDigitsData(1)->Fill(dig->Charge());
187   }
188
189
190 //__________________________________________________________________
191 void AliMUONQADataMakerSim::MakeDigits(TTree* digitsTree)         
192 {
193    /// makes data from Digits
194   if (!fDigitStore)
195     fDigitStore = AliMUONVDigitStore::Create(*digitsTree);
196   fDigitStore->Connect(*digitsTree, false);
197   digitsTree->GetEvent(0);
198     
199   TIter next(fDigitStore->CreateIterator());
200
201   AliMUONVDigit* dig = 0x0;
202
203   while ( ( dig = static_cast<AliMUONVDigit*>(next()) ) )
204   {
205     GetDigitsData(0)->Fill(dig->DetElemId());
206     GetDigitsData(1)->Fill(dig->ADC());
207   }
208 }
209       
210 //____________________________________________________________________________ 
211 void AliMUONQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
212 {
213     ///Detector specific actions at end of cycle
214     // do the QA checking
215     AliQAChecker::Instance()->Run(AliQAv1::kMUON, task, list) ;  
216 }
217
218
219 //____________________________________________________________________________ 
220 void AliMUONQADataMakerSim::StartOfDetectorCycle()
221 {
222     /// Detector specific actions at start of cycle
223   
224 }