1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //-----------------------------------------------------------------------------
19 /// \namespace AliMUONCDB
21 /// Helper functions to experience the OCDB
23 /// They allow to read magnetic field, mapping and recoParam from OCDB
25 /// And also to generate dummy (but complete) containers for all the
26 /// calibration data types we have for tracker and trigger, and to write
29 /// For more information, please see READMEcalib
31 /// \author Laurent Aphecetche
32 //-----------------------------------------------------------------------------
34 #include "AliMUONCDB.h"
36 #include "AliMUON1DArray.h"
37 #include "AliMUON1DMap.h"
38 #include "AliMUON2DMap.h"
39 #include "AliMUON2DStoreValidator.h"
40 #include "AliMUONCalibParamND.h"
41 #include "AliMUONCalibParamNF.h"
42 #include "AliMUONCalibParamNI.h"
43 #include "AliMUONCalibrationData.h"
44 #include "AliMUONConstants.h"
45 #include "AliMUONGlobalCrateConfig.h"
46 #include "AliMUONPadStatusMaker.h"
47 #include "AliMUONPadStatusMapMaker.h"
48 #include "AliMUONRecoParam.h"
49 #include "AliMUONRegionalTriggerConfig.h"
50 #include "AliMUONRejectList.h"
51 #include "AliMUONTrackerIO.h"
52 #include "AliMUONTriggerEfficiencyCells.h"
53 #include "AliMUONTriggerLut.h"
54 #include "AliMUONVCalibParam.h"
55 #include "AliMUONVCalibParam.h"
56 #include "AliMUONVStore.h"
60 #include "AliMpConstants.h"
61 #include "AliMpDEStore.h"
62 #include "AliMpDDLStore.h"
63 #include "AliMpManuStore.h"
64 #include "AliMpDEManager.h"
65 #include "AliMpDetElement.h"
66 #include "AliMpFiles.h"
67 #include "AliMpDCSNamer.h"
68 #include "AliMpManuIterator.h"
69 #include "AliMpSegmentation.h"
70 #include "AliMpStationType.h"
71 #include "AliMpVSegmentation.h"
73 #include "AliCodeTimer.h"
74 #include "AliCDBEntry.h"
75 #include "AliCDBManager.h"
76 #include "AliGRPManager.h"
77 #include "AliDCSValue.h"
79 #include "AliMpBusPatch.h"
81 #include <Riostream.h>
87 #include <TObjString.h>
90 #include <TStopwatch.h>
93 #include <TGeoGlobalMagField.h>
94 #include <TClonesArray.h>
99 //_____________________________________________________________________________
100 AliMUONVStore* Create2DMap()
102 return new AliMUON2DMap(true);
105 //_____________________________________________________________________________
106 void getBoundaries(const AliMUONVStore& store, Int_t dim,
107 Float_t* xmin, Float_t* xmax)
109 /// Assuming the store contains AliMUONVCalibParam objects, compute the
110 /// limits of the value contained in the VCalibParam, for each of its dimensions
111 /// xmin and xmax must be of dimension dim
113 for ( Int_t i = 0; i < dim; ++i )
119 TIter next(store.CreateIterator());
120 AliMUONVCalibParam* value;
122 while ( ( value = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
124 Int_t detElemId = value->ID0();
125 Int_t manuId = value->ID1();
127 const AliMpVSegmentation* seg =
128 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
130 for ( Int_t manuChannel = 0; manuChannel < value->Size(); ++manuChannel )
132 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
133 if (!pad.IsValid()) continue;
135 for ( Int_t i = 0; i < dim; ++i )
137 Float_t x0 = value->ValueAsFloat(manuChannel,i);
139 xmin[i] = TMath::Min(xmin[i],x0);
140 xmax[i] = TMath::Max(xmax[i],x0);
145 for ( Int_t i = 0; i < dim; ++i )
147 if ( TMath::Abs(xmin[i]-xmax[i]) < 1E-3 )
155 //_____________________________________________________________________________
156 Double_t GetRandom(Double_t mean, Double_t sigma, Bool_t mustBePositive)
159 if ( mustBePositive )
163 x = gRandom->Gaus(mean,sigma);
168 x = gRandom->Gaus(mean,sigma);
175 //_____________________________________________________________________________
176 Bool_t AliMUONCDB::CheckOCDB(Bool_t pathOnly)
178 /// Check that OCDB path and run number are properly set
180 AliCDBManager* man = AliCDBManager::Instance();
183 if (!man->IsDefaultStorageSet()) {
184 AliErrorGeneral("AliMUONCDB", "OCDB path must be properly set");
188 // then run number if required
189 if (pathOnly) return kTRUE;
190 if (man->GetRun() < 0) {
191 AliErrorGeneral("AliMUONCDB", "Run number must be properly set");
199 //_____________________________________________________________________________
200 Bool_t AliMUONCDB::CheckMapping(Bool_t segmentationOnly)
202 /// Check that the mapping has been loaded
204 // first the segmentation
205 if (!AliMpSegmentation::Instance(false)) {
206 AliErrorGeneral("AliMUONCDB", "Mapping segmentation must be loaded first");
210 // then the others if required
211 if (segmentationOnly) return kTRUE;
212 if (!AliMpDDLStore::Instance(false) || !AliMpDEStore::Instance(false) || !AliMpManuStore::Instance(false)) {
213 AliErrorGeneral("AliMUONCDB", "Full mapping must be loaded first");
221 //_____________________________________________________________________________
222 Bool_t AliMUONCDB::LoadField()
224 /// Load magnetic field (existing field will be deleted).
225 /// OCDB path and run number are supposed to be set.
227 AliInfoGeneral("AliMUONCDB","Loading field map from GRP...");
229 if (!AliMUONCDB::CheckOCDB()) return kFALSE;
231 AliGRPManager grpMan;
233 // make sure the old field is deleted even if it is locked
234 if(TGeoGlobalMagField::Instance()->IsLocked()) delete TGeoGlobalMagField::Instance();
236 if (!grpMan.ReadGRPEntry() || !grpMan.SetMagField()) {
237 AliErrorGeneral("AliMUONCDB", "failed to load magnetic field from OCDB");
245 //_____________________________________________________________________________
246 Bool_t AliMUONCDB::LoadMapping(Bool_t segmentationOnly)
248 /// Load mapping (existing mapping will be unloaded).
249 /// OCDB path and run number are supposed to be set.
251 AliInfoGeneral("AliMUONCDB", "Loading mapping from OCDB...");
253 if (!AliMUONCDB::CheckOCDB()) return kFALSE;
255 // in case it has already been set
256 AliMpCDB::UnloadAll();
258 if (segmentationOnly) {
260 if (!AliMpCDB::LoadMpSegmentation(kTRUE)){
261 AliErrorGeneral("AliMUONCDB", "failed to load segmentation from OCDB");
267 if (!AliMpCDB::LoadAll(kTRUE)) {
268 AliErrorGeneral("AliMUONCDB", "failed to load mapping from OCDB");
278 //_____________________________________________________________________________
279 AliMUONRecoParam* AliMUONCDB::LoadRecoParam()
281 /// Load and return reconstruction parameters.
282 /// OCDB path is supposed to be set.
284 AliInfoGeneral("AliMUONCDB", "Loading RecoParam from OCDB...");
286 if (!AliMUONCDB::CheckOCDB()) return kFALSE;
288 AliMUONRecoParam* recoParam = 0x0;
289 AliCDBEntry* entry = AliCDBManager::Instance()->Get("MUON/Calib/RecoParam");
293 // load recoParam according OCDB content (single or array)
294 if (!(recoParam = dynamic_cast<AliMUONRecoParam*>(entry->GetObject()))) {
296 TObjArray* recoParamArray = static_cast<TObjArray*>(entry->GetObject());
298 for(Int_t i = 0; i < recoParamArray->GetEntriesFast(); i++) {
299 recoParam = static_cast<AliMUONRecoParam*>(recoParamArray->UncheckedAt(i));
300 if (recoParam->IsDefault()) break;
308 if (!recoParam) AliErrorGeneral("AliMUONCDB", "failed to load RecoParam from OCDB");
314 //_____________________________________________________________________________
315 TClonesArray* AliMUONCDB::LoadAlignmentData()
317 /// Load and return the array of alignment objects.
319 AliInfoGeneral("AliMUONCDB", "Loading Alignemnt from OCDB...");
321 if (!AliMUONCDB::CheckOCDB()) return kFALSE;
323 TClonesArray* alignmentArray = 0x0;
324 AliCDBEntry* entry = AliCDBManager::Instance()->Get("MUON/Align/Data");
327 // load alignement array
328 alignmentArray = dynamic_cast<TClonesArray*>(entry->GetObject());
331 if (!alignmentArray) {
332 AliErrorGeneral("AliMUONCDB", "failed to load Alignemnt from OCDB");
335 return alignmentArray;
338 //_____________________________________________________________________________
340 AliMUONCDB::Diff(AliMUONVStore& store1, AliMUONVStore& store2,
343 /// creates a store which contains store1-store2
344 /// if opt="abs" the difference is absolute one,
345 /// if opt="rel" then what is stored is (store1-store2)/store1
346 /// if opt="percent" then what is stored is rel*100
348 /// WARNING Works only for stores which holds AliMUONVCalibParam objects
353 if ( !sopt.Contains("ABS") && !sopt.Contains("REL") && !sopt.Contains("PERCENT") )
355 AliErrorGeneral("AliMUONCDB", Form("opt %s not supported. Only ABS, REL, PERCENT are",opt));
359 AliMUONVStore* d = static_cast<AliMUONVStore*>(store1.Clone());
361 TIter next(d->CreateIterator());
363 AliMUONVCalibParam* param;
365 while ( ( param = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
367 Int_t detElemId = param->ID0();
368 Int_t manuId = param->ID1();
370 AliMUONVCalibParam* param2 = dynamic_cast<AliMUONVCalibParam*>(store2.FindObject(detElemId,manuId));
371 //FIXME: this might happen. Handle it.
374 cerr << "param2 is null : FIXME : this might happen !" << endl;
379 for ( Int_t i = 0; i < param->Size(); ++i )
381 for ( Int_t j = 0; j < param->Dimension(); ++j )
384 if ( sopt.Contains("ABS") )
386 value = param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j);
388 else if ( sopt.Contains("REL") || sopt.Contains("PERCENT") )
390 if ( param->ValueAsFloat(i,j) )
392 value = (param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j))/param->ValueAsFloat(i,j);
398 if ( sopt.Contains("PERCENT") ) value *= 100.0;
400 param->SetValueAsFloat(i,j,value);
407 //_____________________________________________________________________________
409 AliMUONCDB::Plot(const AliMUONVStore& store, const char* name, Int_t nbins)
411 /// Make histograms of each dimension of the AliMUONVCalibParam
412 /// contained inside store.
413 /// It produces histograms named name_0, name_1, etc...
415 if (!AliMUONCDB::CheckMapping(kTRUE)) return;
417 TIter next(store.CreateIterator());
418 AliMUONVCalibParam* param;
420 const Int_t kNStations = AliMpConstants::NofTrackingChambers()/2;
421 Int_t* nPerStation = new Int_t[kNStations];
424 for ( Int_t i = 0; i < kNStations; ++i ) nPerStation[i]=0;
426 while ( ( param = static_cast<AliMUONVCalibParam*>(next()) ) )
430 Int_t dim = param->Dimension();
432 Float_t* xmin = new Float_t[dim];
433 Float_t* xmax = new Float_t[dim];
434 getBoundaries(store,dim,xmin,xmax);
436 for ( Int_t i = 0; i < dim; ++i )
438 h[i] = new TH1F(Form("%s_%d",name,i),Form("%s_%d",name,i),
439 nbins,xmin[i],xmax[i]);
440 AliInfoGeneral("AliMUONCDB", Form("Created histogram %s",h[i]->GetName()));
444 Int_t detElemId = param->ID0();
445 Int_t manuId = param->ID1();
446 Int_t station = AliMpDEManager::GetChamberId(detElemId)/2;
448 const AliMpVSegmentation* seg =
449 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
451 for ( Int_t manuChannel = 0; manuChannel < param->Size(); ++manuChannel )
453 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
454 if (!pad.IsValid()) continue;
457 ++nPerStation[station];
459 for ( Int_t dim = 0; dim < param->Dimension(); ++dim )
461 h[dim]->Fill(param->ValueAsFloat(manuChannel,dim));
466 for ( Int_t i = 0; i < kNStations; ++i )
468 AliInfoGeneral("AliMUONCDB", Form("Station %d %d ",(i+1),nPerStation[i]));
471 AliInfoGeneral("AliMUONCDB", Form("Number of channels = %d",n));
473 delete[] nPerStation;
476 //_____________________________________________________________________________
478 AliMUONCDB::MakeHVStore(TMap& aliasMap, Bool_t defaultValues)
480 /// Create a HV store
482 if (!AliMUONCDB::CheckMapping()) return 0;
484 AliMpDCSNamer hvNamer("TRACKER");
486 TObjArray* aliases = hvNamer.GenerateAliases();
491 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
493 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
494 TString& aliasName = alias->String();
495 if ( aliasName.Contains("sw") )
497 // HV Switch (St345 only)
498 TObjArray* valueSet = new TObjArray;
499 valueSet->SetOwner(kTRUE);
501 Bool_t value = kTRUE;
505 Float_t r = gRandom->Uniform();
506 if ( r < 0.007 ) value = kFALSE;
509 for ( UInt_t timeStamp = 0; timeStamp < 60*3; timeStamp += 60 )
511 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
512 valueSet->Add(dcsValue);
514 aliasMap.Add(new TObjString(*alias),valueSet);
519 TObjArray* valueSet = new TObjArray;
520 valueSet->SetOwner(kTRUE);
521 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
523 Float_t value = 1500;
524 if (!defaultValues) value = GetRandom(1750,62.5,true);
525 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
526 valueSet->Add(dcsValue);
528 aliasMap.Add(new TObjString(*alias),valueSet);
535 AliInfoGeneral("AliMUONCDB", Form("%d HV channels and %d switches",nChannels,nSwitch));
537 return nChannels+nSwitch;
540 //_____________________________________________________________________________
542 AliMUONCDB::MakeTriggerDCSStore(TMap& aliasMap, Bool_t defaultValues)
544 /// Create a Trigger HV and Currents store
546 if (!AliMUONCDB::CheckMapping()) return 0;
548 AliMpDCSNamer triggerDCSNamer("TRIGGER");
550 TObjArray* aliases = triggerDCSNamer.GenerateAliases();
552 Int_t nChannels[2] = {0, 0};
554 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
556 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
557 TString& aliasName = alias->String();
559 TObjArray* valueSet = new TObjArray;
560 valueSet->SetOwner(kTRUE);
562 Int_t measureType = triggerDCSNamer.DCSvariableFromDCSAlias(aliasName.Data());
563 if ( measureType < 0 ) {
564 AliErrorGeneralStream("AliMUONCDB")
565 << "Failed to get DCS variable from an alias (trigger): "
566 << aliasName.Data() << endl;
570 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
573 (measureType == AliMpDCSNamer::kDCSI) ? 2. : 8000.;
574 if (!defaultValues) {
575 switch (measureType){
576 case AliMpDCSNamer::kDCSI:
577 value = GetRandom(2.,0.4,true);
579 case AliMpDCSNamer::kDCSHV:
580 value = GetRandom(8000.,16.,true);
584 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
585 valueSet->Add(dcsValue);
587 aliasMap.Add(new TObjString(*alias),valueSet);
588 ++nChannels[measureType];
593 AliInfoGeneral("AliMUONCDB", Form("Trigger channels I -> %i HV -> %i",nChannels[0], nChannels[1]));
595 return nChannels[0] + nChannels[1];
598 //_____________________________________________________________________________
600 AliMUONCDB::MakePedestalStore(AliMUONVStore& pedestalStore, Bool_t defaultValues)
602 /// Create a pedestal store. if defaultValues=true, ped.mean=ped.sigma=1,
603 /// otherwise mean and sigma are from a gaussian (with parameters
604 /// defined below by the kPedestal* constants)
606 AliCodeTimerAutoGeneral("",0);
608 if (!AliMUONCDB::CheckMapping()) return 0;
613 const Int_t kChannels(AliMpConstants::ManuNofChannels());
616 const Float_t kPedestalMeanMeanB(200.);
617 const Float_t kPedestalMeanSigmaB(10.);
618 const Float_t kPedestalSigmaMeanB(1.);
619 const Float_t kPedestalSigmaSigmaB(0.2);
622 const Float_t kPedestalMeanMeanNB(200.);
623 const Float_t kPedestalMeanSigmaNB(10.);
624 const Float_t kPedestalSigmaMeanNB(1.);
625 const Float_t kPedestalSigmaSigmaNB(0.2);
627 const Float_t kFractionOfDeadManu(0.); // within [0.,1.]
632 AliMpManuIterator it;
634 while ( it.Next(detElemId,manuId) )
636 // skip a given fraction of manus
637 if (kFractionOfDeadManu > 0. && gRandom->Uniform() < kFractionOfDeadManu) continue;
641 AliMUONVCalibParam* ped =
642 new AliMUONCalibParamNF(2,kChannels,detElemId,manuId,AliMUONVCalibParam::InvalidFloatValue());
644 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
646 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
648 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
652 Float_t meanPedestal;
653 Float_t sigmaPedestal;
662 Bool_t positive(kTRUE);
665 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) ) { // manu in non bending plane
667 while ( meanPedestal == 0.0 ) // avoid strict zero
669 meanPedestal = GetRandom(kPedestalMeanMeanNB,kPedestalMeanSigmaNB,positive);
671 sigmaPedestal = GetRandom(kPedestalSigmaMeanNB,kPedestalSigmaSigmaNB,positive);
673 } else { // manu in bending plane
675 while ( meanPedestal == 0.0 ) // avoid strict zero
677 meanPedestal = GetRandom(kPedestalMeanMeanB,kPedestalMeanSigmaB,positive);
679 sigmaPedestal = GetRandom(kPedestalSigmaMeanB,kPedestalSigmaSigmaB,positive);
685 ped->SetValueAsFloat(manuChannel,0,meanPedestal);
686 ped->SetValueAsFloat(manuChannel,1,sigmaPedestal);
689 Bool_t ok = pedestalStore.Add(ped);
692 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
696 AliInfoGeneral("AliMUONCDB", Form("%d Manus and %d channels.",nmanus,nchannels));
700 //_____________________________________________________________________________
702 AliMUONCDB::MakeRejectListStore(Bool_t defaultValues)
704 /// Create a reject list
706 AliCodeTimerAutoGeneral("",0);
708 AliMUONRejectList* rl = new AliMUONRejectList;
712 rl->SetDetectionElementProbability(510);
713 rl->SetDetectionElementProbability(508);
720 //_____________________________________________________________________________
722 AliMUONCDB::MakeOccupancyMapStore(AliMUONVStore& occupancyMapStore, Bool_t defaultValues)
724 /// Create an occupancy map.
726 AliCodeTimerAutoGeneral("",0);
728 if (!AliMUONCDB::CheckMapping()) return 0;
735 AliMpManuIterator it;
739 while ( it.Next(detElemId,manuId) )
743 AliMUONVCalibParam* occupancy = new AliMUONCalibParamND(5,1,detElemId,manuId,0);
747 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
749 Int_t numberOfChannelsInManu = de->NofChannelsInManu(manuId);
751 if (!defaultValues) occ = gRandom->Rndm(1);
753 Double_t sumn = occ*nevents;
755 occupancy->SetValueAsFloat(0,0,sumn);
756 occupancy->SetValueAsFloat(0,1,sumn);
757 occupancy->SetValueAsFloat(0,2,sumn);
758 occupancy->SetValueAsInt(0,3,numberOfChannelsInManu);
759 occupancy->SetValueAsInt(0,4,nevents);
761 Bool_t ok = occupancyMapStore.Add(occupancy);
764 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
771 //_____________________________________________________________________________
773 AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, const char* file)
775 /// Read the capacitance values from file and append them to the capaStore
777 if (!AliMUONCDB::CheckMapping()) return 0;
779 return AliMUONTrackerIO::ReadCapacitances(file,capaStore);
782 //_____________________________________________________________________________
784 AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, Bool_t defaultValues)
786 /// Create a capacitance store. if defaultValues=true, all capa are 1.0,
787 /// otherwise they are from a gaussian with parameters defined in the
788 /// kCapa* constants below.
790 AliCodeTimerAutoGeneral("",0);
792 if (!AliMUONCDB::CheckMapping()) return 0;
796 Int_t nmanusOK(0); // manus for which we got the serial number
798 const Float_t kCapaMean(0.3);
799 const Float_t kCapaSigma(0.1);
800 const Float_t kInjectionGainMean(3);
801 const Float_t kInjectionGainSigma(1);
806 AliMpManuIterator it;
808 while ( it.Next(detElemId,manuId) )
812 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
813 Int_t serialNumber = AliMpManuStore::Instance()->GetManuSerial(detElemId, manuId);
815 if ( serialNumber <= 0 ) continue;
819 AliMUONVCalibParam* capa = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
823 capa = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
824 Bool_t ok = capaStore.Add(capa);
827 AliErrorGeneral("AliMUONCDB", Form("Could not set serialNumber=%d manuId=%d",serialNumber,manuId));
831 for ( Int_t manuChannel = 0; manuChannel < capa->Size(); ++manuChannel )
833 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
838 Float_t injectionGain;
847 capaValue = GetRandom(kCapaMean,kCapaSigma,kTRUE);
848 injectionGain = GetRandom(kInjectionGainMean,kInjectionGainSigma,kTRUE);
850 capa->SetValueAsFloat(manuChannel,0,capaValue);
851 capa->SetValueAsFloat(manuChannel,1,injectionGain);
856 if ( nmanus ) percent = 100*nmanusOK/nmanus;
857 AliInfoGeneral("AliMUONCDB", Form("%5d manus with serial number (out of %5d manus = %3.0f%%)",
858 nmanusOK,nmanus,percent));
859 AliInfoGeneral("AliMUONCDB", Form("%5d channels",nchannels));
862 AliWarningGeneral("AliMUONCDB", "Did not get all serial numbers. capaStore is incomplete !!!!");
868 //_____________________________________________________________________________
870 AliMUONCDB::MakeGainStore(AliMUONVStore& gainStore, Bool_t defaultValues)
872 /// Create a gain store. if defaultValues=true, all gains set so that
873 /// charge = (adc-ped)
875 /// otherwise parameters are taken from gaussians with parameters
876 /// defined in the k* constants below.
878 AliCodeTimerAutoGeneral("",0);
880 if (!AliMUONCDB::CheckMapping()) return 0;
885 const Int_t kSaturation(3000);
886 const Double_t kA0Mean(1.2);
887 const Double_t kA0Sigma(0.1);
888 const Double_t kA1Mean(1E-5);
889 const Double_t kA1Sigma(1E-6);
890 const Double_t kQualMean(0xFF);
891 const Double_t kQualSigma(0x10);
892 const Int_t kThresMean(1600);
893 const Int_t kThresSigma(100);
898 AliMpManuIterator it;
900 while ( it.Next(detElemId,manuId) )
904 AliMUONVCalibParam* gain =
905 new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),
908 AliMUONVCalibParam::InvalidFloatValue());
910 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
912 for ( Int_t manuChannel = 0; manuChannel < gain->Size(); ++manuChannel )
914 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
920 gain->SetValueAsFloat(manuChannel,0,1.0);
921 gain->SetValueAsFloat(manuChannel,1,0.0);
922 gain->SetValueAsInt(manuChannel,2,4095);
923 gain->SetValueAsInt(manuChannel,3,1);
924 gain->SetValueAsInt(manuChannel,4,kSaturation);
928 Bool_t positive(kTRUE);
929 gain->SetValueAsFloat(manuChannel,0,GetRandom(kA0Mean,kA0Sigma,positive));
930 gain->SetValueAsFloat(manuChannel,1,GetRandom(kA1Mean,kA1Sigma,!positive));
931 gain->SetValueAsInt(manuChannel,2,(Int_t)TMath::Nint(GetRandom(kThresMean,kThresSigma,positive)));
932 gain->SetValueAsInt(manuChannel,3,(Int_t)TMath::Nint(GetRandom(kQualMean,kQualSigma,positive)));
933 gain->SetValueAsInt(manuChannel,4,kSaturation);
937 Bool_t ok = gainStore.Add(gain);
940 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
944 AliInfoGeneral("AliMUONCDB", Form("%d Manus and %d channels.",nmanus,nchannels));
948 //_____________________________________________________________________________
950 AliMUONCDB::MakeLocalTriggerMaskStore(AliMUONVStore& localBoardMasks)
952 /// Generate local trigger masks store. All masks are set to FFFF
954 AliCodeTimerAutoGeneral("",0);
957 // Generate fake mask values for all localboards and put that into
958 // one single container (localBoardMasks)
959 for ( Int_t i = 1; i <= AliMpConstants::TotalNofLocalBoards(); ++i )
961 AliMUONVCalibParam* localBoard = new AliMUONCalibParamNI(1,8,i,0,0);
962 for ( Int_t x = 0; x < 2; ++x )
964 for ( Int_t y = 0; y < 4; ++y )
967 localBoard->SetValueAsInt(index,0,0xFFFF);
971 localBoardMasks.Add(localBoard);
976 //_____________________________________________________________________________
978 AliMUONCDB::MakeRegionalTriggerConfigStore(AliMUONRegionalTriggerConfig& rtm)
980 /// Make a regional trigger config store. Mask is set to FFFF for each local board (Ch.F.)
982 AliCodeTimerAutoGeneral("",0);
984 if ( ! rtm.ReadData(AliMpFiles::LocalTriggerBoardMapping()) ) {
985 AliErrorGeneral("AliMUONCDB", "Error when reading from mapping file");
989 return rtm.GetNofTriggerCrates();
993 //_____________________________________________________________________________
995 AliMUONCDB::MakeGlobalTriggerConfigStore(AliMUONGlobalCrateConfig& gtm)
997 /// Make a global trigger config store. All masks (disable) set to 0x00 for each Darc board (Ch.F.)
999 AliCodeTimerAutoGeneral("",0);
1001 return gtm.ReadData(AliMpFiles::GlobalTriggerBoardMapping());
1005 //_____________________________________________________________________________
1007 AliMUONCDB::MakeTriggerLUT(const char* file)
1009 /// Make a triggerlut object, from a file.
1011 AliCodeTimerAutoGeneral("",0);
1013 AliMUONTriggerLut* lut = new AliMUONTriggerLut;
1014 lut->ReadFromFile(file);
1018 //_____________________________________________________________________________
1019 AliMUONTriggerEfficiencyCells*
1020 AliMUONCDB::MakeTriggerEfficiency(const char* file)
1022 /// Make a trigger efficiency object from a file.
1024 AliCodeTimerAutoGeneral("",0);
1026 return new AliMUONTriggerEfficiencyCells(file);
1029 //_____________________________________________________________________________
1031 AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
1032 Int_t startRun, Int_t endRun,
1033 const char* filename)
1035 /// Write a given object to OCDB
1037 TString comment(gSystem->ExpandPathName(filename));
1039 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
1042 //_____________________________________________________________________________
1044 AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
1045 Int_t startRun, Int_t endRun, Bool_t defaultValues)
1047 /// Write a given object to OCDB
1050 if ( defaultValues ) comment += "Test with default values";
1051 else comment += "Test with random values";
1053 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
1056 //_____________________________________________________________________________
1058 AliMUONCDB::WriteToCDB(TObject* object, const char* calibpath, Int_t startRun, Int_t endRun,
1059 const char* comment, const char* responsible)
1061 /// Write a given object to OCDB
1063 if (!AliMUONCDB::CheckOCDB(kTRUE)) return;
1065 AliCDBId id(calibpath,startRun,endRun);
1067 md.SetAliRootVersion(gROOT->GetVersion());
1068 md.SetComment(comment);
1069 md.SetResponsible(responsible);
1070 AliCDBManager::Instance()->Put(object,id,&md);
1073 //_____________________________________________________________________________
1075 AliMUONCDB::MakeNeighbourStore(AliMUONVStore& neighbourStore)
1077 /// Fill the neighbours store with, for each channel, a TObjArray of its
1078 /// neighbouring pads (including itself)
1080 AliCodeTimerAutoGeneral("",0);
1082 if (!AliMUONCDB::CheckMapping()) return 0;
1084 AliInfoGeneral("AliMUONCDB", "Generating NeighbourStore. This will take a while. Please be patient.");
1093 AliMpManuIterator it;
1095 while ( it.Next(detElemId,manuId) )
1097 const AliMpVSegmentation* seg =
1098 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
1100 AliMUONVCalibParam* calibParam = static_cast<AliMUONVCalibParam*>(neighbourStore.FindObject(detElemId,manuId));
1103 Int_t dimension(11);
1104 Int_t size(AliMpConstants::ManuNofChannels());
1105 Int_t defaultValue(-1);
1106 Int_t packingFactor(size);
1108 calibParam = new AliMUONCalibParamNI(dimension,size,detElemId,manuId,defaultValue,packingFactor);
1109 Bool_t ok = neighbourStore.Add(calibParam);
1112 AliErrorGeneral("AliMUONCDB", Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
1117 for ( Int_t manuChannel = 0; manuChannel < AliMpConstants::ManuNofChannels(); ++manuChannel )
1119 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
1125 seg->GetNeighbours(pad,tmp,true,true);
1126 Int_t nofPadNeighbours = tmp.GetEntriesFast();
1128 for ( Int_t i = 0; i < nofPadNeighbours; ++i )
1130 AliMpPad* p = static_cast<AliMpPad*>(tmp.UncheckedAt(i));
1133 calibParam->PackValues(p->GetManuId(),p->GetManuChannel(),x);
1136 // AliError("Could not pack value. Something is seriously wrong. Please check");
1137 // StdoutToAliError(pad->Print(););
1140 calibParam->SetValueAsInt(manuChannel,i,x);
1149 //_____________________________________________________________________________
1151 AliMUONCDB::WriteLocalTriggerMasks(Int_t startRun, Int_t endRun)
1153 /// Write local trigger masks to OCDB
1155 AliMUONVStore* ltm = new AliMUON1DArray(AliMpConstants::TotalNofLocalBoards()+1);
1156 Int_t ngenerated = MakeLocalTriggerMaskStore(*ltm);
1157 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1160 WriteToCDB("MUON/Calib/LocalTriggerBoardMasks",ltm,startRun,endRun,true);
1165 //_____________________________________________________________________________
1167 AliMUONCDB::WriteRegionalTriggerConfig(Int_t startRun, Int_t endRun)
1169 /// Write regional trigger masks to OCDB
1171 AliMUONRegionalTriggerConfig* rtm = new AliMUONRegionalTriggerConfig();
1172 Int_t ngenerated = MakeRegionalTriggerConfigStore(*rtm);
1173 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1176 WriteToCDB("MUON/Calib/RegionalTriggerConfig",rtm,startRun,endRun,true);
1182 //_____________________________________________________________________________
1184 AliMUONCDB::WriteGlobalTriggerConfig(Int_t startRun, Int_t endRun)
1186 /// Write global trigger masks to OCDB
1188 AliMUONGlobalCrateConfig* gtm = new AliMUONGlobalCrateConfig();
1190 Int_t ngenerated = MakeGlobalTriggerConfigStore(*gtm);
1191 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1194 WriteToCDB("MUON/Calib/GlobalTriggerCrateConfig",gtm,startRun,endRun,true);
1200 //_____________________________________________________________________________
1202 AliMUONCDB::WriteTriggerLut(Int_t startRun, Int_t endRun)
1204 /// Write trigger LUT to OCDB
1206 AliMUONTriggerLut* lut = MakeTriggerLUT();
1209 WriteToCDB("MUON/Calib/TriggerLut",lut,startRun,endRun,true);
1214 //_____________________________________________________________________________
1216 AliMUONCDB::WriteTriggerEfficiency(Int_t startRun, Int_t endRun)
1218 /// Write trigger efficiency to OCDB
1220 AliMUONTriggerEfficiencyCells* eff = MakeTriggerEfficiency();
1223 WriteToCDB("MUON/Calib/TriggerEfficiency",eff,startRun,endRun,true);
1228 //_____________________________________________________________________________
1230 AliMUONCDB::WriteNeighbours(Int_t startRun, Int_t endRun)
1232 /// Write neighbours to OCDB
1234 AliMUONVStore* neighbours = Create2DMap();
1235 Int_t ngenerated = MakeNeighbourStore(*neighbours);
1236 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1239 WriteToCDB("MUON/Calib/Neighbours",neighbours,startRun,endRun,true);
1244 //_____________________________________________________________________________
1246 AliMUONCDB::WriteHV(Bool_t defaultValues,
1247 Int_t startRun, Int_t endRun)
1249 /// generate HV values (either cste = 1500 V) if defaultValues=true or random
1250 /// if defaultValues=false, see makeHVStore) and
1251 /// store them into CDB located at cdbpath, with a validity period
1252 /// ranging from startRun to endRun
1254 TMap* hvStore = new TMap;
1255 Int_t ngenerated = MakeHVStore(*hvStore,defaultValues);
1256 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1259 WriteToCDB("MUON/Calib/HV",hvStore,startRun,endRun,defaultValues);
1264 //_____________________________________________________________________________
1266 AliMUONCDB::WriteTriggerDCS(Bool_t defaultValues,
1267 Int_t startRun, Int_t endRun)
1269 /// generate Trigger HV and current values (either const if defaultValues=true or random
1270 /// if defaultValues=false, see makeTriggerDCSStore) and
1271 /// store them into CDB located at cdbpath, with a validity period
1272 /// ranging from startRun to endRun
1274 TMap* triggerDCSStore = new TMap;
1275 Int_t ngenerated = MakeTriggerDCSStore(*triggerDCSStore,defaultValues);
1276 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1279 WriteToCDB("MUON/Calib/TriggerDCS",triggerDCSStore,startRun,endRun,defaultValues);
1281 delete triggerDCSStore;
1284 //_____________________________________________________________________________
1286 AliMUONCDB::WritePedestals(Bool_t defaultValues,
1287 Int_t startRun, Int_t endRun)
1289 /// generate pedestal values (either 0 if defaultValues=true or random
1290 /// if defaultValues=false, see makePedestalStore) and
1291 /// store them into CDB located at cdbpath, with a validity period
1292 /// ranging from startRun to endRun
1294 AliMUONVStore* pedestalStore = Create2DMap();
1295 Int_t ngenerated = MakePedestalStore(*pedestalStore,defaultValues);
1296 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1297 WriteToCDB("MUON/Calib/Pedestals",pedestalStore,startRun,endRun,defaultValues);
1298 delete pedestalStore;
1301 //_____________________________________________________________________________
1303 AliMUONCDB::WriteOccupancyMap(Bool_t defaultValues,
1304 Int_t startRun, Int_t endRun)
1306 /// generate occupancy map values (either empty one if defaultValues=true, or
1307 /// random one, see MakeOccupancyMapStore) and
1308 /// store them into CDB located at cdbpath, with a validity period
1309 /// ranging from startRun to endRun
1311 AliMUONVStore* occupancyMapStore = Create2DMap();
1312 Int_t ngenerated = MakeOccupancyMapStore(*occupancyMapStore,defaultValues);
1313 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1314 WriteToCDB("MUON/Calib/OccupancyMap",occupancyMapStore,startRun,endRun,defaultValues);
1315 delete occupancyMapStore;
1318 //_____________________________________________________________________________
1320 AliMUONCDB::WriteRejectList(Bool_t defaultValues,
1321 Int_t startRun, Int_t endRun)
1323 /// generate reject list values (either empty one if defaultValues=true, or
1324 /// random one, see MakeRejectListStore) and
1325 /// store them into CDB located at cdbpath, with a validity period
1326 /// ranging from startRun to endRun
1328 AliMUONRejectList* rl = MakeRejectListStore(defaultValues);
1329 WriteToCDB("MUON/Calib/RejectList",rl,startRun,endRun,defaultValues);
1334 //_____________________________________________________________________________
1336 AliMUONCDB::WriteGains(Bool_t defaultValues,
1337 Int_t startRun, Int_t endRun)
1339 /// generate gain values (either 1 if defaultValues=true or random
1340 /// if defaultValues=false, see makeGainStore) and
1341 /// store them into CDB located at cdbpath, with a validity period
1342 /// ranging from startRun to endRun
1344 AliMUONVStore* gainStore = Create2DMap();
1345 Int_t ngenerated = MakeGainStore(*gainStore,defaultValues);
1346 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1347 WriteToCDB("MUON/Calib/Gains",gainStore,startRun,endRun,defaultValues);
1351 //_____________________________________________________________________________
1353 AliMUONCDB::WriteCapacitances(const char* filename,
1354 Int_t startRun, Int_t endRun)
1356 /// read manu capacitance and injection gain values from file
1357 /// and store them into CDB located at cdbpath, with a validity period
1358 /// ranging from startRun to endRun
1360 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1361 Int_t ngenerated = MakeCapacitanceStore(*capaStore,filename);
1362 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1363 if ( ngenerated > 0 )
1365 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,filename);
1370 //_____________________________________________________________________________
1372 AliMUONCDB::WriteCapacitances(Bool_t defaultValues,
1373 Int_t startRun, Int_t endRun)
1375 /// generate manu capacitance values (either 1 if defaultValues=true or random
1376 /// if defaultValues=false, see makeCapacitanceStore) and
1377 /// store them into CDB located at cdbpath, with a validity period
1378 /// ranging from startRun to endRun
1380 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1381 Int_t ngenerated = MakeCapacitanceStore(*capaStore,defaultValues);
1382 AliInfoGeneral("AliMUONCDB", Form("Ngenerated = %d",ngenerated));
1383 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,defaultValues);
1387 //_____________________________________________________________________________
1389 AliMUONCDB::WriteTrigger(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1391 /// Writes all Trigger related calibration to CDB
1392 WriteTriggerDCS(defaultValues,startRun,endRun);
1393 WriteLocalTriggerMasks(startRun,endRun);
1394 WriteRegionalTriggerConfig(startRun,endRun);
1395 WriteGlobalTriggerConfig(startRun,endRun);
1396 WriteTriggerLut(startRun,endRun);
1397 WriteTriggerEfficiency(startRun,endRun);
1400 //_____________________________________________________________________________
1402 AliMUONCDB::WriteConfig(Int_t startRun, Int_t endRun)
1404 /// Write complete tracker configuration to OCDB
1405 ostringstream lines;
1406 TIter next(AliMpDDLStore::Instance()->CreateBusPatchIterator());
1408 while ( ( bp = static_cast<AliMpBusPatch*>(next()) ) )
1410 for (Int_t imanu = 0; imanu < bp->GetNofManus(); ++imanu)
1412 lines << bp->GetId() << " " << bp->GetManuId(imanu) << endl;
1416 AliMUON2DMap config(kTRUE);
1418 AliMUONTrackerIO::DecodeConfig(lines.str().c_str(),config);
1420 WriteToCDB("MUON/Calib/Config",&config,startRun,endRun,kTRUE);
1423 //_____________________________________________________________________________
1425 AliMUONCDB::WriteTracker(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1427 /// Writes all Tracker related calibration to CDB
1428 WriteHV(defaultValues,startRun,endRun);
1429 WritePedestals(defaultValues,startRun,endRun);
1430 WriteGains(defaultValues,startRun,endRun);
1431 WriteCapacitances(defaultValues,startRun,endRun);
1432 WriteNeighbours(startRun,endRun);
1433 WriteOccupancyMap(defaultValues,startRun,endRun);
1434 WriteRejectList(defaultValues,startRun,endRun);
1435 WriteConfig(startRun,endRun);
1438 //_____________________________________________________________________________
1440 AliMUONCDB::ShowCapacitances()
1442 /// Show briefly the number of capa values we have in the OCDB,
1443 /// and the list of manu that are actually in the config and for which
1444 /// we miss the capa (if any).
1446 if (!AliMUONCDB::CheckOCDB()) return;
1448 AliMUONCDB::LoadMapping();
1450 if (!AliMUONCDB::CheckMapping()) return;
1452 AliCDBEntry* e = AliCDBManager::Instance()->Get("MUON/Calib/Config");
1456 AliMUONVStore* config = static_cast<AliMUONVStore*>(e->GetObject());
1458 e = AliCDBManager::Instance()->Get("MUON/Calib/Capacitances");
1462 AliMUONVStore* capacitances = static_cast<AliMUONVStore*>(e->GetObject());
1464 AliInfoGeneral("ShowCapacitances",Form("%d capacitances are in OCDB",capacitances->GetSize()));
1466 TIter nextManu(config->CreateIterator());
1467 AliMUONVCalibParam* param;
1469 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1471 Int_t detElemId = param->ID0();
1472 Int_t manuId = param->ID1();
1475 = AliMpManuStore::Instance()->GetManuSerial(detElemId, manuId);
1479 AliErrorGeneral("ShowCapacitances",Form("Did not find serial for DE %04d MANUID %04d",detElemId,manuId));
1483 AliMUONVCalibParam* capa = static_cast<AliMUONVCalibParam*>(capacitances->FindObject(serialNumber));
1486 AliErrorGeneral("ShowCapacitances",Form("Did not find capacitance for DE %04d MANUID %04d SERIAL %d",detElemId,manuId,serialNumber));
1493 //_____________________________________________________________________________
1495 AliMUONCDB::ShowConfig(Bool_t withStatusMap)
1497 /// Dumps the current tracker configuration, i.e. number and identity of missing buspatches
1498 /// If statusMap is true, will also take into account the status map to report the number
1499 /// of good channels
1501 if (!AliMUONCDB::CheckOCDB()) return;
1503 AliMUONCDB::LoadMapping();
1505 if (!AliMUONCDB::CheckMapping()) return;
1507 AliCDBEntry* e = AliCDBManager::Instance()->Get("MUON/Calib/Config");
1511 AliMUONVStore* config = static_cast<AliMUONVStore*>(e->GetObject());
1513 AliMUONPadStatusMapMaker* statusMapMaker(0x0);
1514 AliMUONCalibrationData* cd(0x0);
1515 AliMUONPadStatusMaker* statusMaker(0x0);
1517 if ( withStatusMap )
1519 cd = new AliMUONCalibrationData(AliCDBManager::Instance()->GetRun());
1521 statusMaker = new AliMUONPadStatusMaker(*cd);
1523 AliMUONRecoParam* recoParam = AliMUONCDB::LoadRecoParam();
1525 statusMaker->SetLimits(*recoParam);
1529 UInt_t mask = recoParam->PadGoodnessMask();
1531 const Bool_t deferredInitialization = kFALSE;
1533 statusMapMaker = new AliMUONPadStatusMapMaker(*cd,mask,deferredInitialization);
1536 TIter nextManu(config->CreateIterator());
1537 AliMUONVCalibParam* param;
1539 AliMpExMap buspatches;
1541 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1543 Int_t detElemId = param->ID0();
1544 Int_t manuId = param->ID1();
1545 Int_t busPatchId = AliMpDDLStore::Instance()->GetBusPatchId(detElemId,manuId);
1546 if ( buspatches.GetValue(busPatchId) == 0x0 )
1548 buspatches.Add(busPatchId,new TObjString(Form("BP%04d",busPatchId)));
1552 TArrayI removed(buspatches.GetSize());
1554 TIter next(AliMpDDLStore::Instance()->CreateBusPatchIterator());
1560 // accounting of bus patches first
1562 while ( ( bp = static_cast<AliMpBusPatch*>(next())))
1564 if ( buspatches.GetValue(bp->GetId()) )
1570 removed.SetAt(bp->GetId(),nremoved++);
1574 // accounting of channels
1576 AliMpManuIterator it;
1578 Int_t totalNumberOfChannels(0);
1579 Int_t removedChannels(0);
1580 Int_t badChannels(0);
1581 Int_t badAndRemovedChannels(0);
1582 Int_t badOrRemovedChannels(0);
1584 Int_t detElemId, manuId;
1586 while ( it.Next(detElemId,manuId) )
1588 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
1589 for ( Int_t i = 0; i < AliMpConstants::ManuNofChannels(); ++i )
1591 Int_t busPatchId = AliMpDDLStore::Instance()->GetBusPatchId(detElemId,manuId);
1593 if ( de->IsConnectedChannel(manuId,i) )
1595 ++totalNumberOfChannels;
1596 Bool_t badBusPatch = ( buspatches.GetValue(busPatchId) == 0x0 );
1598 if ( withStatusMap )
1600 Bool_t badChannel = ( ( statusMapMaker->StatusMap(detElemId,manuId,i) & AliMUONPadStatusMapMaker::SelfDeadMask() ) != 0);
1601 if ( badChannel ) ++badChannels;
1602 if ( badBusPatch && badChannel ) ++badAndRemovedChannels;
1603 if ( badBusPatch || badChannel ) ++badOrRemovedChannels;
1606 if ( badBusPatch) ++removedChannels;
1612 Int_t* indices = new Int_t[nremoved];
1614 TMath::Sort(nremoved,removed.GetArray(),indices,kFALSE);
1616 for ( Int_t i = 0; i < nremoved; ++i )
1618 Int_t busPatchId = removed[indices[i]];
1619 bp = AliMpDDLStore::Instance()->GetBusPatch(busPatchId);
1626 cout << Form("Bus patches n=%3d nok=%3d nremoved=%3d",n,nok,nremoved) << endl;
1628 cout << Form("Channels n=%6d nremoved=%6d bad=%6d bad and removed=%6d bad or removed=%6d",
1629 totalNumberOfChannels,removedChannels,badChannels,badAndRemovedChannels,badOrRemovedChannels) << endl;
1631 if (totalNumberOfChannels>0)
1633 cout << Form("Percentage of readout channels %5.1f %%",removedChannels*100.0/totalNumberOfChannels) << endl;
1634 if ( withStatusMap )
1636 cout << Form("Percentage of non useable channels (bad or removed) %5.1f %%",
1637 badOrRemovedChannels*100.0/totalNumberOfChannels) << endl;
1642 delete statusMapMaker;