macosx is already taken care by root

[u/mrichter/AliRoot.git] / PHOS / AliPHOSCalibData.cxx

diff --git a/PHOS/AliPHOSCalibData.cxx b/PHOS/AliPHOSCalibData.cxx
index 4b9bf08e51f40aef9b5737ac95fd9742c12263b9..d18cae2dc286c1270c91a6ad5bda3ac55a8cc5bd 100644 (file)
--- a/PHOS/AliPHOSCalibData.cxx
+++ b/PHOS/AliPHOSCalibData.cxx
@@ -22,6 +22,7 @@
 ///////////////////////////////////////////////////////////////////////////////
 
 #include "TRandom.h"
+#include "AliLog.h"
 #include "AliPHOSCalibData.h"
 #include "AliCDBManager.h"
 #include "AliCDBStorage.h"
@@ -44,16 +45,24 @@ ClassImp(AliPHOSCalibData)
     fCpvDataPath("PHOS/Calib/CpvGainPedestals"),
     fEmcBadChannelsMapPath("PHOS/Calib/EmcBadChannels")
 {
-  // Default constructor
+  // 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
   
    AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data());
   if(entryEmc)
     fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
-  
+
+  if(!fCalibDataEmc)
+    AliWarning("Calibration parameters for PHOS EMC not found. Create a new set.\n");
+    
   AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data());
   if(entryCpv)
     fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
 
+  if(!fCalibDataCpv)
+    AliWarning("Calibration parameters for PHOS CPV not found. Create a new set.\n");
+
   AliCDBEntry* entryEmcBadMap = AliCDBManager::Instance()->Get(fEmcBadChannelsMapPath.Data());
   if(entryEmcBadMap)
     fEmcBadChannelsMap = (AliPHOSEmcBadChannelsMap*)entryEmcBadMap->GetObject(); 
@@ -69,14 +78,23 @@ AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
   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)
@@ -122,10 +140,12 @@ AliPHOSCalibData & AliPHOSCalibData::operator = (const AliPHOSCalibData & rhs)
 //________________________________________________________________
 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();
 }
 
 //________________________________________________________________
@@ -233,6 +253,17 @@ Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row
     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
 {
@@ -249,27 +280,97 @@ Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t ro
 }
 
 //________________________________________________________________
-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->SetADCpedestalEmc(module,column,row,value);
+ 
+  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 1.0; // default width of one EMC ADC channel in GeV
+}
+ 
+//________________________________________________________________
+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);
+}
+//________________________________________________________________
+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
 {
@@ -398,3 +499,13 @@ Int_t AliPHOSCalibData::GetNumOfEmcBadChannels() const
   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);
+}