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