]>
Commit | Line | Data |
---|---|---|
5253c153 | 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 "AliMUONPedestal.h" | |
19 | #include "AliMUONErrorCounter.h" | |
20 | #include "AliMUONVStore.h" | |
21 | #include "AliMUON2DMap.h" | |
22 | #include "AliMUONCalibParamND.h" | |
23 | ||
24 | #include <TString.h> | |
5253c153 | 25 | #include <TTimeStamp.h> |
26 | #include <TMath.h> | |
27 | #include <TTree.h> | |
28 | #include <TFile.h> | |
29 | #include <TH1F.h> | |
30 | #include <Riostream.h> | |
31 | ||
32 | #include <sstream> | |
33 | ||
34 | //----------------------------------------------------------------------------- | |
35 | /// \class AliMUONPedestal | |
36 | /// | |
37 | /// Implementation of the pedestal computing | |
38 | /// | |
39 | /// add | |
40 | /// | |
41 | /// | |
42 | /// \author Alberto Baldisseri, JL Charvet (05/05/2009) | |
43 | //----------------------------------------------------------------------------- | |
44 | ||
45 | /// \cond CLASSIMP | |
46 | ClassImp(AliMUONPedestal) | |
47 | /// \endcond | |
48 | ||
49 | //______________________________________________________________________________ | |
50 | AliMUONPedestal::AliMUONPedestal() | |
51 | : TObject(), | |
5431405e | 52 | fN(0), |
5253c153 | 53 | fNEvents(0), |
54 | fRunNumber(0), | |
55 | fNChannel(0), | |
56 | fNManu(0), | |
a54854fc | 57 | fErrorBuspatchTable(new AliMUON2DMap(kFALSE)), |
58 | fManuBuspatchTable(new AliMUON2DMap(kFALSE)), | |
5253c153 | 59 | fDate(new TTimeStamp()), |
60 | fFilcout(0), | |
a54854fc | 61 | fPedestalStore(new AliMUON2DMap(kTRUE)), |
5253c153 | 62 | fIndex(-1) |
63 | { | |
64 | /// Default constructor | |
521ca393 | 65 | sprintf(fHistoFileName," "); |
66 | sprintf(fprefixDA," "); | |
5253c153 | 67 | } |
68 | // AliMUONPedestal& operator=(const AliMUONPedestal& other); Copy ctor | |
69 | ||
70 | //______________________________________________________________________________ | |
71 | AliMUONPedestal::~AliMUONPedestal() | |
72 | { | |
73 | /// Destructor | |
a54854fc | 74 | delete fErrorBuspatchTable; |
75 | delete fManuBuspatchTable; | |
dbbb2c64 | 76 | delete fPedestalStore; |
5253c153 | 77 | } |
78 | ||
79 | //______________________________________________________________________________ | |
80 | void AliMUONPedestal::MakePed(Int_t busPatchId, Int_t manuId, Int_t channelId, Int_t charge) | |
81 | { | |
82 | /// Compute pedestals values | |
83 | AliMUONVCalibParam* ped = | |
84 | static_cast<AliMUONVCalibParam*>(fPedestalStore ->FindObject(busPatchId, manuId)); | |
85 | ||
86 | if (!ped) { | |
87 | fNManu++; | |
88 | ped = new AliMUONCalibParamND(2, kNChannels,busPatchId, manuId, -1.); // put default wise -1, not connected channel | |
89 | fPedestalStore ->Add(ped); | |
90 | } | |
91 | ||
92 | // Initialization for the first value | |
93 | if (ped->ValueAsDouble(channelId, 0) == -1) ped->SetValueAsDouble(channelId, 0, 0.); | |
94 | if (ped->ValueAsDouble(channelId, 1) == -1) ped->SetValueAsDouble(channelId, 1, 0.); | |
95 | ||
96 | Double_t pedMean = ped->ValueAsDouble(channelId, 0) + (Double_t) charge; | |
97 | Double_t pedSigma = ped->ValueAsDouble(channelId, 1) + (Double_t) charge*charge; | |
98 | ||
99 | ped->SetValueAsDouble(channelId, 0, pedMean); | |
100 | ped->SetValueAsDouble(channelId, 1, pedSigma); | |
5431405e | 101 | |
5431405e | 102 | AliMUONErrorCounter* manuCounter; |
a54854fc | 103 | if (!(manuCounter = static_cast<AliMUONErrorCounter*>(fManuBuspatchTable->FindObject(busPatchId,manuId)))) |
5431405e | 104 | { |
105 | // New (buspatch,manu) | |
106 | manuCounter = new AliMUONErrorCounter(busPatchId,manuId); | |
a54854fc | 107 | fManuBuspatchTable->Add(manuCounter); |
5431405e | 108 | } |
109 | else | |
110 | { | |
111 | // Existing buspatch | |
112 | manuCounter->Increment(); | |
113 | } | |
5253c153 | 114 | } |
115 | ||
116 | //______________________________________________________________________________ | |
dbbb2c64 | 117 | void AliMUONPedestal::Finalize() |
5253c153 | 118 | { |
5253c153 | 119 | Double_t pedMean; |
120 | Double_t pedSigma; | |
5253c153 | 121 | Int_t busPatchId; |
122 | Int_t manuId; | |
123 | Int_t channelId; | |
124 | ||
5253c153 | 125 | // print in logfile |
126 | if (fErrorBuspatchTable->GetSize()) | |
127 | { | |
128 | cout<<"\n* Buspatches with less statistics (due to parity errors)"<<endl; | |
129 | (*fFilcout)<<"\n* Buspatches with less statistics (due to parity errors)"<<endl; | |
a54854fc | 130 | TIter nextParityError(fErrorBuspatchTable->CreateIterator()); |
5253c153 | 131 | AliMUONErrorCounter* parityerror; |
a54854fc | 132 | while((parityerror = static_cast<AliMUONErrorCounter*>(nextParityError()))) |
5253c153 | 133 | { |
134 | cout<<" bp "<<parityerror->BusPatch()<<": events used = "<<fNEvents-parityerror->Events()<<endl; | |
135 | (*fFilcout)<<" bp "<<parityerror->BusPatch()<<": events used = "<<fNEvents-parityerror->Events()<<endl; | |
136 | } | |
5253c153 | 137 | } |
138 | ||
5253c153 | 139 | // iterator over pedestal |
140 | TIter next(fPedestalStore ->CreateIterator()); | |
141 | AliMUONVCalibParam* ped; | |
142 | ||
143 | while ( ( ped = dynamic_cast<AliMUONVCalibParam*>(next() ) ) ) | |
144 | { | |
145 | busPatchId = ped->ID0(); | |
146 | manuId = ped->ID1(); | |
5431405e | 147 | if(manuId==0) |
148 | { | |
149 | cout << " !!! BIG WARNING: ManuId = " << manuId << " !!! in BP = " << busPatchId << endl; | |
150 | (*fFilcout) << " !!! BIG WARNING: ManuId = " << manuId << " !!! in BP = " << busPatchId << endl; | |
151 | } | |
5253c153 | 152 | Int_t eventCounter; |
5253c153 | 153 | // Correct the number of events for buspatch with errors |
5253c153 | 154 | AliMUONErrorCounter* errorCounter; |
a54854fc | 155 | if ((errorCounter = (AliMUONErrorCounter*)fErrorBuspatchTable->FindObject(busPatchId))) |
5253c153 | 156 | { |
157 | eventCounter = fNEvents - errorCounter->Events(); | |
158 | } | |
159 | else | |
160 | { | |
161 | eventCounter = fNEvents; | |
162 | } | |
163 | ||
5431405e | 164 | Int_t occupancy; |
165 | // value of (buspatch, manu) occupancy | |
5431405e | 166 | AliMUONErrorCounter* manuCounter; |
a54854fc | 167 | manuCounter = static_cast<AliMUONErrorCounter*>(fManuBuspatchTable->FindObject(busPatchId,manuId)); |
5431405e | 168 | occupancy = manuCounter->Events()/64/eventCounter; |
169 | if(occupancy>1) | |
dbbb2c64 | 170 | { |
a54854fc | 171 | cout << " !!! BIG WARNING: ManuId = " << manuId << " !!! in BP = " << busPatchId << " occupancy (>1) = " << occupancy << endl; |
172 | (*fFilcout) << " !!! BIG WARNING: ManuId = " << manuId << " !!! in BP = " << busPatchId << " occupancy (>1) = " << occupancy <<endl; | |
dbbb2c64 | 173 | } |
5431405e | 174 | |
dbbb2c64 | 175 | for (channelId = 0; channelId < ped->Size() ; ++channelId) |
176 | { | |
5253c153 | 177 | pedMean = ped->ValueAsDouble(channelId, 0); |
178 | ||
5431405e | 179 | if (pedMean > 0) // connected channels |
dbbb2c64 | 180 | { |
5431405e | 181 | ped->SetValueAsDouble(channelId, 0, pedMean/(Double_t)eventCounter); |
182 | pedMean = ped->ValueAsDouble(channelId, 0); | |
183 | pedSigma = ped->ValueAsDouble(channelId, 1); | |
184 | ped->SetValueAsDouble(channelId, 1, TMath::Sqrt(TMath::Abs(pedSigma/(Double_t)eventCounter - pedMean*pedMean))); | |
185 | if(manuId == 0 || occupancy>1) | |
dbbb2c64 | 186 | { |
5431405e | 187 | ped->SetValueAsDouble(channelId, 0, kADCMax); |
188 | ped->SetValueAsDouble(channelId, 1, kADCMax); | |
dbbb2c64 | 189 | } |
190 | } | |
5431405e | 191 | else |
dbbb2c64 | 192 | { |
5431405e | 193 | ped->SetValueAsDouble(channelId, 0, kADCMax); |
194 | ped->SetValueAsDouble(channelId, 1, kADCMax); | |
dbbb2c64 | 195 | } |
196 | } | |
197 | } | |
198 | } | |
5253c153 | 199 | |
dbbb2c64 | 200 | //______________________________________________________________________________ |
201 | void AliMUONPedestal::MakeASCIIoutput(ostream& out) const | |
202 | { | |
203 | /// put pedestal store in the output stream | |
204 | ||
205 | out<<"//===========================================================================" << endl; | |
206 | out<<"// Pedestal file calculated by MUONTRKda"<<endl; | |
207 | out<<"//===========================================================================" << endl; | |
208 | out<<"// * Run : " << fRunNumber << endl; | |
209 | out<<"// * Date : " << fDate->AsString("l") <<endl; | |
210 | out<<"// * Statictics : " << fNEvents << endl; | |
211 | out<<"// * # of MANUS : " << fNManu << endl; | |
212 | out<<"// * # of channels : " << fNChannel << endl; | |
213 | if (fErrorBuspatchTable->GetSize()) | |
214 | { | |
215 | out<<"//"<<endl; | |
216 | out<<"// * Buspatches with less statistics (due to parity errors)"<<endl; | |
a54854fc | 217 | TIter next(fErrorBuspatchTable->CreateIterator()); |
dbbb2c64 | 218 | AliMUONErrorCounter* parityerror; |
a54854fc | 219 | while((parityerror = static_cast<AliMUONErrorCounter*>(next()))) |
dbbb2c64 | 220 | { |
221 | out<<"// bp "<<parityerror->BusPatch()<<" events used "<<fNEvents-parityerror->Events()<<endl; | |
222 | } | |
223 | } | |
224 | out<<"//"<<endl; | |
225 | out<<"//---------------------------------------------------------------------------" << endl; | |
226 | out<<"//---------------------------------------------------------------------------" << endl; | |
227 | out<<"// BP MANU CH. MEAN SIGMA"<<endl; | |
228 | out<<"//---------------------------------------------------------------------------" << endl; | |
5253c153 | 229 | |
dbbb2c64 | 230 | // iterator over pedestal |
231 | TIter next(fPedestalStore ->CreateIterator()); | |
232 | AliMUONVCalibParam* ped; | |
233 | ||
234 | while ( ( ped = dynamic_cast<AliMUONVCalibParam*>(next() ) ) ) | |
235 | { | |
236 | Int_t busPatchId = ped->ID0(); | |
237 | Int_t manuId = ped->ID1(); | |
238 | ||
239 | for ( Int_t channelId = 0; channelId < ped->Size(); ++channelId ) | |
240 | { | |
241 | Double_t pedMean = ped->ValueAsDouble(channelId, 0); | |
242 | Double_t pedSigma = ped->ValueAsDouble(channelId, 1); | |
243 | ||
244 | out << "\t" << busPatchId << "\t" << manuId <<"\t"<< channelId << "\t" << pedMean <<"\t"<< pedSigma << endl; | |
5253c153 | 245 | } |
246 | } | |
dbbb2c64 | 247 | } |
248 | ||
249 | //______________________________________________________________________________ | |
250 | void AliMUONPedestal::MakeControlHistos() | |
251 | { | |
5253c153 | 252 | |
dbbb2c64 | 253 | if (fIndex>=0) return; // Pedestal run (fIndex=-1) |
254 | ||
255 | Double_t pedMean; | |
256 | Double_t pedSigma; | |
257 | Int_t busPatchId; | |
258 | Int_t manuId; | |
259 | Int_t channelId; | |
260 | ||
261 | // histo | |
262 | TFile* histoFile = 0; | |
263 | TTree* tree = 0; | |
264 | TH1F* pedMeanHisto = 0; | |
265 | TH1F* pedSigmaHisto = 0; | |
266 | ||
267 | sprintf(fHistoFileName,"%s.root",fprefixDA); | |
268 | histoFile = new TFile(fHistoFileName,"RECREATE","MUON Tracking pedestals"); | |
269 | ||
270 | Char_t name[255]; | |
271 | Char_t title[255]; | |
272 | sprintf(name,"pedmean_allch"); | |
273 | sprintf(title,"Pedestal mean all channels"); | |
274 | Int_t nx = kADCMax+1; | |
275 | Int_t xmin = 0; | |
276 | Int_t xmax = kADCMax; | |
277 | pedMeanHisto = new TH1F(name,title,nx,xmin,xmax); | |
278 | pedMeanHisto->SetDirectory(histoFile); | |
279 | ||
280 | sprintf(name,"pedsigma_allch"); | |
281 | sprintf(title,"Pedestal sigma all channels"); | |
282 | nx = 201; | |
283 | xmin = 0; | |
284 | xmax = 200; | |
285 | pedSigmaHisto = new TH1F(name,title,nx,xmin,xmax); | |
286 | pedSigmaHisto->SetDirectory(histoFile); | |
287 | ||
288 | tree = new TTree("t","Pedestal tree"); | |
289 | tree->Branch("bp",&busPatchId,"bp/I"); | |
290 | tree->Branch("manu",&manuId,",manu/I"); | |
291 | tree->Branch("channel",&channelId,",channel/I"); | |
292 | tree->Branch("pedMean",&pedMean,",pedMean/D"); | |
293 | tree->Branch("pedSigma",&pedSigma,",pedSigma/D"); | |
294 | ||
295 | // iterator over pedestal | |
296 | TIter next(fPedestalStore ->CreateIterator()); | |
297 | AliMUONVCalibParam* ped; | |
298 | ||
299 | while ( ( ped = dynamic_cast<AliMUONVCalibParam*>(next() ) ) ) | |
300 | { | |
301 | busPatchId = ped->ID0(); | |
302 | manuId = ped->ID1(); | |
303 | ||
304 | for ( channelId = 0; channelId < ped->Size(); ++channelId ) | |
305 | { | |
306 | pedMean = ped->ValueAsDouble(channelId, 0); | |
307 | pedSigma = ped->ValueAsDouble(channelId, 1); | |
308 | ||
309 | pedMeanHisto->Fill(pedMean); | |
310 | pedSigmaHisto->Fill(pedSigma); | |
311 | tree->Fill(); | |
312 | } | |
5253c153 | 313 | } |
dbbb2c64 | 314 | |
315 | histoFile->Write(); | |
316 | histoFile->Close(); | |
317 | ||
4af2c34a | 318 | } |