Removed hardcoded restriction on the minimal digit energy.
[u/mrichter/AliRoot.git] / PHOS / AliPHOSCalibHistoProducer.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 ///////////////////////////////////////////////////////////////////////////////
19 // Class AliPHOSCalibHistoProducer accumulating histograms
20 // with amplitudes per PHOS channel
21 // It is intended to run at DAQ computers (LDC, GDC, HLT or MOOD)
22 // and it fills the histograms with amplitudes per channel.
23 // Usage example see in PHOS/macros/Shuttle/AliPHOSCalibHistoProducer.C
24 //
25 // Author: Boris Polichtchouk, 4 October 2006
26 ///////////////////////////////////////////////////////////////////////////////
27
28 #include "AliLog.h"
29 #include "AliPHOSCalibHistoProducer.h"
30 #include "TH1.h"
31 #include "TH2F.h"
32 #include "TFile.h"
33 #include "AliPHOSRawDecoder.h"
34
35 ClassImp(AliPHOSCalibHistoProducer)
36
37 //-----------------------------------------------------------------------------
38 AliPHOSCalibHistoProducer::AliPHOSCalibHistoProducer() : 
39   fRawDecoder(0),fHistoFile(0),fUpdatingRate(100),fIsOldRCUFormat(kFALSE),
40   fEvents(0),fNbins(100),fXlow(0.),fXup(1000.)
41 {
42   // Constructor: initializes data members
43   // Checks existence of histograms which might have been left
44   // from the previous runs to continues their filling
45
46   fHistoFile =  new TFile("calibHisto.root","update");
47   
48   for(Int_t module=0; module<5; module++) {
49     for(Int_t column=0; column<56; column++) {
50       for(Int_t row=0; row<64; row++) {
51         char hname[128];
52         sprintf(hname,"mod%dcol%drow%d",module,column,row);
53         TH1F* hist = (TH1F*)fHistoFile->Get(hname);
54         if(hist) 
55           fAmpHisto[module][column][row]=hist;
56         else 
57           fAmpHisto[module][column][row] = 0;     
58       }
59     }
60   }
61 }
62
63 //-----------------------------------------------------------------------------            
64 AliPHOSCalibHistoProducer::AliPHOSCalibHistoProducer(Int_t nbinsx, Double_t xlow, Double_t xup) :
65   fRawDecoder(0),fHistoFile(0),fUpdatingRate(100),fIsOldRCUFormat(kFALSE),
66   fEvents(0),fNbins(nbinsx),fXlow(xlow),fXup(xup)
67 {
68   // Constructor: initializes data members.
69   // Checks existence of histograms which might have been left
70   // from the previous runs to continues their filling.
71   // In addition sets number of bins, low and upper limits common for all histograms.
72
73   fHistoFile =  new TFile("calibHisto.root","update");
74
75   for(Int_t module=0; module<5; module++) {
76     for(Int_t column=0; column<56; column++) {
77       for(Int_t row=0; row<64; row++) {
78         char hname[128];
79         sprintf(hname,"mod%dcol%drow%d",module,column,row);
80         TH1F* hist = (TH1F*)fHistoFile->Get(hname);
81         if(hist)
82           fAmpHisto[module][column][row]=hist;
83         else
84           fAmpHisto[module][column][row] = 0;
85       }
86     }
87   }
88 }
89
90 //-----------------------------------------------------------------------------
91 AliPHOSCalibHistoProducer::~AliPHOSCalibHistoProducer()
92 {
93   // Destructor
94   
95   UpdateHistoFile();
96   if(fHistoFile) delete fHistoFile;
97
98 }
99
100 //-----------------------------------------------------------------------------
101 AliPHOSCalibHistoProducer::AliPHOSCalibHistoProducer(const AliPHOSCalibHistoProducer &histoproducer) :
102   TObject(histoproducer),fRawDecoder(histoproducer.fRawDecoder),fHistoFile(histoproducer.fHistoFile),
103   fUpdatingRate(histoproducer.fUpdatingRate),fIsOldRCUFormat(histoproducer.fIsOldRCUFormat),
104   fEvents(histoproducer.fEvents),fNbins(histoproducer.fNbins),fXlow(histoproducer.fXlow),fXup(histoproducer.fXup)
105 {
106   //Copy constructor.
107
108   for(Int_t module=0; module<5; module++) {
109     for(Int_t column=0; column<56; column++) {
110       for(Int_t row=0; row<64; row++) {
111         char hname[128];
112         sprintf(hname,"mod%dcol%drow%d",module,column,row);
113         TH1F* hist = (TH1F*)histoproducer.fHistoFile->Get(hname);
114         if(hist) 
115           fAmpHisto[module][column][row]= new TH1F(*hist);
116         else
117           fAmpHisto[module][column][row]=0;
118       }
119     }
120   }
121 }
122
123 //-----------------------------------------------------------------------------
124 AliPHOSCalibHistoProducer& AliPHOSCalibHistoProducer::operator= 
125 (const AliPHOSCalibHistoProducer &histoproducer)
126 {
127   //Assignment operator.
128
129   if(this != &histoproducer) {
130
131     fRawDecoder = histoproducer.fRawDecoder;
132     fHistoFile = histoproducer.fHistoFile;
133     fUpdatingRate = histoproducer.fUpdatingRate;
134     fIsOldRCUFormat = histoproducer.fIsOldRCUFormat;
135     fEvents = histoproducer.fEvents;
136     fEvents = histoproducer.fEvents;
137     fNbins = histoproducer.fNbins;
138     fXlow = histoproducer.fXlow;
139     fXup = histoproducer.fXup;   
140     
141     for(Int_t module=0; module<5; module++) {
142       for(Int_t column=0; column<56; column++) {
143         for(Int_t row=0; row<64; row++) {
144           if(fAmpHisto[module][column][row]){
145             delete fAmpHisto[module][column][row];
146             fAmpHisto[module][column][row] = histoproducer.fAmpHisto[module][column][row];
147           }
148           else
149           fAmpHisto[module][column][row] = histoproducer.fAmpHisto[module][column][row];
150         }
151       }
152     }
153
154
155   }
156
157   return *this;
158 }
159 //-----------------------------------------------------------------------------
160 void AliPHOSCalibHistoProducer::Run()
161 {
162   // Reads raw data of current event and fills amplitude histograms
163   // The histograms are written to file every fUpdatingRate events
164
165   if(!fRawDecoder) AliFatal("Raw decoder not set!");
166   
167   Double_t energy;
168   Int_t mod,col,row;
169   
170   if(fIsOldRCUFormat)
171     fRawDecoder->SetOldRCUFormat(kTRUE);
172
173   while(fRawDecoder->NextDigit()) {
174     
175     if(fRawDecoder->IsLowGain()) continue; 
176
177     energy = fRawDecoder->GetEnergy();
178     
179     mod = fRawDecoder->GetModule()-1;
180     col = fRawDecoder->GetColumn()-1;
181     row = fRawDecoder->GetRow()-1;
182     
183     if(fAmpHisto[mod][col][row]) {
184       fAmpHisto[mod][col][row]->Fill(energy);
185     }
186     else {
187       char hname[128];
188       sprintf(hname,"mod%dcol%drow%d",mod,col,row);
189       fAmpHisto[mod][col][row] = new TH1F(hname,hname,fNbins,fXlow,fXup);
190       fAmpHisto[mod][col][row]->Fill(energy);
191     }
192     
193     // update histograms in local file every 100th event
194     if(fEvents%fUpdatingRate == 0) {
195       AliInfo(Form("Updating histo file, event %d, run %d\n",
196                    fEvents,fRawDecoder->GetRawReader()->GetRunNumber()));
197       UpdateHistoFile();
198     }
199     
200     //   UpdateHistoFile();
201     //   AliInfo(Form("%d events of run %d processed.",iEvent,runNum));
202   }
203   
204   fEvents++;
205   
206 }
207
208 //-----------------------------------------------------------------------------
209 void AliPHOSCalibHistoProducer::UpdateHistoFile()
210 {
211   // Write histograms to file
212
213   if(!fHistoFile) return;
214   if(!fHistoFile->IsOpen()) return;
215
216   TH1F* hist=0;
217   char hname[128];
218   char htitle[128];
219
220   for(Int_t module=0; module<5; module++) {
221     sprintf(hname,"hMeanE%d",module);
222     sprintf(htitle,"Mean energies in module %d",module);
223     TH2F hMeanE(hname,htitle,56,0.,56.,64,0.,64);
224
225     for(Int_t column=0; column<56; column++) {
226       for(Int_t row=0; row<64; row++) {
227         hist = fAmpHisto[module][column][row]; 
228         if(hist) hist->Write(hist->GetName(),TObject::kWriteDelete);
229         if(hist) hMeanE.SetBinContent(column,row,hist->GetMean());
230       }
231     }
232     hMeanE.Write(hMeanE.GetName(),TObject::kWriteDelete);
233   }
234
235 }