including task to produce final histograms on isolated photon analysis
[u/mrichter/AliRoot.git] / MUON / AliMUONCDB.cxx
CommitLineData
de01cdf0 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$ */
ee523da4 17
3d1463c8 18//-----------------------------------------------------------------------------
a99c3449 19/// \namespace AliMUONCDB
de01cdf0 20///
a99c3449 21/// Helper functions to experience the OCDB
22///
23/// They allow to read magnetic field, mapping and recoParam from OCDB
24///
25/// And also to generate dummy (but complete) containers for all the
de01cdf0 26/// calibration data types we have for tracker and trigger, and to write
27/// them into OCDB.
28///
a99c3449 29/// For more information, please see READMEcalib
de01cdf0 30///
baf1a7a0 31/// \author Laurent Aphecetche
3d1463c8 32//-----------------------------------------------------------------------------
de01cdf0 33
34#include "AliMUONCDB.h"
35
de01cdf0 36#include "AliMUON1DArray.h"
37#include "AliMUON1DMap.h"
38#include "AliMUON2DMap.h"
39#include "AliMUON2DStoreValidator.h"
7eafe398 40#include "AliMUONCalibParamND.h"
de01cdf0 41#include "AliMUONCalibParamNF.h"
42#include "AliMUONCalibParamNI.h"
207947f7 43#include "AliMUONCalibrationData.h"
de01cdf0 44#include "AliMUONConstants.h"
207947f7 45#include "AliMUONGlobalCrateConfig.h"
b37b9546 46#include "AliMUONLogger.h"
207947f7 47#include "AliMUONPadStatusMaker.h"
48#include "AliMUONPadStatusMapMaker.h"
49#include "AliMUONRecoParam.h"
50#include "AliMUONRegionalTriggerConfig.h"
0045b488 51#include "AliMUONRejectList.h"
9e7c4517 52#include "AliMUONTrackerData.h"
49419555 53#include "AliMUONTrackerIO.h"
de01cdf0 54#include "AliMUONTriggerEfficiencyCells.h"
55#include "AliMUONTriggerLut.h"
de01cdf0 56#include "AliMUONVCalibParam.h"
57#include "AliMUONVCalibParam.h"
207947f7 58#include "AliMUONVStore.h"
59
25e1df3e 60#include "AliMpCDB.h"
a0eca509 61#include "AliMpConstants.h"
a99c3449 62#include "AliMpDEStore.h"
de01cdf0 63#include "AliMpDDLStore.h"
ab167304 64#include "AliMpManuStore.h"
de01cdf0 65#include "AliMpDEManager.h"
66#include "AliMpDetElement.h"
92c23b09 67#include "AliMpFiles.h"
49e110ec 68#include "AliMpDCSNamer.h"
409de10f 69#include "AliMpManuIterator.h"
de01cdf0 70#include "AliMpSegmentation.h"
71#include "AliMpStationType.h"
72#include "AliMpVSegmentation.h"
92e36663 73
74#include "AliCodeTimer.h"
75#include "AliCDBEntry.h"
76#include "AliCDBManager.h"
a99c3449 77#include "AliGRPManager.h"
92e36663 78#include "AliDCSValue.h"
79#include "AliLog.h"
a31df857 80#include "AliMpBusPatch.h"
92e36663 81
de01cdf0 82#include <Riostream.h>
83#include <TArrayI.h>
84#include <TClass.h>
9e7c4517 85#include <TFile.h>
de01cdf0 86#include <TH1F.h>
87#include <TList.h>
88#include <TMap.h>
89#include <TObjString.h>
90#include <TROOT.h>
91#include <TRandom.h>
92#include <TStopwatch.h>
93#include <TSystem.h>
cf27231a 94#include <TMath.h>
a99c3449 95#include <TGeoGlobalMagField.h>
73d3aca5 96#include <TClonesArray.h>
a31df857 97#include <sstream>
9e7c4517 98#include <set>
c05673c9 99
de01cdf0 100namespace
101{
409de10f 102 //_____________________________________________________________________________
103AliMUONVStore* Create2DMap()
104{
105 return new AliMUON2DMap(true);
106}
107
108 //_____________________________________________________________________________
046e5fd4 109void getBoundaries(const AliMUONVStore& store, Int_t dim,
110 Float_t* xmin, Float_t* xmax)
de01cdf0 111{
046e5fd4 112 /// Assuming the store contains AliMUONVCalibParam objects, compute the
113 /// limits of the value contained in the VCalibParam, for each of its dimensions
114 /// xmin and xmax must be of dimension dim
115
116 for ( Int_t i = 0; i < dim; ++i )
117 {
118 xmin[i]=1E30;
119 xmax[i]=-1E30;
120 }
de01cdf0 121
a0eca509 122 TIter next(store.CreateIterator());
123 AliMUONVCalibParam* value;
de01cdf0 124
a0eca509 125 while ( ( value = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
de01cdf0 126 {
a0eca509 127 Int_t detElemId = value->ID0();
128 Int_t manuId = value->ID1();
de01cdf0 129
130 const AliMpVSegmentation* seg =
131 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
132
b37b9546 133 if (!seg) continue;
134
de01cdf0 135 for ( Int_t manuChannel = 0; manuChannel < value->Size(); ++manuChannel )
136 {
168e9c4d 137 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
de01cdf0 138 if (!pad.IsValid()) continue;
139
046e5fd4 140 for ( Int_t i = 0; i < dim; ++i )
de01cdf0 141 {
046e5fd4 142 Float_t x0 = value->ValueAsFloat(manuChannel,i);
143
144 xmin[i] = TMath::Min(xmin[i],x0);
145 xmax[i] = TMath::Max(xmax[i],x0);
de01cdf0 146 }
147 }
046e5fd4 148 }
149
150 for ( Int_t i = 0; i < dim; ++i )
151 {
152 if ( TMath::Abs(xmin[i]-xmax[i]) < 1E-3 )
153 {
154 xmin[i] -= 1;
155 xmax[i] += 1;
156 }
157 }
de01cdf0 158}
159
cf27231a 160//_____________________________________________________________________________
161Double_t GetRandom(Double_t mean, Double_t sigma, Bool_t mustBePositive)
162{
163 Double_t x(-1);
164 if ( mustBePositive )
165 {
166 while ( x < 0 )
167 {
168 x = gRandom->Gaus(mean,sigma);
169 }
170 }
171 else
172 {
173 x = gRandom->Gaus(mean,sigma);
174 }
175 return x;
176}
177
de01cdf0 178}
179
180//_____________________________________________________________________________
a99c3449 181Bool_t AliMUONCDB::CheckOCDB(Bool_t pathOnly)
de01cdf0 182{
a99c3449 183 /// Check that OCDB path and run number are properly set
184
185 AliCDBManager* man = AliCDBManager::Instance();
186
187 // first OCDB path
188 if (!man->IsDefaultStorageSet()) {
189 AliErrorGeneral("AliMUONCDB", "OCDB path must be properly set");
190 return kFALSE;
191 }
192
193 // then run number if required
194 if (pathOnly) return kTRUE;
195 if (man->GetRun() < 0) {
196 AliErrorGeneral("AliMUONCDB", "Run number must be properly set");
197 return kFALSE;
198 }
199
200 return kTRUE;
201
202}
203
204//_____________________________________________________________________________
205Bool_t AliMUONCDB::CheckMapping(Bool_t segmentationOnly)
206{
207 /// Check that the mapping has been loaded
208
209 // first the segmentation
210 if (!AliMpSegmentation::Instance(false)) {
211 AliErrorGeneral("AliMUONCDB", "Mapping segmentation must be loaded first");
212 return kFALSE;
213 }
214
215 // then the others if required
216 if (segmentationOnly) return kTRUE;
217 if (!AliMpDDLStore::Instance(false) || !AliMpDEStore::Instance(false) || !AliMpManuStore::Instance(false)) {
218 AliErrorGeneral("AliMUONCDB", "Full mapping must be loaded first");
219 return kFALSE;
220 }
221
222 return kTRUE;
223
224}
225
226//_____________________________________________________________________________
227Bool_t AliMUONCDB::LoadField()
228{
229 /// Load magnetic field (existing field will be deleted).
230 /// OCDB path and run number are supposed to be set.
231
232 AliInfoGeneral("AliMUONCDB","Loading field map from GRP...");
233
234 if (!AliMUONCDB::CheckOCDB()) return kFALSE;
235
236 AliGRPManager grpMan;
237
238 // make sure the old field is deleted even if it is locked
239 if(TGeoGlobalMagField::Instance()->IsLocked()) delete TGeoGlobalMagField::Instance();
240
241 if (!grpMan.ReadGRPEntry() || !grpMan.SetMagField()) {
242 AliErrorGeneral("AliMUONCDB", "failed to load magnetic field from OCDB");
243 return kFALSE;
244 }
245
246 return kTRUE;
247
248}
ab167304 249
a99c3449 250//_____________________________________________________________________________
251Bool_t AliMUONCDB::LoadMapping(Bool_t segmentationOnly)
252{
253 /// Load mapping (existing mapping will be unloaded).
254 /// OCDB path and run number are supposed to be set.
255
256 AliInfoGeneral("AliMUONCDB", "Loading mapping from OCDB...");
257
258 if (!AliMUONCDB::CheckOCDB()) return kFALSE;
259
260 // in case it has already been set
261 AliMpCDB::UnloadAll();
262
263 if (segmentationOnly) {
264
265 if (!AliMpCDB::LoadMpSegmentation(kTRUE)){
266 AliErrorGeneral("AliMUONCDB", "failed to load segmentation from OCDB");
267 return kFALSE;
268 }
269
270 } else {
271
272 if (!AliMpCDB::LoadAll(kTRUE)) {
273 AliErrorGeneral("AliMUONCDB", "failed to load mapping from OCDB");
274 return kFALSE;
ab167304 275 }
a99c3449 276
277 }
278
279 return kTRUE;
280
de01cdf0 281}
282
283//_____________________________________________________________________________
a99c3449 284AliMUONRecoParam* AliMUONCDB::LoadRecoParam()
de01cdf0 285{
a99c3449 286 /// Load and return reconstruction parameters.
287 /// OCDB path is supposed to be set.
288
289 AliInfoGeneral("AliMUONCDB", "Loading RecoParam from OCDB...");
290
2ff3de0a 291 if (!AliMUONCDB::CheckOCDB()) return 0x0;
a99c3449 292
293 AliMUONRecoParam* recoParam = 0x0;
294 AliCDBEntry* entry = AliCDBManager::Instance()->Get("MUON/Calib/RecoParam");
295
296 if(entry) {
297
298 // load recoParam according OCDB content (single or array)
299 if (!(recoParam = dynamic_cast<AliMUONRecoParam*>(entry->GetObject()))) {
300
301 TObjArray* recoParamArray = static_cast<TObjArray*>(entry->GetObject());
cdffeaea 302// recoParamArray->SetOwner(kTRUE); // FIXME: this should be done, but is causing a problem at the end of the reco... investigate why...
a99c3449 303
304 for(Int_t i = 0; i < recoParamArray->GetEntriesFast(); i++) {
305 recoParam = static_cast<AliMUONRecoParam*>(recoParamArray->UncheckedAt(i));
306 if (recoParam->IsDefault()) break;
307 recoParam = 0x0;
308 }
309
310 }
311
312 }
313
314 if (!recoParam) AliErrorGeneral("AliMUONCDB", "failed to load RecoParam from OCDB");
315
316 return recoParam;
317
a0eca509 318}
319
320//_____________________________________________________________________________
73d3aca5 321TClonesArray* AliMUONCDB::LoadAlignmentData()
322{
323 /// Load and return the array of alignment objects.
324
325 AliInfoGeneral("AliMUONCDB", "Loading Alignemnt from OCDB...");
326
2ff3de0a 327 if (!AliMUONCDB::CheckOCDB()) return 0x0;
73d3aca5 328
329 TClonesArray* alignmentArray = 0x0;
330 AliCDBEntry* entry = AliCDBManager::Instance()->Get("MUON/Align/Data");
331
332 if (entry) {
333 // load alignement array
334 alignmentArray = dynamic_cast<TClonesArray*>(entry->GetObject());
335 }
336
337 if (!alignmentArray) {
338 AliErrorGeneral("AliMUONCDB", "failed to load Alignemnt from OCDB");
339 }
340
341 return alignmentArray;
342}
343
344//_____________________________________________________________________________
a0eca509 345AliMUONVStore*
346AliMUONCDB::Diff(AliMUONVStore& store1, AliMUONVStore& store2,
de01cdf0 347 const char* opt)
348{
349 /// creates a store which contains store1-store2
350 /// if opt="abs" the difference is absolute one,
351 /// if opt="rel" then what is stored is (store1-store2)/store1
046e5fd4 352 /// if opt="percent" then what is stored is rel*100
353 ///
de01cdf0 354 /// WARNING Works only for stores which holds AliMUONVCalibParam objects
355
356 TString sopt(opt);
357 sopt.ToUpper();
358
046e5fd4 359 if ( !sopt.Contains("ABS") && !sopt.Contains("REL") && !sopt.Contains("PERCENT") )
de01cdf0 360 {
a99c3449 361 AliErrorGeneral("AliMUONCDB", Form("opt %s not supported. Only ABS, REL, PERCENT are",opt));
de01cdf0 362 return 0x0;
363 }
364
a0eca509 365 AliMUONVStore* d = static_cast<AliMUONVStore*>(store1.Clone());
de01cdf0 366
a0eca509 367 TIter next(d->CreateIterator());
de01cdf0 368
a0eca509 369 AliMUONVCalibParam* param;
de01cdf0 370
a0eca509 371 while ( ( param = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
de01cdf0 372 {
a0eca509 373 Int_t detElemId = param->ID0();
374 Int_t manuId = param->ID1();
de01cdf0 375
a0eca509 376 AliMUONVCalibParam* param2 = dynamic_cast<AliMUONVCalibParam*>(store2.FindObject(detElemId,manuId));
de01cdf0 377 //FIXME: this might happen. Handle it.
378 if (!param2)
379 {
380 cerr << "param2 is null : FIXME : this might happen !" << endl;
381 delete d;
382 return 0;
383 }
384
385 for ( Int_t i = 0; i < param->Size(); ++i )
386 {
387 for ( Int_t j = 0; j < param->Dimension(); ++j )
388 {
389 Float_t value(0);
390 if ( sopt.Contains("ABS") )
391 {
392 value = param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j);
393 }
046e5fd4 394 else if ( sopt.Contains("REL") || sopt.Contains("PERCENT") )
de01cdf0 395 {
396 if ( param->ValueAsFloat(i,j) )
397 {
398 value = (param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j))/param->ValueAsFloat(i,j);
399 }
400 else
401 {
402 continue;
403 }
046e5fd4 404 if ( sopt.Contains("PERCENT") ) value *= 100.0;
de01cdf0 405 }
406 param->SetValueAsFloat(i,j,value);
407 }
408 }
de01cdf0 409 }
410 return d;
411}
412
413//_____________________________________________________________________________
b37b9546 414TH1**
a0eca509 415AliMUONCDB::Plot(const AliMUONVStore& store, const char* name, Int_t nbins)
de01cdf0 416{
046e5fd4 417 /// Make histograms of each dimension of the AliMUONVCalibParam
de01cdf0 418 /// contained inside store.
419 /// It produces histograms named name_0, name_1, etc...
420
b37b9546 421 if (!AliMUONCDB::CheckMapping(kTRUE)) return 0x0;
a99c3449 422
046e5fd4 423 TIter next(store.CreateIterator());
424 AliMUONVCalibParam* param;
de01cdf0 425 Int_t n(0);
046e5fd4 426 const Int_t kNStations = AliMpConstants::NofTrackingChambers()/2;
427 Int_t* nPerStation = new Int_t[kNStations];
428 TH1** h(0x0);
de01cdf0 429
046e5fd4 430 for ( Int_t i = 0; i < kNStations; ++i ) nPerStation[i]=0;
de01cdf0 431
046e5fd4 432 while ( ( param = static_cast<AliMUONVCalibParam*>(next()) ) )
de01cdf0 433 {
046e5fd4 434 if (!h)
435 {
436 Int_t dim = param->Dimension();
437 h = new TH1*[dim];
438 Float_t* xmin = new Float_t[dim];
439 Float_t* xmax = new Float_t[dim];
440 getBoundaries(store,dim,xmin,xmax);
441
442 for ( Int_t i = 0; i < dim; ++i )
443 {
444 h[i] = new TH1F(Form("%s_%d",name,i),Form("%s_%d",name,i),
445 nbins,xmin[i],xmax[i]);
a99c3449 446 AliInfoGeneral("AliMUONCDB", Form("Created histogram %s",h[i]->GetName()));
046e5fd4 447 }
f1bf7efd 448 delete [] xmin;
449 delete [] xmax;
046e5fd4 450 }
451
452 Int_t detElemId = param->ID0();
453 Int_t manuId = param->ID1();
454 Int_t station = AliMpDEManager::GetChamberId(detElemId)/2;
de01cdf0 455
456 const AliMpVSegmentation* seg =
457 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
458
b37b9546 459 if (!seg) continue;
460
046e5fd4 461 for ( Int_t manuChannel = 0; manuChannel < param->Size(); ++manuChannel )
de01cdf0 462 {
168e9c4d 463 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
046e5fd4 464 if (!pad.IsValid()) continue;
de01cdf0 465
046e5fd4 466 ++n;
467 ++nPerStation[station];
468
469 for ( Int_t dim = 0; dim < param->Dimension(); ++dim )
470 {
471 h[dim]->Fill(param->ValueAsFloat(manuChannel,dim));
de01cdf0 472 }
473 }
046e5fd4 474 }
de01cdf0 475
046e5fd4 476 for ( Int_t i = 0; i < kNStations; ++i )
de01cdf0 477 {
a99c3449 478 AliInfoGeneral("AliMUONCDB", Form("Station %d %d ",(i+1),nPerStation[i]));
de01cdf0 479 }
046e5fd4 480
a99c3449 481 AliInfoGeneral("AliMUONCDB", Form("Number of channels = %d",n));
046e5fd4 482
483 delete[] nPerStation;
b37b9546 484
485 return h;
de01cdf0 486}
487
488//_____________________________________________________________________________
489Int_t
490AliMUONCDB::MakeHVStore(TMap& aliasMap, Bool_t defaultValues)
491{
492 /// Create a HV store
493
a99c3449 494 if (!AliMUONCDB::CheckMapping()) return 0;
495
49e110ec 496 AliMpDCSNamer hvNamer("TRACKER");
de01cdf0 497
498 TObjArray* aliases = hvNamer.GenerateAliases();
499
500 Int_t nSwitch(0);
501 Int_t nChannels(0);
502
503 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
504 {
505 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
506 TString& aliasName = alias->String();
507 if ( aliasName.Contains("sw") )
508 {
509 // HV Switch (St345 only)
510 TObjArray* valueSet = new TObjArray;
511 valueSet->SetOwner(kTRUE);
512
513 Bool_t value = kTRUE;
514
515 if (!defaultValues)
516 {
a0eca509 517 Float_t r = gRandom->Uniform();
de01cdf0 518 if ( r < 0.007 ) value = kFALSE;
519 }
520
521 for ( UInt_t timeStamp = 0; timeStamp < 60*3; timeStamp += 60 )
522 {
523 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
524 valueSet->Add(dcsValue);
525 }
526 aliasMap.Add(new TObjString(*alias),valueSet);
527 ++nSwitch;
528 }
529 else
530 {
531 TObjArray* valueSet = new TObjArray;
532 valueSet->SetOwner(kTRUE);
533 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
534 {
535 Float_t value = 1500;
cf27231a 536 if (!defaultValues) value = GetRandom(1750,62.5,true);
de01cdf0 537 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
538 valueSet->Add(dcsValue);
539 }
540 aliasMap.Add(new TObjString(*alias),valueSet);
541 ++nChannels;
542 }
543 }
544
545 delete aliases;
546
a99c3449 547 AliInfoGeneral("AliMUONCDB", Form("%d HV channels and %d switches",nChannels,nSwitch));
de01cdf0 548
549 return nChannels+nSwitch;
550}
551
552//_____________________________________________________________________________
553Int_t
49e110ec 554AliMUONCDB::MakeTriggerDCSStore(TMap& aliasMap, Bool_t defaultValues)
555{
556 /// Create a Trigger HV and Currents store
557
a99c3449 558 if (!AliMUONCDB::CheckMapping()) return 0;
559
49e110ec 560 AliMpDCSNamer triggerDCSNamer("TRIGGER");
561
562 TObjArray* aliases = triggerDCSNamer.GenerateAliases();
563
564 Int_t nChannels[2] = {0, 0};
565
566 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
567 {
568 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
569 TString& aliasName = alias->String();
570
571 TObjArray* valueSet = new TObjArray;
572 valueSet->SetOwner(kTRUE);
ad30b53f 573
49e110ec 574 Int_t measureType = triggerDCSNamer.DCSvariableFromDCSAlias(aliasName.Data());
ad30b53f 575 if ( measureType < 0 ) {
576 AliErrorGeneralStream("AliMUONCDB")
577 << "Failed to get DCS variable from an alias (trigger): "
578 << aliasName.Data() << endl;
579 return 0;
580 }
581
49e110ec 582 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
583 {
584 Float_t value =
9b1e069f 585 (measureType == AliMpDCSNamer::kDCSI) ? 2. : 8000.;
49e110ec 586 if (!defaultValues) {
587 switch (measureType){
588 case AliMpDCSNamer::kDCSI:
9b1e069f 589 value = GetRandom(2.,0.4,true);
49e110ec 590 break;
591 case AliMpDCSNamer::kDCSHV:
9b1e069f 592 value = GetRandom(8000.,16.,true);
49e110ec 593 break;
594 }
595 }
596 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
597 valueSet->Add(dcsValue);
598 }
599 aliasMap.Add(new TObjString(*alias),valueSet);
600 ++nChannels[measureType];
601 }
602
603 delete aliases;
604
a99c3449 605 AliInfoGeneral("AliMUONCDB", Form("Trigger channels I -> %i HV -> %i",nChannels[0], nChannels[1]));
49e110ec 606
607 return nChannels[0] + nChannels[1];
608}
609
610//_____________________________________________________________________________
611Int_t
a0eca509 612AliMUONCDB::MakePedestalStore(AliMUONVStore& pedestalStore, Bool_t defaultValues)
de01cdf0 613{
614 /// Create a pedestal store. if defaultValues=true, ped.mean=ped.sigma=1,
615 /// otherwise mean and sigma are from a gaussian (with parameters
616 /// defined below by the kPedestal* constants)
409de10f 617
a99c3449 618 AliCodeTimerAutoGeneral("",0);
619
620 if (!AliMUONCDB::CheckMapping()) return 0;
de01cdf0 621
622 Int_t nchannels(0);
623 Int_t nmanus(0);
624
a0eca509 625 const Int_t kChannels(AliMpConstants::ManuNofChannels());
e7ee8d13 626
627 // bending
628 const Float_t kPedestalMeanMeanB(200.);
629 const Float_t kPedestalMeanSigmaB(10.);
630 const Float_t kPedestalSigmaMeanB(1.);
631 const Float_t kPedestalSigmaSigmaB(0.2);
632
633 // non bending
634 const Float_t kPedestalMeanMeanNB(200.);
635 const Float_t kPedestalMeanSigmaNB(10.);
636 const Float_t kPedestalSigmaMeanNB(1.);
637 const Float_t kPedestalSigmaSigmaNB(0.2);
638
639 const Float_t kFractionOfDeadManu(0.); // within [0.,1.]
409de10f 640
641 Int_t detElemId;
642 Int_t manuId;
643
644 AliMpManuIterator it;
de01cdf0 645
409de10f 646 while ( it.Next(detElemId,manuId) )
de01cdf0 647 {
e7ee8d13 648 // skip a given fraction of manus
649 if (kFractionOfDeadManu > 0. && gRandom->Uniform() < kFractionOfDeadManu) continue;
650
de01cdf0 651 ++nmanus;
a0eca509 652
046e5fd4 653 AliMUONVCalibParam* ped =
654 new AliMUONCalibParamNF(2,kChannels,detElemId,manuId,AliMUONVCalibParam::InvalidFloatValue());
409de10f 655
656 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
de01cdf0 657
658 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
659 {
409de10f 660 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
de01cdf0 661
662 ++nchannels;
663
664 Float_t meanPedestal;
665 Float_t sigmaPedestal;
666
667 if ( defaultValues )
668 {
669 meanPedestal = 0.0;
670 sigmaPedestal = 1.0;
671 }
672 else
673 {
cf27231a 674 Bool_t positive(kTRUE);
046e5fd4 675 meanPedestal = 0.0;
e7ee8d13 676
677 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) ) { // manu in non bending plane
678
679 while ( meanPedestal == 0.0 ) // avoid strict zero
680 {
681 meanPedestal = GetRandom(kPedestalMeanMeanNB,kPedestalMeanSigmaNB,positive);
682 }
683 sigmaPedestal = GetRandom(kPedestalSigmaMeanNB,kPedestalSigmaSigmaNB,positive);
684
685 } else { // manu in bending plane
686
687 while ( meanPedestal == 0.0 ) // avoid strict zero
688 {
689 meanPedestal = GetRandom(kPedestalMeanMeanB,kPedestalMeanSigmaB,positive);
690 }
691 sigmaPedestal = GetRandom(kPedestalSigmaMeanB,kPedestalSigmaSigmaB,positive);
692
693 }
694
de01cdf0 695 }
e7ee8d13 696
de01cdf0 697 ped->SetValueAsFloat(manuChannel,0,meanPedestal);
698 ped->SetValueAsFloat(manuChannel,1,sigmaPedestal);
699
700 }
a0eca509 701 Bool_t ok = pedestalStore.Add(ped);
de01cdf0 702 if (!ok)
703 {
a99c3449 704 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
de01cdf0 705 }
706 }
707
a99c3449 708 AliInfoGeneral("AliMUONCDB", Form("%d Manus and %d channels.",nmanus,nchannels));
de01cdf0 709 return nchannels;
de01cdf0 710}
711
712//_____________________________________________________________________________
0045b488 713AliMUONRejectList*
714AliMUONCDB::MakeRejectListStore(Bool_t defaultValues)
715{
716 /// Create a reject list
717
a99c3449 718 AliCodeTimerAutoGeneral("",0);
0045b488 719
720 AliMUONRejectList* rl = new AliMUONRejectList;
721
722 if (!defaultValues)
723 {
724 rl->SetDetectionElementProbability(510);
725 rl->SetDetectionElementProbability(508);
726 return rl;
727 }
728
729 return rl;
730}
731
732//_____________________________________________________________________________
2b8a1212 733Int_t
7eafe398 734AliMUONCDB::MakeOccupancyMapStore(AliMUONVStore& occupancyMapStore, Bool_t defaultValues)
2b8a1212 735{
7eafe398 736 /// Create an occupancy map.
2b8a1212 737
a99c3449 738 AliCodeTimerAutoGeneral("",0);
739
740 if (!AliMUONCDB::CheckMapping()) return 0;
2b8a1212 741
2b8a1212 742 Int_t nmanus(0);
743
2b8a1212 744 Int_t detElemId;
745 Int_t manuId;
746
747 AliMpManuIterator it;
748
7eafe398 749 Int_t nevents(1000);
750
2b8a1212 751 while ( it.Next(detElemId,manuId) )
752 {
2b8a1212 753 ++nmanus;
754
7eafe398 755 AliMUONVCalibParam* occupancy = new AliMUONCalibParamND(5,1,detElemId,manuId,0);
2b8a1212 756
7eafe398 757 Double_t occ = 0.0;
758
2b8a1212 759 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
760
7eafe398 761 Int_t numberOfChannelsInManu = de->NofChannelsInManu(manuId);
762
a99c3449 763 if (!defaultValues) occ = gRandom->Rndm(1);
7eafe398 764
765 Double_t sumn = occ*nevents;
766
767 occupancy->SetValueAsFloat(0,0,sumn);
768 occupancy->SetValueAsFloat(0,1,sumn);
769 occupancy->SetValueAsFloat(0,2,sumn);
770 occupancy->SetValueAsInt(0,3,numberOfChannelsInManu);
771 occupancy->SetValueAsInt(0,4,nevents);
772
773 Bool_t ok = occupancyMapStore.Add(occupancy);
2b8a1212 774 if (!ok)
775 {
a99c3449 776 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
2b8a1212 777 }
2b8a1212 778 }
779
7eafe398 780 return nmanus;
2b8a1212 781}
782
783//_____________________________________________________________________________
25e1df3e 784Int_t
785AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, const char* file)
786{
6687de6d 787 /// Read the capacitance values from file and append them to the capaStore
788
a99c3449 789 if (!AliMUONCDB::CheckMapping()) return 0;
790
49419555 791 return AliMUONTrackerIO::ReadCapacitances(file,capaStore);
25e1df3e 792}
793
794//_____________________________________________________________________________
de01cdf0 795Int_t
a0eca509 796AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, Bool_t defaultValues)
de01cdf0 797{
798 /// Create a capacitance store. if defaultValues=true, all capa are 1.0,
799 /// otherwise they are from a gaussian with parameters defined in the
800 /// kCapa* constants below.
409de10f 801
a99c3449 802 AliCodeTimerAutoGeneral("",0);
803
804 if (!AliMUONCDB::CheckMapping()) return 0;
de01cdf0 805
806 Int_t nchannels(0);
807 Int_t nmanus(0);
808 Int_t nmanusOK(0); // manus for which we got the serial number
a0eca509 809
25e1df3e 810 const Float_t kCapaMean(0.3);
811 const Float_t kCapaSigma(0.1);
812 const Float_t kInjectionGainMean(3);
813 const Float_t kInjectionGainSigma(1);
409de10f 814
815 Int_t detElemId;
816 Int_t manuId;
817
818 AliMpManuIterator it;
de01cdf0 819
409de10f 820 while ( it.Next(detElemId,manuId) )
de01cdf0 821 {
822 ++nmanus;
823
de01cdf0 824 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
ab167304 825 Int_t serialNumber = AliMpManuStore::Instance()->GetManuSerial(detElemId, manuId);
de01cdf0 826
a0eca509 827 if ( serialNumber <= 0 ) continue;
828
829 ++nmanusOK;
de01cdf0 830
a0eca509 831 AliMUONVCalibParam* capa = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
de01cdf0 832
833 if (!capa)
834 {
25e1df3e 835 capa = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
a0eca509 836 Bool_t ok = capaStore.Add(capa);
de01cdf0 837 if (!ok)
838 {
a99c3449 839 AliErrorGeneral("AliMUONCDB", Form("Could not set serialNumber=%d manuId=%d",serialNumber,manuId));
de01cdf0 840 }
841 }
842
843 for ( Int_t manuChannel = 0; manuChannel < capa->Size(); ++manuChannel )
844 {
409de10f 845 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
de01cdf0 846
847 ++nchannels;
848
849 Float_t capaValue;
25e1df3e 850 Float_t injectionGain;
de01cdf0 851
852 if ( defaultValues )
853 {
854 capaValue = 1.0;
25e1df3e 855 injectionGain = 1.0;
de01cdf0 856 }
857 else
858 {
cf27231a 859 capaValue = GetRandom(kCapaMean,kCapaSigma,kTRUE);
25e1df3e 860 injectionGain = GetRandom(kInjectionGainMean,kInjectionGainSigma,kTRUE);
de01cdf0 861 }
862 capa->SetValueAsFloat(manuChannel,0,capaValue);
25e1df3e 863 capa->SetValueAsFloat(manuChannel,1,injectionGain);
de01cdf0 864 }
865 }
866
867 Float_t percent = 0;
868 if ( nmanus ) percent = 100*nmanusOK/nmanus;
a99c3449 869 AliInfoGeneral("AliMUONCDB", Form("%5d manus with serial number (out of %5d manus = %3.0f%%)",
de01cdf0 870 nmanusOK,nmanus,percent));
a99c3449 871 AliInfoGeneral("AliMUONCDB", Form("%5d channels",nchannels));
de01cdf0 872 if ( percent < 100 )
873 {
a99c3449 874 AliWarningGeneral("AliMUONCDB", "Did not get all serial numbers. capaStore is incomplete !!!!");
de01cdf0 875 }
876 return nchannels;
877
878}
879
880//_____________________________________________________________________________
881Int_t
a0eca509 882AliMUONCDB::MakeGainStore(AliMUONVStore& gainStore, Bool_t defaultValues)
de01cdf0 883{
cf27231a 884 /// Create a gain store. if defaultValues=true, all gains set so that
885 /// charge = (adc-ped)
886 ///
887 /// otherwise parameters are taken from gaussians with parameters
888 /// defined in the k* constants below.
de01cdf0 889
a99c3449 890 AliCodeTimerAutoGeneral("",0);
891
892 if (!AliMUONCDB::CheckMapping()) return 0;
de01cdf0 893
894 Int_t nchannels(0);
895 Int_t nmanus(0);
a0eca509 896
cf27231a 897 const Int_t kSaturation(3000);
898 const Double_t kA0Mean(1.2);
409de10f 899 const Double_t kA0Sigma(0.1);
900 const Double_t kA1Mean(1E-5);
901 const Double_t kA1Sigma(1E-6);
cf27231a 902 const Double_t kQualMean(0xFF);
903 const Double_t kQualSigma(0x10);
904 const Int_t kThresMean(1600);
409de10f 905 const Int_t kThresSigma(100);
906
907 Int_t detElemId;
908 Int_t manuId;
de01cdf0 909
409de10f 910 AliMpManuIterator it;
911
912 while ( it.Next(detElemId,manuId) )
de01cdf0 913 {
914 ++nmanus;
de01cdf0 915
a0eca509 916 AliMUONVCalibParam* gain =
cf27231a 917 new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),
a0eca509 918 detElemId,
919 manuId,
920 AliMUONVCalibParam::InvalidFloatValue());
921
409de10f 922 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
de01cdf0 923
924 for ( Int_t manuChannel = 0; manuChannel < gain->Size(); ++manuChannel )
925 {
409de10f 926 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
de01cdf0 927
928 ++nchannels;
929
de01cdf0 930 if ( defaultValues )
931 {
cf27231a 932 gain->SetValueAsFloat(manuChannel,0,1.0);
933 gain->SetValueAsFloat(manuChannel,1,0.0);
934 gain->SetValueAsInt(manuChannel,2,4095);
935 gain->SetValueAsInt(manuChannel,3,1);
936 gain->SetValueAsInt(manuChannel,4,kSaturation);
de01cdf0 937 }
938 else
939 {
cf27231a 940 Bool_t positive(kTRUE);
941 gain->SetValueAsFloat(manuChannel,0,GetRandom(kA0Mean,kA0Sigma,positive));
942 gain->SetValueAsFloat(manuChannel,1,GetRandom(kA1Mean,kA1Sigma,!positive));
943 gain->SetValueAsInt(manuChannel,2,(Int_t)TMath::Nint(GetRandom(kThresMean,kThresSigma,positive)));
944 gain->SetValueAsInt(manuChannel,3,(Int_t)TMath::Nint(GetRandom(kQualMean,kQualSigma,positive)));
945 gain->SetValueAsInt(manuChannel,4,kSaturation);
de01cdf0 946 }
de01cdf0 947
948 }
a0eca509 949 Bool_t ok = gainStore.Add(gain);
de01cdf0 950 if (!ok)
951 {
a99c3449 952 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
de01cdf0 953 }
954 }
955
a99c3449 956 AliInfoGeneral("AliMUONCDB", Form("%d Manus and %d channels.",nmanus,nchannels));
de01cdf0 957 return nchannels;
958}
959
960//_____________________________________________________________________________
961Int_t
a99c3449 962AliMUONCDB::MakeLocalTriggerMaskStore(AliMUONVStore& localBoardMasks)
de01cdf0 963{
964 /// Generate local trigger masks store. All masks are set to FFFF
965
a99c3449 966 AliCodeTimerAutoGeneral("",0);
409de10f 967
de01cdf0 968 Int_t ngenerated(0);
630711ed 969 // Generate fake mask values for all localboards and put that into
de01cdf0 970 // one single container (localBoardMasks)
0145e89a 971 for ( Int_t i = 1; i <= AliMpConstants::TotalNofLocalBoards(); ++i )
de01cdf0 972 {
a0eca509 973 AliMUONVCalibParam* localBoard = new AliMUONCalibParamNI(1,8,i,0,0);
de01cdf0 974 for ( Int_t x = 0; x < 2; ++x )
975 {
976 for ( Int_t y = 0; y < 4; ++y )
977 {
978 Int_t index = x*4+y;
979 localBoard->SetValueAsInt(index,0,0xFFFF);
980 ++ngenerated;
981 }
982 }
a0eca509 983 localBoardMasks.Add(localBoard);
de01cdf0 984 }
985 return ngenerated;
986}
987
988//_____________________________________________________________________________
989Int_t
a99c3449 990AliMUONCDB::MakeRegionalTriggerConfigStore(AliMUONRegionalTriggerConfig& rtm)
de01cdf0 991{
92c23b09 992 /// Make a regional trigger config store. Mask is set to FFFF for each local board (Ch.F.)
de01cdf0 993
a99c3449 994 AliCodeTimerAutoGeneral("",0);
409de10f 995
5cc125b2 996 if ( ! rtm.ReadData(AliMpFiles::LocalTriggerBoardMapping()) ) {
a99c3449 997 AliErrorGeneral("AliMUONCDB", "Error when reading from mapping file");
5cc125b2 998 return 0;
999 }
1000
1001 return rtm.GetNofTriggerCrates();
de01cdf0 1002}
1003
92c23b09 1004
de01cdf0 1005//_____________________________________________________________________________
1006Int_t
a99c3449 1007AliMUONCDB::MakeGlobalTriggerConfigStore(AliMUONGlobalCrateConfig& gtm)
de01cdf0 1008{
92c23b09 1009 /// Make a global trigger config store. All masks (disable) set to 0x00 for each Darc board (Ch.F.)
de01cdf0 1010
a99c3449 1011 AliCodeTimerAutoGeneral("",0);
913f1b43 1012
1013 return gtm.ReadData(AliMpFiles::GlobalTriggerBoardMapping());
de01cdf0 1014}
1015
92c23b09 1016
de01cdf0 1017//_____________________________________________________________________________
1018AliMUONTriggerLut*
a99c3449 1019AliMUONCDB::MakeTriggerLUT(const char* file)
de01cdf0 1020{
1021 /// Make a triggerlut object, from a file.
1022
a99c3449 1023 AliCodeTimerAutoGeneral("",0);
409de10f 1024
de01cdf0 1025 AliMUONTriggerLut* lut = new AliMUONTriggerLut;
1026 lut->ReadFromFile(file);
1027 return lut;
1028}
1029
1030//_____________________________________________________________________________
1031AliMUONTriggerEfficiencyCells*
a99c3449 1032AliMUONCDB::MakeTriggerEfficiency(const char* file)
de01cdf0 1033{
1034 /// Make a trigger efficiency object from a file.
1035
a99c3449 1036 AliCodeTimerAutoGeneral("",0);
409de10f 1037
de01cdf0 1038 return new AliMUONTriggerEfficiencyCells(file);
1039}
1040
1041//_____________________________________________________________________________
1042void
1043AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
25e1df3e 1044 Int_t startRun, Int_t endRun,
1045 const char* filename)
1046{
1047 /// Write a given object to OCDB
1048
913f1b43 1049 TString comment(gSystem->ExpandPathName(filename));
1050
1051 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
25e1df3e 1052}
1053
1054//_____________________________________________________________________________
1055void
1056AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
de01cdf0 1057 Int_t startRun, Int_t endRun, Bool_t defaultValues)
1058{
1059 /// Write a given object to OCDB
1060
913f1b43 1061 TString comment;
1062 if ( defaultValues ) comment += "Test with default values";
1063 else comment += "Test with random values";
1064
1065 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
1066}
1067
1068//_____________________________________________________________________________
1069void
1070AliMUONCDB::WriteToCDB(TObject* object, const char* calibpath, Int_t startRun, Int_t endRun,
1071 const char* comment, const char* responsible)
1072{
1073 /// Write a given object to OCDB
1074
a99c3449 1075 if (!AliMUONCDB::CheckOCDB(kTRUE)) return;
1076
de01cdf0 1077 AliCDBId id(calibpath,startRun,endRun);
1078 AliCDBMetaData md;
1079 md.SetAliRootVersion(gROOT->GetVersion());
913f1b43 1080 md.SetComment(comment);
1081 md.SetResponsible(responsible);
a99c3449 1082 AliCDBManager::Instance()->Put(object,id,&md);
de01cdf0 1083}
1084
1085//_____________________________________________________________________________
1086Int_t
a0eca509 1087AliMUONCDB::MakeNeighbourStore(AliMUONVStore& neighbourStore)
de01cdf0 1088{
1089 /// Fill the neighbours store with, for each channel, a TObjArray of its
1090 /// neighbouring pads (including itself)
1091
a99c3449 1092 AliCodeTimerAutoGeneral("",0);
1093
1094 if (!AliMUONCDB::CheckMapping()) return 0;
de01cdf0 1095
a99c3449 1096 AliInfoGeneral("AliMUONCDB", "Generating NeighbourStore. This will take a while. Please be patient.");
de01cdf0 1097
1098 Int_t nchannels(0);
1099
1100 TObjArray tmp;
409de10f 1101
1102 Int_t detElemId;
1103 Int_t manuId;
de01cdf0 1104
409de10f 1105 AliMpManuIterator it;
1106
1107 while ( it.Next(detElemId,manuId) )
de01cdf0 1108 {
de01cdf0 1109 const AliMpVSegmentation* seg =
1110 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
1111
a0eca509 1112 AliMUONVCalibParam* calibParam = static_cast<AliMUONVCalibParam*>(neighbourStore.FindObject(detElemId,manuId));
de01cdf0 1113 if (!calibParam)
1114 {
1115 Int_t dimension(11);
a0eca509 1116 Int_t size(AliMpConstants::ManuNofChannels());
de01cdf0 1117 Int_t defaultValue(-1);
a0eca509 1118 Int_t packingFactor(size);
de01cdf0 1119
a0eca509 1120 calibParam = new AliMUONCalibParamNI(dimension,size,detElemId,manuId,defaultValue,packingFactor);
1121 Bool_t ok = neighbourStore.Add(calibParam);
de01cdf0 1122 if (!ok)
1123 {
a99c3449 1124 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
de01cdf0 1125 return -1;
1126 }
1127 }
1128
a0eca509 1129 for ( Int_t manuChannel = 0; manuChannel < AliMpConstants::ManuNofChannels(); ++manuChannel )
de01cdf0 1130 {
168e9c4d 1131 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
de01cdf0 1132
1133 if (pad.IsValid())
1134 {
1135 ++nchannels;
1136
1137 seg->GetNeighbours(pad,tmp,true,true);
1138 Int_t nofPadNeighbours = tmp.GetEntriesFast();
1139
1140 for ( Int_t i = 0; i < nofPadNeighbours; ++i )
1141 {
bf0d3528 1142 AliMpPad* p = static_cast<AliMpPad*>(tmp.UncheckedAt(i));
de01cdf0 1143 Int_t x;
409de10f 1144// Bool_t ok =
6e97fbb8 1145 calibParam->PackValues(p->GetManuId(),p->GetManuChannel(),x);
409de10f 1146// if (!ok)
1147// {
1148// AliError("Could not pack value. Something is seriously wrong. Please check");
1149// StdoutToAliError(pad->Print(););
1150// return -1;
1151// }
de01cdf0 1152 calibParam->SetValueAsInt(manuChannel,i,x);
1153 }
1154 }
1155 }
1156 }
1157
de01cdf0 1158 return nchannels;
1159}
1160
1161//_____________________________________________________________________________
1162void
1163AliMUONCDB::WriteLocalTriggerMasks(Int_t startRun, Int_t endRun)
1164{
1165 /// Write local trigger masks to OCDB
1166
630711ed 1167 AliMUONVStore* ltm = new AliMUON1DArray(AliMpConstants::TotalNofLocalBoards()+1);
de01cdf0 1168 Int_t ngenerated = MakeLocalTriggerMaskStore(*ltm);
a99c3449 1169 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1170 if (ngenerated>0)
1171 {
1172 WriteToCDB("MUON/Calib/LocalTriggerBoardMasks",ltm,startRun,endRun,true);
1173 }
1174 delete ltm;
1175}
1176
1177//_____________________________________________________________________________
1178void
92c23b09 1179AliMUONCDB::WriteRegionalTriggerConfig(Int_t startRun, Int_t endRun)
de01cdf0 1180{
1181 /// Write regional trigger masks to OCDB
1182
92c23b09 1183 AliMUONRegionalTriggerConfig* rtm = new AliMUONRegionalTriggerConfig();
1184 Int_t ngenerated = MakeRegionalTriggerConfigStore(*rtm);
a99c3449 1185 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1186 if (ngenerated>0)
1187 {
92c23b09 1188 WriteToCDB("MUON/Calib/RegionalTriggerConfig",rtm,startRun,endRun,true);
de01cdf0 1189 }
1190 delete rtm;
1191}
1192
92c23b09 1193
de01cdf0 1194//_____________________________________________________________________________
1195void
92c23b09 1196AliMUONCDB::WriteGlobalTriggerConfig(Int_t startRun, Int_t endRun)
de01cdf0 1197{
1198 /// Write global trigger masks to OCDB
1199
92c23b09 1200 AliMUONGlobalCrateConfig* gtm = new AliMUONGlobalCrateConfig();
de01cdf0 1201
92c23b09 1202 Int_t ngenerated = MakeGlobalTriggerConfigStore(*gtm);
a99c3449 1203 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1204 if (ngenerated>0)
1205 {
92c23b09 1206 WriteToCDB("MUON/Calib/GlobalTriggerCrateConfig",gtm,startRun,endRun,true);
de01cdf0 1207 }
1208 delete gtm;
1209}
1210
92c23b09 1211
de01cdf0 1212//_____________________________________________________________________________
1213void
1214AliMUONCDB::WriteTriggerLut(Int_t startRun, Int_t endRun)
1215{
1216 /// Write trigger LUT to OCDB
1217
1218 AliMUONTriggerLut* lut = MakeTriggerLUT();
1219 if (lut)
1220 {
1221 WriteToCDB("MUON/Calib/TriggerLut",lut,startRun,endRun,true);
1222 }
1223 delete lut;
1224}
1225
1226//_____________________________________________________________________________
1227void
1228AliMUONCDB::WriteTriggerEfficiency(Int_t startRun, Int_t endRun)
1229{
1230 /// Write trigger efficiency to OCDB
1231
1232 AliMUONTriggerEfficiencyCells* eff = MakeTriggerEfficiency();
1233 if (eff)
1234 {
1235 WriteToCDB("MUON/Calib/TriggerEfficiency",eff,startRun,endRun,true);
1236 }
1237 delete eff;
1238}
1239
1240//_____________________________________________________________________________
1241void
1242AliMUONCDB::WriteNeighbours(Int_t startRun, Int_t endRun)
1243{
1244 /// Write neighbours to OCDB
1245
409de10f 1246 AliMUONVStore* neighbours = Create2DMap();
de01cdf0 1247 Int_t ngenerated = MakeNeighbourStore(*neighbours);
a99c3449 1248 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1249 if (ngenerated>0)
1250 {
1251 WriteToCDB("MUON/Calib/Neighbours",neighbours,startRun,endRun,true);
1252 }
1253 delete neighbours;
1254}
1255
1256//_____________________________________________________________________________
1257void
1258AliMUONCDB::WriteHV(Bool_t defaultValues,
1259 Int_t startRun, Int_t endRun)
1260{
1261 /// generate HV values (either cste = 1500 V) if defaultValues=true or random
1262 /// if defaultValues=false, see makeHVStore) and
1263 /// store them into CDB located at cdbpath, with a validity period
1264 /// ranging from startRun to endRun
1265
1266 TMap* hvStore = new TMap;
1267 Int_t ngenerated = MakeHVStore(*hvStore,defaultValues);
a99c3449 1268 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1269 if (ngenerated>0)
1270 {
1271 WriteToCDB("MUON/Calib/HV",hvStore,startRun,endRun,defaultValues);
1272 }
1273 delete hvStore;
1274}
1275
1276//_____________________________________________________________________________
1277void
49e110ec 1278AliMUONCDB::WriteTriggerDCS(Bool_t defaultValues,
1279 Int_t startRun, Int_t endRun)
1280{
1281 /// generate Trigger HV and current values (either const if defaultValues=true or random
1282 /// if defaultValues=false, see makeTriggerDCSStore) and
1283 /// store them into CDB located at cdbpath, with a validity period
1284 /// ranging from startRun to endRun
1285
1286 TMap* triggerDCSStore = new TMap;
1287 Int_t ngenerated = MakeTriggerDCSStore(*triggerDCSStore,defaultValues);
a99c3449 1288 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
49e110ec 1289 if (ngenerated>0)
1290 {
1291 WriteToCDB("MUON/Calib/TriggerDCS",triggerDCSStore,startRun,endRun,defaultValues);
1292 }
1293 delete triggerDCSStore;
1294}
1295
1296//_____________________________________________________________________________
1297void
de01cdf0 1298AliMUONCDB::WritePedestals(Bool_t defaultValues,
1299 Int_t startRun, Int_t endRun)
1300{
1301 /// generate pedestal values (either 0 if defaultValues=true or random
1302 /// if defaultValues=false, see makePedestalStore) and
1303 /// store them into CDB located at cdbpath, with a validity period
1304 /// ranging from startRun to endRun
1305
409de10f 1306 AliMUONVStore* pedestalStore = Create2DMap();
de01cdf0 1307 Int_t ngenerated = MakePedestalStore(*pedestalStore,defaultValues);
a99c3449 1308 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1309 WriteToCDB("MUON/Calib/Pedestals",pedestalStore,startRun,endRun,defaultValues);
1310 delete pedestalStore;
1311}
1312
2b8a1212 1313//_____________________________________________________________________________
1314void
7eafe398 1315AliMUONCDB::WriteOccupancyMap(Bool_t defaultValues,
1316 Int_t startRun, Int_t endRun)
2b8a1212 1317{
7eafe398 1318 /// generate occupancy map values (either empty one if defaultValues=true, or
1319 /// random one, see MakeOccupancyMapStore) and
2b8a1212 1320 /// store them into CDB located at cdbpath, with a validity period
1321 /// ranging from startRun to endRun
1322
7eafe398 1323 AliMUONVStore* occupancyMapStore = Create2DMap();
1324 Int_t ngenerated = MakeOccupancyMapStore(*occupancyMapStore,defaultValues);
a99c3449 1325 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
7eafe398 1326 WriteToCDB("MUON/Calib/OccupancyMap",occupancyMapStore,startRun,endRun,defaultValues);
1327 delete occupancyMapStore;
2b8a1212 1328}
1329
0045b488 1330//_____________________________________________________________________________
1331void
1332AliMUONCDB::WriteRejectList(Bool_t defaultValues,
1333 Int_t startRun, Int_t endRun)
1334{
1335 /// generate reject list values (either empty one if defaultValues=true, or
1336 /// random one, see MakeRejectListStore) and
1337 /// store them into CDB located at cdbpath, with a validity period
1338 /// ranging from startRun to endRun
1339
1340 AliMUONRejectList* rl = MakeRejectListStore(defaultValues);
1341 WriteToCDB("MUON/Calib/RejectList",rl,startRun,endRun,defaultValues);
1342 delete rl;
1343}
1344
de01cdf0 1345
1346//_____________________________________________________________________________
1347void
1348AliMUONCDB::WriteGains(Bool_t defaultValues,
1349 Int_t startRun, Int_t endRun)
1350{
1351 /// generate gain values (either 1 if defaultValues=true or random
1352 /// if defaultValues=false, see makeGainStore) and
1353 /// store them into CDB located at cdbpath, with a validity period
1354 /// ranging from startRun to endRun
1355
409de10f 1356 AliMUONVStore* gainStore = Create2DMap();
de01cdf0 1357 Int_t ngenerated = MakeGainStore(*gainStore,defaultValues);
a99c3449 1358 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1359 WriteToCDB("MUON/Calib/Gains",gainStore,startRun,endRun,defaultValues);
1360 delete gainStore;
1361}
1362
1363//_____________________________________________________________________________
1364void
25e1df3e 1365AliMUONCDB::WriteCapacitances(const char* filename,
1366 Int_t startRun, Int_t endRun)
1367{
1368 /// read manu capacitance and injection gain values from file
1369 /// and store them into CDB located at cdbpath, with a validity period
1370 /// ranging from startRun to endRun
1371
1372 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1373 Int_t ngenerated = MakeCapacitanceStore(*capaStore,filename);
a99c3449 1374 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
25e1df3e 1375 if ( ngenerated > 0 )
1376 {
1377 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,filename);
1378 }
1379 delete capaStore;
1380}
1381
1382//_____________________________________________________________________________
1383void
de01cdf0 1384AliMUONCDB::WriteCapacitances(Bool_t defaultValues,
1385 Int_t startRun, Int_t endRun)
1386{
1387 /// generate manu capacitance values (either 1 if defaultValues=true or random
1388 /// if defaultValues=false, see makeCapacitanceStore) and
1389 /// store them into CDB located at cdbpath, with a validity period
1390 /// ranging from startRun to endRun
1391
a0eca509 1392 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
de01cdf0 1393 Int_t ngenerated = MakeCapacitanceStore(*capaStore,defaultValues);
a99c3449 1394 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
de01cdf0 1395 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,defaultValues);
1396 delete capaStore;
1397}
1398
1399//_____________________________________________________________________________
1400void
49e110ec 1401AliMUONCDB::WriteTrigger(Bool_t defaultValues, Int_t startRun, Int_t endRun)
de01cdf0 1402{
1403 /// Writes all Trigger related calibration to CDB
49e110ec 1404 WriteTriggerDCS(defaultValues,startRun,endRun);
de01cdf0 1405 WriteLocalTriggerMasks(startRun,endRun);
92c23b09 1406 WriteRegionalTriggerConfig(startRun,endRun);
1407 WriteGlobalTriggerConfig(startRun,endRun);
de01cdf0 1408 WriteTriggerLut(startRun,endRun);
1409 WriteTriggerEfficiency(startRun,endRun);
1410}
1411
1412//_____________________________________________________________________________
1413void
a31df857 1414AliMUONCDB::WriteConfig(Int_t startRun, Int_t endRun)
1415{
1416 /// Write complete tracker configuration to OCDB
1417 ostringstream lines;
1418 TIter next(AliMpDDLStore::Instance()->CreateBusPatchIterator());
1419 AliMpBusPatch* bp;
1420 while ( ( bp = static_cast<AliMpBusPatch*>(next()) ) )
1421 {
1422 for (Int_t imanu = 0; imanu < bp->GetNofManus(); ++imanu)
1423 {
1424 lines << bp->GetId() << " " << bp->GetManuId(imanu) << endl;
1425 }
1426 }
1427
1428 AliMUON2DMap config(kTRUE);
a31df857 1429
042cd64e 1430 AliMUONTrackerIO::DecodeConfig(lines.str().c_str(),config);
a31df857 1431
1432 WriteToCDB("MUON/Calib/Config",&config,startRun,endRun,kTRUE);
1433}
1434
1435//_____________________________________________________________________________
1436void
de01cdf0 1437AliMUONCDB::WriteTracker(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1438{
1439 /// Writes all Tracker related calibration to CDB
1440 WriteHV(defaultValues,startRun,endRun);
1441 WritePedestals(defaultValues,startRun,endRun);
1442 WriteGains(defaultValues,startRun,endRun);
1443 WriteCapacitances(defaultValues,startRun,endRun);
1444 WriteNeighbours(startRun,endRun);
c21c5694 1445 WriteOccupancyMap(defaultValues,startRun,endRun);
1446 WriteRejectList(defaultValues,startRun,endRun);
a31df857 1447 WriteConfig(startRun,endRun);
de01cdf0 1448}
913f1b43 1449
5148da98 1450//_____________________________________________________________________________
1451void
c241bfaa 1452AliMUONCDB::ShowCapacitances()
1453{
1454 /// Show briefly the number of capa values we have in the OCDB,
1455 /// and the list of manu that are actually in the config and for which
1456 /// we miss the capa (if any).
1457
1458 if (!AliMUONCDB::CheckOCDB()) return;
1459
1460 AliMUONCDB::LoadMapping();
1461
1462 if (!AliMUONCDB::CheckMapping()) return;
1463
1464 AliCDBEntry* e = AliCDBManager::Instance()->Get("MUON/Calib/Config");
1465
1466 if (!e) return ;
1467
1468 AliMUONVStore* config = static_cast<AliMUONVStore*>(e->GetObject());
1469
1470 e = AliCDBManager::Instance()->Get("MUON/Calib/Capacitances");
1471
1472 if (!e) return;
1473
1474 AliMUONVStore* capacitances = static_cast<AliMUONVStore*>(e->GetObject());
1475
1476 AliInfoGeneral("ShowCapacitances",Form("%d capacitances are in OCDB",capacitances->GetSize()));
1477
1478 TIter nextManu(config->CreateIterator());
1479 AliMUONVCalibParam* param;
1480
1481 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1482 {
1483 Int_t detElemId = param->ID0();
1484 Int_t manuId = param->ID1();
1485
1486 Int_t serialNumber
1487 = AliMpManuStore::Instance()->GetManuSerial(detElemId, manuId);
1488
1489 if (serialNumber<0)
1490 {
1491 AliErrorGeneral("ShowCapacitances",Form("Did not find serial for DE %04d MANUID %04d",detElemId,manuId));
1492 }
1493 else
1494 {
1495 AliMUONVCalibParam* capa = static_cast<AliMUONVCalibParam*>(capacitances->FindObject(serialNumber));
1496 if (!capa)
1497 {
1498 AliErrorGeneral("ShowCapacitances",Form("Did not find capacitance for DE %04d MANUID %04d SERIAL %d",detElemId,manuId,serialNumber));
1499 }
1500 }
1501 }
1502
1503}
1504
1505//_____________________________________________________________________________
1506void
207947f7 1507AliMUONCDB::ShowConfig(Bool_t withStatusMap)
5148da98 1508{
1509 /// Dumps the current tracker configuration, i.e. number and identity of missing buspatches
207947f7 1510 /// If statusMap is true, will also take into account the status map to report the number
1511 /// of good channels
1512
5148da98 1513 if (!AliMUONCDB::CheckOCDB()) return;
1514
1515 AliMUONCDB::LoadMapping();
1516
1517 if (!AliMUONCDB::CheckMapping()) return;
1518
1519 AliCDBEntry* e = AliCDBManager::Instance()->Get("MUON/Calib/Config");
1520
1521 if (!e) return ;
1522
1523 AliMUONVStore* config = static_cast<AliMUONVStore*>(e->GetObject());
1524
207947f7 1525 AliMUONPadStatusMapMaker* statusMapMaker(0x0);
1526 AliMUONCalibrationData* cd(0x0);
1527 AliMUONPadStatusMaker* statusMaker(0x0);
1528
1529 if ( withStatusMap )
1530 {
1531 cd = new AliMUONCalibrationData(AliCDBManager::Instance()->GetRun());
1532
1533 statusMaker = new AliMUONPadStatusMaker(*cd);
1534
1535 AliMUONRecoParam* recoParam = AliMUONCDB::LoadRecoParam();
1536
0c53c46d 1537 if (!recoParam)
1538 {
1539 AliErrorGeneral("AliMUONCDB::ShowConfig","Cannot get recoParams from OCDB !");
1540 return;
1541 }
1542
207947f7 1543 statusMaker->SetLimits(*recoParam);
1544
207947f7 1545 UInt_t mask = recoParam->PadGoodnessMask();
9b780418 1546
1547 delete recoParam;
1548
207947f7 1549 const Bool_t deferredInitialization = kFALSE;
1550
1551 statusMapMaker = new AliMUONPadStatusMapMaker(*cd,mask,deferredInitialization);
1552 }
1553
5148da98 1554 TIter nextManu(config->CreateIterator());
1555 AliMUONVCalibParam* param;
1556
1557 AliMpExMap buspatches;
1558
1559 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1560 {
1561 Int_t detElemId = param->ID0();
1562 Int_t manuId = param->ID1();
1563 Int_t busPatchId = AliMpDDLStore::Instance()->GetBusPatchId(detElemId,manuId);
1564 if ( buspatches.GetValue(busPatchId) == 0x0 )
1565 {
1566 buspatches.Add(busPatchId,new TObjString(Form("BP%04d",busPatchId)));
1567 }
1568 }
1569
1570 TArrayI removed(buspatches.GetSize());
1571
1572 TIter next(AliMpDDLStore::Instance()->CreateBusPatchIterator());
1573 AliMpBusPatch* bp;
1574 Int_t n(0);
1575 Int_t nok(0);
1576 Int_t nremoved(0);
1577
207947f7 1578 // accounting of bus patches first
1579
5148da98 1580 while ( ( bp = static_cast<AliMpBusPatch*>(next())))
1581 {
1582 if ( buspatches.GetValue(bp->GetId()) )
1583 {
1584 ++nok;
1585 }
1586 else
1587 {
1588 removed.SetAt(bp->GetId(),nremoved++);
1589 }
5148da98 1590 }
1591
207947f7 1592 // accounting of channels
1593
1594 AliMpManuIterator it;
1595
1596 Int_t totalNumberOfChannels(0);
1597 Int_t removedChannels(0);
1598 Int_t badChannels(0);
1599 Int_t badAndRemovedChannels(0);
1600 Int_t badOrRemovedChannels(0);
1601
1602 Int_t detElemId, manuId;
1603
1604 while ( it.Next(detElemId,manuId) )
1605 {
1606 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
1607 for ( Int_t i = 0; i < AliMpConstants::ManuNofChannels(); ++i )
1608 {
1609 Int_t busPatchId = AliMpDDLStore::Instance()->GetBusPatchId(detElemId,manuId);
1610
1611 if ( de->IsConnectedChannel(manuId,i) )
1612 {
1613 ++totalNumberOfChannels;
1614 Bool_t badBusPatch = ( buspatches.GetValue(busPatchId) == 0x0 );
1615
1616 if ( withStatusMap )
1617 {
47274b13 1618 Bool_t badChannel = ( ( statusMapMaker->StatusMap(detElemId,manuId,i) & AliMUONPadStatusMapMaker::SelfDeadMask() ) != 0);
207947f7 1619 if ( badChannel ) ++badChannels;
1620 if ( badBusPatch && badChannel ) ++badAndRemovedChannels;
1621 if ( badBusPatch || badChannel ) ++badOrRemovedChannels;
1622 }
1623
1624 if ( badBusPatch) ++removedChannels;
1625 }
1626 }
1627 }
1628
5148da98 1629
1630 Int_t* indices = new Int_t[nremoved];
1631
1632 TMath::Sort(nremoved,removed.GetArray(),indices,kFALSE);
1633
1634 for ( Int_t i = 0; i < nremoved; ++i )
1635 {
1636 Int_t busPatchId = removed[indices[i]];
1637 bp = AliMpDDLStore::Instance()->GetBusPatch(busPatchId);
1638 bp->Print();
1639 }
1640
1641 delete[] indices;
207947f7 1642
1643 cout << endl;
1644 cout << Form("Bus patches n=%3d nok=%3d nremoved=%3d",n,nok,nremoved) << endl;
1645
1646 cout << Form("Channels n=%6d nremoved=%6d bad=%6d bad and removed=%6d bad or removed=%6d",
1647 totalNumberOfChannels,removedChannels,badChannels,badAndRemovedChannels,badOrRemovedChannels) << endl;
1648
1649 if (totalNumberOfChannels>0)
1650 {
1651 cout << Form("Percentage of readout channels %5.1f %%",removedChannels*100.0/totalNumberOfChannels) << endl;
1652 if ( withStatusMap )
1653 {
1654 cout << Form("Percentage of non useable channels (bad or removed) %5.1f %%",
1655 badOrRemovedChannels*100.0/totalNumberOfChannels) << endl;
1656 }
1657 }
1658
1659
1660 delete statusMapMaker;
1661 delete cd;
1662 delete statusMaker;
5148da98 1663}
9e7c4517 1664
1665//______________________________________________________________________________
1666void AliMUONCDB::ReadIntegers(const char* filename, std::vector<int>& integers)
1667{
1668 /// Read integers from filename, where integers are either
1669 /// separated by "," or by return carriage
1670 ifstream in(gSystem->ExpandPathName(filename));
1671 int i;
1672
1673 std::set<int> runset;
1674
1675 char line[10000];
1676
1677 in.getline(line,10000,'\n');
1678
1679 TString sline(line);
1680
1681 if (sline.Contains(","))
1682 {
1683 TObjArray* a = sline.Tokenize(",");
1684 TIter next(a);
1685 TObjString* s;
1686 while ( ( s = static_cast<TObjString*>(next()) ) )
1687 {
1688 runset.insert(s->String().Atoi());
1689 }
1690 delete a;
1691 }
1692 else
1693 {
1694 runset.insert(sline.Atoi());
1695
1696 while ( in >> i )
1697 {
1698 runset.insert(i);
1699 }
1700 }
1701
1702 for ( std::set<int>::const_iterator it = runset.begin(); it != runset.end(); ++it )
1703 {
1704 integers.push_back((*it));
1705 }
1706
1707 std::sort(integers.begin(),integers.end());
1708}
1709
1710//______________________________________________________________________________
1711void AliMUONCDB::ShowFaultyBusPatches(const char* runlist, double occLimit,
1712 const char* outputBaseName,
1713 const char* ocdbPath)
1714{
1715 /// Shows the list of bus patches above a given occupancy limit,
1716 /// for each run in the runlist
1717
1718 AliLog::GetRootLogger()->SetGlobalLogLevel(AliLog::kError);
1719
1720 // AliLog::SetPrintType(AliLog::kInfo,kFALSE);
1721 // AliLog::SetPrintType(AliLog::kWarning,kFALSE);
1722 // gErrorIgnoreLevel=kError; // to avoid all the TAlienFile::Open messages...
1723
1724 AliCDBManager* man = AliCDBManager::Instance();
1725
1726 man->SetDefaultStorage(ocdbPath);
1727
1728 Bool_t first(kTRUE);
1729
1730 std::vector<int> runnumbers;
1731
1732 ReadIntegers(runlist,runnumbers);
1733
1734 AliMUON2DMap bpValues(kFALSE);
1735
1736 std::ofstream outfile(Form("%s.txt",outputBaseName));
1737
1738 for ( unsigned int i = 0 ; i < runnumbers.size(); ++i )
1739 {
1740 int runNumber = runnumbers[i];
1741
1742 man->SetRun(runNumber);
1743
1744 if ( first )
1745 {
1746 AliMpCDB::LoadAll();
1747 first = kFALSE;
1748 }
1749
1750 AliCDBEntry* e = man->Get("MUON/Calib/OccupancyMap",runNumber);
1751
b37b9546 1752 if (!e)
1753 {
1754 AliErrorGeneral("AliMUONCDB::ShowFaultyBusPatches",
1755 Form("Could not get OccupancyMap for run %09d",runNumber));
1756 continue;
1757 }
1758
9e7c4517 1759 AliMUONVStore* occmap = static_cast<AliMUONVStore*>(e->GetObject());
1760
1761 AliMUONTrackerData td("occ","occ",*occmap);
1762
1763 TIter nextBP(AliMpDDLStore::Instance()->CreateBusPatchIterator());
1764 AliMpBusPatch* bp;
1765 std::set<int> buspatches;
1766 Double_t sumn = 1000.0;
1767
1768 while ( ( bp = static_cast<AliMpBusPatch*>(nextBP()) ) )
1769 {
1770 Double_t occ = td.BusPatch(bp->GetId(),2);
1771
1772 if (occ>occLimit)
1773 {
1774 buspatches.insert(bp->GetId());
1775
1776 AliMUONVCalibParam* param = static_cast<AliMUONVCalibParam*>(bpValues.FindObject(bp->GetId()));
1777
1778 if (!param)
1779 {
1780 param = new AliMUONCalibParamND(5, 1, bp->GetId(), 0);
1781 bpValues.Add(param);
1782
1783 Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(bp->GetId());
1784 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
1785
1786 Int_t nchannels(0);
1787
1788 for ( Int_t imanu = 0; imanu < bp->GetNofManus(); ++imanu )
1789 {
1790 Int_t manuId = bp->GetManuId(imanu);
1791 nchannels += de->NofChannelsInManu(manuId);
1792 }
1793
1794 param->SetValueAsDouble(0,2,sumn);
1795 param->SetValueAsDouble(0,3,nchannels);
1796 param->SetValueAsDouble(0,4,1);
1797 }
1798
1799 Double_t sumw = sumn*(param->ValueAsDouble(0)/sumn+1.0/runnumbers.size());
1800 Double_t sumw2 = 0.0; //(sumn-1)*ey*ey+sumw*sumw/sumn;
1801
1802 param->SetValueAsDouble(0,0,sumw);
1803 param->SetValueAsDouble(0,1,sumw2);
1804
1805 }
1806 }
1807
1808 outfile << Form("RUN %09d",runNumber);
1809
1810 for ( std::set<int>::const_iterator bit = buspatches.begin(); bit != buspatches.end(); ++bit )
1811 {
1812 outfile << Form(" %4d",*bit);
1813 }
1814 outfile << endl;
1815 }
1816
1817 outfile.close();
1818
1819 const char* name = "BPfailureRate";
1820
1821 AliMUONTrackerData* mpData = new AliMUONTrackerData(name,name,bpValues,2);
1822 mpData->SetDimensionName(0,name);
1823
1824 TFile f(Form("%s.root",outputBaseName),"recreate");
1825 mpData->Write();
1826 f.Close();
1827
1828 cout << Form("Results are in %s.txt and %s.root",outputBaseName,outputBaseName) << endl;
1829
1830 gSystem->Exec(Form("cat %s.txt",outputBaseName));
1831
1832}
1833
b37b9546 1834//______________________________________________________________________________
1835void AliMUONCDB::CheckHV(Int_t runNumber, Int_t verbose)
1836{
1837 /// Check the HV values in OCDB for a given run
1838
1839 TList messages;
1840 messages.SetOwner(kTRUE);
1841
1842 Bool_t patched(kTRUE);
1843
1844 if (!AliCDBManager::Instance()->IsDefaultStorageSet())
1845 {
1846 AliCDBManager::Instance()->SetDefaultStorage("raw://");
1847 }
1848
1849 AliCDBManager::Instance()->SetRun(runNumber);
1850
1851 LoadMapping();
1852
1853 AliMUONCalibrationData::CreateHV(runNumber,0,patched,&messages);
1854
1855 AliMUONCalibrationData cd(runNumber,true);
1856
1857 AliMUONPadStatusMaker statusMaker(cd);
1858
1859 AliMUONRecoParam* rp = AliMUONCDB::LoadRecoParam();
1860
1861 if (!rp)
1862 {
1863 AliErrorGeneral("AliMUONCDB::CheckHV","Could not get RecoParam !!!");
1864 return;
1865 }
1866
1867 statusMaker.SetLimits(*rp);
1868
1869 TIter next(&messages);
1870 TObjString* s;
1871 AliMpDCSNamer hvNamer("TRACKER");
1872 AliMUONLogger log;
1873
1874 while ( ( s = static_cast<TObjString*>(next()) ) )
1875 {
1876 TObjArray* a = s->String().Tokenize(":");
1877
1878 TString name(static_cast<TObjString*>(a->At(0))->String());
1879
1880 if ( name.Contains("sw") || name.Contains("SUMMARY") ) continue;
1881
1882 Int_t index = hvNamer.DCSIndexFromDCSAlias(name.Data());
1883
1884 Int_t detElemId = hvNamer.DetElemIdFromDCSAlias(name.Data());
1885
1886 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
1887
1888 if (!de)
1889 {
1890 AliErrorGeneral("AliMUONCDB::CheckHV",Form("Could not get detElemId from dcsAlias %s",name.Data()));
1891 continue;
1892 }
1893
1894 Int_t manuId;
1895
1896 if ( index >= 0 )
1897 {
1898 const AliMpArrayI* array = de->ManusForHV(index);
1899 manuId = array->GetValue(0);
1900 }
1901 else
1902
1903 {
1904 AliMpBusPatch* bp = AliMpDDLStore::Instance()->GetBusPatch(de->GetBusPatchId(0));
1905
1906 manuId = bp->GetManuId(0);
1907 }
1908
1909 Int_t status = statusMaker.HVStatus(detElemId,manuId);
1910
1911 log.Log(AliMUONPadStatusMaker::AsString(status).Data());
1912
1913 s->String() += Form(" (DE %4d) ",detElemId);
1914 s->String() += AliMUONPadStatusMaker::AsString(status).Data();
1915
1916 delete a;
1917 }
1918
1919 TIter nextMessage(&messages);
1920 TObjString* msg;
1921
1922 while ( ( msg = static_cast<TObjString*>(nextMessage()) ) )
1923 {
1924 if ( verbose > 0 || msg->String().Contains("SUMMARY") )
1925 {
1926 AliInfoGeneral("AliMUONCDB::CheckHV",Form("RUN %09d HVchannel %s",runNumber,msg->String().Data()));
1927 }
1928 }
1929
1930 TString lmsg;
1931 Int_t occurance;
1932 TString totalLog;
1933
1934 while (log.Next(lmsg,occurance))
1935 {
1936 totalLog += Form("%s(%d)",lmsg.Data(),occurance);
1937 totalLog += " | ";
1938 }
1939
1940 AliInfoGeneral("AliMUONCDB::CheckHV",Form("RUN %09d %s",runNumber,totalLog.Data()));
1941
1942 // one last loop to get the list of problematic HV channels
1943 nextMessage.Reset();
1944
1945 while ( ( msg = static_cast<TObjString*>(nextMessage()) ) )
1946 {
1947 if ( msg->String().Contains("HV ") )
1948 {
1949 AliInfoGeneral("AliMUONCDB::CheckHV",Form(" Problem at %s",msg->String().Data()));
1950 }
1951 }
1952
1953 AliCDBManager::Instance()->ClearCache();
1954}
1955
1956