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 "AliMUONCalibParamND.h"
38 #include "AliMUONCalibParamNF.h"
39 #include "AliMUONCalibParamNI.h"
40 #include "AliMUONConstants.h"
41 #include "AliMUONRejectList.h"
42 #include "AliMUONTrackerIO.h"
43 #include "AliMUONTriggerEfficiencyCells.h"
44 #include "AliMUONTriggerLut.h"
45 #include "AliMUONVStore.h"
46 #include "AliMUONVCalibParam.h"
47 #include "AliMUONVCalibParam.h"
48 #include "AliMUONGlobalCrateConfig.h"
49 #include "AliMUONRegionalTriggerConfig.h"
52 #include "AliMpConstants.h"
53 #include "AliMpDDLStore.h"
54 #include "AliMpManuStore.h"
55 #include "AliMpDEManager.h"
56 #include "AliMpDetElement.h"
57 #include "AliMpFiles.h"
58 #include "AliMpDCSNamer.h"
59 #include "AliMpManuIterator.h"
60 #include "AliMpSegmentation.h"
61 #include "AliMpStationType.h"
62 #include "AliMpVSegmentation.h"
64 #include "AliCodeTimer.h"
65 #include "AliCDBEntry.h"
66 #include "AliCDBManager.h"
67 #include "AliDCSValue.h"
70 #include <Riostream.h>
76 #include <TObjString.h>
79 #include <TStopwatch.h>
90 //_____________________________________________________________________________
91 AliMUONVStore* Create2DMap()
93 return new AliMUON2DMap(true);
96 //_____________________________________________________________________________
97 void getBoundaries(const AliMUONVStore& store, Int_t dim,
98 Float_t* xmin, Float_t* xmax)
100 /// Assuming the store contains AliMUONVCalibParam objects, compute the
101 /// limits of the value contained in the VCalibParam, for each of its dimensions
102 /// xmin and xmax must be of dimension dim
104 for ( Int_t i = 0; i < dim; ++i )
110 TIter next(store.CreateIterator());
111 AliMUONVCalibParam* value;
113 while ( ( value = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
115 Int_t detElemId = value->ID0();
116 Int_t manuId = value->ID1();
118 const AliMpVSegmentation* seg =
119 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
121 for ( Int_t manuChannel = 0; manuChannel < value->Size(); ++manuChannel )
123 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
124 if (!pad.IsValid()) continue;
126 for ( Int_t i = 0; i < dim; ++i )
128 Float_t x0 = value->ValueAsFloat(manuChannel,i);
130 xmin[i] = TMath::Min(xmin[i],x0);
131 xmax[i] = TMath::Max(xmax[i],x0);
136 for ( Int_t i = 0; i < dim; ++i )
138 if ( TMath::Abs(xmin[i]-xmax[i]) < 1E-3 )
146 //_____________________________________________________________________________
147 Double_t GetRandom(Double_t mean, Double_t sigma, Bool_t mustBePositive)
150 if ( mustBePositive )
154 x = gRandom->Gaus(mean,sigma);
159 x = gRandom->Gaus(mean,sigma);
166 //_____________________________________________________________________________
167 AliMUONCDB::AliMUONCDB(const char* cdbpath)
170 fMaxNofChannelsToGenerate(-1)
174 if ( ! AliMpCDB::LoadDDLStore() ) {
175 AliFatal("Could not access mapping from OCDB !");
178 if ( ! AliMpCDB::LoadManuStore() ) {
179 AliFatal("Could not access run-dependent mapping from OCDB !");
183 //_____________________________________________________________________________
184 AliMUONCDB::~AliMUONCDB()
189 //_____________________________________________________________________________
191 AliMUONCDB::Diff(AliMUONVStore& store1, AliMUONVStore& store2,
194 /// creates a store which contains store1-store2
195 /// if opt="abs" the difference is absolute one,
196 /// if opt="rel" then what is stored is (store1-store2)/store1
197 /// if opt="percent" then what is stored is rel*100
199 /// WARNING Works only for stores which holds AliMUONVCalibParam objects
204 if ( !sopt.Contains("ABS") && !sopt.Contains("REL") && !sopt.Contains("PERCENT") )
206 AliErrorClass(Form("opt %s not supported. Only ABS, REL, PERCENT are",opt));
210 AliMUONVStore* d = static_cast<AliMUONVStore*>(store1.Clone());
212 TIter next(d->CreateIterator());
214 AliMUONVCalibParam* param;
216 while ( ( param = dynamic_cast<AliMUONVCalibParam*>(next() ) ) )
218 Int_t detElemId = param->ID0();
219 Int_t manuId = param->ID1();
221 AliMUONVCalibParam* param2 = dynamic_cast<AliMUONVCalibParam*>(store2.FindObject(detElemId,manuId));
222 //FIXME: this might happen. Handle it.
225 cerr << "param2 is null : FIXME : this might happen !" << endl;
230 for ( Int_t i = 0; i < param->Size(); ++i )
232 for ( Int_t j = 0; j < param->Dimension(); ++j )
235 if ( sopt.Contains("ABS") )
237 value = param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j);
239 else if ( sopt.Contains("REL") || sopt.Contains("PERCENT") )
241 if ( param->ValueAsFloat(i,j) )
243 value = (param->ValueAsFloat(i,j) - param2->ValueAsFloat(i,j))/param->ValueAsFloat(i,j);
249 if ( sopt.Contains("PERCENT") ) value *= 100.0;
251 param->SetValueAsFloat(i,j,value);
258 //_____________________________________________________________________________
260 AliMUONCDB::Plot(const AliMUONVStore& store, const char* name, Int_t nbins)
262 /// Make histograms of each dimension of the AliMUONVCalibParam
263 /// contained inside store.
264 /// It produces histograms named name_0, name_1, etc...
266 TIter next(store.CreateIterator());
267 AliMUONVCalibParam* param;
269 const Int_t kNStations = AliMpConstants::NofTrackingChambers()/2;
270 Int_t* nPerStation = new Int_t[kNStations];
273 for ( Int_t i = 0; i < kNStations; ++i ) nPerStation[i]=0;
275 while ( ( param = static_cast<AliMUONVCalibParam*>(next()) ) )
279 Int_t dim = param->Dimension();
281 Float_t* xmin = new Float_t[dim];
282 Float_t* xmax = new Float_t[dim];
283 getBoundaries(store,dim,xmin,xmax);
285 for ( Int_t i = 0; i < dim; ++i )
287 h[i] = new TH1F(Form("%s_%d",name,i),Form("%s_%d",name,i),
288 nbins,xmin[i],xmax[i]);
289 AliInfo(Form("Created histogram %s",h[i]->GetName()));
293 Int_t detElemId = param->ID0();
294 Int_t manuId = param->ID1();
295 Int_t station = AliMpDEManager::GetChamberId(detElemId)/2;
297 const AliMpVSegmentation* seg =
298 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
300 for ( Int_t manuChannel = 0; manuChannel < param->Size(); ++manuChannel )
302 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
303 if (!pad.IsValid()) continue;
306 ++nPerStation[station];
308 for ( Int_t dim = 0; dim < param->Dimension(); ++dim )
310 h[dim]->Fill(param->ValueAsFloat(manuChannel,dim));
315 for ( Int_t i = 0; i < kNStations; ++i )
317 AliInfo(Form("Station %d %d ",(i+1),nPerStation[i]));
320 AliInfo(Form("Number of channels = %d",n));
322 delete[] nPerStation;
325 //_____________________________________________________________________________
327 AliMUONCDB::MakeHVStore(TMap& aliasMap, Bool_t defaultValues)
329 /// Create a HV store
331 AliMpDCSNamer hvNamer("TRACKER");
333 TObjArray* aliases = hvNamer.GenerateAliases();
338 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
340 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
341 TString& aliasName = alias->String();
342 if ( aliasName.Contains("sw") )
344 // HV Switch (St345 only)
345 TObjArray* valueSet = new TObjArray;
346 valueSet->SetOwner(kTRUE);
348 Bool_t value = kTRUE;
352 Float_t r = gRandom->Uniform();
353 if ( r < 0.007 ) value = kFALSE;
356 for ( UInt_t timeStamp = 0; timeStamp < 60*3; timeStamp += 60 )
358 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
359 valueSet->Add(dcsValue);
361 aliasMap.Add(new TObjString(*alias),valueSet);
366 TObjArray* valueSet = new TObjArray;
367 valueSet->SetOwner(kTRUE);
368 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
370 Float_t value = 1500;
371 if (!defaultValues) value = GetRandom(1750,62.5,true);
372 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
373 valueSet->Add(dcsValue);
375 aliasMap.Add(new TObjString(*alias),valueSet);
382 AliInfo(Form("%d HV channels and %d switches",nChannels,nSwitch));
384 return nChannels+nSwitch;
387 //_____________________________________________________________________________
389 AliMUONCDB::MakeTriggerDCSStore(TMap& aliasMap, Bool_t defaultValues)
391 /// Create a Trigger HV and Currents store
393 AliMpDCSNamer triggerDCSNamer("TRIGGER");
395 TObjArray* aliases = triggerDCSNamer.GenerateAliases();
397 Int_t nChannels[2] = {0, 0};
399 for ( Int_t i = 0; i < aliases->GetEntries(); ++i )
401 TObjString* alias = static_cast<TObjString*>(aliases->At(i));
402 TString& aliasName = alias->String();
404 TObjArray* valueSet = new TObjArray;
405 valueSet->SetOwner(kTRUE);
406 Int_t measureType = triggerDCSNamer.DCSvariableFromDCSAlias(aliasName.Data());
407 for ( UInt_t timeStamp = 0; timeStamp < 60*15; timeStamp += 120 )
410 (measureType == AliMpDCSNamer::kDCSI) ? 2. : 8000.;
411 if (!defaultValues) {
412 switch (measureType){
413 case AliMpDCSNamer::kDCSI:
414 value = GetRandom(2.,0.4,true);
416 case AliMpDCSNamer::kDCSHV:
417 value = GetRandom(8000.,16.,true);
421 AliDCSValue* dcsValue = new AliDCSValue(value,timeStamp);
422 valueSet->Add(dcsValue);
424 aliasMap.Add(new TObjString(*alias),valueSet);
425 ++nChannels[measureType];
430 AliInfo(Form("Trigger channels I -> %i HV -> %i",nChannels[0], nChannels[1]));
432 return nChannels[0] + nChannels[1];
435 //_____________________________________________________________________________
437 AliMUONCDB::MakePedestalStore(AliMUONVStore& pedestalStore, Bool_t defaultValues)
439 /// Create a pedestal store. if defaultValues=true, ped.mean=ped.sigma=1,
440 /// otherwise mean and sigma are from a gaussian (with parameters
441 /// defined below by the kPedestal* constants)
443 AliCodeTimerAuto("");
448 const Int_t kChannels(AliMpConstants::ManuNofChannels());
451 const Float_t kPedestalMeanMeanB(200.);
452 const Float_t kPedestalMeanSigmaB(10.);
453 const Float_t kPedestalSigmaMeanB(1.);
454 const Float_t kPedestalSigmaSigmaB(0.2);
457 const Float_t kPedestalMeanMeanNB(200.);
458 const Float_t kPedestalMeanSigmaNB(10.);
459 const Float_t kPedestalSigmaMeanNB(1.);
460 const Float_t kPedestalSigmaSigmaNB(0.2);
462 const Float_t kFractionOfDeadManu(0.); // within [0.,1.]
467 AliMpManuIterator it;
469 while ( it.Next(detElemId,manuId) )
471 // skip a given fraction of manus
472 if (kFractionOfDeadManu > 0. && gRandom->Uniform() < kFractionOfDeadManu) continue;
476 AliMUONVCalibParam* ped =
477 new AliMUONCalibParamNF(2,kChannels,detElemId,manuId,AliMUONVCalibParam::InvalidFloatValue());
479 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
481 for ( Int_t manuChannel = 0; manuChannel < kChannels; ++manuChannel )
483 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
487 Float_t meanPedestal;
488 Float_t sigmaPedestal;
497 Bool_t positive(kTRUE);
500 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) ) { // manu in non bending plane
502 while ( meanPedestal == 0.0 ) // avoid strict zero
504 meanPedestal = GetRandom(kPedestalMeanMeanNB,kPedestalMeanSigmaNB,positive);
506 sigmaPedestal = GetRandom(kPedestalSigmaMeanNB,kPedestalSigmaSigmaNB,positive);
508 } else { // manu in bending plane
510 while ( meanPedestal == 0.0 ) // avoid strict zero
512 meanPedestal = GetRandom(kPedestalMeanMeanB,kPedestalMeanSigmaB,positive);
514 sigmaPedestal = GetRandom(kPedestalSigmaMeanB,kPedestalSigmaSigmaB,positive);
520 ped->SetValueAsFloat(manuChannel,0,meanPedestal);
521 ped->SetValueAsFloat(manuChannel,1,sigmaPedestal);
524 Bool_t ok = pedestalStore.Add(ped);
527 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
529 if ( fMaxNofChannelsToGenerate > 0 && nchannels >= fMaxNofChannelsToGenerate ) break;
532 AliInfo(Form("%d Manus and %d channels.",nmanus,nchannels));
536 //_____________________________________________________________________________
538 AliMUONCDB::MakeRejectListStore(Bool_t defaultValues)
540 /// Create a reject list
542 AliCodeTimerAuto("");
544 AliMUONRejectList* rl = new AliMUONRejectList;
548 rl->SetDetectionElementProbability(510);
549 rl->SetDetectionElementProbability(508);
556 //_____________________________________________________________________________
558 AliMUONCDB::MakeOccupancyMapStore(AliMUONVStore& occupancyMapStore, Bool_t defaultValues)
560 /// Create an occupancy map.
562 AliCodeTimerAuto("");
569 AliMpManuIterator it;
573 while ( it.Next(detElemId,manuId) )
577 AliMUONVCalibParam* occupancy = new AliMUONCalibParamND(5,1,detElemId,manuId,0);
581 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
583 Int_t numberOfChannelsInManu = de->NofChannelsInManu(manuId);
585 if (!defaultValues) occ = gRandom->Rndm(1.0);
587 Double_t sumn = occ*nevents;
589 occupancy->SetValueAsFloat(0,0,sumn);
590 occupancy->SetValueAsFloat(0,1,sumn);
591 occupancy->SetValueAsFloat(0,2,sumn);
592 occupancy->SetValueAsInt(0,3,numberOfChannelsInManu);
593 occupancy->SetValueAsInt(0,4,nevents);
595 Bool_t ok = occupancyMapStore.Add(occupancy);
598 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
605 //_____________________________________________________________________________
607 AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, const char* file)
609 /// Read the capacitance values from file and append them to the capaStore
611 return AliMUONTrackerIO::ReadCapacitances(file,capaStore);
614 //_____________________________________________________________________________
616 AliMUONCDB::MakeCapacitanceStore(AliMUONVStore& capaStore, Bool_t defaultValues)
618 /// Create a capacitance store. if defaultValues=true, all capa are 1.0,
619 /// otherwise they are from a gaussian with parameters defined in the
620 /// kCapa* constants below.
622 AliCodeTimerAuto("");
626 Int_t nmanusOK(0); // manus for which we got the serial number
628 const Float_t kCapaMean(0.3);
629 const Float_t kCapaSigma(0.1);
630 const Float_t kInjectionGainMean(3);
631 const Float_t kInjectionGainSigma(1);
636 AliMpManuIterator it;
638 while ( it.Next(detElemId,manuId) )
642 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
643 Int_t serialNumber = AliMpManuStore::Instance()->GetManuSerial(detElemId, manuId);
645 if ( serialNumber <= 0 ) continue;
649 AliMUONVCalibParam* capa = static_cast<AliMUONVCalibParam*>(capaStore.FindObject(serialNumber));
653 capa = new AliMUONCalibParamNF(2,AliMpConstants::ManuNofChannels(),serialNumber,0,1.0);
654 Bool_t ok = capaStore.Add(capa);
657 AliError(Form("Could not set serialNumber=%d manuId=%d",serialNumber,manuId));
661 for ( Int_t manuChannel = 0; manuChannel < capa->Size(); ++manuChannel )
663 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
668 Float_t injectionGain;
677 capaValue = GetRandom(kCapaMean,kCapaSigma,kTRUE);
678 injectionGain = GetRandom(kInjectionGainMean,kInjectionGainSigma,kTRUE);
680 capa->SetValueAsFloat(manuChannel,0,capaValue);
681 capa->SetValueAsFloat(manuChannel,1,injectionGain);
686 if ( nmanus ) percent = 100*nmanusOK/nmanus;
687 AliInfo(Form("%5d manus with serial number (out of %5d manus = %3.0f%%)",
688 nmanusOK,nmanus,percent));
689 AliInfo(Form("%5d channels",nchannels));
692 AliWarning("Did not get all serial numbers. capaStore is incomplete !!!!");
698 //_____________________________________________________________________________
700 AliMUONCDB::MakeGainStore(AliMUONVStore& gainStore, Bool_t defaultValues)
702 /// Create a gain store. if defaultValues=true, all gains set so that
703 /// charge = (adc-ped)
705 /// otherwise parameters are taken from gaussians with parameters
706 /// defined in the k* constants below.
708 AliCodeTimerAuto("");
713 const Int_t kSaturation(3000);
714 const Double_t kA0Mean(1.2);
715 const Double_t kA0Sigma(0.1);
716 const Double_t kA1Mean(1E-5);
717 const Double_t kA1Sigma(1E-6);
718 const Double_t kQualMean(0xFF);
719 const Double_t kQualSigma(0x10);
720 const Int_t kThresMean(1600);
721 const Int_t kThresSigma(100);
726 AliMpManuIterator it;
728 while ( it.Next(detElemId,manuId) )
732 AliMUONVCalibParam* gain =
733 new AliMUONCalibParamNF(5,AliMpConstants::ManuNofChannels(),
736 AliMUONVCalibParam::InvalidFloatValue());
738 AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
740 for ( Int_t manuChannel = 0; manuChannel < gain->Size(); ++manuChannel )
742 if ( ! de->IsConnectedChannel(manuId,manuChannel) ) continue;
748 gain->SetValueAsFloat(manuChannel,0,1.0);
749 gain->SetValueAsFloat(manuChannel,1,0.0);
750 gain->SetValueAsInt(manuChannel,2,4095);
751 gain->SetValueAsInt(manuChannel,3,1);
752 gain->SetValueAsInt(manuChannel,4,kSaturation);
756 Bool_t positive(kTRUE);
757 gain->SetValueAsFloat(manuChannel,0,GetRandom(kA0Mean,kA0Sigma,positive));
758 gain->SetValueAsFloat(manuChannel,1,GetRandom(kA1Mean,kA1Sigma,!positive));
759 gain->SetValueAsInt(manuChannel,2,(Int_t)TMath::Nint(GetRandom(kThresMean,kThresSigma,positive)));
760 gain->SetValueAsInt(manuChannel,3,(Int_t)TMath::Nint(GetRandom(kQualMean,kQualSigma,positive)));
761 gain->SetValueAsInt(manuChannel,4,kSaturation);
765 Bool_t ok = gainStore.Add(gain);
768 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
770 if ( fMaxNofChannelsToGenerate > 0 && nchannels >= fMaxNofChannelsToGenerate ) break;
773 AliInfo(Form("%d Manus and %d channels.",nmanus,nchannels));
777 //_____________________________________________________________________________
779 AliMUONCDB::MakeLocalTriggerMaskStore(AliMUONVStore& localBoardMasks) const
781 /// Generate local trigger masks store. All masks are set to FFFF
783 AliCodeTimerAuto("");
786 // Generate fake mask values for all localboards and put that into
787 // one single container (localBoardMasks)
788 for ( Int_t i = 1; i <= AliMpConstants::TotalNofLocalBoards(); ++i )
790 AliMUONVCalibParam* localBoard = new AliMUONCalibParamNI(1,8,i,0,0);
791 for ( Int_t x = 0; x < 2; ++x )
793 for ( Int_t y = 0; y < 4; ++y )
796 localBoard->SetValueAsInt(index,0,0xFFFF);
800 localBoardMasks.Add(localBoard);
805 //_____________________________________________________________________________
807 AliMUONCDB::MakeRegionalTriggerConfigStore(AliMUONRegionalTriggerConfig& rtm) const
809 /// Make a regional trigger config store. Mask is set to FFFF for each local board (Ch.F.)
811 AliCodeTimerAuto("");
813 if ( ! rtm.ReadData(AliMpFiles::LocalTriggerBoardMapping()) ) {
814 AliErrorStream() << "Error when reading from mapping file" << endl;
818 return rtm.GetNofTriggerCrates();
822 //_____________________________________________________________________________
824 AliMUONCDB::MakeGlobalTriggerConfigStore(AliMUONGlobalCrateConfig& gtm) const
826 /// Make a global trigger config store. All masks (disable) set to 0x00 for each Darc board (Ch.F.)
828 AliCodeTimerAuto("");
830 return gtm.ReadData(AliMpFiles::GlobalTriggerBoardMapping());
834 //_____________________________________________________________________________
836 AliMUONCDB::MakeTriggerLUT(const char* file) const
838 /// Make a triggerlut object, from a file.
840 AliCodeTimerAuto("");
842 AliMUONTriggerLut* lut = new AliMUONTriggerLut;
843 lut->ReadFromFile(file);
847 //_____________________________________________________________________________
848 AliMUONTriggerEfficiencyCells*
849 AliMUONCDB::MakeTriggerEfficiency(const char* file) const
851 /// Make a trigger efficiency object from a file.
853 AliCodeTimerAuto("");
855 return new AliMUONTriggerEfficiencyCells(file);
858 //_____________________________________________________________________________
860 AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
861 Int_t startRun, Int_t endRun,
862 const char* filename)
864 /// Write a given object to OCDB
866 TString comment(gSystem->ExpandPathName(filename));
868 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
871 //_____________________________________________________________________________
873 AliMUONCDB::WriteToCDB(const char* calibpath, TObject* object,
874 Int_t startRun, Int_t endRun, Bool_t defaultValues)
876 /// Write a given object to OCDB
879 if ( defaultValues ) comment += "Test with default values";
880 else comment += "Test with random values";
882 WriteToCDB(object, calibpath, startRun, endRun, comment.Data());
885 //_____________________________________________________________________________
887 AliMUONCDB::WriteToCDB(TObject* object, const char* calibpath, Int_t startRun, Int_t endRun,
888 const char* comment, const char* responsible)
890 /// Write a given object to OCDB
892 AliCDBId id(calibpath,startRun,endRun);
894 md.SetAliRootVersion(gROOT->GetVersion());
895 md.SetComment(comment);
896 md.SetResponsible(responsible);
897 AliCDBManager* man = AliCDBManager::Instance();
898 if (!man->IsDefaultStorageSet()) man->SetDefaultStorage(fCDBPath);
899 man->Put(object,id,&md);
902 //_____________________________________________________________________________
904 AliMUONCDB::MakeNeighbourStore(AliMUONVStore& neighbourStore)
906 /// Fill the neighbours store with, for each channel, a TObjArray of its
907 /// neighbouring pads (including itself)
909 AliCodeTimerAuto("");
911 AliInfo("Generating NeighbourStore. This will take a while. Please be patient.");
920 AliMpManuIterator it;
922 while ( it.Next(detElemId,manuId) )
924 const AliMpVSegmentation* seg =
925 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
927 AliMUONVCalibParam* calibParam = static_cast<AliMUONVCalibParam*>(neighbourStore.FindObject(detElemId,manuId));
931 Int_t size(AliMpConstants::ManuNofChannels());
932 Int_t defaultValue(-1);
933 Int_t packingFactor(size);
935 calibParam = new AliMUONCalibParamNI(dimension,size,detElemId,manuId,defaultValue,packingFactor);
936 Bool_t ok = neighbourStore.Add(calibParam);
939 AliError(Form("Could not set DetElemId=%d manuId=%d",detElemId,manuId));
944 for ( Int_t manuChannel = 0; manuChannel < AliMpConstants::ManuNofChannels(); ++manuChannel )
946 AliMpPad pad = seg->PadByLocation(manuId,manuChannel,kFALSE);
952 seg->GetNeighbours(pad,tmp,true,true);
953 Int_t nofPadNeighbours = tmp.GetEntriesFast();
955 for ( Int_t i = 0; i < nofPadNeighbours; ++i )
957 AliMpPad* p = static_cast<AliMpPad*>(tmp.UncheckedAt(i));
960 calibParam->PackValues(p->GetManuId(),p->GetManuChannel(),x);
963 // AliError("Could not pack value. Something is seriously wrong. Please check");
964 // StdoutToAliError(pad->Print(););
967 calibParam->SetValueAsInt(manuChannel,i,x);
976 //_____________________________________________________________________________
978 AliMUONCDB::SetMaxNofChannelsToGenerate(Int_t n)
980 /// Set the maximum number of channels to generate (used for testing only)
981 /// n < 0 means no limit
982 fMaxNofChannelsToGenerate = n;
985 //_____________________________________________________________________________
987 AliMUONCDB::WriteLocalTriggerMasks(Int_t startRun, Int_t endRun)
989 /// Write local trigger masks to OCDB
991 AliMUONVStore* ltm = new AliMUON1DArray(AliMpConstants::TotalNofLocalBoards()+1);
992 Int_t ngenerated = MakeLocalTriggerMaskStore(*ltm);
993 AliInfo(Form("Ngenerated = %d",ngenerated));
996 WriteToCDB("MUON/Calib/LocalTriggerBoardMasks",ltm,startRun,endRun,true);
1001 //_____________________________________________________________________________
1003 AliMUONCDB::WriteRegionalTriggerConfig(Int_t startRun, Int_t endRun)
1005 /// Write regional trigger masks to OCDB
1007 AliMUONRegionalTriggerConfig* rtm = new AliMUONRegionalTriggerConfig();
1008 Int_t ngenerated = MakeRegionalTriggerConfigStore(*rtm);
1009 AliInfo(Form("Ngenerated = %d",ngenerated));
1012 WriteToCDB("MUON/Calib/RegionalTriggerConfig",rtm,startRun,endRun,true);
1018 //_____________________________________________________________________________
1020 AliMUONCDB::WriteGlobalTriggerConfig(Int_t startRun, Int_t endRun)
1022 /// Write global trigger masks to OCDB
1024 AliMUONGlobalCrateConfig* gtm = new AliMUONGlobalCrateConfig();
1026 Int_t ngenerated = MakeGlobalTriggerConfigStore(*gtm);
1027 AliInfo(Form("Ngenerated = %d",ngenerated));
1030 WriteToCDB("MUON/Calib/GlobalTriggerCrateConfig",gtm,startRun,endRun,true);
1036 //_____________________________________________________________________________
1038 AliMUONCDB::WriteTriggerLut(Int_t startRun, Int_t endRun)
1040 /// Write trigger LUT to OCDB
1042 AliMUONTriggerLut* lut = MakeTriggerLUT();
1045 WriteToCDB("MUON/Calib/TriggerLut",lut,startRun,endRun,true);
1050 //_____________________________________________________________________________
1052 AliMUONCDB::WriteTriggerEfficiency(Int_t startRun, Int_t endRun)
1054 /// Write trigger efficiency to OCDB
1056 AliMUONTriggerEfficiencyCells* eff = MakeTriggerEfficiency();
1059 WriteToCDB("MUON/Calib/TriggerEfficiency",eff,startRun,endRun,true);
1064 //_____________________________________________________________________________
1066 AliMUONCDB::WriteNeighbours(Int_t startRun, Int_t endRun)
1068 /// Write neighbours to OCDB
1070 AliMUONVStore* neighbours = Create2DMap();
1071 Int_t ngenerated = MakeNeighbourStore(*neighbours);
1072 AliInfo(Form("Ngenerated = %d",ngenerated));
1075 WriteToCDB("MUON/Calib/Neighbours",neighbours,startRun,endRun,true);
1080 //_____________________________________________________________________________
1082 AliMUONCDB::WriteHV(Bool_t defaultValues,
1083 Int_t startRun, Int_t endRun)
1085 /// generate HV values (either cste = 1500 V) if defaultValues=true or random
1086 /// if defaultValues=false, see makeHVStore) and
1087 /// store them into CDB located at cdbpath, with a validity period
1088 /// ranging from startRun to endRun
1090 TMap* hvStore = new TMap;
1091 Int_t ngenerated = MakeHVStore(*hvStore,defaultValues);
1092 AliInfo(Form("Ngenerated = %d",ngenerated));
1095 WriteToCDB("MUON/Calib/HV",hvStore,startRun,endRun,defaultValues);
1100 //_____________________________________________________________________________
1102 AliMUONCDB::WriteTriggerDCS(Bool_t defaultValues,
1103 Int_t startRun, Int_t endRun)
1105 /// generate Trigger HV and current values (either const if defaultValues=true or random
1106 /// if defaultValues=false, see makeTriggerDCSStore) and
1107 /// store them into CDB located at cdbpath, with a validity period
1108 /// ranging from startRun to endRun
1110 TMap* triggerDCSStore = new TMap;
1111 Int_t ngenerated = MakeTriggerDCSStore(*triggerDCSStore,defaultValues);
1112 AliInfo(Form("Ngenerated = %d",ngenerated));
1115 WriteToCDB("MUON/Calib/TriggerDCS",triggerDCSStore,startRun,endRun,defaultValues);
1117 delete triggerDCSStore;
1120 //_____________________________________________________________________________
1122 AliMUONCDB::WritePedestals(Bool_t defaultValues,
1123 Int_t startRun, Int_t endRun)
1125 /// generate pedestal values (either 0 if defaultValues=true or random
1126 /// if defaultValues=false, see makePedestalStore) and
1127 /// store them into CDB located at cdbpath, with a validity period
1128 /// ranging from startRun to endRun
1130 AliMUONVStore* pedestalStore = Create2DMap();
1131 Int_t ngenerated = MakePedestalStore(*pedestalStore,defaultValues);
1132 AliInfo(Form("Ngenerated = %d",ngenerated));
1133 WriteToCDB("MUON/Calib/Pedestals",pedestalStore,startRun,endRun,defaultValues);
1134 delete pedestalStore;
1137 //_____________________________________________________________________________
1139 AliMUONCDB::WriteOccupancyMap(Bool_t defaultValues,
1140 Int_t startRun, Int_t endRun)
1142 /// generate occupancy map values (either empty one if defaultValues=true, or
1143 /// random one, see MakeOccupancyMapStore) and
1144 /// store them into CDB located at cdbpath, with a validity period
1145 /// ranging from startRun to endRun
1147 AliMUONVStore* occupancyMapStore = Create2DMap();
1148 Int_t ngenerated = MakeOccupancyMapStore(*occupancyMapStore,defaultValues);
1149 AliInfo(Form("Ngenerated = %d",ngenerated));
1150 WriteToCDB("MUON/Calib/OccupancyMap",occupancyMapStore,startRun,endRun,defaultValues);
1151 delete occupancyMapStore;
1154 //_____________________________________________________________________________
1156 AliMUONCDB::WriteRejectList(Bool_t defaultValues,
1157 Int_t startRun, Int_t endRun)
1159 /// generate reject list values (either empty one if defaultValues=true, or
1160 /// random one, see MakeRejectListStore) and
1161 /// store them into CDB located at cdbpath, with a validity period
1162 /// ranging from startRun to endRun
1164 AliMUONRejectList* rl = MakeRejectListStore(defaultValues);
1165 WriteToCDB("MUON/Calib/RejectList",rl,startRun,endRun,defaultValues);
1170 //_____________________________________________________________________________
1172 AliMUONCDB::WriteGains(Bool_t defaultValues,
1173 Int_t startRun, Int_t endRun)
1175 /// generate gain values (either 1 if defaultValues=true or random
1176 /// if defaultValues=false, see makeGainStore) and
1177 /// store them into CDB located at cdbpath, with a validity period
1178 /// ranging from startRun to endRun
1180 AliMUONVStore* gainStore = Create2DMap();
1181 Int_t ngenerated = MakeGainStore(*gainStore,defaultValues);
1182 AliInfo(Form("Ngenerated = %d",ngenerated));
1183 WriteToCDB("MUON/Calib/Gains",gainStore,startRun,endRun,defaultValues);
1187 //_____________________________________________________________________________
1189 AliMUONCDB::WriteCapacitances(const char* filename,
1190 Int_t startRun, Int_t endRun)
1192 /// read manu capacitance and injection gain values from file
1193 /// and store them into CDB located at cdbpath, with a validity period
1194 /// ranging from startRun to endRun
1196 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1197 Int_t ngenerated = MakeCapacitanceStore(*capaStore,filename);
1198 AliInfo(Form("Ngenerated = %d",ngenerated));
1199 if ( ngenerated > 0 )
1201 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,filename);
1206 //_____________________________________________________________________________
1208 AliMUONCDB::WriteCapacitances(Bool_t defaultValues,
1209 Int_t startRun, Int_t endRun)
1211 /// generate manu capacitance values (either 1 if defaultValues=true or random
1212 /// if defaultValues=false, see makeCapacitanceStore) and
1213 /// store them into CDB located at cdbpath, with a validity period
1214 /// ranging from startRun to endRun
1216 AliMUONVStore* capaStore = new AliMUON1DMap(16828);
1217 Int_t ngenerated = MakeCapacitanceStore(*capaStore,defaultValues);
1218 AliInfo(Form("Ngenerated = %d",ngenerated));
1219 WriteToCDB("MUON/Calib/Capacitances",capaStore,startRun,endRun,defaultValues);
1223 //_____________________________________________________________________________
1225 AliMUONCDB::WriteTrigger(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1227 /// Writes all Trigger related calibration to CDB
1228 WriteTriggerDCS(defaultValues,startRun,endRun);
1229 WriteLocalTriggerMasks(startRun,endRun);
1230 WriteRegionalTriggerConfig(startRun,endRun);
1231 WriteGlobalTriggerConfig(startRun,endRun);
1232 WriteTriggerLut(startRun,endRun);
1233 WriteTriggerEfficiency(startRun,endRun);
1236 //_____________________________________________________________________________
1238 AliMUONCDB::WriteTracker(Bool_t defaultValues, Int_t startRun, Int_t endRun)
1240 /// Writes all Tracker related calibration to CDB
1241 WriteHV(defaultValues,startRun,endRun);
1242 WritePedestals(defaultValues,startRun,endRun);
1243 WriteGains(defaultValues,startRun,endRun);
1244 WriteCapacitances(defaultValues,startRun,endRun);
1245 WriteNeighbours(startRun,endRun);
1246 WriteOccupancyMap(startRun,endRun,defaultValues);
1247 WriteRejectList(startRun,endRun,defaultValues);