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