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 //-----------------------------------------------------------------------------
21 /// Helper class to experience the OCDB
22 /// It allows to generate dummy (but complete) containers for all the
23 /// calibration data types we have for tracker and trigger, and to write
26 /// For more information, please see READMECDB
28 // \author Laurent Aphecetche
29 //-----------------------------------------------------------------------------
31 #include "AliMUONCDB.h"
33 #include "AliMUON1DArray.h"
34 #include "AliMUON1DMap.h"
35 #include "AliMUON2DMap.h"
36 #include "AliMUON2DStoreValidator.h"
37 #include "AliMUONCalibParamNF.h"
38 #include "AliMUONCalibParamNI.h"
39 #include "AliMUONConstants.h"
40 #include "AliMUONTrackerIO.h"
41 #include "AliMUONTriggerEfficiencyCells.h"
42 #include "AliMUONTriggerLut.h"
43 #include "AliMUONVStore.h"
44 #include "AliMUONVCalibParam.h"
45 #include "AliMUONVCalibParam.h"
46 #include "AliMUONGlobalCrateConfig.h"
47 #include "AliMUONRegionalTriggerConfig.h"
50 #include "AliMpConstants.h"
51 #include "AliMpDDLStore.h"
52 #include "AliMpManuStore.h"
53 #include "AliMpDEManager.h"
54 #include "AliMpDetElement.h"
55 #include "AliMpFiles.h"
56 #include "AliMpDCSNamer.h"
57 #include "AliMpManuIterator.h"
58 #include "AliMpSegmentation.h"
59 #include "AliMpStationType.h"
60 #include "AliMpVSegmentation.h"
62 #include "AliCodeTimer.h"
63 #include "AliCDBEntry.h"
64 #include "AliCDBManager.h"
65 #include "AliDCSValue.h"
68 #include <Riostream.h>
74 #include <TObjString.h>
77 #include <TStopwatch.h>
88 //_____________________________________________________________________________
89 AliMUONVStore* Create2DMap()
91 return new AliMUON2DMap(true);
94 //_____________________________________________________________________________
95 void getBoundaries(const AliMUONVStore& store, Int_t dim,
96 Float_t* xmin, Float_t* xmax)
98 /// Assuming the store contains AliMUONVCalibParam objects, compute the
99 /// limits of the value contained in the VCalibParam, for each of its dimensions
100 /// xmin and xmax must be of dimension dim
102 for ( Int_t i = 0; i < dim; ++i )
108 TIter next(store.CreateIterator());
109 AliMUONVCalibParam* value;
111 while ( ( value = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
113 Int_t detElemId = value->ID0();
114 Int_t manuId = value->ID1();
116 const AliMpVSegmentation* seg =
117 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
119 for ( Int_t manuChannel = 0; manuChannel < value->Size(); ++manuChannel )
121 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
122 if (!pad.IsValid()) continue;
124 for ( Int_t i = 0; i < dim; ++i )
126 Float_t x0 = value->ValueAsFloat(manuChannel,i);
128 xmin[i] = TMath::Min(xmin[i],x0);
129 xmax[i] = TMath::Max(xmax[i],x0);
134 for ( Int_t i = 0; i < dim; ++i )
136 if ( TMath::Abs(xmin[i]-xmax[i]) < 1E-3 )
144 //_____________________________________________________________________________
145 Double_t GetRandom(Double_t mean, Double_t sigma, Bool_t mustBePositive)
148 if ( mustBePositive )
152 x = gRandom->Gaus(mean,sigma);
157 x = gRandom->Gaus(mean,sigma);
164 //_____________________________________________________________________________
165 AliMUONCDB::AliMUONCDB(const char* cdbpath)
168 fMaxNofChannelsToGenerate(-1)
172 if ( ! AliMpCDB::LoadDDLStore() ) {
173 AliFatal("Could not access mapping from OCDB !");
176 if ( ! AliMpCDB::LoadManuStore() ) {
177 AliFatal("Could not access run-dependent mapping from OCDB !");
181 //_____________________________________________________________________________
182 AliMUONCDB::~AliMUONCDB()
187 //_____________________________________________________________________________
189 AliMUONCDB::Diff(AliMUONVStore& store1, AliMUONVStore& store2,
192 /// creates a store which contains store1-store2
193 /// if opt="abs" the difference is absolute one,
194 /// if opt="rel" then what is stored is (store1-store2)/store1
195 /// if opt="percent" then what is stored is rel*100
197 /// WARNING Works only for stores which holds AliMUONVCalibParam objects
202 if ( !sopt.Contains("ABS") && !sopt.Contains("REL") && !sopt.Contains("PERCENT") )
204 AliErrorClass(Form("opt %s not supported. Only ABS, REL, PERCENT are",opt));
208 AliMUONVStore* d = static_cast<AliMUONVStore*>(store1.Clone());
210 TIter next(d->CreateIterator());
212 AliMUONVCalibParam* param;
214 while ( ( param = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
216 Int_t detElemId = param->ID0();
217 Int_t manuId = param->ID1();
219 AliMUONVCalibParam* param2 = dynamic_cast<AliMUONVCalibParam*>(store2.FindObject(detElemId,manuId));
220 //FIXME: this might happen. Handle it.
223 cerr << "param2 is null : FIXME : this might happen !" << endl;
228 for ( Int_t i = 0; i < param->Size(); ++i )
230 for ( Int_t j = 0; j < param->Dimension(); ++j )
233 if ( sopt.Contains("ABS") )
235 value = param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j);
237 else if ( sopt.Contains("REL") || sopt.Contains("PERCENT") )
239 if ( param->ValueAsFloat(i,j) )
241 value = (param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j))/param->ValueAsFloat(i,j);
247 if ( sopt.Contains("PERCENT") ) value *= 100.0;
249 param->SetValueAsFloat(i,j,value);
256 //_____________________________________________________________________________
258 AliMUONCDB::Plot(const AliMUONVStore& store, const char* name, Int_t nbins)
260 /// Make histograms of each dimension of the AliMUONVCalibParam
261 /// contained inside store.
262 /// It produces histograms named name_0, name_1, etc...
264 TIter next(store.CreateIterator());
265 AliMUONVCalibParam* param;
267 const Int_t kNStations = AliMpConstants::NofTrackingChambers()/2;
268 Int_t* nPerStation = new Int_t[kNStations];
271 for ( Int_t i = 0; i < kNStations; ++i ) nPerStation[i]=0;
273 while ( ( param = static_cast<AliMUONVCalibParam*>(next()) ) )
277 Int_t dim = param->Dimension();
279 Float_t* xmin = new Float_t[dim];
280 Float_t* xmax = new Float_t[dim];
281 getBoundaries(store,dim,xmin,xmax);
283 for ( Int_t i = 0; i < dim; ++i )
285 h[i] = new TH1F(Form("%s_%d",name,i),Form("%s_%d",name,i),
286 nbins,xmin[i],xmax[i]);
287 AliInfo(Form("Created histogram %s",h[i]->GetName()));
291 Int_t detElemId = param->ID0();
292 Int_t manuId = param->ID1();
293 Int_t station = AliMpDEManager::GetChamberId(detElemId)/2;
295 const AliMpVSegmentation* seg =
296 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
298 for ( Int_t manuChannel = 0; manuChannel < param->Size(); ++manuChannel )
300 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
301 if (!pad.IsValid()) continue;
304 ++nPerStation[station];
306 for ( Int_t dim = 0; dim < param->Dimension(); ++dim )
308 h[dim]->Fill(param->ValueAsFloat(manuChannel,dim));
313 for ( Int_t i = 0; i < kNStations; ++i )
315 AliInfo(Form("Station %d %d ",(i+1),nPerStation[i]));
318 AliInfo(Form("Number of channels = %d",n));
320 delete[] nPerStation;
323 //_____________________________________________________________________________
325 AliMUONCDB::MakeHVStore(TMap& aliasMap, Bool_t defaultValues)
327 /// Create a HV store
329 AliMpDCSNamer hvNamer("TRACKER");
331 TObjArray* aliases = hvNamer.GenerateAliases();
336 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
338 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
339 TString& aliasName = alias->String();
340 if ( aliasName.Contains("sw") )
342 // HV Switch (St345 only)
343 TObjArray* valueSet = new TObjArray;
344 valueSet->SetOwner(kTRUE);
346 Bool_t value = kTRUE;
350 Float_t r = gRandom->Uniform();
351 if ( r < 0.007 ) value = kFALSE;
354 for ( UInt_t timeStamp = 0; timeStamp < 60*3; timeStamp += 60 )
356 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
357 valueSet->Add(dcsValue);
359 aliasMap.Add(new TObjString(*alias),valueSet);
364 TObjArray* valueSet = new TObjArray;
365 valueSet->SetOwner(kTRUE);
366 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
368 Float_t value = 1500;
369 if (!defaultValues) value = GetRandom(1750,62.5,true);
370 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
371 valueSet->Add(dcsValue);
373 aliasMap.Add(new TObjString(*alias),valueSet);
380 AliInfo(Form("%d HV channels and %d switches",nChannels,nSwitch));
382 return nChannels+nSwitch;
385 //_____________________________________________________________________________
387 AliMUONCDB::MakeTriggerDCSStore(TMap& aliasMap, Bool_t defaultValues)
389 /// Create a Trigger HV and Currents store
391 AliMpDCSNamer triggerDCSNamer("TRIGGER");
393 TObjArray* aliases = triggerDCSNamer.GenerateAliases();
395 Int_t nChannels[2] = {0, 0};
397 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
399 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
400 TString& aliasName = alias->String();
402 TObjArray* valueSet = new TObjArray;
403 valueSet->SetOwner(kTRUE);
404 Int_t measureType = triggerDCSNamer.DCSvariableFromDCSAlias(aliasName.Data());
405 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
408 (measureType == AliMpDCSNamer::kDCSI) ? 2. : 8000.;
409 if (!defaultValues) {
410 switch (measureType){
411 case AliMpDCSNamer::kDCSI:
412 value = GetRandom(2.,0.4,true);
414 case AliMpDCSNamer::kDCSHV:
415 value = GetRandom(8000.,16.,true);
419 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
420 valueSet->Add(dcsValue);
422 aliasMap.Add(new TObjString(*alias),valueSet);
423 ++nChannels[measureType];
428 AliInfo(Form("Trigger channels I -> %i HV -> %i",nChannels[0], nChannels[1]));
430 return nChannels[0] + nChannels[1];
433 //_____________________________________________________________________________
435 AliMUONCDB::MakePedestalStore(AliMUONVStore& pedestalStore, Bool_t defaultValues)
437 /// Create a pedestal store. if defaultValues=true, ped.mean=ped.sigma=1,
438 /// otherwise mean and sigma are from a gaussian (with parameters
439 /// defined below by the kPedestal* constants)
441 AliCodeTimerAuto("");
446 const Int_t kChannels(AliMpConstants::ManuNofChannels());
449 const Float_t kPedestalMeanMeanB(200.);
450 const Float_t kPedestalMeanSigmaB(10.);
451 const Float_t kPedestalSigmaMeanB(1.);
452 const Float_t kPedestalSigmaSigmaB(0.2);
455 const Float_t kPedestalMeanMeanNB(200.);
456 const Float_t kPedestalMeanSigmaNB(10.);
457 const Float_t kPedestalSigmaMeanNB(1.);
458 const Float_t kPedestalSigmaSigmaNB(0.2);
460 const Float_t kFractionOfDeadManu(0.); // within [0.,1.]
465 AliMpManuIterator it;
467 while ( it.Next(detElemId,manuId) )
469 // skip a given fraction of manus
470 if (kFractionOfDeadManu > 0. && gRandom->Uniform() < kFractionOfDeadManu) continue;
474 AliMUONVCalibParam* ped =
475 new AliMUONCalibParamNF(2,kChannels,detElemId,manuId,AliMUONVCalibParam::InvalidFloatValue());
477 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
479 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
481 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
485 Float_t meanPedestal;
486 Float_t sigmaPedestal;
495 Bool_t positive(kTRUE);
498 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) ) { // manu in non bending plane
500 while ( meanPedestal == 0.0 ) // avoid strict zero
502 meanPedestal = GetRandom(kPedestalMeanMeanNB,kPedestalMeanSigmaNB,positive);
504 sigmaPedestal = GetRandom(kPedestalSigmaMeanNB,kPedestalSigmaSigmaNB,positive);
506 } else { // manu in bending plane
508 while ( meanPedestal == 0.0 ) // avoid strict zero
510 meanPedestal = GetRandom(kPedestalMeanMeanB,kPedestalMeanSigmaB,positive);
512 sigmaPedestal = GetRandom(kPedestalSigmaMeanB,kPedestalSigmaSigmaB,positive);
518 ped->SetValueAsFloat(manuChannel,0,meanPedestal);
519 ped->SetValueAsFloat(manuChannel,1,sigmaPedestal);
522 Bool_t ok = pedestalStore.Add(ped);
525 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
527 if ( fMaxNofChannelsToGenerate > 0 && nchannels >= fMaxNofChannelsToGenerate ) break;
530 AliInfo(Form("%d Manus and %d channels.",nmanus,nchannels));
534 //_____________________________________________________________________________
536 AliMUONCDB::MakeKillMapStore(AliMUONVStore& killMapStore)
538 /// Create a kill map.
540 AliCodeTimerAuto("");
545 const Int_t kChannels(AliMpConstants::ManuNofChannels());
547 const Float_t kFractionOfDeadManu(0.1); // within [0.,1.]
549 if ( kFractionOfDeadManu == 0.0 ) return 0.0;
554 AliMpManuIterator it;
556 while ( it.Next(detElemId,manuId) )
558 // skip a given fraction of manus
559 if ( gRandom->Uniform() > kFractionOfDeadManu) continue;
563 AliMUONVCalibParam* kill = new AliMUONCalibParamNI(1,kChannels,detElemId,manuId,0);
565 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
567 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
569 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
573 kill->SetValueAsInt(manuChannel,0,1);
576 Bool_t ok = killMapStore.Add(kill);
579 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
581 if ( fMaxNofChannelsToGenerate > 0 && nchannels >= fMaxNofChannelsToGenerate ) break;
584 AliInfo(Form("%d Manus and %d channels.",nmanus,nchannels));
588 //_____________________________________________________________________________
590 AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, const char* file)
592 /// Read the capacitance values from file and append them to the capaStore
594 return AliMUONTrackerIO::ReadCapacitances(file,capaStore);
597 //_____________________________________________________________________________
599 AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, Bool_t defaultValues)
601 /// Create a capacitance store. if defaultValues=true, all capa are 1.0,
602 /// otherwise they are from a gaussian with parameters defined in the
603 /// kCapa* constants below.
605 AliCodeTimerAuto("");
609 Int_t nmanusOK(0); // manus for which we got the serial number
611 const Float_t kCapaMean(0.3);
612 const Float_t kCapaSigma(0.1);
613 const Float_t kInjectionGainMean(3);
614 const Float_t kInjectionGainSigma(1);
619 AliMpManuIterator it;
621 while ( it.Next(detElemId,manuId) )
625 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
626 Int_t serialNumber = AliMpManuStore::Instance()->GetManuSerial(detElemId, manuId);
628 if ( serialNumber <= 0 ) continue;
632 AliMUONVCalibParam* capa = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
636 capa = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
637 Bool_t ok = capaStore.Add(capa);
640 AliError(Form("Could not set serialNumber=%d manuId=%d",serialNumber,manuId));
644 for ( Int_t manuChannel = 0; manuChannel < capa->Size(); ++manuChannel )
646 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
651 Float_t injectionGain;
660 capaValue = GetRandom(kCapaMean,kCapaSigma,kTRUE);
661 injectionGain = GetRandom(kInjectionGainMean,kInjectionGainSigma,kTRUE);
663 capa->SetValueAsFloat(manuChannel,0,capaValue);
664 capa->SetValueAsFloat(manuChannel,1,injectionGain);
669 if ( nmanus ) percent = 100*nmanusOK/nmanus;
670 AliInfo(Form("%5d manus with serial number (out of %5d manus = %3.0f%%)",
671 nmanusOK,nmanus,percent));
672 AliInfo(Form("%5d channels",nchannels));
675 AliWarning("Did not get all serial numbers. capaStore is incomplete !!!!");
681 //_____________________________________________________________________________
683 AliMUONCDB::MakeGainStore(AliMUONVStore& gainStore, Bool_t defaultValues)
685 /// Create a gain store. if defaultValues=true, all gains set so that
686 /// charge = (adc-ped)
688 /// otherwise parameters are taken from gaussians with parameters
689 /// defined in the k* constants below.
691 AliCodeTimerAuto("");
696 const Int_t kSaturation(3000);
697 const Double_t kA0Mean(1.2);
698 const Double_t kA0Sigma(0.1);
699 const Double_t kA1Mean(1E-5);
700 const Double_t kA1Sigma(1E-6);
701 const Double_t kQualMean(0xFF);
702 const Double_t kQualSigma(0x10);
703 const Int_t kThresMean(1600);
704 const Int_t kThresSigma(100);
709 AliMpManuIterator it;
711 while ( it.Next(detElemId,manuId) )
715 AliMUONVCalibParam* gain =
716 new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),
719 AliMUONVCalibParam::InvalidFloatValue());
721 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
723 for ( Int_t manuChannel = 0; manuChannel < gain->Size(); ++manuChannel )
725 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
731 gain->SetValueAsFloat(manuChannel,0,1.0);
732 gain->SetValueAsFloat(manuChannel,1,0.0);
733 gain->SetValueAsInt(manuChannel,2,4095);
734 gain->SetValueAsInt(manuChannel,3,1);
735 gain->SetValueAsInt(manuChannel,4,kSaturation);
739 Bool_t positive(kTRUE);
740 gain->SetValueAsFloat(manuChannel,0,GetRandom(kA0Mean,kA0Sigma,positive));
741 gain->SetValueAsFloat(manuChannel,1,GetRandom(kA1Mean,kA1Sigma,!positive));
742 gain->SetValueAsInt(manuChannel,2,(Int_t)TMath::Nint(GetRandom(kThresMean,kThresSigma,positive)));
743 gain->SetValueAsInt(manuChannel,3,(Int_t)TMath::Nint(GetRandom(kQualMean,kQualSigma,positive)));
744 gain->SetValueAsInt(manuChannel,4,kSaturation);
748 Bool_t ok = gainStore.Add(gain);
751 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
753 if ( fMaxNofChannelsToGenerate > 0 && nchannels >= fMaxNofChannelsToGenerate ) break;
756 AliInfo(Form("%d Manus and %d channels.",nmanus,nchannels));
760 //_____________________________________________________________________________
762 AliMUONCDB::MakeLocalTriggerMaskStore(AliMUONVStore& localBoardMasks) const
764 /// Generate local trigger masks store. All masks are set to FFFF
766 AliCodeTimerAuto("");
769 // Generate fake mask values for all localboards and put that into
770 // one single container (localBoardMasks)
771 for ( Int_t i = 1; i <= AliMpConstants::TotalNofLocalBoards(); ++i )
773 AliMUONVCalibParam* localBoard = new AliMUONCalibParamNI(1,8,i,0,0);
774 for ( Int_t x = 0; x < 2; ++x )
776 for ( Int_t y = 0; y < 4; ++y )
779 localBoard->SetValueAsInt(index,0,0xFFFF);
783 localBoardMasks.Add(localBoard);
788 //_____________________________________________________________________________
790 AliMUONCDB::MakeRegionalTriggerConfigStore(AliMUONRegionalTriggerConfig& rtm) const
792 /// Make a regional trigger config store. Mask is set to FFFF for each local board (Ch.F.)
794 AliCodeTimerAuto("");
796 if ( ! rtm.ReadData(AliMpFiles::LocalTriggerBoardMapping()) ) {
797 AliErrorStream() << "Error when reading from mapping file" << endl;
801 return rtm.GetNofTriggerCrates();
805 //_____________________________________________________________________________
807 AliMUONCDB::MakeGlobalTriggerConfigStore(AliMUONGlobalCrateConfig& gtm) const
809 /// Make a global trigger config store. All masks (disable) set to 0x00 for each Darc board (Ch.F.)
811 AliCodeTimerAuto("");
813 return gtm.ReadData(AliMpFiles::GlobalTriggerBoardMapping());
817 //_____________________________________________________________________________
819 AliMUONCDB::MakeTriggerLUT(const char* file) const
821 /// Make a triggerlut object, from a file.
823 AliCodeTimerAuto("");
825 AliMUONTriggerLut* lut = new AliMUONTriggerLut;
826 lut->ReadFromFile(file);
830 //_____________________________________________________________________________
831 AliMUONTriggerEfficiencyCells*
832 AliMUONCDB::MakeTriggerEfficiency(const char* file) const
834 /// Make a trigger efficiency object from a file.
836 AliCodeTimerAuto("");
838 return new AliMUONTriggerEfficiencyCells(file);
841 //_____________________________________________________________________________
843 AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
844 Int_t startRun, Int_t endRun,
845 const char* filename)
847 /// Write a given object to OCDB
849 TString comment(gSystem->ExpandPathName(filename));
851 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
854 //_____________________________________________________________________________
856 AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
857 Int_t startRun, Int_t endRun, Bool_t defaultValues)
859 /// Write a given object to OCDB
862 if ( defaultValues ) comment += "Test with default values";
863 else comment += "Test with random values";
865 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
868 //_____________________________________________________________________________
870 AliMUONCDB::WriteToCDB(TObject* object, const char* calibpath, Int_t startRun, Int_t endRun,
871 const char* comment, const char* responsible)
873 /// Write a given object to OCDB
875 AliCDBId id(calibpath,startRun,endRun);
877 md.SetAliRootVersion(gROOT->GetVersion());
878 md.SetComment(comment);
879 md.SetResponsible(responsible);
880 AliCDBManager* man = AliCDBManager::Instance();
881 if (!man->IsDefaultStorageSet()) man->SetDefaultStorage(fCDBPath);
882 man->Put(object,id,&md);
885 //_____________________________________________________________________________
887 AliMUONCDB::MakeNeighbourStore(AliMUONVStore& neighbourStore)
889 /// Fill the neighbours store with, for each channel, a TObjArray of its
890 /// neighbouring pads (including itself)
892 AliCodeTimerAuto("");
894 AliInfo("Generating NeighbourStore. This will take a while. Please be patient.");
903 AliMpManuIterator it;
905 while ( it.Next(detElemId,manuId) )
907 const AliMpVSegmentation* seg =
908 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
910 AliMUONVCalibParam* calibParam = static_cast<AliMUONVCalibParam*>(neighbourStore.FindObject(detElemId,manuId));
914 Int_t size(AliMpConstants::ManuNofChannels());
915 Int_t defaultValue(-1);
916 Int_t packingFactor(size);
918 calibParam = new AliMUONCalibParamNI(dimension,size,detElemId,manuId,defaultValue,packingFactor);
919 Bool_t ok = neighbourStore.Add(calibParam);
922 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
927 for ( Int_t manuChannel = 0; manuChannel < AliMpConstants::ManuNofChannels(); ++manuChannel )
929 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
935 seg->GetNeighbours(pad,tmp,true,true);
936 Int_t nofPadNeighbours = tmp.GetEntriesFast();
938 for ( Int_t i = 0; i < nofPadNeighbours; ++i )
940 AliMpPad* p = static_cast<AliMpPad*>(tmp.UncheckedAt(i));
943 calibParam->PackValues(p->GetManuId(),p->GetManuChannel(),x);
946 // AliError("Could not pack value. Something is seriously wrong. Please check");
947 // StdoutToAliError(pad->Print(););
950 calibParam->SetValueAsInt(manuChannel,i,x);
959 //_____________________________________________________________________________
961 AliMUONCDB::SetMaxNofChannelsToGenerate(Int_t n)
963 /// Set the maximum number of channels to generate (used for testing only)
964 /// n < 0 means no limit
965 fMaxNofChannelsToGenerate = n;
968 //_____________________________________________________________________________
970 AliMUONCDB::WriteLocalTriggerMasks(Int_t startRun, Int_t endRun)
972 /// Write local trigger masks to OCDB
974 AliMUONVStore* ltm = new AliMUON1DArray(AliMpConstants::TotalNofLocalBoards()+1);
975 Int_t ngenerated = MakeLocalTriggerMaskStore(*ltm);
976 AliInfo(Form("Ngenerated = %d",ngenerated));
979 WriteToCDB("MUON/Calib/LocalTriggerBoardMasks",ltm,startRun,endRun,true);
984 //_____________________________________________________________________________
986 AliMUONCDB::WriteRegionalTriggerConfig(Int_t startRun, Int_t endRun)
988 /// Write regional trigger masks to OCDB
990 AliMUONRegionalTriggerConfig* rtm = new AliMUONRegionalTriggerConfig();
991 Int_t ngenerated = MakeRegionalTriggerConfigStore(*rtm);
992 AliInfo(Form("Ngenerated = %d",ngenerated));
995 WriteToCDB("MUON/Calib/RegionalTriggerConfig",rtm,startRun,endRun,true);
1001 //_____________________________________________________________________________
1003 AliMUONCDB::WriteGlobalTriggerConfig(Int_t startRun, Int_t endRun)
1005 /// Write global trigger masks to OCDB
1007 AliMUONGlobalCrateConfig* gtm = new AliMUONGlobalCrateConfig();
1009 Int_t ngenerated = MakeGlobalTriggerConfigStore(*gtm);
1010 AliInfo(Form("Ngenerated = %d",ngenerated));
1013 WriteToCDB("MUON/Calib/GlobalTriggerCrateConfig",gtm,startRun,endRun,true);
1019 //_____________________________________________________________________________
1021 AliMUONCDB::WriteTriggerLut(Int_t startRun, Int_t endRun)
1023 /// Write trigger LUT to OCDB
1025 AliMUONTriggerLut* lut = MakeTriggerLUT();
1028 WriteToCDB("MUON/Calib/TriggerLut",lut,startRun,endRun,true);
1033 //_____________________________________________________________________________
1035 AliMUONCDB::WriteTriggerEfficiency(Int_t startRun, Int_t endRun)
1037 /// Write trigger efficiency to OCDB
1039 AliMUONTriggerEfficiencyCells* eff = MakeTriggerEfficiency();
1042 WriteToCDB("MUON/Calib/TriggerEfficiency",eff,startRun,endRun,true);
1047 //_____________________________________________________________________________
1049 AliMUONCDB::WriteNeighbours(Int_t startRun, Int_t endRun)
1051 /// Write neighbours to OCDB
1053 AliMUONVStore* neighbours = Create2DMap();
1054 Int_t ngenerated = MakeNeighbourStore(*neighbours);
1055 AliInfo(Form("Ngenerated = %d",ngenerated));
1058 WriteToCDB("MUON/Calib/Neighbours",neighbours,startRun,endRun,true);
1063 //_____________________________________________________________________________
1065 AliMUONCDB::WriteHV(Bool_t defaultValues,
1066 Int_t startRun, Int_t endRun)
1068 /// generate HV values (either cste = 1500 V) if defaultValues=true or random
1069 /// if defaultValues=false, see makeHVStore) and
1070 /// store them into CDB located at cdbpath, with a validity period
1071 /// ranging from startRun to endRun
1073 TMap* hvStore = new TMap;
1074 Int_t ngenerated = MakeHVStore(*hvStore,defaultValues);
1075 AliInfo(Form("Ngenerated = %d",ngenerated));
1078 WriteToCDB("MUON/Calib/HV",hvStore,startRun,endRun,defaultValues);
1083 //_____________________________________________________________________________
1085 AliMUONCDB::WriteTriggerDCS(Bool_t defaultValues,
1086 Int_t startRun, Int_t endRun)
1088 /// generate Trigger HV and current values (either const if defaultValues=true or random
1089 /// if defaultValues=false, see makeTriggerDCSStore) and
1090 /// store them into CDB located at cdbpath, with a validity period
1091 /// ranging from startRun to endRun
1093 TMap* triggerDCSStore = new TMap;
1094 Int_t ngenerated = MakeTriggerDCSStore(*triggerDCSStore,defaultValues);
1095 AliInfo(Form("Ngenerated = %d",ngenerated));
1098 WriteToCDB("MUON/Calib/TriggerDCS",triggerDCSStore,startRun,endRun,defaultValues);
1100 delete triggerDCSStore;
1103 //_____________________________________________________________________________
1105 AliMUONCDB::WritePedestals(Bool_t defaultValues,
1106 Int_t startRun, Int_t endRun)
1108 /// generate pedestal values (either 0 if defaultValues=true or random
1109 /// if defaultValues=false, see makePedestalStore) and
1110 /// store them into CDB located at cdbpath, with a validity period
1111 /// ranging from startRun to endRun
1113 AliMUONVStore* pedestalStore = Create2DMap();
1114 Int_t ngenerated = MakePedestalStore(*pedestalStore,defaultValues);
1115 AliInfo(Form("Ngenerated = %d",ngenerated));
1116 WriteToCDB("MUON/Calib/Pedestals",pedestalStore,startRun,endRun,defaultValues);
1117 delete pedestalStore;
1120 //_____________________________________________________________________________
1122 AliMUONCDB::WriteKillMap(Bool_t defaultValues,
1123 Int_t startRun, Int_t endRun)
1125 /// generate kill map values (either empty one if defaultValues=true, or
1126 /// random one, see MakeKillMapStore) and
1127 /// store them into CDB located at cdbpath, with a validity period
1128 /// ranging from startRun to endRun
1130 AliMUONVStore* killMapStore = Create2DMap();
1131 if ( !defaultValues )
1133 Int_t ngenerated = MakeKillMapStore(*killMapStore);
1134 AliInfo(Form("Ngenerated = %d",ngenerated));
1136 WriteToCDB("MUON/Calib/KillMap",killMapStore,startRun,endRun,defaultValues);
1137 delete killMapStore;
1141 //_____________________________________________________________________________
1143 AliMUONCDB::WriteGains(Bool_t defaultValues,
1144 Int_t startRun, Int_t endRun)
1146 /// generate gain values (either 1 if defaultValues=true or random
1147 /// if defaultValues=false, see makeGainStore) and
1148 /// store them into CDB located at cdbpath, with a validity period
1149 /// ranging from startRun to endRun
1151 AliMUONVStore* gainStore = Create2DMap();
1152 Int_t ngenerated = MakeGainStore(*gainStore,defaultValues);
1153 AliInfo(Form("Ngenerated = %d",ngenerated));
1154 WriteToCDB("MUON/Calib/Gains",gainStore,startRun,endRun,defaultValues);
1158 //_____________________________________________________________________________
1160 AliMUONCDB::WriteCapacitances(const char* filename,
1161 Int_t startRun, Int_t endRun)
1163 /// read manu capacitance and injection gain values from file
1164 /// and store them into CDB located at cdbpath, with a validity period
1165 /// ranging from startRun to endRun
1167 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1168 Int_t ngenerated = MakeCapacitanceStore(*capaStore,filename);
1169 AliInfo(Form("Ngenerated = %d",ngenerated));
1170 if ( ngenerated > 0 )
1172 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,filename);
1177 //_____________________________________________________________________________
1179 AliMUONCDB::WriteCapacitances(Bool_t defaultValues,
1180 Int_t startRun, Int_t endRun)
1182 /// generate manu capacitance values (either 1 if defaultValues=true or random
1183 /// if defaultValues=false, see makeCapacitanceStore) and
1184 /// store them into CDB located at cdbpath, with a validity period
1185 /// ranging from startRun to endRun
1187 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1188 Int_t ngenerated = MakeCapacitanceStore(*capaStore,defaultValues);
1189 AliInfo(Form("Ngenerated = %d",ngenerated));
1190 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,defaultValues);
1194 //_____________________________________________________________________________
1196 AliMUONCDB::WriteTrigger(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1198 /// Writes all Trigger related calibration to CDB
1199 WriteTriggerDCS(defaultValues,startRun,endRun);
1200 WriteLocalTriggerMasks(startRun,endRun);
1201 WriteRegionalTriggerConfig(startRun,endRun);
1202 WriteGlobalTriggerConfig(startRun,endRun);
1203 WriteTriggerLut(startRun,endRun);
1204 WriteTriggerEfficiency(startRun,endRun);
1207 //_____________________________________________________________________________
1209 AliMUONCDB::WriteTracker(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1211 /// Writes all Tracker related calibration to CDB
1212 WriteHV(defaultValues,startRun,endRun);
1213 WritePedestals(defaultValues,startRun,endRun);
1214 WriteGains(defaultValues,startRun,endRun);
1215 WriteCapacitances(defaultValues,startRun,endRun);
1216 WriteNeighbours(startRun,endRun);
1217 WriteKillMap(startRun,endRun,defaultValues);