///////////////////////////////////////////////////////////////////////////////
#include "TRandom.h"
+#include "AliLog.h"
#include "AliPHOSCalibData.h"
#include "AliCDBManager.h"
#include "AliCDBStorage.h"
#include "AliCDBEntry.h"
#include "AliPHOSEmcCalibData.h"
#include "AliPHOSCpvCalibData.h"
+#include "AliPHOSEmcBadChannelsMap.h"
#include "AliCDBMetaData.h"
ClassImp(AliPHOSCalibData)
TNamed(),
fCalibDataEmc(0x0),
fCalibDataCpv(0x0),
+ fEmcBadChannelsMap(0x0),
fEmcDataPath("PHOS/Calib/EmcGainPedestals"),
- fCpvDataPath("PHOS/Calib/CpvGainPedestals")
+ fCpvDataPath("PHOS/Calib/CpvGainPedestals"),
+ fEmcBadChannelsMapPath("PHOS/Calib/EmcBadChannels")
{
- // Default constructor
-
- 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();
-
+ // Default constructor.
+ // Open CDB entry, get EMC and CPV calibration data and bad channel map.
+ // If EMC or CPV calibration data does not exist, stop the run
+
}
//________________________________________________________________
AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
TNamed("phosCalib","PHOS Calibration Data Manager"),
- fCalibDataEmc(0x0), fCalibDataCpv(0x0),
+ fCalibDataEmc(0x0), fCalibDataCpv(0x0), fEmcBadChannelsMap(0x0),
fEmcDataPath("PHOS/Calib/EmcGainPedestals"),
- fCpvDataPath("PHOS/Calib/CpvGainPedestals")
+ fCpvDataPath("PHOS/Calib/CpvGainPedestals"),
+ fEmcBadChannelsMapPath("PHOS/Calib/EmcBadChannels")
{
// Constructor
+ // Open CDB entry, get EMC and CPV calibration data and bad channel map.
+ // If EMC or CPV calibration data does not exist, stop the run
+
AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
if(entryEmc)
fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
+
+ if(!fCalibDataEmc)
+ AliFatal("Calibration parameters for PHOS EMC not found. Stop reconstruction!\n");
AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
if(entryCpv)
fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
+ if(!fCalibDataCpv)
+ AliFatal("Calibration parameters for PHOS CPV not found. Stop reconstruction!\n");
+
+ AliCDBEntry* entryEmcBadMap = AliCDBManager::Instance()->
+ Get(fEmcBadChannelsMapPath.Data(),runNumber);
+ if(entryEmcBadMap)
+ fEmcBadChannelsMap = (AliPHOSEmcBadChannelsMap*)entryEmcBadMap->GetObject();
+
}
//________________________________________________________________
TNamed(phosCDB),
fCalibDataEmc(phosCDB.fCalibDataEmc),
fCalibDataCpv(phosCDB.fCalibDataCpv),
+ fEmcBadChannelsMap(phosCDB.fEmcBadChannelsMap),
fEmcDataPath(phosCDB.fEmcDataPath),
- fCpvDataPath(phosCDB.fCpvDataPath)
+ fCpvDataPath(phosCDB.fCpvDataPath),
+ fEmcBadChannelsMapPath(phosCDB.fEmcBadChannelsMapPath)
{
// Copy constructor
}
}
+//________________________________________________________________
AliPHOSCalibData & AliPHOSCalibData::operator = (const AliPHOSCalibData & rhs)
{
//Copy-assignment. Does not delete anything (see destructor)
//compiler generated is ok, but ... because -Weffc++ and pointer
//members we have to define it explicitly.
TNamed::operator=(rhs);
- fCalibDataEmc = rhs.fCalibDataEmc;
- fCalibDataCpv = rhs.fCalibDataCpv;
- fEmcDataPath = rhs.fEmcDataPath;
- fCpvDataPath = rhs.fCpvDataPath;
-
+ if (this != &rhs) {
+ fCalibDataEmc = rhs.fCalibDataEmc;
+ fCalibDataCpv = rhs.fCalibDataCpv;
+ fEmcBadChannelsMap = rhs.fEmcBadChannelsMap;
+ fEmcDataPath = rhs.fEmcDataPath;
+ fCpvDataPath = rhs.fCpvDataPath;
+ fEmcBadChannelsMapPath = rhs.fEmcBadChannelsMapPath;
+ }
+ else {
+ AliFatal("Self assignment!");
+ }
return *this;
}
//________________________________________________________________
void AliPHOSCalibData::Reset()
{
- // Set all pedestals to 0 and all ADC channels to 1
+ // Set all pedestals to 0 and all ADC channels to 1,
+ // and all channels are good (alive)
- fCalibDataEmc->Reset();
- fCalibDataCpv->Reset();
+ fCalibDataEmc ->Reset();
+ fCalibDataCpv ->Reset();
+ fEmcBadChannelsMap->Reset();
}
//________________________________________________________________
if(fCalibDataCpv) delete fCalibDataCpv;
fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
+ if(fEmcBadChannelsMap) delete fEmcBadChannelsMap;
+ fEmcBadChannelsMap = new AliPHOSEmcBadChannelsMap();
+
}
//________________________________________________________________
}
+
+//________________________________________________________________
+Bool_t AliPHOSCalibData::WriteEmcBadChannelsMap(Int_t firstRun,Int_t lastRun,AliCDBMetaData *md)
+{
+ //Write EMC bad channels map into CDB.
+
+ if(!fEmcBadChannelsMap) return kFALSE;
+
+ AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS/*");
+ if(!storage)
+ storage = AliCDBManager::Instance()->GetDefaultStorage();
+
+ if(storage) {
+ AliCDBId id(fEmcBadChannelsMapPath.Data(),firstRun,lastRun);
+ storage->Put(fEmcBadChannelsMap,id, md);
+ return kTRUE;
+ }
+ else
+ return kFALSE;
+}
+
//________________________________________________________________
Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
{
return 1.0; // default width of one EMC ADC channel in GeV
}
+//________________________________________________________________
+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);
+}
+
//________________________________________________________________
Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
{
}
//________________________________________________________________
-void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
+void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
{
- // Set EMC calibration coefficient for (module,column,row)
+ // Set EMC pedestal for (module,column,row)
if(!fCalibDataEmc)
fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
- fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
+ fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
}
//________________________________________________________________
-void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
+Float_t AliPHOSCalibData::GetHighLowRatioEmc(Int_t module, Int_t column, Int_t row) const
{
- // Set EMC pedestal for (module,column,row)
-
+ // 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->GetHighLowRatioEmc(module,column,row);
+ else
+ return 1.0; // default width of one EMC ADC channel in GeV
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetHighLowRatioEmc(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->SetHighLowRatioEmc(module,column,row,value);
+}
+
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetTimeShiftEmc(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->GetTimeShiftEmc(module,column,row);
+ else
+ return 0.0; // by default no offset
+}
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetLGTimeShiftEmc(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->GetLGTimeShiftEmc(module,column,row);
+ else
+ return 0.0; // no offset by default
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetLGTimeShiftEmc(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->SetLGTimeShiftEmc(module,column,row,value);
+}
+//________________________________________________________________
+void AliPHOSCalibData::SetTimeShiftEmc(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->SetTimeShiftEmc(module,column,row,value);
+}
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetSampleTimeStep() const
+{
+ //Get conversion coeff. from sample time step to seconds.
+ //Negative value means that it is not used in reconstruction
+ //but only in simulation of raw.
+ if(fCalibDataEmc)
+ return fCalibDataEmc->GetSampleTimeStep();
+ else
+ return 0.0; // default width of one EMC ADC channel in GeV
+}
+//________________________________________________________________
+void AliPHOSCalibData::SetSampleTimeStep(Float_t step)
+{
+ //Set conversion coeff. from sample time step to seconds.
+ //Negative value means that it is not used in reconstruction
+ //but only in simulation of raw.
if(!fCalibDataEmc)
fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
- fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
+ fCalibDataEmc->SetSampleTimeStep(step) ;
+}
+//________________________________________________________________
+Int_t AliPHOSCalibData::GetAltroOffsetEmc(Int_t module, Int_t column, Int_t row) const
+{
+ // Return ALTRO pedestal 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->GetAltroOffsetEmc(module,column,row);
+ else
+ return 0; // default width of one EMC ADC channel in GeV
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetAltroOffsetEmc(Int_t module, Int_t column, Int_t row, Int_t value)
+{
+ // Set altro offset for (module,column,row)
+
+ if(!fCalibDataEmc)
+ fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
+
+ fCalibDataEmc->SetAltroOffsetEmc(module,column,row,value);
}
+
//________________________________________________________________
Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
{
}
}
}
+//________________________________________________________________
+Bool_t AliPHOSCalibData::IsBadChannelEmc(Int_t module, Int_t col, Int_t row) const
+{
+ //If no bad channels map found, channel considered good
+
+ if(fEmcBadChannelsMap)
+ return fEmcBadChannelsMap->IsBadChannel(module,col,row);
+ else
+ return kFALSE;
+}
+
+//________________________________________________________________
+Int_t AliPHOSCalibData::GetNumOfEmcBadChannels() const
+{
+ if(fEmcBadChannelsMap)
+ return fEmcBadChannelsMap->GetNumOfBadChannels();
+ else
+ return 0;
+}
+//________________________________________________________________
+void AliPHOSCalibData::EmcBadChannelIds(Int_t *badIds)
+{
+ //Fill array badIds by the Ids of EMC bad channels.
+ //Array badIds of length GetNumOfBadChannels() should be prepared in advance.
+
+ if(fEmcBadChannelsMap)
+ fEmcBadChannelsMap->BadChannelIds(badIds);
+}
+
+//________________________________________________________________
+Float_t AliPHOSCalibData::GetADCchannelEmcDecalib(Int_t module, Int_t column, Int_t row) const
+{
+ // Return random EMC (de)calibration factor O(1) for channel defined by (module,column,row).
+ // Used in simulation.
+
+ // 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 (no decalibration).
+
+ if(fCalibDataEmc)
+ return fCalibDataEmc->GetADCchannelEmcDecalib(module,column,row);
+ else
+ return 1.0; // no decalibration by default
+}
+
+//________________________________________________________________
+void AliPHOSCalibData::SetADCchannelEmcDecalib(Int_t module, Int_t column, Int_t row, Float_t value)
+{
+ // Set EMC (de)calibration factor for (module,column,row).
+ // Used in simulation.
+
+ if(!fCalibDataEmc)
+ fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
+
+ fCalibDataEmc->SetADCchannelEmcDecalib(module,column,row,value);
+}