// //
///////////////////////////////////////////////////////////////////////////////
+#include "TRandom.h"
#include "AliPHOSCalibData.h"
+#include "AliCDBManager.h"
+#include "AliCDBStorage.h"
+#include "AliCDBId.h"
+#include "AliCDBEntry.h"
+#include "AliPHOSEmcCalibData.h"
+#include "AliPHOSCpvCalibData.h"
+#include "AliCDBMetaData.h"
ClassImp(AliPHOSCalibData)
//________________________________________________________________
-AliPHOSCalibData::AliPHOSCalibData()
+ AliPHOSCalibData::AliPHOSCalibData():
+ TNamed(), fCalibDataEmc(0x0), fCalibDataCpv(0x0)
{
- // Default constructor
- Reset();
+ // Default constructor
+
+ fEmcDataPath="PHOS/Calib/EmcGainPedestals";
+ fCpvDataPath="PHOS/Calib/CpvGainPedestals";
+
+ AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data());
+ if(entryEmc)
+ fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
+
+ AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data());
+ if(entryCpv)
+ fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
+
}
//________________________________________________________________
-AliPHOSCalibData::AliPHOSCalibData(const char* name)
+AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
+ TNamed("phosCalib","PHOS Calibration Data Manager"),
+ fCalibDataEmc(0x0), fCalibDataCpv(0x0)
{
// Constructor
- TString namst = "Calib_";
- namst += name;
- SetName(namst.Data());
- SetTitle(namst.Data());
- Reset();
-}
-
-//________________________________________________________________
-AliPHOSCalibData::AliPHOSCalibData(const AliPHOSCalibData& calibda) :
- TNamed(calibda)
-{
- // copy constructor
- SetName(calibda.GetName());
- SetTitle(calibda.GetName());
- Reset();
- for(Int_t module=0; module<5; module++) {
- for(Int_t column=0; column<64; column++) {
- for(Int_t row=0; row<64; row++) {
- fADCchannelEmc[module][column][row] = calibda.GetADCchannelEmc(module,column,row);
- fADCpedestalEmc[module][column][row] = calibda.GetADCpedestalEmc(module,column,row);
- }
- }
- }
+
+ fEmcDataPath="PHOS/Calib/EmcGainPedestals";
+ fCpvDataPath="PHOS/Calib/CpvGainPedestals";
+
+ AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
+ if(entryEmc)
+ fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
+
+ AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
+ if(entryCpv)
+ fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
+
}
//________________________________________________________________
-AliPHOSCalibData &AliPHOSCalibData::operator =(const AliPHOSCalibData& calibda)
+AliPHOSCalibData::AliPHOSCalibData(AliPHOSCalibData & phosCDB) :
+ TNamed(phosCDB)
{
- // assignment operator
- SetName(calibda.GetName());
- SetTitle(calibda.GetName());
- Reset();
- for(Int_t module=0; module<5; module++) {
- for(Int_t column=0; column<64; column++) {
- for(Int_t row=0; row<64; row++) {
- fADCchannelEmc[module][column][row] = calibda.GetADCchannelEmc(module,column,row);
- fADCpedestalEmc[module][column][row] = calibda.GetADCpedestalEmc(module,column,row);
- }
- }
- }
- return *this;
-}
+ // Copy constructor
+ fCalibDataEmc = phosCDB.fCalibDataEmc;
+ fCalibDataCpv = phosCDB.fCalibDataCpv;
+
+ fEmcDataPath = phosCDB.fEmcDataPath;
+ fCpvDataPath = phosCDB.fCpvDataPath;
+}
//________________________________________________________________
AliPHOSCalibData::~AliPHOSCalibData()
{
// Destructor
+
}
//________________________________________________________________
void AliPHOSCalibData::Reset()
{
// Set all pedestals to 0 and all ADC channels to 1
- memset(fADCchannelEmc ,0,5*64*56*sizeof(Float_t));
- memset(fADCpedestalEmc,1,5*64*56*sizeof(Float_t));
+
+ fCalibDataEmc->Reset();
+ fCalibDataCpv->Reset();
}
//________________________________________________________________
void AliPHOSCalibData::Print(Option_t *option) const
{
- // Print tables of pedestals and ADC channels
-
- if (strstr(option,"ped") == 0) {
- printf("\n ---- Pedestal values ----\n\n");
- for (Int_t module=0; module<5; module++){
- printf("============== Module %d\n",module+1);
- for (Int_t column=0; column<64; column++){
- for (Int_t row=0; row<64; row++){
- printf("%4.1f",fADCpedestalEmc[module][column][row]);
- }
- printf("\n");
+ // Print EMC and CPV calibration containers
+ // Input: option="ped" to print pedestals
+ // option="gain" to print calibration coefficients
+ if (fCalibDataEmc) fCalibDataEmc->Print(option);
+ if (fCalibDataCpv) fCalibDataCpv->Print(option);
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::CreateNew()
+{
+ // Create new EMC and CPV calibration containers with ideal coefficients
+
+ if(fCalibDataEmc) delete fCalibDataEmc;
+ fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
+
+ if(fCalibDataCpv) delete fCalibDataCpv;
+ fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
+
+}
+
+//________________________________________________________________
+Bool_t AliPHOSCalibData::WriteEmc(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
+{
+ // Write EMC calibration container to CDB
+
+ if(!fCalibDataEmc) return kFALSE;
+
+ AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
+ if(!storage)
+ storage = AliCDBManager::Instance()->GetDefaultStorage();
+
+ if(storage) {
+ AliCDBId id(fEmcDataPath.Data(),firstRun,lastRun);
+ storage->Put(fCalibDataEmc,id, md);
+ return kTRUE;
+ }
+ else
+ return kFALSE;
+
+}
+
+//________________________________________________________________
+Bool_t AliPHOSCalibData::WriteCpv(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
+{
+ // Write CPV calibration container to CDB
+
+ if(!fCalibDataCpv) return kFALSE;
+
+ AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
+ if(!storage)
+ storage = AliCDBManager::Instance()->GetDefaultStorage();
+
+ if(storage) {
+ AliCDBId id(fCpvDataPath.Data(),firstRun,lastRun);
+ storage->Put(fCalibDataCpv,id, md);
+ return kTRUE;
+ }
+ else
+ return kFALSE;
+
+}
+
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
+{
+ // Return EMC calibration coefficient
+ // for channel defined by (module,column,row)
+ // module, column,raw should follow the internal PHOS convention:
+ // module 1:5, column 1:56, row 1:64
+ // if CBD instance exists, the value is taken from CDB.
+ // Otherwise it is an ideal one
+
+ if(fCalibDataEmc)
+ return fCalibDataEmc->GetADCchannelEmc(module,column,row);
+ else
+ return 0.0015; // default width of one EMC ADC channel in GeV
+}
+
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
+{
+ // Return EMC pedestal for channel defined by (module,column,row)
+ // module, column,raw should follow the internal PHOS convention:
+ // module 1:5, column 1:56, row 1:64
+ // if CBD instance exists, the value is taken from CDB.
+ // Otherwise it is an ideal one
+
+ if(fCalibDataEmc)
+ return fCalibDataEmc->GetADCpedestalEmc(module,column,row);
+ else
+ return 0.005; // default EMC ADC pedestal
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
+{
+ // Set EMC calibration coefficient for (module,column,row)
+
+ if(!fCalibDataEmc)
+ fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
+
+ fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
+{
+ // Set EMC pedestal for (module,column,row)
+
+ if(!fCalibDataEmc)
+ fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
+
+ fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
+}
+
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
+{
+ // Return CPV calibration coefficient
+ // for channel defined by (module,column,row)
+ // module, column,raw should follow the internal CPV convention:
+ // module 1:5, column 1:56, row 1:128
+ // if CBD instance exists, the value is taken from CDB.
+ // Otherwise it is an ideal one
+
+ if(fCalibDataCpv)
+ return fCalibDataCpv->GetADCchannelCpv(module,column,row);
+ else
+ return 0.0012; // default width of one ADC channel in CPV 'popugais'
+}
+
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetADCpedestalCpv(Int_t module, Int_t column, Int_t row) const
+{
+ // Return CPV pedestal
+ // for channel defined by (module,column,row)
+ // module, column,raw should follow the internal CPV convention:
+ // module 1:5, column 1:56, row 1:128
+ // if CBD instance exists, the value is taken from CDB.
+ // Otherwise it is an ideal one
+
+ if(fCalibDataCpv)
+ return fCalibDataCpv->GetADCpedestalCpv(module,column,row);
+ else
+ return 0.012; // default CPV ADC pedestal
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetADCchannelCpv(Int_t module, Int_t column, Int_t row, Float_t value)
+{
+ // Set CPV calibration coefficient for (module,column,row)
+
+ if(!fCalibDataCpv)
+ fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
+
+ fCalibDataCpv->SetADCchannelCpv(module,column,row,value);
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetADCpedestalCpv(Int_t module, Int_t column, Int_t row, Float_t value)
+{
+ // Set CPV pedestal for (module,column,row)
+
+ if(!fCalibDataCpv)
+ fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
+
+ fCalibDataCpv->SetADCpedestalCpv(module,column,row,value);
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::RandomEmc()
+{
+ // Create decalibrated EMC with calibration coefficients and pedestals
+ // randomly distributed within hard-coded limits
+
+ if(fCalibDataEmc) delete fCalibDataEmc;
+ fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
+
+ TRandom rn;
+ rn.SetSeed(0); //the seed is set to the current machine clock
+
+ Float_t adcChannelEmc,adcPedestalEmc;
+
+ for(Int_t module=1; module<6; module++) {
+ for(Int_t column=1; column<57; column++) {
+ for(Int_t row=1; row<65; row++) {
+ adcChannelEmc=rn.Uniform(0.00075,0.00375); // Cmax/Cmin = 5, (Cmax-Cmin)/2 = 0.0015
+ adcPedestalEmc=rn.Uniform(0.0045,0.0055); //+-10% spread of pedestals from 0.005
+ fCalibDataEmc->SetADCchannelEmc(module,column,row,adcChannelEmc);
+ fCalibDataEmc->SetADCpedestalEmc(module,column,row,adcPedestalEmc);
}
}
}
- if (strstr(option,"gain") == 0) {
- printf("\n ---- ADC channel values ----\n\n");
- for (Int_t module=0; module<5; module++){
- printf("============== Module %d\n",module+1);
- for (Int_t column=0; column<64; column++){
- for (Int_t row=0; row<64; row++){
- printf("%4.1f",fADCchannelEmc[module][column][row]);
- }
- printf("\n");
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::RandomCpv()
+{
+ // Create decalibrated CPV with calibration coefficients and pedestals
+ // randomly distributed within hard-coded limits
+
+ if(fCalibDataCpv) delete fCalibDataCpv;
+ fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
+
+ TRandom rn;
+ rn.SetSeed(0); //the seed is set to the current machine clock
+
+ Float_t adcChannelCpv,adcPedestalCpv;
+
+ for(Int_t module=1; module<6; module++) {
+ for(Int_t column=1; column<65; column++) {
+ for(Int_t row=1; row<129; row++) {
+ adcChannelCpv =rn.Uniform(0.0009,0.0015); // 0.0012 +- 25%
+ adcPedestalCpv=rn.Uniform(0.0048,0.0192); // Ped[max]/Ped[min] = 4, <Ped> = 0.012
+ fCalibDataCpv->SetADCchannelCpv(module,column,row,adcChannelCpv);
+ fCalibDataCpv->SetADCpedestalCpv(module,column,row,adcPedestalCpv);
}
}
}