74 #include "AliCodeTimer.h"
75 #include "AliCDBEntry.h"
76 #include "AliCDBManager.h"
77 #include "AliGRPManager.h"
78 #include "AliDCSValue.h"
82 #include <Riostream.h>
89 #include <TObjString.h>
92 #include <TStopwatch.h>
95 #include <TGeoGlobalMagField.h>
96 #include <TClonesArray.h>
103 using std::ostringstream;
116 Float_t* xmin, Float_t* xmax)
122 for ( Int_t i = 0; i < dim; ++i )
131 while ( ( value = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
133 Int_t detElemId = value->
ID0();
134 Int_t manuId = value->ID1();
141 for ( Int_t manuChannel = 0; manuChannel < value->Size(); ++manuChannel )
146 for ( Int_t i = 0; i < dim; ++i )
148 Float_t x0 = value->ValueAsFloat(manuChannel,i);
150 xmin[i] = TMath::Min(xmin[i],x0);
151 xmax[i] = TMath::Max(xmax[i],x0);
156 for ( Int_t i = 0; i < dim; ++i )
158 if ( TMath::Abs(xmin[i]-xmax[i]) < 1E-3 )
167 Double_t GetRandom(Double_t mean, Double_t sigma, Bool_t mustBePositive)
170 if ( mustBePositive )
174 x = gRandom->Gaus(mean,sigma);
179 x = gRandom->Gaus(mean,sigma);
191 AliCDBManager* man = AliCDBManager::Instance();
194 if (!man->IsDefaultStorageSet()) {
195 AliErrorGeneral(
"AliMUONCDB",
"OCDB path must be properly set");
200 if (pathOnly)
return kTRUE;
201 if (man->GetRun() < 0) {
202 AliErrorGeneral(
"AliMUONCDB",
"Run number must be properly set");
217 AliErrorGeneral(
"AliMUONCDB",
"Mapping segmentation must be loaded first");
222 if (segmentationOnly)
return kTRUE;
224 AliErrorGeneral(
"AliMUONCDB",
"Full mapping must be loaded first");
238 AliInfoGeneral(
"AliMUONCDB",
"Loading field map from GRP...");
242 AliGRPManager grpMan;
245 if(TGeoGlobalMagField::Instance()->IsLocked())
delete TGeoGlobalMagField::Instance();
247 if (!grpMan.ReadGRPEntry() || !grpMan.SetMagField()) {
248 AliErrorGeneral(
"AliMUONCDB",
"failed to load magnetic field from OCDB");
262 AliInfoGeneral(
"AliMUONCDB",
"Loading mapping from OCDB...");
269 if (segmentationOnly) {
272 AliErrorGeneral(
"AliMUONCDB",
"failed to load segmentation from OCDB");
279 AliErrorGeneral(
"AliMUONCDB",
"failed to load mapping from OCDB");
295 AliInfoGeneral(
"AliMUONCDB",
"Loading RecoParam from OCDB...");
300 AliCDBEntry* entry = AliCDBManager::Instance()->Get(
"MUON/Calib/RecoParam");
305 if (!(recoParam = dynamic_cast<AliMUONRecoParam*>(entry->GetObject()))) {
310 for(Int_t i = 0; i < recoParamArray->GetEntriesFast(); i++) {
312 if (recoParam->IsDefault())
break;
320 if (!recoParam) AliErrorGeneral(
"AliMUONCDB",
"failed to load RecoParam from OCDB");
331 AliInfoGeneral(
"AliMUONCDB",
"Loading Alignemnt from OCDB...");
335 TClonesArray* alignmentArray = 0x0;
336 AliCDBEntry* entry = AliCDBManager::Instance()->Get(
"MUON/Align/Data");
340 alignmentArray =
dynamic_cast<TClonesArray*
>(entry->GetObject());
343 if (!alignmentArray) {
344 AliErrorGeneral(
"AliMUONCDB",
"failed to load Alignemnt from OCDB");
347 return alignmentArray;
365 if ( !sopt.Contains(
"ABS") && !sopt.Contains(
"REL") && !sopt.Contains(
"PERCENT") )
367 AliErrorGeneral(
"AliMUONCDB", Form(
"opt %s not supported. Only ABS, REL, PERCENT are",opt));
377 while ( ( param = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
379 Int_t detElemId = param->
ID0();
380 Int_t manuId = param->ID1();
386 cerr <<
"param2 is null : FIXME : this might happen !" << endl;
391 for ( Int_t i = 0; i < param->Size(); ++i )
393 for ( Int_t j = 0; j < param->Dimension(); ++j )
396 if ( sopt.Contains(
"ABS") )
398 value = param->ValueAsFloat(i,j) - param2->
ValueAsFloat(i,j);
400 else if ( sopt.Contains(
"REL") || sopt.Contains(
"PERCENT") )
402 if ( param->ValueAsFloat(i,j) )
404 value = (param->ValueAsFloat(i,j) - param2->
ValueAsFloat(i,j))/param->ValueAsFloat(i,j);
410 if ( sopt.Contains(
"PERCENT") ) value *= 100.0;
412 param->SetValueAsFloat(i,j,value);
433 Int_t* nPerStation =
new Int_t[kNStations];
436 for ( Int_t i = 0; i < kNStations; ++i ) nPerStation[i]=0;
438 while ( ( param = static_cast<AliMUONVCalibParam*>(next()) ) )
442 Int_t dim = param->Dimension();
444 Float_t* xmin =
new Float_t[dim];
445 Float_t* xmax =
new Float_t[dim];
446 getBoundaries(store,dim,xmin,xmax);
448 for ( Int_t i = 0; i < dim; ++i )
450 h[i] =
new TH1F(Form(
"%s_%d",name,i),Form(
"%s_%d",name,i),
451 nbins,xmin[i],xmax[i]);
452 AliInfoGeneral(
"AliMUONCDB", Form(
"Created histogram %s",h[i]->GetName()));
458 Int_t detElemId = param->ID0();
459 Int_t manuId = param->ID1();
467 for ( Int_t manuChannel = 0; manuChannel < param->Size(); ++manuChannel )
473 ++nPerStation[station];
475 for ( Int_t dim = 0; dim < param->Dimension(); ++dim )
477 h[dim]->Fill(param->ValueAsFloat(manuChannel,dim));
482 for ( Int_t i = 0; i < kNStations; ++i )
484 AliInfoGeneral(
"AliMUONCDB", Form(
"Station %d %d ",(i+1),nPerStation[i]));
487 AliInfoGeneral(
"AliMUONCDB", Form(
"Number of channels = %d",n));
489 delete[] nPerStation;
509 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
511 TObjString* alias =
static_cast<TObjString*
>(aliases->At(i));
512 TString& aliasName = alias->String();
513 if ( aliasName.Contains(
"sw") )
517 valueSet->SetOwner(kTRUE);
519 Bool_t value = kTRUE;
523 Float_t r = gRandom->Uniform();
524 if ( r < 0.007 ) value = kFALSE;
527 for ( UInt_t timeStamp = 0; timeStamp < 60*3; timeStamp += 60 )
529 AliDCSValue* dcsValue =
new AliDCSValue(value,timeStamp);
530 valueSet->Add(dcsValue);
532 aliasMap.Add(
new TObjString(*alias),valueSet);
538 valueSet->SetOwner(kTRUE);
539 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
541 Float_t value = 1500;
542 if (!defaultValues) value = GetRandom(1750,62.5,
true);
543 AliDCSValue* dcsValue =
new AliDCSValue(value,timeStamp);
544 valueSet->Add(dcsValue);
546 aliasMap.Add(
new TObjString(*alias),valueSet);
553 AliInfoGeneral(
"AliMUONCDB", Form(
"%d HV channels and %d switches",nChannels,nSwitch));
555 return nChannels+nSwitch;
562 TString alias = Form(
"MTR_%s_%s_RPC%i_%s",sInOut,smt,rpc,sMeasure.Data());
564 valueSet->SetOwner(kTRUE);
566 for ( UInt_t timeStamp = 0; timeStamp < 60*2; timeStamp += 60 ) {
567 AliDCSValue* dcsValue =
new AliDCSValue(value,timeStamp);
568 valueSet->Add(dcsValue);
571 aliasMap.Add(
new TObjString(alias),valueSet);
582 Int_t nChannels[2] = {0, 0};
666 TString chName[4] = {
"MT11",
"MT12",
"MT21",
"MT22"};
668 for ( Int_t ich=0; ich<4; ich++ ) {
669 for ( Int_t iside=0; iside<2; iside++ ) {
670 TString sInOut = ( iside == 0 ) ?
"INSIDE" :
"OUTSIDE";
671 for ( Int_t irpc=1; irpc<=9; irpc++ ) {
678 AliInfoGeneral(
"AliMUONCDB", Form(
"Trigger channels I -> %i HV -> %i",nChannels[0], nChannels[1]));
680 return nChannels[0] + nChannels[1];
691 AliCodeTimerAutoGeneral(
"",0);
701 const Float_t kPedestalMeanMeanB(200.);
702 const Float_t kPedestalMeanSigmaB(10.);
703 const Float_t kPedestalSigmaMeanB(1.);
704 const Float_t kPedestalSigmaSigmaB(0.2);
707 const Float_t kPedestalMeanMeanNB(200.);
708 const Float_t kPedestalMeanSigmaNB(10.);
709 const Float_t kPedestalSigmaMeanNB(1.);
710 const Float_t kPedestalSigmaSigmaNB(0.2);
712 const Float_t kFractionOfDeadManu(0.);
719 while ( it.
Next(detElemId,manuId) )
722 if (kFractionOfDeadManu > 0. && gRandom->Uniform() < kFractionOfDeadManu)
continue;
731 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
737 Float_t meanPedestal;
738 Float_t sigmaPedestal;
747 Bool_t positive(kTRUE);
752 while ( meanPedestal == 0.0 )
754 meanPedestal = GetRandom(kPedestalMeanMeanNB,kPedestalMeanSigmaNB,positive);
756 sigmaPedestal = GetRandom(kPedestalSigmaMeanNB,kPedestalSigmaSigmaNB,positive);
760 while ( meanPedestal == 0.0 )
762 meanPedestal = GetRandom(kPedestalMeanMeanB,kPedestalMeanSigmaB,positive);
764 sigmaPedestal = GetRandom(kPedestalSigmaMeanB,kPedestalSigmaSigmaB,positive);
774 Bool_t ok = pedestalStore.
Add(ped);
777 AliErrorGeneral(
"AliMUONCDB", Form(
"Could not set DetElemId=%d manuId=%d",detElemId,manuId));
781 AliInfoGeneral(
"AliMUONCDB", Form(
"%d Manus and %d channels.",nmanus,nchannels));
791 AliCodeTimerAutoGeneral(
"",0);
811 AliCodeTimerAutoGeneral(
"",0);
824 while ( it.
Next(detElemId,manuId) )
836 if (!defaultValues) occ = gRandom->Rndm(1);
838 Double_t sumn = occ*nevents;
846 Bool_t ok = occupancyMapStore.
Add(occupancy);
849 AliErrorGeneral(
"AliMUONCDB", Form(
"Could not set DetElemId=%d manuId=%d",detElemId,manuId));
875 AliCodeTimerAutoGeneral(
"",0);
883 const Float_t kCapaMean(0.3);
884 const Float_t kCapaSigma(0.1);
885 const Float_t kInjectionGainMean(3);
886 const Float_t kInjectionGainSigma(1);
893 while ( it.
Next(detElemId,manuId) )
900 if ( serialNumber <= 0 )
continue;
909 Bool_t ok = capaStore.
Add(capa);
912 AliErrorGeneral(
"AliMUONCDB", Form(
"Could not set serialNumber=%d manuId=%d",serialNumber,manuId));
916 for ( Int_t manuChannel = 0; manuChannel < capa->
Size(); ++manuChannel )
923 Float_t injectionGain;
932 capaValue = GetRandom(kCapaMean,kCapaSigma,kTRUE);
933 injectionGain = GetRandom(kInjectionGainMean,kInjectionGainSigma,kTRUE);
941 if ( nmanus ) percent = 100*nmanusOK/nmanus;
942 AliInfoGeneral(
"AliMUONCDB", Form(
"%5d manus with serial number (out of %5d manus = %3.0f%%)",
943 nmanusOK,nmanus,percent));
944 AliInfoGeneral(
"AliMUONCDB", Form(
"%5d channels",nchannels));
947 AliWarningGeneral(
"AliMUONCDB",
"Did not get all serial numbers. capaStore is incomplete !!!!");
963 AliCodeTimerAutoGeneral(
"",0);
970 const Int_t kSaturation(3000);
971 const Double_t kA0Mean(1.2);
972 const Double_t kA0Sigma(0.1);
973 const Double_t kA1Mean(1E-5);
974 const Double_t kA1Sigma(1E-6);
975 const Double_t kQualMean(0xFF);
976 const Double_t kQualSigma(0x10);
977 const Int_t kThresMean(1600);
978 const Int_t kThresSigma(100);
985 while ( it.
Next(detElemId,manuId) )
997 for ( Int_t manuChannel = 0; manuChannel < gain->
Size(); ++manuChannel )
1003 if ( defaultValues )
1013 Bool_t positive(kTRUE);
1014 gain->
SetValueAsFloat(manuChannel,0,GetRandom(kA0Mean,kA0Sigma,positive));
1015 gain->
SetValueAsFloat(manuChannel,1,GetRandom(kA1Mean,kA1Sigma,!positive));
1016 gain->
SetValueAsInt(manuChannel,2,(Int_t)TMath::Nint(GetRandom(kThresMean,kThresSigma,positive)));
1017 gain->
SetValueAsInt(manuChannel,3,(Int_t)TMath::Nint(GetRandom(kQualMean,kQualSigma,positive)));
1022 Bool_t ok = gainStore.
Add(gain);
1025 AliErrorGeneral(
"AliMUONCDB", Form(
"Could not set DetElemId=%d manuId=%d",detElemId,manuId));
1029 AliInfoGeneral(
"AliMUONCDB", Form(
"%d Manus and %d channels.",nmanus,nchannels));
1039 AliCodeTimerAutoGeneral(
"",0);
1041 Int_t ngenerated(0);
1047 for ( Int_t x = 0; x < 2; ++x )
1049 for ( Int_t y = 0; y < 4; ++y )
1051 Int_t index = x*4+y;
1056 localBoardMasks.
Add(localBoard);
1067 AliCodeTimerAutoGeneral(
"",0);
1070 AliErrorGeneral(
"AliMUONCDB",
"Error when reading from mapping file");
1084 AliCodeTimerAutoGeneral(
"",0);
1096 AliCodeTimerAutoGeneral(
"",0);
1109 AliCodeTimerAutoGeneral(
"",0);
1117 Int_t startRun, Int_t endRun,
1118 const char* filename)
1122 TString
comment(gSystem->ExpandPathName(filename));
1130 Int_t startRun, Int_t endRun, Bool_t defaultValues)
1135 if ( defaultValues ) comment +=
"Test with default values";
1136 else comment +=
"Test with random values";
1138 WriteToCDB(
object, calibpath, startRun, endRun, comment.Data());
1144 const char*
comment,
const char* responsible)
1150 AliCDBId id(calibpath,startRun,endRun);
1152 md.SetAliRootVersion(
gROOT->GetVersion());
1153 md.SetComment(comment);
1154 md.SetResponsible(responsible);
1155 AliCDBManager::Instance()->Put(
object,
id,&md);
1166 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1169 WriteToCDB(
"MUON/Calib/LocalTriggerBoardMasks",ltm,startRun,endRun,
true);
1182 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1185 WriteToCDB(
"MUON/Calib/RegionalTriggerConfig",rtm,startRun,endRun,
true);
1200 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1203 WriteToCDB(
"MUON/Calib/GlobalTriggerCrateConfig",gtm,startRun,endRun,
true);
1218 WriteToCDB(
"MUON/Calib/TriggerLut",lut,startRun,endRun,
true);
1232 WriteToCDB(
"MUON/Calib/TriggerEfficiency",eff,startRun,endRun,
true);
1245 TFile*
f = TFile::Open(gSystem->ExpandPathName(inputFile));
1247 if (!f->IsOpen())
return;
1249 TMap* hvStore =
static_cast<TMap*
>(f->Get(
"map"));
1251 WriteToCDB(
"MUON/Calib/HV",hvStore,runNumber,runNumber,kFALSE);
1258 Int_t startRun, Int_t endRun)
1265 TMap* hvStore =
new TMap;
1266 Int_t ngenerated =
MakeHVStore(*hvStore,defaultValues);
1267 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1270 WriteToCDB(
"MUON/Calib/HV",hvStore,startRun,endRun,defaultValues);
1283 TMap* triggerDCSStore =
new TMap;
1285 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1288 WriteToCDB(
"MUON/Calib/TriggerDCS",triggerDCSStore,startRun,endRun,
true);
1290 delete triggerDCSStore;
1296 Int_t startRun, Int_t endRun)
1305 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1306 WriteToCDB(
"MUON/Calib/Pedestals",pedestalStore,startRun,endRun,defaultValues);
1307 delete pedestalStore;
1313 Int_t startRun, Int_t endRun)
1322 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1323 WriteToCDB(
"MUON/Calib/OccupancyMap",occupancyMapStore,startRun,endRun,defaultValues);
1324 delete occupancyMapStore;
1330 Int_t startRun, Int_t endRun)
1338 WriteToCDB(
"MUON/Calib/RejectList",rl,startRun,endRun,defaultValues);
1346 Int_t startRun, Int_t endRun)
1355 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1356 WriteToCDB(
"MUON/Calib/Gains",gainStore,startRun,endRun,defaultValues);
1363 Int_t startRun, Int_t endRun)
1371 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1372 if ( ngenerated > 0 )
1374 WriteToCDB(
"MUON/Calib/Capacitances",capaStore,startRun,endRun,filename);
1382 Int_t startRun, Int_t endRun)
1391 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1392 WriteToCDB(
"MUON/Calib/Capacitances",capaStore,startRun,endRun,defaultValues);
1399 gSystem->Setenv(
"MINSTALL",gSystem->ExpandPathName(
"$ALICE_ROOT/MUON/mapping"));
1422 ostringstream lines;
1425 while ( ( bp = static_cast<AliMpBusPatch*>(next()) ) )
1427 for (Int_t imanu = 0; imanu < bp->
GetNofManus(); ++imanu)
1437 WriteToCDB(
"MUON/Calib/Config",&config,startRun,endRun,kTRUE);
1446 WriteHV(defaultValues,startRun,endRun);
1469 AliCDBEntry* e = AliCDBManager::Instance()->Get(
"MUON/Calib/Config");
1475 e = AliCDBManager::Instance()->Get(
"MUON/Calib/Capacitances");
1481 AliInfoGeneral(
"ShowCapacitances",Form(
"%d capacitances are in OCDB",capacitances->
GetSize()));
1486 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1488 Int_t detElemId = param->
ID0();
1489 Int_t manuId = param->ID1();
1496 AliErrorGeneral(
"ShowCapacitances",Form(
"Did not find serial for DE %04d MANUID %04d",detElemId,manuId));
1503 AliErrorGeneral(
"ShowCapacitances",Form(
"Did not find capacitance for DE %04d MANUID %04d SERIAL %d",detElemId,manuId,serialNumber));
1524 AliCDBEntry* e = AliCDBManager::Instance()->Get(
"MUON/Calib/Config");
1534 if ( withStatusMap )
1544 AliErrorGeneral(
"AliMUONCDB::ShowConfig",
"Cannot get recoParams from OCDB !");
1554 const Bool_t deferredInitialization = kFALSE;
1564 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1566 Int_t detElemId = param->ID0();
1567 Int_t manuId = param->ID1();
1569 if ( buspatches.
GetValue(busPatchId) == 0x0 )
1571 buspatches.
Add(busPatchId,
new TObjString(Form(
"BP%04d",busPatchId)));
1575 TArrayI removed(buspatches.
GetSize());
1585 while ( ( bp = static_cast<AliMpBusPatch*>(next())))
1593 removed.SetAt(bp->
GetId(),nremoved++);
1601 Int_t totalNumberOfChannels(0);
1602 Int_t removedChannels(0);
1603 Int_t badChannels(0);
1604 Int_t badAndRemovedChannels(0);
1605 Int_t badOrRemovedChannels(0);
1607 Int_t detElemId, manuId;
1609 while ( it.
Next(detElemId,manuId) )
1618 ++totalNumberOfChannels;
1619 Bool_t badBusPatch = ( buspatches.
GetValue(busPatchId) == 0x0 );
1621 if ( withStatusMap )
1624 if ( badChannel ) ++badChannels;
1625 if ( badBusPatch && badChannel ) ++badAndRemovedChannels;
1626 if ( badBusPatch || badChannel ) ++badOrRemovedChannels;
1629 if ( badBusPatch) ++removedChannels;
1635 Int_t* indices =
new Int_t[nremoved];
1637 TMath::Sort(nremoved,removed.GetArray(),indices,kFALSE);
1639 for ( Int_t i = 0; i < nremoved; ++i )
1641 Int_t busPatchId = removed[indices[i]];
1649 cout << Form(
"Bus patches n=%3d nok=%3d nremoved=%3d",n,nok,nremoved) << endl;
1651 cout << Form(
"Channels n=%6d nremoved=%6d bad=%6d bad and removed=%6d bad or removed=%6d",
1652 totalNumberOfChannels,removedChannels,badChannels,badAndRemovedChannels,badOrRemovedChannels) << endl;
1654 if (totalNumberOfChannels>0)
1656 cout << Form(
"Percentage of readout channels %5.1f %%",removedChannels*100.0/totalNumberOfChannels) << endl;
1657 if ( withStatusMap )
1659 cout << Form(
"Percentage of non useable channels (bad or removed) %5.1f %%",
1660 badOrRemovedChannels*100.0/totalNumberOfChannels) << endl;
1665 delete statusMapMaker;
1675 ifstream in(gSystem->ExpandPathName(filename));
1678 std::set<int> runset;
1682 in.getline(line,10000,
'\n');
1684 TString sline(line);
1686 if (sline.Contains(
","))
1691 while ( ( s = static_cast<TObjString*>(next()) ) )
1693 runset.insert(s->String().Atoi());
1699 runset.insert(sline.Atoi());
1707 for ( std::set<int>::const_iterator it = runset.begin(); it != runset.end(); ++it )
1709 integers.push_back((*it));
1712 std::sort(integers.begin(),integers.end());
1717 const char* outputBaseName,
1718 const char* ocdbPath)
1723 AliLog::GetRootLogger()->SetGlobalLogLevel(AliLog::kError);
1729 AliCDBManager* man = AliCDBManager::Instance();
1731 man->SetDefaultStorage(ocdbPath);
1733 Bool_t first(kTRUE);
1735 std::vector<int> runnumbers;
1741 std::ofstream outfile(Form(
"%s.txt",outputBaseName));
1743 for (
unsigned int i = 0 ; i < runnumbers.size(); ++i )
1745 int runNumber = runnumbers[i];
1747 man->SetRun(runNumber);
1755 AliCDBEntry* e = man->Get(
"MUON/Calib/OccupancyMap",runNumber);
1759 AliErrorGeneral(
"AliMUONCDB::ShowFaultyBusPatches",
1760 Form(
"Could not get OccupancyMap for run %09d",runNumber));
1770 std::set<int> buspatches;
1771 Double_t sumn = 1000.0;
1773 while ( ( bp = static_cast<AliMpBusPatch*>(nextBP()) ) )
1779 buspatches.insert(bp->
GetId());
1786 bpValues.
Add(param);
1793 for ( Int_t imanu = 0; imanu < bp->
GetNofManus(); ++imanu )
1804 Double_t sumw = sumn*(param->
ValueAsDouble(0)/sumn+1.0/runnumbers.size());
1805 Double_t sumw2 = 0.0;
1813 outfile << Form(
"RUN %09d",runNumber);
1815 for ( std::set<int>::const_iterator bit = buspatches.begin(); bit != buspatches.end(); ++bit )
1817 outfile << Form(
" %4d",*bit);
1824 if ( bpValues.
GetSize() == 0 )
1826 cout << Form(
"Great. No faulty bus patch (at the %g occupancy limit) found.",occLimit) << endl;
1827 gSystem->Exec(Form(
"rm %s.txt",outputBaseName));
1831 const char* name =
"BPfailureRate";
1836 TFile
f(Form(
"%s.root",outputBaseName),
"recreate");
1840 cout << Form(
"Results are in %s.txt and %s.root",outputBaseName,outputBaseName) << endl;
1842 gSystem->Exec(Form(
"cat %s.txt",outputBaseName));
1850 TPair* hvPair =
static_cast<TPair*
>(hvMap.FindObject(hvChannel));
1853 AliErrorGeneral(
"AliMUONCDB::MeanHVValueForDCSAlias",Form(
"Did not find expected alias (%s)",hvChannel));
1861 AliErrorGeneral(
"AliMUONCDB::MeanHVValueForDCSAlias",Form(
"Could not get values for alias %s",hvChannel));
1872 while ( ( val = static_cast<AliDCSValue*>(next()) ) )
1874 hv += val->GetFloat();
1894 messages.SetOwner(kTRUE);
1896 Bool_t patched(kTRUE);
1898 if (!AliCDBManager::Instance()->IsDefaultStorageSet())
1900 AliCDBManager::Instance()->SetDefaultStorage(
"raw://");
1903 AliCDBManager::Instance()->SetRun(runNumber);
1917 AliErrorGeneral(
"AliMUONCDB::CheckHV",
"Could not get RecoParam !!!");
1923 TIter next(&messages);
1927 Double_t meanHVValue(0.0);
1928 Double_t nofHVValues(0.0);
1930 while ( ( s = static_cast<TObjString*>(next()) ) )
1932 TObjArray* a = s->String().Tokenize(
":");
1934 TString name(static_cast<TObjString*>(a->At(0))->String());
1938 name =
static_cast<TObjString*
>(b->At(0))->String();
1943 if ( name.Contains(
"sw") || name.Contains(
"SUMMARY") ) {
continue;}
1953 AliErrorGeneral(
"AliMUONCDB::CheckHV",Form(
"Could not get detElemId from dcsAlias %s",name.Data()));
1971 Int_t status = statusMaker.
HVStatus(detElemId,manuId);
1979 s->String() += Form(
" (DE %4d) ",detElemId);
1983 TIter nextMessage(&messages);
1986 while ( ( msg = static_cast<TObjString*>(nextMessage()) ) )
1988 if ( verbose > 0 || msg->String().Contains(
"SUMMARY") )
1990 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
"RUN %09d HVchannel %s",runNumber,msg->String().Data()));
1998 while (log.
Next(lmsg,occurance))
2000 totalLog += Form(
"%s(%d)",lmsg.Data(),occurance);
2004 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
"RUN %09d %s",runNumber,totalLog.Data()));
2007 nextMessage.Reset();
2009 while ( ( msg = static_cast<TObjString*>(nextMessage()) ) )
2011 if ( msg->String().Contains(
"HV ") )
2013 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
" Problem at %s",msg->String().Data()));
2019 meanHVValue /= nofHVValues;
2020 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
"Mean HV for run %09d was %7.2f",runNumber,meanHVValue));
2023 AliCDBManager::Instance()->ClearCache();
Int_t GetNofManus() const
Convert a pad status container into a pad status map container.
virtual TIterator * CreateIterator() const =0
Return an iterator to loop over the whole store.
Implementation of AliMUONVTrackerData.
The class defines the configuration of global crate.
static Float_t InvalidFloatValue()
Return 1E38 as invalid float value.
void WritePedestals(Bool_t defaultValues, Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
void WriteHV(Bool_t defaultValues, Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
AliMUONRejectList * MakeRejectListStore(Bool_t defaultValues)
static AliMpDEStore * Instance(Bool_t warn=true)
static Bool_t WriteMpData()
void WriteTriggerEfficiency(Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
void WriteTriggerLut(Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
AliMUONTriggerLut * MakeTriggerLUT(const char *file="$(ALICE_ROOT)/MUON/data/lutAptLpt1Hpt1p7.root")
static TString AsString(Int_t status)
static AliMpSegmentation * Instance(Bool_t warn=true)
virtual Float_t ValueAsFloat(Int_t i, Int_t j=0) const =0
TMap * HV(Bool_t patched=kTRUE) const
Get the HV values. Use patched=kFALSE to get unprocessed (i.e. "raw") values as they are in the OCDB...
Implementation of AliMUONVCalibParam for tuples of floats.
virtual Int_t Size() const =0
The number of channels handled by this object.
virtual void SetValueAsFloat(Int_t i, Int_t j, Float_t value)=0
Set one value, for channel i, dimension j. Consider value is a float.
void WriteGains(Bool_t defaultValues, Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
void CheckHV(Int_t runNumber, Int_t verbose=0)
TClonesArray * LoadAlignmentData()
TFile f("CalibObjects.root")
virtual Double_t BusPatch(Int_t busPatchId, Int_t dim=0) const
Get the value for a given buspatch and given dimension.
Int_t Log(const char *message)
Int_t GetNofTriggerCrates() const
virtual Int_t GetSize() const =0
The number of objects stored.
static Bool_t WriteMpRunData()
Int_t GetDEfromBus(Int_t busPatchId) const
Store and give access to the trigger chamber efficiency.
Int_t MakeLocalTriggerMaskStore(AliMUONVStore <m)
Collection of methods usefull to DCS handling for MUON TRK and TRG.
void WriteGlobalTriggerConfig(Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
Int_t MakePedestalStore(AliMUONVStore &pedestalStore, Bool_t defaultValues)
const AliMpVSegmentation * GetMpSegmentationByElectronics(Int_t detElemId, Int_t elCardID, Bool_t warn=true) const
virtual void SetDimensionName(Int_t index, const char *value)
Set the name of a given dimension.
static Int_t DecodeConfig(const char *data, AliMUONVStore &confStore)
virtual void SetValueAsInt(Int_t i, Int_t j, Int_t value)=0
Set one value, for channel i, dimension j. Consider value is an integer.
Class with MUON reconstruction parameters.
AliMpDetElement * GetDetElement(Int_t detElemId, Bool_t warn=true) const
virtual Double_t ValueAsDouble(Int_t i, Int_t j=0) const
Int_t GetManuId(Int_t index) const
static AliMpManuStore * Instance(Bool_t warn=true)
The class defines the electronics properties of detection element.
TObjArray * GenerateAliases() const
Int_t MakeTriggerDCSStore(TMap &aliasMap)
Bool_t Next(TString &msg, Int_t &occurance)
void WriteTriggerDCS(Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
Bool_t Next(Int_t &detElemId, Int_t &manuId)
Class to loop over all manus of MUON Tracker.
static TString LocalTriggerBoardMapping()
void ShowFaultyBusPatches(const char *runlist, double occLimit=0.1, const char *outputBaseName="faulty.buspatches", const char *ocdbPath="raw://")
static Int_t ReadCapacitances(const char *filename, AliMUONVStore &capaStore)
void SetDetectionElementProbability(Int_t detElemId, Float_t proba=1.0)
A logger that keeps track of the number of times a message appeared.
void AddDCSValue(TMap &aliasMap, Int_t imeas, const char *smt, const char *sInOut, Int_t rpc, Float_t value)
Container of calibration values for a given number of channels.
static AliMpDDLStore * Instance(Bool_t warn=true)
static Int_t GetChamberId(Int_t detElemId, Bool_t warn=true)
AliMUONVStore * StatusMap() const
static Int_t ManuNofChannels()
Max number of channels per manu.
Int_t MakeGainStore(AliMUONVStore &gainStore, Bool_t defaultValues)
AliMUONTriggerEfficiencyCells * MakeTriggerEfficiency(const char *file="$ALICE_ROOT/MUON/data/efficiencyCells.dat")
void ReadFromFile(const char *filename)
Int_t GetBusPatchId(Int_t detElemId, Int_t manuId) const
Int_t GetManuSerial(Int_t detElemId, Int_t manuId) const
virtual TObject * FindObject(Int_t i, Int_t j) const
Find an object using 2 ids.
void WriteOccupancyMap(Bool_t defaultValues, Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
Implementation of AliMUONVCalibParam for tuples of ints.
void WriteConfig(Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
Int_t MakeRegionalTriggerConfigStore(AliMUONRegionalTriggerConfig &rtm)
Int_t MakeHVStore(TMap &aliasMap, Bool_t defaultValues)
AliMUONVStore * Diff(AliMUONVStore &store1, AliMUONVStore &store2, const char *opt="abs")
Int_t MakeCapacitanceStore(AliMUONVStore &capaStore, Bool_t defaultValues)
virtual Int_t ID0() const
First id of this object.
Implementation of AliMUONVCalibParam for tuples of double.
void ReadIntegers(const char *filename, std::vector< int > &integers)
void ShowConfig(Bool_t withStatusMap=kFALSE)
AliMpBusPatch * GetBusPatch(Int_t busPatchId, Bool_t warn=true) const
TObject * GetValue(Int_t keyFirst, Int_t keySecond) const
TH1 ** Plot(const AliMUONVStore &store, const char *name, Int_t nbins=512)
Int_t GetId() const
Return the unique Id.
void WriteTracker(Bool_t defaultValues=kTRUE, Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
Int_t DCSIndexFromDCSAlias(const char *dcsAlias) const
static Int_t TotalNofLocalBoards()
Return total number of trigger local boards.
Int_t MakeGlobalTriggerConfigStore(AliMUONGlobalCrateConfig >m)
Int_t MakeOccupancyMapStore(AliMUONVStore &occupancyMap, Bool_t defaultValues)
void Add(Int_t keyFirst, Int_t keySecond, TObject *object)
Helper class for sorted integer array.
Implementation of AliMUONVStore.
static Int_t SelfDeadMask()
Return status bit map to tell a pad is bad.
void WriteTrigger(Bool_t defaultValues=kTRUE, Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
Int_t NofChannelsInManu(Int_t manuId) const
void WriteMapping(Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
Bool_t CheckMapping(Bool_t segmentationOnly=kFALSE)
Basic implementation of AliMUONVStore container using AliMpExMap internally.
void WriteRejectList(Bool_t defaultValues, Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
Bool_t IsValid() const
Return validity.
virtual Bool_t Add(TObject *object)=0
Add an object to the store.
The class defines the properties of BusPatch.
Single entry point to access MUON calibration data.
virtual Bool_t Add(TObject *object)
Add an object to the store.
The abstract base class for the segmentation.
Implementation of AliMUONVStore.
Int_t GetBusPatchId(Int_t index) const
static TString GlobalTriggerBoardMapping()
Object to hold the list of elements we want to reject from the reconstruction.
void WriteLocalTriggerMasks(Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
void ReadIntegers(const char *filename, std::vector< int > &integers)
Class which encapsuate all information about a pad.
Bool_t CheckOCDB(Bool_t pathOnly=kFALSE)
Int_t HVStatus(Int_t detElemId, Int_t manuId) const
static TMap * CreateHV(Int_t runNumber, Int_t *startOfValidity=0, Bool_t patched=kTRUE, TList *messages=0x0, Bool_t dryRun=kFALSE)
Create a hv map (which must be deleted) from OCDB for the given run.
const AliMpArrayI * ManusForHV(Int_t hvIndex) const
Double_t MeanHVValueForDCSAlias(TMap &hvMap, const char *hvChannel)
UInt_t PadGoodnessMask() const
Get the goodness mask.
Int_t GetValue(Int_t index) const
Bool_t IsConnectedChannel(Int_t manuId, Int_t manuChannel) const
virtual AliMpPad PadByLocation(Int_t manuId, Int_t manuChannel, Bool_t warning=true) const =0
Find pad by location.
Int_t DetElemIdFromDCSAlias(const char *dcsAlias) const
Base class for MUON data stores.
Bool_t LoadMapping(Bool_t segmentationOnly=kFALSE)
void WriteCapacitances(Bool_t defaultValues, Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
static Bool_t LoadAll(Bool_t warn=false)
Int_t ReadData(const TString &fileName="")
virtual Int_t GetSize() const
The number of objects stored.
void WriteRegionalTriggerConfig(Int_t startRun=0, Int_t endRun=AliCDBRunRange::Infinity())
MUON trigger look up table class.
static Bool_t LoadMpSegmentation(Bool_t warn=false)
Int_t ReadData(const TString &fileName="")
static Int_t NofTrackingChambers()
Return number of tracking chambers.
The class defines the properties of regional trigger crate.
Make a 2DStore of pad statuses, using different sources of information.
void WriteToCDB(const char *calibpath, TObject *object, Int_t startRun, Int_t endRun, Bool_t defaultValues)
virtual TObject * FindObject(const char *name) const
Find an object by name.
Helper class making Root persistent TExMap.
static Int_t ManuMask(AliMp::PlaneType planeType)
AliMUONRecoParam * LoadRecoParam()
virtual void SetValueAsDouble(Int_t i, Int_t j, Double_t value)
void SetLimits(const AliMUONRecoParam &recoParams)
virtual void Print(Option_t *opt="") const