]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PHOS/AliPHOSRcuDA1.cxx
Starting a collection of QA/Comparison macros
[u/mrichter/AliRoot.git] / PHOS / AliPHOSRcuDA1.cxx
CommitLineData
354c9df9 1#include "AliPHOSRcuDA1.h"
2#include "TString.h"
3
4ClassImp(AliPHOSRcuDA1)
5
6//----------------------------------------------------------------
7AliPHOSRcuDA1::AliPHOSRcuDA1(Int_t module, Int_t rcu) : TNamed(),
8 fHistoFile(0),fMod(module),fRCU(rcu),fWriteToFile(kTRUE),fHistoArray()
9
10{
11 // Create DA1 ("Calibration DA") object.
12 // module is the PHOS module number (0..4).
13 // Checks existence of histograms which might have been left
14 // from the previous runs to continue their filling.
15 // Histogram names: module_iX_iZ_gain for TH2 and module_iX_iZ for TH1.
4056a2e7 16 // Root file name: PHOS_ModuleX_RCUY_Calib.root, where X - module number,
17 // Y - RCU number. If no RCU specified (rcu<0), file name is simply
18 // PHOS_ModuleX_Calib.root.
354c9df9 19
6ba1dd81 20 char name[128]; TString sname;
21
22 if(rcu<0) {
23 sname="PHOS_Module%d_Calib";
24 snprintf(name,sname.Length(),sname.Data(),fMod);
25 }
26 else {
27 sname="PHOS_Module%d_RCU%d_Calib";
28 snprintf(name,sname.Length(),sname.Data(),fMod,fRCU);
29 }
354c9df9 30
6ba1dd81 31 SetName(name);
354c9df9 32 SetTitle("Calibration Detector Algorithm");
33
6ba1dd81 34 char rootname[128];
35 TString srootname="%s.root";
36 snprintf(rootname,srootname.Length(),srootname.Data(),GetName());
354c9df9 37
38 fHistoFile = new TFile(rootname,"update");
39 fHistoArray.SetName("histo_container");
40
354c9df9 41 char hname[128];
42 TH1F* hist1=0;
43 TH2F* hist2=0;
6ba1dd81 44 TString shname;
354c9df9 45
46 for(Int_t iX=0; iX<64; iX++) {
47 for(Int_t iZ=0; iZ<56; iZ++) {
6ba1dd81 48
49 shname = "%d_%d_%d";
50 snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ);
354c9df9 51
52fccb4e 52 hist1 = (TH1F*)fHistoFile->Get(hname);
354c9df9 53 if(hist1) {
54 fHgLgRatio[iX][iZ] = hist1;
55 fHistoArray.Add(hist1);
56 }
57 else
58 fHgLgRatio[iX][iZ] = 0;
59
60 for(Int_t iGain=0; iGain<2; iGain++) {
6ba1dd81 61
62 shname = "%d_%d_%d_%d";
63 snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ,iGain);
64
52fccb4e 65 hist2 = (TH2F*)fHistoFile->Get(hname);
354c9df9 66 if(hist2) {
67 fTimeEnergy[iX][iZ][iGain] = hist2;
68 fHistoArray.Add(hist2);
69 }
70 else
71 fTimeEnergy[iX][iZ][iGain] = 0;
72 }
73
74 }
75 }
76
77}
78
79//-------------------------------------------------------------------
80AliPHOSRcuDA1::AliPHOSRcuDA1(Int_t module, Int_t rcu, TObjArray* oldTimeEnergy) :
81 TNamed(),fHistoFile(0),fMod(module),fRCU(rcu),fWriteToFile(kFALSE),
82 fHistoArray()
83{
84 // Constructor.
85 // oldTimeEnergy is an array of histograms kept from the previous run.
86 // By default the final histograms will not be saved to the root file.
87
6ba1dd81 88 char name[128]; TString sname;
89
90 if(rcu<0) {
91 sname="PHOS_Module%d_Calib";
92 snprintf(name,sname.Length(),sname.Data(),fMod);
93 }
94 else {
95 sname="PHOS_Module%d_RCU%d_Calib";
96 snprintf(name,sname.Length(),sname.Data(),fMod,fRCU);
97 }
354c9df9 98
6ba1dd81 99 SetName(name);
354c9df9 100 SetTitle("Calibration Detector Algorithm");
6ba1dd81 101
354c9df9 102 fHistoArray.SetName("histo_container");
103
104 char hname[128];
105 TH1F* hist1=0;
106 TH2F* hist2=0;
6ba1dd81 107 TString shname;
354c9df9 108
109 for(Int_t iX=0; iX<64; iX++) {
110 for(Int_t iZ=0; iZ<56; iZ++) {
111
6ba1dd81 112 shname = "%d_%d_%d";
113 snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ);
114
354c9df9 115 if(oldTimeEnergy)
116 hist1 = (TH1F*)oldTimeEnergy->FindObject(hname);
117 if(hist1) {
118 fHgLgRatio[iX][iZ] = hist1;
119 fHistoArray.Add(hist1);
120 }
121 else
122 fHgLgRatio[iX][iZ] = 0;
123
124 for(Int_t iGain=0; iGain<2; iGain++) {
6ba1dd81 125
126 shname = "%d_%d_%d_%d";
127 snprintf(hname,shname.Length(),shname.Data(),fMod,iX,iZ,iGain);
128
354c9df9 129 if(oldTimeEnergy)
130 hist2 = (TH2F*)oldTimeEnergy->FindObject(hname);
131 if(hist2) {
132 fTimeEnergy[iX][iZ][iGain] = hist2;
133 fHistoArray.Add(hist2);
134 }
135 else
136 fTimeEnergy[iX][iZ][iGain] = 0;
137 }
138
139 }
140 }
141
142}
143
144//-------------------------------------------------------------------
145AliPHOSRcuDA1::~AliPHOSRcuDA1()
146{
147 // Destructor
148
149 UpdateHistoFile();
150 if(fHistoFile) delete fHistoFile;
151
152}
153
154//-------------------------------------------------------------------
155void AliPHOSRcuDA1::FillHistograms(Float_t e[64][56][2], Float_t t[64][56][2])
156{
157 // Fill energy vs time-of-flight and HG/LG ratio histograms in one event.
158 // Energy and TOF data are encoded as e[X][Z][gain] and t[X][Z][gain],
159 // where X(0..63) and Z(0..55) - crystal position in the module,
160 // gain=0 - low gain, gain=1 - high gain.
161 // Energy in ADC counts, time in samples.
162 // If no energy or time read for particular channel,
163 // the correspondent array entry should be filled by zero.
164 // WARNING: this function should be called once per event!
165
7e6cd815 166 TString hname;
167 TString htitle;
354c9df9 168
169 for(Int_t iX=0; iX<64; iX++) {
170 for (Int_t iZ=0; iZ<56; iZ++) {
171
172 // HG/LG
2d03aabf 173 if(e[iX][iZ][0]>10. && e[iX][iZ][1]>10. && e[iX][iZ][1]<900.) {
354c9df9 174
175 if(fHgLgRatio[iX][iZ]) {
176 //printf("iX=%d iZ=%d,e[iX][iZ][1]=%.3f,e[iX][iZ][0]=%.3f, t1=%.3f,t0=%.3f\n",
177 // iX,iZ,e[iX][iZ][1],e[iX][iZ][0], t[iX][iZ][1],t[iX][iZ][0]);
178 fHgLgRatio[iX][iZ]->Fill(e[iX][iZ][1]/e[iX][iZ][0]);
179 }
180 else
181 {
7e6cd815 182 hname.Clear(); htitle.Clear();
183 hname += fMod; hname += "_"; hname += iX; hname += "_"; hname += iZ;
184 htitle += "HG/LG ratio for crystal "; htitle += hname;
185
186 fHgLgRatio[iX][iZ] = new TH1F(hname,htitle,400,14.,18.);
354c9df9 187// printf("iX=%d iZ=%d,e[iX][iZ][1]=%.3f,e[iX][iZ][0]=%.3f\n",
188// iX,iZ,e[iX][iZ][1],e[iX][iZ][0]);
189 fHgLgRatio[iX][iZ]->Fill(e[iX][iZ][1]/e[iX][iZ][0]);
190 fHistoArray.Add(fHgLgRatio[iX][iZ]);
191 }
192 }
193
194 // Energy vs TOF
195 for(Int_t iGain=0; iGain<2; iGain++) {
196
197 if(e[iX][iZ][iGain]<10) continue;
198 if(!(t[iX][iZ][iGain]>0)) continue;
199
200 if(fTimeEnergy[iX][iZ][iGain])
201 fTimeEnergy[iX][iZ][iGain]->Fill(e[iX][iZ][iGain],t[iX][iZ][iGain]);
7e6cd815 202 else
203 {
204 hname.Clear(); htitle.Clear();
205 hname += fMod; hname += "_"; hname += iX; hname += "_"; hname += iZ;
206 htitle = "Energy vs TOF for crystal "; htitle += hname; htitle += " and gain "; htitle += iGain;
207 hname += "_"; hname += iGain;
208
209 fTimeEnergy[iX][iZ][iGain] = new TH2F(hname,htitle,100,0.,1024.,50,0.,50.);
210 fTimeEnergy[iX][iZ][iGain]->Fill(e[iX][iZ][iGain],t[iX][iZ][iGain]);
211 fHistoArray.Add(fTimeEnergy[iX][iZ][iGain]);
212 }
354c9df9 213 }
7e6cd815 214
354c9df9 215 }
216 }
7e6cd815 217
354c9df9 218}
219
220//-------------------------------------------------------------------
221void AliPHOSRcuDA1::UpdateHistoFile()
222{
223 // Write histograms to file
224
225 if(!fWriteToFile) return;
226 if(!fHistoFile) return;
227 if(!fHistoFile->IsOpen()) return;
228
f39fcb60 229 fHistoArray.Write();
354c9df9 230 fHistoFile->Purge();
231
232}
233
234//-----------------------------------------------------------------------
235void AliPHOSRcuDA1::SetWriteToFile(Bool_t write)
236{
237 if(!write) {
238 fWriteToFile = write;
239 return;
240 }
241
242 if(!fHistoFile) {
7e6cd815 243 TString rootname(GetName());
244 rootname += ".root";
245 fHistoFile = new TFile(rootname.Data(),"update");
354c9df9 246 }
247
248 fWriteToFile = write;
249}