74 #include "AliCodeTimer.h"
75 #include "AliCDBEntry.h"
76 #include "AliCDBManager.h"
77 #include "AliGRPManager.h"
78 #include "AliDCSValue.h"
81 #include "AliMergeableCollection.h"
83 #include <Riostream.h>
91 #include <TObjString.h>
94 #include <TStopwatch.h>
97 #include <TGeoGlobalMagField.h>
98 #include <TClonesArray.h>
105 using std::ostringstream;
118 Float_t* xmin, Float_t* xmax)
124 for ( Int_t i = 0; i < dim; ++i )
133 while ( ( value = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
135 Int_t detElemId = value->
ID0();
136 Int_t manuId = value->ID1();
143 for ( Int_t manuChannel = 0; manuChannel < value->Size(); ++manuChannel )
148 for ( Int_t i = 0; i < dim; ++i )
150 Float_t x0 = value->ValueAsFloat(manuChannel,i);
152 xmin[i] = TMath::Min(xmin[i],x0);
153 xmax[i] = TMath::Max(xmax[i],x0);
158 for ( Int_t i = 0; i < dim; ++i )
160 if ( TMath::Abs(xmin[i]-xmax[i]) < 1E-3 )
169 Double_t GetRandom(Double_t mean, Double_t sigma, Bool_t mustBePositive)
172 if ( mustBePositive )
176 x = gRandom->Gaus(mean,sigma);
181 x = gRandom->Gaus(mean,sigma);
193 AliCDBManager* man = AliCDBManager::Instance();
196 if (!man->IsDefaultStorageSet()) {
197 AliErrorGeneral(
"AliMUONCDB",
"OCDB path must be properly set");
202 if (pathOnly)
return kTRUE;
203 if (man->GetRun() < 0) {
204 AliErrorGeneral(
"AliMUONCDB",
"Run number must be properly set");
218 messages.SetOwner(kTRUE);
220 Bool_t patched(kTRUE);
222 if (!AliCDBManager::Instance()->IsDefaultStorageSet())
224 AliCDBManager::Instance()->SetDefaultStorage(
"raw://");
227 AliCDBManager::Instance()->SetRun(runNumber);
241 AliErrorGeneral(
"AliMUONCDB::CheckHV",
"Could not get RecoParam !!!");
247 TIter next(&messages);
251 Double_t meanHVValue(0.0);
252 Double_t nofHVValues(0.0);
254 while ( ( s = static_cast<TObjString*>(next()) ) )
256 TObjArray* a = s->String().Tokenize(
":");
258 TString name(static_cast<TObjString*>(a->At(0))->String());
262 name =
static_cast<TObjString*
>(b->At(0))->String();
267 if ( name.Contains(
"sw") || name.Contains(
"SUMMARY") ) {
continue;}
277 AliErrorGeneral(
"AliMUONCDB::CheckHV",Form(
"Could not get detElemId from dcsAlias %s",name.Data()));
295 Int_t status = statusMaker.
HVStatus(detElemId,manuId);
303 s->String() += Form(
" (DE %4d) ",detElemId);
307 TIter nextMessage(&messages);
310 while ( ( msg = static_cast<TObjString*>(nextMessage()) ) )
312 if ( verbose > 0 || msg->String().Contains(
"SUMMARY") )
314 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
"RUN %09d HVchannel %s",runNumber,msg->String().Data()));
322 while (log.
Next(lmsg,occurance))
324 totalLog += Form(
"%s(%d)",lmsg.Data(),occurance);
328 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
"RUN %09d %s",runNumber,totalLog.Data()));
333 while ( ( msg = static_cast<TObjString*>(nextMessage()) ) )
335 if ( msg->String().Contains(
"HV ") )
337 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
" Problem at %s",msg->String().Data()));
343 meanHVValue /= nofHVValues;
344 AliInfoGeneral(
"AliMUONCDB::CheckHV",Form(
"Mean HV for run %09d was %7.2f",runNumber,meanHVValue));
347 AliCDBManager::Instance()->ClearCache();
357 std::vector<int> runnumbers;
361 std::vector<int> affectedRuns;
363 for (
unsigned int i = 0 ; i < runnumbers.size(); ++i )
365 int runNumber = runnumbers[i];
369 std::cout << Form(
"RUN %09d is potentially affected by bug ALIROOT-6402 : %s",runNumber,affected ?
"YES":
"NO") << std::endl;
373 affectedRuns.push_back(runNumber);
377 std::cout << Form(
"%4lu runs affected in the list of %4lu runs : ",affectedRuns.size(),runnumbers.size());
378 for (
unsigned int i = 0 ; i < affectedRuns.size(); ++i )
380 std::cout << affectedRuns[i] <<
",";
382 std::cout << std::endl;
395 AliLog::GetRootLogger()->SetGlobalLogLevel(AliLog::kError);
398 messages.SetOwner(kTRUE);
400 if (!AliCDBManager::Instance()->IsDefaultStorageSet())
402 AliCDBManager::Instance()->SetDefaultStorage(
"raw://");
405 AliCDBManager::Instance()->SetRun(runNumber);
411 PatchHV(*hvMap,&messages,kTRUE);
415 TIter next(&messages);
417 while ( ( msg = static_cast<TObjString*>(next())))
419 std::cout << Form(
"RUN %09d %s",runNumber,msg->String().Data()) << std::endl;
424 TObjString* hvChannelName;
426 Bool_t affected(kFALSE);
428 while ( ( hvChannelName = static_cast<TObjString*>(next()) ) )
430 TString name(hvChannelName->String());
434 if ( hvvalue < 1590.0 )
439 std::cout << Form(
"RUN %09d %40s HV VALUE %7.2f",runNumber,name.Data(),hvvalue) << std::endl;
445 AliCDBManager::Instance()->ClearCache();
457 AliErrorGeneral(
"AliMUONCDB",
"Mapping segmentation must be loaded first");
462 if (segmentationOnly)
return kTRUE;
464 AliErrorGeneral(
"AliMUONCDB",
"Full mapping must be loaded first");
477 if ( name.Contains(
"Chamber00Left") )
479 if (name.Contains(
"Quad1Sect0")) isit = kTRUE;
481 if (name.Contains(
"Quad1Sect1")) isit = kTRUE;
482 if (name.Contains(
"Quad1Sect2")) isit = kTRUE;
484 if (name.Contains(
"Quad2Sect2")) isit = kTRUE;
485 if (name.Contains(
"Quad2Sect1")) isit = kTRUE;
486 if (name.Contains(
"Quad2Sect0")) isit = kTRUE;
488 else if ( name.Contains(
"Chamber01Left"))
490 if (name.Contains(
"Quad2Sect2")) isit = kTRUE;
491 if (name.Contains(
"Quad2Sect0")) isit = kTRUE;
503 AliInfoGeneral(
"AliMUONCDB",
"Loading field map from GRP...");
507 AliGRPManager grpMan;
510 if(TGeoGlobalMagField::Instance()->IsLocked())
delete TGeoGlobalMagField::Instance();
512 if (!grpMan.ReadGRPEntry() || !grpMan.SetMagField()) {
513 AliErrorGeneral(
"AliMUONCDB",
"failed to load magnetic field from OCDB");
527 AliInfoGeneral(
"AliMUONCDB",
"Loading mapping from OCDB...");
534 if (segmentationOnly) {
537 AliErrorGeneral(
"AliMUONCDB",
"failed to load segmentation from OCDB");
544 AliErrorGeneral(
"AliMUONCDB",
"failed to load mapping from OCDB");
560 AliInfoGeneral(
"AliMUONCDB",
"Loading RecoParam from OCDB...");
565 AliCDBEntry* entry = AliCDBManager::Instance()->Get(
"MUON/Calib/RecoParam");
570 if (!(recoParam = dynamic_cast<AliMUONRecoParam*>(entry->GetObject()))) {
575 for(Int_t i = 0; i < recoParamArray->GetEntriesFast(); i++) {
577 if (recoParam->IsDefault())
break;
585 if (!recoParam) AliErrorGeneral(
"AliMUONCDB",
"failed to load RecoParam from OCDB");
596 AliInfoGeneral(
"AliMUONCDB",
"Loading Alignemnt from OCDB...");
600 TClonesArray* alignmentArray = 0x0;
601 AliCDBEntry* entry = AliCDBManager::Instance()->Get(
"MUON/Align/Data");
605 alignmentArray =
dynamic_cast<TClonesArray*
>(entry->GetObject());
608 if (!alignmentArray) {
609 AliErrorGeneral(
"AliMUONCDB",
"failed to load Alignemnt from OCDB");
612 return alignmentArray;
630 if ( !sopt.Contains(
"ABS") && !sopt.Contains(
"REL") && !sopt.Contains(
"PERCENT") )
632 AliErrorGeneral(
"AliMUONCDB", Form(
"opt %s not supported. Only ABS, REL, PERCENT are",opt));
642 while ( ( param = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
644 Int_t detElemId = param->
ID0();
645 Int_t manuId = param->ID1();
651 cerr <<
"param2 is null : FIXME : this might happen !" << endl;
656 for ( Int_t i = 0; i < param->Size(); ++i )
658 for ( Int_t j = 0; j < param->Dimension(); ++j )
661 if ( sopt.Contains(
"ABS") )
663 value = param->ValueAsFloat(i,j) - param2->
ValueAsFloat(i,j);
665 else if ( sopt.Contains(
"REL") || sopt.Contains(
"PERCENT") )
667 if ( param->ValueAsFloat(i,j) )
669 value = (param->ValueAsFloat(i,j) - param2->
ValueAsFloat(i,j))/param->ValueAsFloat(i,j);
675 if ( sopt.Contains(
"PERCENT") ) value *= 100.0;
677 param->SetValueAsFloat(i,j,value);
698 Int_t* nPerStation =
new Int_t[kNStations];
701 for ( Int_t i = 0; i < kNStations; ++i ) nPerStation[i]=0;
703 while ( ( param = static_cast<AliMUONVCalibParam*>(next()) ) )
707 Int_t dim = param->Dimension();
709 Float_t* xmin =
new Float_t[dim];
710 Float_t* xmax =
new Float_t[dim];
711 getBoundaries(store,dim,xmin,xmax);
713 for ( Int_t i = 0; i < dim; ++i )
715 h[i] =
new TH1F(Form(
"%s_%d",name,i),Form(
"%s_%d",name,i),
716 nbins,xmin[i],xmax[i]);
717 AliInfoGeneral(
"AliMUONCDB", Form(
"Created histogram %s",h[i]->GetName()));
723 Int_t detElemId = param->ID0();
724 Int_t manuId = param->ID1();
732 for ( Int_t manuChannel = 0; manuChannel < param->Size(); ++manuChannel )
738 ++nPerStation[station];
740 for ( Int_t dim = 0; dim < param->Dimension(); ++dim )
742 h[dim]->Fill(param->ValueAsFloat(manuChannel,dim));
747 for ( Int_t i = 0; i < kNStations; ++i )
749 AliInfoGeneral(
"AliMUONCDB", Form(
"Station %d %d ",(i+1),nPerStation[i]));
752 AliInfoGeneral(
"AliMUONCDB", Form(
"Number of channels = %d",n));
754 delete[] nPerStation;
770 double xmax = xmin + 3600;
772 int nbins = TMath::Nint((xmax-xmin)/timeResolution);
779 TF1 f1(
"f1",
"pol0",xmin,xmax);
781 while ((bp = static_cast<AliMpBusPatch*>(next())))
784 TH1* h =
new TH1F(Form(
"BP%04d",bp->
GetId()),Form(
"Number of hits in %d s bins",timeResolution),nbins,xmin,xmax);
785 f1.SetParameter(0,bp->
GetId());
786 h->FillRandom(
"f1",10000);
787 h->GetXaxis()->SetTimeDisplay(1);
788 h->GetXaxis()->SetTimeFormat(
"%d/%m/%y %H:%M");
789 h->GetXaxis()->SetTimeOffset(origin.Convert());
790 hc.Adopt(Form(
"/BUSPATCH/HITS/%ds",timeResolution),h);
795 TH1* h =
new TH1F(Form(
"Nevents%ds",timeResolution),Form(
"Number of events %d s bins",timeResolution),nbins,xmin,xmax);
797 f1.SetParameter(0,4200);
799 h->FillRandom(
"f1",10000);
801 h->GetXaxis()->SetTimeDisplay(1);
802 h->GetXaxis()->SetTimeFormat(
"%d/%m/%y %H:%M");
803 h->GetXaxis()->SetTimeOffset(origin.Convert());
807 return (total == 888);
825 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
827 TObjString* alias =
static_cast<TObjString*
>(aliases->At(i));
828 TString& aliasName = alias->String();
829 if ( aliasName.Contains(
"sw") )
833 valueSet->SetOwner(kTRUE);
835 Bool_t value = kTRUE;
839 Float_t r = gRandom->Uniform();
840 if ( r < 0.007 ) value = kFALSE;
843 for ( UInt_t timeStamp = 0; timeStamp < 60*3; timeStamp += 60 )
845 AliDCSValue* dcsValue =
new AliDCSValue(value,timeStamp);
846 valueSet->Add(dcsValue);
848 aliasMap.Add(
new TObjString(*alias),valueSet);
851 else if (aliasName.Contains(
"Mon") )
854 valueSet->SetOwner(kTRUE);
855 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
857 Float_t value = 1500;
858 if (!defaultValues) value = GetRandom(1750,62.5,
true);
859 AliDCSValue* dcsValue =
new AliDCSValue(value,timeStamp);
860 valueSet->Add(dcsValue);
862 aliasMap.Add(
new TObjString(*alias),valueSet);
869 AliInfoGeneral(
"AliMUONCDB", Form(
"%d HV channels and %d switches",nChannels,nSwitch));
871 return nChannels+nSwitch;
886 Int_t npos(0), nneg(0), ndig(0);
888 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
890 TObjString* alias =
static_cast<TObjString*
>(aliases->At(i));
891 TString& aliasName = alias->String();
892 Float_t refValue = 0;
894 if ( aliasName.Contains(
"anp") )
899 else if ( aliasName.Contains(
"dig") )
904 else if ( aliasName.Contains(
"ann") )
912 AliErrorGeneral(
"AliMUONCDB",
"Should not be here ! CHECK ME !");
917 valueSet->SetOwner(kTRUE);
919 Float_t value = refValue;
921 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
923 if (!defaultValues) value = GetRandom(refValue,0.05,
false);
924 AliDCSValue* dcsValue =
new AliDCSValue(value,refTime+timeStamp);
925 valueSet->Add(dcsValue);
927 aliasMap.Add(
new TObjString(*alias),valueSet);
930 Bool_t ok = (npos==nneg) && (npos==ndig) && (ndig==nneg);
933 AliErrorGeneral(
"AliMUONCDB",Form(
"Wrong number of LV channels : npos=%d nneg=%d ndig=%d",npos,nneg,ndig));
936 AliInfoGeneral(
"AliMUONCDB", Form(
"%d LV groups - %d aliases",npos,aliasMap.GetEntries()));
948 TString alias = Form(
"MTR_%s_%s_RPC%i_%s",sInOut,smt,rpc,sMeasure.Data());
950 valueSet->SetOwner(kTRUE);
952 for ( UInt_t timeStamp = 0; timeStamp < 60*2; timeStamp += 60 ) {
953 AliDCSValue* dcsValue =
new AliDCSValue(value,timeStamp);
954 valueSet->Add(dcsValue);
957 aliasMap.Add(
new TObjString(alias),valueSet);
968 Int_t nChannels[2] = {0, 0};
1052 TString chName[4] = {
"MT11",
"MT12",
"MT21",
"MT22"};
1054 for ( Int_t ich=0; ich<4; ich++ ) {
1055 for ( Int_t iside=0; iside<2; iside++ ) {
1056 TString sInOut = ( iside == 0 ) ?
"INSIDE" :
"OUTSIDE";
1057 for ( Int_t irpc=1; irpc<=9; irpc++ ) {
1064 AliInfoGeneral(
"AliMUONCDB", Form(
"Trigger channels I -> %i HV -> %i",nChannels[0], nChannels[1]));
1066 return nChannels[0] + nChannels[1];
1077 AliCodeTimerAutoGeneral(
"",0);
1087 const Float_t kPedestalMeanMeanB(200.);
1088 const Float_t kPedestalMeanSigmaB(10.);
1089 const Float_t kPedestalSigmaMeanB(1.);
1090 const Float_t kPedestalSigmaSigmaB(0.2);
1093 const Float_t kPedestalMeanMeanNB(200.);
1094 const Float_t kPedestalMeanSigmaNB(10.);
1095 const Float_t kPedestalSigmaMeanNB(1.);
1096 const Float_t kPedestalSigmaSigmaNB(0.2);
1098 const Float_t kFractionOfDeadManu(0.);
1105 while ( it.
Next(detElemId,manuId) )
1108 if (kFractionOfDeadManu > 0. && gRandom->Uniform() < kFractionOfDeadManu)
continue;
1117 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
1123 Float_t meanPedestal;
1124 Float_t sigmaPedestal;
1126 if ( defaultValues )
1129 sigmaPedestal = 1.0;
1133 Bool_t positive(kTRUE);
1138 while ( meanPedestal == 0.0 )
1140 meanPedestal = GetRandom(kPedestalMeanMeanNB,kPedestalMeanSigmaNB,positive);
1142 sigmaPedestal = GetRandom(kPedestalSigmaMeanNB,kPedestalSigmaSigmaNB,positive);
1146 while ( meanPedestal == 0.0 )
1148 meanPedestal = GetRandom(kPedestalMeanMeanB,kPedestalMeanSigmaB,positive);
1150 sigmaPedestal = GetRandom(kPedestalSigmaMeanB,kPedestalSigmaSigmaB,positive);
1160 Bool_t ok = pedestalStore.
Add(ped);
1163 AliErrorGeneral(
"AliMUONCDB", Form(
"Could not set DetElemId=%d manuId=%d",detElemId,manuId));
1167 AliInfoGeneral(
"AliMUONCDB", Form(
"%d Manus and %d channels.",nmanus,nchannels));
1177 AliCodeTimerAutoGeneral(
"",0);
1197 AliCodeTimerAutoGeneral(
"",0);
1208 Int_t nevents(1000);
1210 while ( it.
Next(detElemId,manuId) )
1222 if (!defaultValues) occ = gRandom->Rndm(1);
1224 Double_t sumn = occ*nevents;
1232 Bool_t ok = occupancyMapStore.
Add(occupancy);
1235 AliErrorGeneral(
"AliMUONCDB", Form(
"Could not set DetElemId=%d manuId=%d",detElemId,manuId));
1248 AliCodeTimerAutoGeneral(
"",0);
1250 Int_t ngenerated(0);
1256 for ( Int_t x = 0; x < 2; ++x )
1258 for ( Int_t y = 0; y < 4; ++y )
1260 Int_t index = x*4+y;
1265 localBoardMasks.
Add(localBoard);
1276 AliCodeTimerAutoGeneral(
"",0);
1279 AliErrorGeneral(
"AliMUONCDB",
"Error when reading from mapping file");
1293 AliCodeTimerAutoGeneral(
"",0);
1305 AliCodeTimerAutoGeneral(
"",0);
1318 AliCodeTimerAutoGeneral(
"",0);
1327 TObjString* hvChannelName;
1329 while ( ( hvChannelName = static_cast<TObjString*>(next()) ) )
1331 TString name(hvChannelName->String());
1333 if ( name.Contains(
"sw") )
continue;
1335 if ( name.Contains(
"iMon") )
continue;
1337 if ( onlySt1remapped )
1345 TPair* hvPair =
static_cast<TPair*
>(hvMap.FindObject(name.Data()));
1349 AliErrorGeneral(
"PatchHV",Form(
"Could not get values for alias %s",name.Data()));
1359 messages->Add(
new TObjString(Form(
"%s %s",hvChannelName->String().Data(),msg.Data())));
1364 AliErrorGeneral(
"PatchHV",Form(
"PatchHVValue was not successfull ! This is serious ! "
1365 "You'll have to check the logic for channel %s",hvChannelName->String().Data()));
1375 Int_t startRun, Int_t endRun,
1376 const char* filename)
1380 TString
comment(gSystem->ExpandPathName(filename));
1388 Int_t startRun, Int_t endRun, Bool_t defaultValues)
1393 if ( defaultValues ) comment +=
"Test with default values";
1394 else comment +=
"Test with random values";
1396 WriteToCDB(
object, calibpath, startRun, endRun, comment.Data());
1402 const char*
comment,
const char* responsible)
1408 AliCDBId id(calibpath,startRun,endRun);
1410 md.SetAliRootVersion(
gROOT->GetVersion());
1411 md.SetComment(comment);
1412 md.SetResponsible(responsible);
1413 AliCDBManager::Instance()->Put(
object,
id,&md);
1424 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1427 WriteToCDB(
"MUON/Calib/LocalTriggerBoardMasks",ltm,startRun,endRun,
true);
1440 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1443 WriteToCDB(
"MUON/Calib/RegionalTriggerConfig",rtm,startRun,endRun,
true);
1458 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1461 WriteToCDB(
"MUON/Calib/GlobalTriggerCrateConfig",gtm,startRun,endRun,
true);
1476 WriteToCDB(
"MUON/Calib/TriggerLut",lut,startRun,endRun,
true);
1490 WriteToCDB(
"MUON/Calib/TriggerEfficiency",eff,startRun,endRun,
true);
1503 TFile*
f =
TFile::Open(gSystem->ExpandPathName(inputFile));
1505 if (!f->IsOpen())
return;
1507 TMap* hvStore =
static_cast<TMap*
>(f->Get(
"map"));
1509 WriteToCDB(
"MUON/Calib/HV",hvStore,runNumber,runNumber,kFALSE);
1517 Int_t startRun, Int_t endRun)
1524 TMap* hvStore =
new TMap;
1525 Int_t ngenerated =
MakeHVStore(*hvStore,defaultValues);
1526 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1529 WriteToCDB(
"MUON/Calib/HV",hvStore,startRun,endRun,defaultValues);
1537 Int_t startRun, Int_t endRun,
1545 TMap* lvStore =
new TMap;
1546 Int_t ngenerated =
MakeLVStore(*lvStore,defaultValues,refTime);
1547 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1550 WriteToCDB(
"MUON/Calib/LV",lvStore,startRun,endRun,defaultValues);
1563 TMap* triggerDCSStore =
new TMap;
1565 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1568 WriteToCDB(
"MUON/Calib/TriggerDCS",triggerDCSStore,startRun,endRun,
true);
1570 delete triggerDCSStore;
1576 Int_t startRun, Int_t endRun)
1585 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1586 WriteToCDB(
"MUON/Calib/Pedestals",pedestalStore,startRun,endRun,defaultValues);
1587 delete pedestalStore;
1593 Int_t startRun, Int_t endRun)
1602 AliInfoGeneral(
"AliMUONCDB", Form(
"Ngenerated = %d",ngenerated));
1603 WriteToCDB(
"MUON/Calib/OccupancyMap",occupancyMapStore,startRun,endRun,defaultValues);
1604 delete occupancyMapStore;
1610 Int_t startRun, Int_t endRun)
1618 WriteToCDB(
"MUON/Calib/RejectList",rl,startRun,endRun,defaultValues);
1625 gSystem->Setenv(
"MINSTALL",gSystem->ExpandPathName(
"$ALICE_ROOT/MUON/mapping"));
1648 ostringstream lines;
1651 while ( ( bp = static_cast<AliMpBusPatch*>(next()) ) )
1653 for (Int_t imanu = 0; imanu < bp->
GetNofManus(); ++imanu)
1663 WriteToCDB(
"MUON/Calib/Config",&config,startRun,endRun,kTRUE);
1672 AliMergeableCollection bpevo(
"BPEVO");
1676 WriteToCDB(
"MUON/Calib/BPEVO",&bpevo,startRun,endRun,kTRUE);
1686 WriteHV(defaultValues,startRun,endRun);
1687 WriteLV(defaultValues,startRun,endRun);
1709 AliCDBEntry* e = AliCDBManager::Instance()->Get(
"MUON/Calib/Config");
1719 if ( withStatusMap )
1729 AliErrorGeneral(
"AliMUONCDB::ShowConfig",
"Cannot get recoParams from OCDB !");
1739 const Bool_t deferredInitialization = kFALSE;
1749 while ( ( param = static_cast<AliMUONVCalibParam*>(nextManu()) ) )
1751 Int_t detElemId = param->ID0();
1752 Int_t manuId = param->ID1();
1754 if ( buspatches.
GetValue(busPatchId) == 0x0 )
1756 buspatches.
Add(busPatchId,
new TObjString(Form(
"BP%04d",busPatchId)));
1760 TArrayI removed(buspatches.
GetSize());
1770 while ( ( bp = static_cast<AliMpBusPatch*>(next())))
1778 removed.SetAt(bp->
GetId(),nremoved++);
1786 Int_t totalNumberOfChannels(0);
1787 Int_t removedChannels(0);
1788 Int_t badChannels(0);
1789 Int_t badAndRemovedChannels(0);
1790 Int_t badOrRemovedChannels(0);
1792 Int_t detElemId, manuId;
1794 while ( it.
Next(detElemId,manuId) )
1803 ++totalNumberOfChannels;
1804 Bool_t badBusPatch = ( buspatches.
GetValue(busPatchId) == 0x0 );
1806 if ( withStatusMap )
1809 if ( badChannel ) ++badChannels;
1810 if ( badBusPatch && badChannel ) ++badAndRemovedChannels;
1811 if ( badBusPatch || badChannel ) ++badOrRemovedChannels;
1814 if ( badBusPatch) ++removedChannels;
1820 Int_t* indices =
new Int_t[nremoved];
1822 TMath::Sort(nremoved,removed.GetArray(),indices,kFALSE);
1824 for ( Int_t i = 0; i < nremoved; ++i )
1826 Int_t busPatchId = removed[indices[i]];
1834 cout << Form(
"Bus patches n=%3d nok=%3d nremoved=%3d",n,nok,nremoved) << endl;
1836 cout << Form(
"Channels n=%6d nremoved=%6d bad=%6d bad and removed=%6d bad or removed=%6d",
1837 totalNumberOfChannels,removedChannels,badChannels,badAndRemovedChannels,badOrRemovedChannels) << endl;
1839 if (totalNumberOfChannels>0)
1841 cout << Form(
"Percentage of readout channels %5.1f %%",removedChannels*100.0/totalNumberOfChannels) << endl;
1842 if ( withStatusMap )
1844 cout << Form(
"Percentage of non useable channels (bad or removed) %5.1f %%",
1845 badOrRemovedChannels*100.0/totalNumberOfChannels) << endl;
1850 delete statusMapMaker;
1860 ifstream in(gSystem->ExpandPathName(filename));
1863 std::set<int> runset;
1867 in.getline(line,10000,
'\n');
1869 TString sline(line);
1871 if (sline.Contains(
","))
1876 while ( ( s = static_cast<TObjString*>(next()) ) )
1878 runset.insert(s->String().Atoi());
1884 runset.insert(sline.Atoi());
1892 for ( std::set<int>::const_iterator it = runset.begin(); it != runset.end(); ++it )
1894 integers.push_back((*it));
1897 std::sort(integers.begin(),integers.end());
1902 const char* outputBaseName,
1903 const char* ocdbPath)
1908 AliLog::GetRootLogger()->SetGlobalLogLevel(AliLog::kError);
1914 AliCDBManager* man = AliCDBManager::Instance();
1916 man->SetDefaultStorage(ocdbPath);
1918 Bool_t first(kTRUE);
1920 std::vector<int> runnumbers;
1926 std::ofstream outfile(Form(
"%s.txt",outputBaseName));
1928 for (
unsigned int i = 0 ; i < runnumbers.size(); ++i )
1930 int runNumber = runnumbers[i];
1932 man->SetRun(runNumber);
1940 AliCDBEntry* e = man->Get(
"MUON/Calib/OccupancyMap",runNumber);
1944 AliErrorGeneral(
"AliMUONCDB::ShowFaultyBusPatches",
1945 Form(
"Could not get OccupancyMap for run %09d",runNumber));
1955 std::set<int> buspatches;
1956 Double_t sumn = 1000.0;
1958 while ( ( bp = static_cast<AliMpBusPatch*>(nextBP()) ) )
1964 buspatches.insert(bp->
GetId());
1971 bpValues.
Add(param);
1978 for ( Int_t imanu = 0; imanu < bp->
GetNofManus(); ++imanu )
1989 Double_t sumw = sumn*(param->
ValueAsDouble(0)/sumn+1.0/runnumbers.size());
1990 Double_t sumw2 = 0.0;
1998 outfile << Form(
"RUN %09d",runNumber);
2000 for ( std::set<int>::const_iterator bit = buspatches.begin(); bit != buspatches.end(); ++bit )
2002 outfile << Form(
" %4d",*bit);
2009 if ( bpValues.
GetSize() == 0 )
2011 cout << Form(
"Great. No faulty bus patch (at the %g occupancy limit) found.",occLimit) << endl;
2012 gSystem->Exec(Form(
"rm %s.txt",outputBaseName));
2016 const char* name =
"BPfailureRate";
2021 TFile
f(Form(
"%s.root",outputBaseName),
"recreate");
2025 cout << Form(
"Results are in %s.txt and %s.root",outputBaseName,outputBaseName) << endl;
2027 gSystem->Exec(Form(
"cat %s.txt",outputBaseName));
2035 TPair* hvPair =
static_cast<TPair*
>(hvMap.FindObject(hvChannel));
2038 AliErrorGeneral(
"AliMUONCDB::MeanHVValueForDCSAlias",Form(
"Did not find expected alias (%s)",hvChannel));
2046 AliErrorGeneral(
"AliMUONCDB::MeanHVValueForDCSAlias",Form(
"Could not get values for alias %s",hvChannel));
2057 while ( ( val = static_cast<AliDCSValue*>(next()) ) )
2059 hv += val->GetFloat();
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.
void WriteLV(Bool_t defaultValues, Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity(), time_t refTime=1449969676)
TFile * Open(const char *filename, Long64_t &nevents)
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.
Bool_t IsSt1DCSAliasRemapped(const TString &name)
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 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 MakeLVStore(TMap &aliasMap, Bool_t defaultValues, time_t refTime)
Int_t Log(const char *message)
Int_t GetNofTriggerCrates() const
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.
void CheckHV_ALIROOT_6402(const char *runlist, Bool_t verbose=kFALSE)
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.
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://")
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.
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
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")
virtual Int_t ID0() const
First id of this object.
static Bool_t PatchHVValues(TObjArray &values, TString *msg=0x0, Bool_t dryRun=kFALSE)
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 MakeBusPatchEvolution(AliMergeableCollection &bpevo, int timeResolution=60)
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.
static TObject * CreateObject(Int_t runNumber, const char *path, Int_t *startOfValidity=0x0)
Int_t GetBusPatchId(Int_t index) const
static TString GlobalTriggerBoardMapping()
TObjArray * GenerateAliases(const char *pattern="") const
Object to hold the list of elements we want to reject from the reconstruction.
void WriteBPEVO(Int_t startRun, Int_t endRun=AliCDBRunRange::Infinity())
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)
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.
void PatchHV(TMap &hvMap, TList *messages, Bool_t onlySt1remapped=kFALSE)
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