]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDcalibDB.cxx
Improved quitting Event Display with close button.
[u/mrichter/AliRoot.git] / TRD / AliTRDcalibDB.cxx
index 4eff278aed89615853b4a44e1c8d6138eca560cf..acd70b6a258417947d65139cdf82b0fab091d56f 100644 (file)
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
-#include <TMath.h>
-#include <TRandom.h>
+#include <TROOT.h>
+#include <TClonesArray.h>
+#include <TObjArray.h>
 
 #include "AliCDBManager.h"
-#include "AliCDBStorage.h"
 #include "AliCDBEntry.h"
 #include "AliLog.h"
 
+#include "AliTRDPIDReference.h"
+#include "AliTRDPIDResponseObject.h"
 #include "AliTRDcalibDB.h"
-#include "AliTRDgeometry.h"
-#include "AliTRDpadPlane.h"
+#include "AliTRDtrapConfig.h"
+#include "AliTRDtrapConfigHandler.h"
 #include "AliTRDCommonParam.h"
+#include "AliTRDgeometry.h"
 
 #include "Cal/AliTRDCalROC.h"
 #include "Cal/AliTRDCalPad.h"
 #include "Cal/AliTRDCalDet.h"
-#include "Cal/AliTRDCalFEE.h"
+#include "Cal/AliTRDCalDCS.h"
+#include "Cal/AliTRDCalDCSFEE.h"
+#include "Cal/AliTRDCalDCSv2.h"
+#include "Cal/AliTRDCalDCSFEEv2.h"
 #include "Cal/AliTRDCalPID.h"
 #include "Cal/AliTRDCalMonitoring.h"
 #include "Cal/AliTRDCalChamberStatus.h"
 #include "Cal/AliTRDCalPadStatus.h"
 #include "Cal/AliTRDCalSingleChamberStatus.h"
+#include "Cal/AliTRDCalTrkAttach.h"
+#include "Cal/AliTRDCalOnlineGainTable.h"
 
 ClassImp(AliTRDcalibDB)
 
@@ -105,6 +113,11 @@ AliTRDcalibDB::AliTRDcalibDB()
   ,fPRFhi(0)
   ,fPRFwid(0)
   ,fPRFpad(0)
+  ,fPIDResponse(NULL)
+  ,fOnlineGainTableID(0)
+  ,fTrapConfig(0x0)
+  ,fTrapConfigName("")
+  ,fTrapConfigVersion("")
 {
   //
   // Default constructor
@@ -134,6 +147,11 @@ AliTRDcalibDB::AliTRDcalibDB(const AliTRDcalibDB &c)
   ,fPRFhi(0)
   ,fPRFwid(0)
   ,fPRFpad(0)
+  ,fPIDResponse(NULL)
+  ,fOnlineGainTableID(0)
+  ,fTrapConfig(0x0)
+  ,fTrapConfigName("")
+  ,fTrapConfigVersion("")
 {
   //
   // Copy constructor (not that it make any sense for a singleton...)
@@ -176,8 +194,13 @@ AliTRDcalibDB::~AliTRDcalibDB()
     fPRFsmp = 0;
   }
 
-  Invalidate();
+  if (fPIDResponse) {
+    delete fPIDResponse;
+    fPIDResponse = 0x0;
+  }
 
+  Invalidate();
+  fgInstance   = 0;
 }
 
 //_caching functions____________________________________________________________
@@ -211,6 +234,9 @@ const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
     case kIDVdriftChamber : 
       return CacheCDBEntry(kIDVdriftChamber     ,"TRD/Calib/ChamberVdrift"); 
       break;
+    case kIDExBChamber : 
+      return CacheCDBEntry(kIDExBChamber        ,"TRD/Calib/ChamberExB"); 
+      break;
     case kIDT0Pad : 
       return CacheCDBEntry(kIDT0Pad             ,"TRD/Calib/LocalT0"); 
       break;
@@ -223,14 +249,58 @@ const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
     case kIDGainFactorChamber : 
       return CacheCDBEntry(kIDGainFactorChamber ,"TRD/Calib/ChamberGainFactor"); 
       break;
+
+    case kIDOnlineGainFactor : 
+      switch(GetOnlineGainTableID()) {
+        case 0:
+          // For testing purposes only !!!
+          AliInfo("No gain table name from OCDB. Use default table!");
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Krypton_2011-01"); 
+          break;
+        case 1:
+         // Online gain table ID 1
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Krypton_2011-01"); 
+          break;
+        case 2:
+         // Online gain table ID 2
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Gaintbl_Uniform_FGAN0_2011-01"); 
+          break;
+        case 3:
+         // Online gain table ID 3
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Gaintbl_Uniform_FGAN8_2011-01"); 
+          break;
+        case 4:
+         // Online gain table ID 4
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Krypton_2011-02"); 
+          break;
+        case 5:
+         // Online gain table ID 5
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Krypton_2011-03"); 
+          break;
+        case 6:
+         // Online gain table ID 6
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Gaintbl_Uniform_FGAN0_2012-01"); 
+          break;
+        case 7:
+         // Online gain table ID 7
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Gaintbl_Uniform_FGAN8_2012-01");
+          break; 
+        case 8:
+         // Online gain table ID 8
+          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Krypton_2012-01"); 
+          break;
+      default:
+       AliError(Form("unknown gaintable requested with ID"));
+      }
+      break;
+
     case kIDNoiseChamber : 
-      return CacheCDBEntry(kIDNoiseChamber ,"TRD/Calib/DetNoise"); 
+      return CacheCDBEntry(kIDNoiseChamber      ,"TRD/Calib/DetNoise"); 
       break;
     case kIDNoisePad : 
-      return CacheCDBEntry(kIDNoisePad ,"TRD/Calib/PadNoise"); 
+      return CacheCDBEntry(kIDNoisePad          ,"TRD/Calib/PadNoise"); 
       break;
 
-
     // Parameters defined per pad
     case kIDPRFWidth : 
       return CacheCDBEntry(kIDPRFWidth          ,"TRD/Calib/PRFWidth"); 
@@ -251,12 +321,30 @@ const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
     case kIDFEE : 
       return CacheCDBEntry(kIDFEE               ,"TRD/Calib/FEE"); 
       break;
+    case kIDTrapConfig :
+      return CacheCDBEntry(kIDTrapConfig        ,"TRD/Calib/TrapConfig"); 
+      break;
+    case kIDDCS :
+      return CacheCDBEntry(kIDDCS               ,"TRD/Calib/DCS");
+      break;
     case kIDPIDNN : 
       return CacheCDBEntry(kIDPIDNN             ,"TRD/Calib/PIDNN");
+      break;
     case kIDPIDLQ : 
       return CacheCDBEntry(kIDPIDLQ             ,"TRD/Calib/PIDLQ"); 
       break;
-
+    case kIDPIDLQ1D:
+      return CacheCDBEntry(kIDPIDLQ1D           ,"TRD/Calib/PIDLQ1D");
+      break;
+    case kIDRecoParam : 
+      return CacheCDBEntry(kIDRecoParam         ,"TRD/Calib/RecoParam"); 
+      break;
+    case kIDAttach : 
+      return CacheCDBEntry(kIDAttach            ,"TRD/Calib/TrkAttach"); 
+      break;
+    case kIDPHQ :
+      return CacheCDBEntry(kIDPHQ               ,"TRD/Calib/PHQ");
+      break;
   }
 
   return 0;
@@ -272,7 +360,7 @@ AliCDBEntry *AliTRDcalibDB::GetCDBEntry(const char *cdbPath)
     
   AliCDBEntry *entry = AliCDBManager::Instance()->Get(cdbPath,fRun);
   if (!entry) { 
-    AliFatal(Form("Failed to get entry: %s",cdbPath));
+    AliError(Form("Failed to get entry: %s",cdbPath));
     return 0; 
   }
   
@@ -286,14 +374,16 @@ const TObject *AliTRDcalibDB::CacheCDBEntry(Int_t id, const char *cdbPath)
   //
   // Caches the entry <id> with cdb path <cdbPath>
   //
-  
+
   if (!fCDBCache[id]) {
     fCDBEntries[id] = GetCDBEntry(cdbPath);
     if (fCDBEntries[id]) {
       fCDBCache[id] = fCDBEntries[id]->GetObject();
+      if (id == kIDOnlineGainFactor)
+       AliInfo(Form("loaded gain table: %s", fCDBEntries[id]->GetId().GetAliCDBPath().GetPath().Data()));
     }
-  }
-
+  } 
+  
   return fCDBCache[id];
 
 }
@@ -337,7 +427,10 @@ void AliTRDcalibDB::Invalidate()
     }
   }
 
+  fOnlineGainTableID = 0;
+
 }
+
 //_____________________________________________________________________________
 Float_t AliTRDcalibDB::GetNoise(Int_t det, Int_t col, Int_t row)
 {
@@ -480,6 +573,16 @@ const AliTRDCalDet *AliTRDcalibDB::GetVdriftDet()
 
 }
 
+//_____________________________________________________________________________
+TObjArray * AliTRDcalibDB::GetPHQ()
+{
+  //
+  //return PHQ calibration object
+  //
+  TObjArray *arr = (TObjArray *) (GetCachedCDBObject(kIDPHQ));
+  return arr;
+}
+
 //_____________________________________________________________________________
 Float_t AliTRDcalibDB::GetVdriftAverage(Int_t det)
 {
@@ -496,7 +599,38 @@ Float_t AliTRDcalibDB::GetVdriftAverage(Int_t det)
   return calDet->GetValue(det);
 
 }
+//_____________________________________________________________________________
+const AliTRDCalDet *AliTRDcalibDB::GetExBDet()
+{
+  //
+  // Returns the exB calibration object
+  // containing one number per detector
+  //
+  
+  const AliTRDCalDet     *calChamber = dynamic_cast<const AliTRDCalDet *> (GetCachedCDBObject(kIDExBChamber));
+  
+  Double_t meanexb = 100.0;
+  if (calChamber) meanexb = calChamber->GetMean();
+  //printf("mean %f\n",meanexb);  
+
+  if ((!calChamber) || (meanexb > 70.0)) {
+    
+    const AliTRDCalDet     *calChambervdrift = dynamic_cast<const AliTRDCalDet *> 
+       (GetCachedCDBObject(kIDVdriftChamber));
+      if (!calChambervdrift) {
+       return 0;
+      }
+      else {
+       AliTRDCalDet *calDetExB = new AliTRDCalDet("lorentz angle tan","lorentz angle tan (detector value)");
+       for(Int_t k = 0; k < 540; k++){
+         calDetExB->SetValue(k,AliTRDCommonParam::Instance()->GetOmegaTau(calChambervdrift->GetValue(k)));
+       }
+       return calDetExB;
+      }
+  }
+  else return calChamber;
 
+}
 //_____________________________________________________________________________
 Float_t AliTRDcalibDB::GetT0(Int_t det, Int_t col, Int_t row)
 {
@@ -592,40 +726,57 @@ Float_t AliTRDcalibDB::GetT0Average(Int_t det)
     return -1;
   }
 
-  Double_t mean = 0.0; 
+  Double_t sum = 0.0; 
   for (Int_t channel = 0; channel < roc->GetNchannels(); ++channel) {
-    mean += (calDet->GetValue(det) + roc->GetValue(channel)) / roc->GetNchannels();
+    sum += roc->GetValue(channel);
   }
-
-  return mean;
+  sum /= roc->GetNchannels();
+  sum += calDet->GetValue(det);
+  return sum;
 
 }
 
 //_____________________________________________________________________________
-Float_t AliTRDcalibDB::GetGainFactor(Int_t det, Int_t col, Int_t row)
+AliTRDCalOnlineGainTableROC* AliTRDcalibDB::GetOnlineGainTableROC(Int_t det)
 {
   //
-  // Returns the gain factor for the given pad.
+  // Returns the online gain factor table for a given ROC.
   //
+
+  if (!HasOnlineFilterGain()) {
+    return 0x0;
+  }
   
-  const AliTRDCalPad *calPad     = dynamic_cast<const AliTRDCalPad *> 
-                                   (GetCachedCDBObject(kIDGainFactorPad));
-  if (!calPad) {
-    return -1;
+  const AliTRDCalOnlineGainTable *calOnline 
+     = dynamic_cast<const AliTRDCalOnlineGainTable *> 
+                                   (GetCachedCDBObject(kIDOnlineGainFactor));
+  if (!calOnline) {
+    return 0x0;
   }
 
-  AliTRDCalROC       *roc        = calPad->GetCalROC(det);
-  if (!roc) {
+  return calOnline->GetGainTableROC(det);
+
+}
+
+//_____________________________________________________________________________
+Float_t AliTRDcalibDB::GetOnlineGainFactor(Int_t det, Int_t col, Int_t row)
+{
+  //
+  // Returns the online gain factor for the given pad.
+  //
+
+  if (!HasOnlineFilterGain()) {
     return -1;
   }
-
-  const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *> 
-                                   (GetCachedCDBObject(kIDGainFactorChamber));
-  if (!calChamber) {
+  
+  const AliTRDCalOnlineGainTable *calOnline 
+     = dynamic_cast<const AliTRDCalOnlineGainTable *> 
+                                   (GetCachedCDBObject(kIDOnlineGainFactor));
+  if (!calOnline) {
     return -1;
   }
 
-  return calChamber->GetValue(det) * roc->GetValue(col,row);
+  return calOnline->GetGainCorrectionFactor(det,row,col);
 
 }
 
@@ -652,6 +803,34 @@ AliTRDCalROC *AliTRDcalibDB::GetGainFactorROC(Int_t det)
 
 }
 
+//_____________________________________________________________________________
+Float_t AliTRDcalibDB::GetGainFactor(Int_t det, Int_t col, Int_t row)
+{
+  //
+  // Returns the gain factor for the given pad.
+  //
+  
+  const AliTRDCalPad *calPad     = dynamic_cast<const AliTRDCalPad *> 
+                                   (GetCachedCDBObject(kIDGainFactorPad));
+  if (!calPad) {
+    return -1;
+  }
+
+  AliTRDCalROC       *roc        = calPad->GetCalROC(det);
+  if (!roc) {
+    return -1;
+  }
+
+  const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *> 
+                                   (GetCachedCDBObject(kIDGainFactorChamber));
+  if (!calChamber) {
+    return -1;
+  }
+
+  return calChamber->GetValue(det) * roc->GetValue(col,row);
+
+}
+
 //_____________________________________________________________________________
 const AliTRDCalDet *AliTRDcalibDB::GetGainFactorDet()
 {
@@ -735,19 +914,667 @@ Float_t AliTRDcalibDB::GetPRFWidth(Int_t det, Int_t col, Int_t row)
 }
   
 //_____________________________________________________________________________
-Int_t AliTRDcalibDB::GetNumberOfTimeBins()
+Int_t AliTRDcalibDB::ExtractTimeBinsFromString(TString tbstr)
+{
+  // Check if there is any content in the string first
+  if (tbstr.Length() == 0) {
+    AliError("Parameter for number of timebins is empty!");
+    return -1;
+  }
+
+  // Check if we have the correct config parameter
+  TString tbident  = "tb";
+  TString tbsubstr = tbstr(0,2);
+  if (!tbsubstr.EqualTo(tbident)) {
+    AliError(Form("Parameter for number of timebins is corrupted (%s)!", tbstr.Data()));
+    return -1;
+  }
+
+  tbstr.Remove(0,2);
+  // check if there is more than a number
+  if (!tbstr.IsDigit()) {
+    AliError(Form("Parameter for number of timebins is corrupted (%s)!", tbstr.Data()));
+    return -1;
+  }
+
+  return tbstr.Atoi();
+
+}
+
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCSBoard(){
+  // This is an old way to extract the number of time bins,
+  // only to be used as a fallback to see whether there's
+  // a patched OCDB object!
+  
+  Int_t nTB=-1;
+  // Get the array with SOR and EOR
+  const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+  if(!dcsArr)
+    return -1; // Error
+  // Check CalDCS version
+  Int_t calver = 0;
+  if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
+  else if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+  // CalDCS version 1
+  if (calver == 1) {
+    // DCS object
+    AliTRDCalDCS *calSOR = dynamic_cast<AliTRDCalDCS *>(dcsArr->At(0));
+    // SOR mandantory
+    if(!calSOR)
+      return -1; // Error
+    else
+      nTB=calSOR->GetGlobalNumberOfTimeBins();
+    AliTRDCalDCS *calEOR = dynamic_cast<AliTRDCalDCS *>(dcsArr->At(1));
+    if(calEOR && calEOR->GetGlobalNumberOfTimeBins()!=nTB)
+      return -2; // Mixed
+  }
+  else if (calver == 2) {
+    // DCS object
+    AliTRDCalDCSv2 *calSOR = dynamic_cast<AliTRDCalDCSv2 *>(dcsArr->At(0));
+    // SOR mandantory
+    if(!calSOR)
+      return -1; // Error
+    else
+      nTB=calSOR->GetGlobalNumberOfTimeBins();
+    AliTRDCalDCSv2 *calEOR = dynamic_cast<AliTRDCalDCSv2 *>(dcsArr->At(1));
+    if(calEOR && calEOR->GetGlobalNumberOfTimeBins()!=nTB)
+      return -2; // Mixed
+  }
+  else{
+    // No version 1 nor version 2 object
+    return -1; // Error
+  }
+
+  // All well
+  return nTB;
+}
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
+{
+  //
+  // Returns number of time bins from the DCS
+  //
+
+  // Initialize with values indicating no information
+  TString cfg="";
+  Int_t nTB = -1;
+  // Extract the global configuration name
+  GetGlobalConfigurationByChamber(cfg,kTimebin);
+  // Extract the number of time bins from
+  // the global configuration 
+  if(cfg.Length()>0){
+    nTB = ExtractTimeBinsFromString(cfg);
+  }
+  if(nTB>0){
+    // All well, we could extract the number
+    // of time bins from the configuration name
+    return nTB;
+  }
+  else{
+    // No number of time bins from config name.
+    // No board responded or similar.
+    // We patched some OCDB entries with 
+    // only the global number of time bins set.
+    // We should get these here
+    nTB=GetNumberOfTimeBinsDCSBoard();
+    if(nTB>0){
+      AliWarning("Using old method for number of time bins."
+                " This is probably a patched OCDB entry");
+      return nTB;
+    }
+    else{
+      // Error
+      AliError("No number of time bins either"
+              " from config name or patched OCDB entry");
+      return -1;
+    }
+  }
+
+  
+  // // Get the corresponding parameter
+  // TString cfgstr = "", cfgname = "";
+  // GetGlobalConfiguration(cfgname);
+  // if(cfgname.Length()==0)
+  //   return -1;
+  // GetDCSConfigParOption(cfgname, kTimebin, 0, cfgstr);
+  // if(cfgstr.Length()==0)
+  //   return -1;
+  // return ExtractTimeBinsFromString(cfgstr);
+
+}
+
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetFilterType(TString &filterType)
+{
+  //
+  // Returns the filter type
+  //
+
+  TString cfgname = "";
+  GetGlobalConfiguration(cfgname);
+  GetDCSConfigParOption(cfgname, kFltrSet, 0, filterType);
+
+}
+
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetOnlineGainTableID()
 {
   //
-  // Returns the number of time bins which are read-out.
+  // Get the gain table ID from the DCS
   //
 
-  const AliTRDCalFEE *calFEE     = dynamic_cast<const AliTRDCalFEE *> 
-                                   (GetCachedCDBObject(kIDFEE));
-  if (!calFEE) {
+  if (fOnlineGainTableID > 0) {
+    return fOnlineGainTableID;
+  }
+
+  const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+  if (!dcsArr){
     return -1;
   }
 
-  return calFEE->GetNumberOfTimeBins();
+  Int_t esor   = 0; // Take SOR
+  Int_t calver = 0; // Check CalDCS version
+  if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS"))   calver = 1;
+  if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+
+  if      (calver == 1) {
+
+    // No data for old DCS object available, anyway
+    return -1;
+
+  } 
+  else if (calver == 2) {
+
+    // DCSv2 object
+    const AliTRDCalDCSv2 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
+    if(!calDCSv2){
+      return -1;
+    }
+
+    TString tableName = "";
+    for (Int_t i = 0; i < 540; i++) {
+      const AliTRDCalDCSFEEv2 *calDCSFEEv2 = calDCSv2->GetCalDCSFEEObj(0);
+      tableName = calDCSFEEv2->GetGainTableName();
+      if (tableName.Length() > 0) {
+        break;
+      }
+    }
+    if (tableName.CompareTo("Krypton_2011-01")               == 0)
+      fOnlineGainTableID = 1;
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2011-01") == 0)
+      fOnlineGainTableID = 2;
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2011-01") == 0)
+      fOnlineGainTableID = 3;
+    else if (tableName.CompareTo("Krypton_2011-02")               == 0)
+      fOnlineGainTableID = 4;
+    else if (tableName.CompareTo("Krypton_2011-03")               == 0)
+      fOnlineGainTableID = 5;
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2012-01") == 0)
+      fOnlineGainTableID = 6;
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2012-01") == 0)
+      fOnlineGainTableID = 7;
+    else if (tableName.CompareTo("Krypton_2012-01")               == 0)
+      fOnlineGainTableID = 8;
+    else
+      AliFatal(Form("unknown gaintable <%s> requested", tableName.Data()));
+
+    AliInfo(Form("looking for gaintable: %s (id %i)",
+                tableName.Data(), fOnlineGainTableID));
+
+    if (fOnlineGainTableID > 0)
+      return fOnlineGainTableID;
+  } 
+  else {
+
+    AliError("NO DCS/DCSv2 OCDB entry found!");
+    return -1;
+
+  }
+
+  return -1;
+
+}
+
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetGlobalConfiguration(TString &config)
+{
+  //
+  // Get Configuration from the DCS
+  //
+
+  const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+  if(!dcsArr){
+    AliError("No DCS CDB Object available!");
+    config = "";
+    return;
+  }
+
+  Int_t idSOR = 0, idEOR=1; // The index of SOR and EOR
+  Bool_t hasSOR = (dcsArr->At(idSOR));
+  Bool_t hasEOR = (dcsArr->At(idEOR));
+  TString cfgSOR = "", cfgEOR = ""; // The configuration at SOR/EOR
+
+  // The SOR object is mandatory
+  if (!hasSOR) {
+    AliError("NO SOR object found in CDB file!");
+    config = "";
+    return;
+  }
+  if (!hasEOR) AliWarning("NO EOR object found in CDB file!");
+
+  // Check CalDCS version
+  Int_t calver = 0;
+  if (!strcmp(dcsArr->At(idSOR)->ClassName(),"AliTRDCalDCS")) calver = 1;
+  else if (!strcmp(dcsArr->At(idSOR)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+
+  // Get the configuration strings
+  if (calver == 1) {
+    // DCS object
+    const AliTRDCalDCS *calSOR = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(idSOR));
+    cfgSOR = calSOR->GetGlobalConfigName();
+    if (hasEOR) {
+      const AliTRDCalDCS *calEOR = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(idEOR));
+      cfgEOR = calEOR->GetGlobalConfigName();
+    }
+  } 
+  else if (calver == 2) {
+    // DCSv2 object
+    const AliTRDCalDCSv2 *calv2SOR = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(idSOR));
+    cfgSOR = calv2SOR->GetGlobalConfigName();
+    if (hasEOR) {
+      const AliTRDCalDCSv2 *calv2EOR = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(idEOR));
+      cfgEOR = calv2EOR->GetGlobalConfigName();
+    }
+  } 
+  else {
+    AliError("NO DCS/DCSv2 OCDB entry found!");
+    config = "";
+    return;
+  }
+
+  // If there is no EOR entry, return the SOR value
+  if (!hasEOR || cfgEOR.Length()==0) {
+    config = cfgSOR;
+    return;
+  }
+
+  // Check if the configuration is the same for both
+  if (cfgSOR.EqualTo(cfgEOR)) {
+    config = cfgSOR;
+    return;
+  }
+
+  // When both SOR and EOR have an entry but are different, the config is not defined
+  AliError("Inconsistent configuration at start and end of run found!");
+  config = "";
+  return;
+
+}
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetGlobalConfigurationByChamber(TString &config,Int_t par, Int_t opt)
+{
+  //
+  // Get Configuration from the DCS
+  //
+
+  // par is the enumeration from kFltrSet = 1 ... kAddOpti 6
+  //   example:
+  //      cf_p_nozs_tb30_csmtrk_ptrg
+  //   par   1   2    3    4     5
+  //      kFltrSet  kTimebin   kTrigSet
+  //          kReadout  kTrkMode
+  // opt is 0 for the value itself and 1,2,3 for the first, second option
+  // opt has standard value of 0 in header file
+
+
+  // The point is that we extract the parameter/option for 
+  // each of the 540 chambers and compare only this 
+  // parameter/option for a global value, not the full
+  // configuration name with all parameters
+
+  // Get the array with SOR and EOR
+  const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+  if(!dcsArr){
+    AliError("No DCS CDB Object available!");
+    config = "";
+    return;
+  }
+
+  // Check CalDCS version
+  Int_t calver = 0;
+  if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
+  else if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+  //
+  // Get the configuration strings
+  //
+  // CalDCS version 1
+  if (calver == 1) {
+
+    // Loop over SOR/EOR
+    for(Int_t iSEOR=0;iSEOR<2;iSEOR++){
+
+      // DCS object
+      AliTRDCalDCS *cal = dynamic_cast<AliTRDCalDCS *>(dcsArr->At(iSEOR));
+      if(!cal){
+       // SOR mandantory
+       if(iSEOR==0){
+         AliError("NO SOR object found in CDB file!");
+         config = "";
+         return;
+       }
+       else{
+         AliWarning("NO EOR object found in CDB file!");
+         continue;
+       }
+      } // Check DCS object is there
+      
+      // Loop over 540 chambers
+      {
+       Int_t iDet=0;
+       // Find the first chamber in SOR
+       if(iSEOR==0){
+         // Loop until we find the first chamber
+         for(;iDet<kNdet;iDet++){
+           const AliTRDCalDCSFEE *DCSFEEObj = cal->GetCalDCSFEEObj(iDet);
+           // Check it's an installed chamber that responded
+           if(DCSFEEObj && !DCSFEEObj->GetStatusBit()) {
+
+             // We ask for a valid configuration name starting with cf_
+             if(DCSFEEObj->GetConfigName().BeginsWith("cf_")){
+               // Set the parameter of the first good ROC as global,
+               // we take only the requested part!
+               GetDCSConfigParOption(DCSFEEObj->GetConfigName(), par, opt, config);
+               // Increase counter to not look at this chamber again
+               // and break loop
+               iDet++;
+               break;
+             } // End: valid config name
+           } // End: first good chamber
+         } // End: find the first chamber
+       } // End: is SOR
+
+       // Now compare the other chambers with the first one
+       for(;iDet<kNdet;iDet++){
+         const AliTRDCalDCSFEE *DCSFEEObj = cal->GetCalDCSFEEObj(iDet);
+         // Check it's an installed chamber that responded
+         if(DCSFEEObj && !DCSFEEObj->GetStatusBit()) {
+           
+           // We ask for a valid configuration name starting with cf_
+           if(DCSFEEObj->GetConfigName().BeginsWith("cf_")){
+
+             // Get the parameter/option of this chamber
+             TString tmpcfg;
+             GetDCSConfigParOption(DCSFEEObj->GetConfigName(), par, opt, tmpcfg);
+             // Compare with the global value
+             if(config.CompareTo(tmpcfg)){
+               // Two cases mixed or changed during run
+               if(iSEOR==0){
+                 AliError("Mixed DCS configuration for different chambers!");
+                 config="mixed";
+                 return;
+               } else {
+                 AliError("Inconsistent configuration at start and end of run found!");
+                 config = "";
+                 return;
+               }
+             }
+           }// Valid config name 
+         }// Good chamber
+       } // Second half loop 
+      } // Loop over 540 chambers
+    } // Loop over SOR / EOR 
+  } // End calver 1
+  //
+  // CalDCS version 2
+  //
+  else if (calver == 2) {
+
+    // Loop over SOR/EOR
+    for(Int_t iSEOR=0;iSEOR<2;iSEOR++){
+
+      // DCS object
+      AliTRDCalDCSv2 *cal = dynamic_cast<AliTRDCalDCSv2 *>(dcsArr->At(iSEOR));
+      if(!cal){
+       // SOR mandantory
+       if(iSEOR==0){
+         AliError("NO SOR object found in CDB file!");
+         config = "";
+         return;
+       }
+       else{
+         AliWarning("NO EOR object found in CDB file!");
+         continue;
+       }
+      } // Check DCS object is there
+      
+      // Loop over 540 chambers
+      {
+       Int_t iDet=0;
+       // Find the first chamber in SOR
+       if(iSEOR==0){
+         // Loop until we find the first chamber
+         for(;iDet<kNdet;iDet++){
+           const AliTRDCalDCSFEEv2 *DCSFEEObj = cal->GetCalDCSFEEObj(iDet);
+           // Check it's an installed chamber that responded
+           if(DCSFEEObj && !DCSFEEObj->GetStatusBit()) {
+             // Check for a valid config name
+             if(DCSFEEObj->GetConfigName().BeginsWith("cf_")){
+
+               // Set the parameter of the first good ROC as global,
+               // we take only the requested part!
+               GetDCSConfigParOption(DCSFEEObj->GetConfigName(), par, opt, config);
+               // Increase counter to not look at this chamber again
+               // and break loop
+               iDet++;
+               break;
+             } // End: valid config name
+           } // End: first good chamber
+         } // End: find the first chamber
+       } // End: is SOR
+
+       // Now compare the other chambers with the first one
+       for(;iDet<kNdet;iDet++){
+         const AliTRDCalDCSFEEv2 *DCSFEEObj = cal->GetCalDCSFEEObj(iDet);
+         // Check it's an installed chamber that responded
+         if(DCSFEEObj && !DCSFEEObj->GetStatusBit()) {
+
+           // Check for a valid config name
+           if(DCSFEEObj->GetConfigName().BeginsWith("cf_")){
+
+             // Get the parameter/option of this chamber
+             TString tmpcfg;
+             GetDCSConfigParOption(DCSFEEObj->GetConfigName(), par, opt, tmpcfg);
+             // Compare with the global value
+             if(config.CompareTo(tmpcfg)){
+               // Two cases mixed or changed during run
+               if(iSEOR==0){
+                 AliError("Mixed DCS configuration for different chambers!");
+                 config="mixed";
+                 return;
+               } else {
+                 AliError("Inconsistent configuration at start and end of run found!");
+                 config = "";
+                 return;
+               }
+             }
+           } // End: valid config name
+         } // End: Good chamber
+       } // End: Second half loop 
+      } // Loop over 540 chambers
+    } // Loop over SOR / EOR 
+  } // End calver 2
+  else {
+    AliError("NO DCS/DCSv2 OCDB entry found!");
+    config = "";
+    return;
+  }
+
+}
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetGlobalConfigurationVersion(TString &version)
+{
+  //
+  // Get Version of Configuration from the DCS
+  //
+
+  const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+  if(!dcsArr){
+    version = "";
+    return;
+  }
+
+  Int_t esor   = 0; // Take SOR
+  Int_t calver = 0; // Check CalDCS version
+  if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS"))   calver = 1;
+  if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+
+  if      (calver == 1) {
+
+    // DCS object
+    const AliTRDCalDCS   *calDCS   = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(esor));
+    if(!calDCS){
+      version = "";
+      return;
+    } 
+    version = calDCS->GetGlobalConfigVersion();
+
+  } 
+  else if (calver == 2) {
+
+    // DCSv2 object
+    const AliTRDCalDCSv2 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
+    if(!calDCSv2){
+      version = "";
+      return;
+    } 
+    version = calDCSv2->GetGlobalConfigVersion();
+
+  } 
+  else {
+
+    AliError("NO DCS/DCSv2 OCDB entry found!");
+
+  }
+
+}
+
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfParsDCS(TString cname, Char_t delimiter)
+{
+  // Get the number of configuration parameters from the DCS config
+
+  //AliInfo(Form("\"%s\" tokenized by \"%c\"", cname.Data(), delimiter));
+  if(!cname.Length()) return -1;  // -1 for the "cf"
+  Int_t nconf(0);
+  for(Int_t ich(1); ich<cname.Length()-1; ich++){ if(cname[ich]==delimiter) nconf++;}
+  return nconf;
+}
+
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfOptsDCS(TString cname, Int_t cfgType)
+{
+  // Get the number of options of a given configuration parameter from DCS
+
+  Char_t cdelim = '_', // define the delimiters
+         odelim = '-';
+  Int_t nconfig = GetNumberOfParsDCS(cname, cdelim);
+
+  // protect
+  if ((nconfig == -1) || (nconfig < cfgType)) {
+    AliError("Not enough parameters in DCS configuration name!");
+    return 0;
+  }
+
+  TObjArray *carr = cname.Tokenize(cdelim);
+  Int_t nopt = GetNumberOfParsDCS(((TObjString*)carr->At(cfgType))->GetString(), odelim);
+  carr->Delete(); delete carr;
+  return nopt;
+}
+
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetDCSConfigParOption(TString cname, Int_t cfgType, Int_t option, TString &cfgo)
+{
+  //
+  // Get a configuration (see enum in header file) or the options of a configuration
+  // option == 0 returns the configuration itself
+  // option >  0 returns the optional parameter Nr. (option) of the configuration (cfgType)
+  //
+
+  Char_t cdelim = '_', // define the delimiters
+         odelim = '-';
+
+  Int_t nconfig = GetNumberOfParsDCS(cname, cdelim);
+  // protect
+  if (nconfig == -1) {
+    AliError("DCS configuration name empty!");
+    cfgo = "";
+    return;
+  } else if (nconfig < cfgType) {
+    AliError(Form("Not enough parameters in DCS configuration name!"
+                 " Name %s",cname.Data()));
+    cfgo = "";
+    return;
+  }
+
+  TObjArray *carr = cname.Tokenize(cdelim);
+  TString cfgString(((TObjString*)carr->At(cfgType))->GetString());
+  Int_t noptions = GetNumberOfParsDCS(cfgString, odelim);
+  // protect
+  if (noptions < option) {
+    AliError(Form("Not enough options in DCS configuration name!"
+                 " Name %s",cname.Data()));
+    cfgo = "";
+    carr->Delete(); delete carr;
+    return;
+  }
+  TObjArray *oarr = cfgString.Tokenize(odelim);
+  cfgo = ((TObjString*)oarr->At(option))->GetString();
+  carr->Delete(); delete carr;
+  oarr->Delete(); delete oarr;
+  return;
+
+}
+
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::HasOnlineFilterPedestal()
+{
+  //
+  // Checks whether pedestal filter was applied online
+  //
+
+  TString filterconfig;
+  GetFilterType(filterconfig);
+
+  return filterconfig.Contains("p");
+
+}
+
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::HasOnlineFilterGain()
+{
+  //
+  // Checks whether online gain filter was applied
+  //
+
+  TString filterconfig;
+  GetFilterType(filterconfig);
+
+  return filterconfig.Contains("g");
+
+}
+
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::HasOnlineTailCancellation()
+{
+  //
+  // Checks whether online tail cancellation was applied
+  //
+
+  TString filterconfig;
+  GetFilterType(filterconfig);
+
+  return filterconfig.Contains("t");
 
 }
 
@@ -813,6 +1640,23 @@ Char_t AliTRDcalibDB::GetChamberStatus(Int_t det)
 
 }
 
+//_____________________________________________________________________________
+AliTRDrecoParam* AliTRDcalibDB::GetRecoParam(Int_t */*eventtype*/)
+{
+  //
+  // Returns the TRD reconstruction parameters from the OCDB
+  //
+
+  const TClonesArray *recos = dynamic_cast<const TClonesArray*>(GetCachedCDBObject(kIDRecoParam));
+  if (!recos) return 0x0;
+
+  // calculate entry based on event type info
+  Int_t n = 0; //f(eventtype[0], eventtype[1], ....)
+
+  return (AliTRDrecoParam *) recos->UncheckedAt(n);
+
+}
+
 //_____________________________________________________________________________
 Bool_t AliTRDcalibDB::IsPadMasked(Int_t det, Int_t col, Int_t row)
 {
@@ -870,6 +1714,7 @@ Bool_t AliTRDcalibDB::IsPadNotConnected(Int_t det, Int_t col, Int_t row)
   //
   // Returns status, see name of functions for details ;-)
   //
+
   const AliTRDCalPadStatus         * cal = dynamic_cast<const AliTRDCalPadStatus *> 
                                            (GetCachedCDBObject(kIDPadStatus));
   if (!cal) {
@@ -881,7 +1726,7 @@ Bool_t AliTRDcalibDB::IsPadNotConnected(Int_t det, Int_t col, Int_t row)
 }
 
 //_____________________________________________________________________________
-Bool_t AliTRDcalibDB::IsChamberInstalled(Int_t det)
+Bool_t AliTRDcalibDB::IsChamberGood(Int_t det)
 {
   //
   // Returns status, see name of functions for details ;-)
@@ -893,12 +1738,12 @@ Bool_t AliTRDcalibDB::IsChamberInstalled(Int_t det)
     return -1;
   }
 
-  return cal->IsInstalled(det);
+  return cal->IsGood(det);
 
 }
 
 //_____________________________________________________________________________
-Bool_t AliTRDcalibDB::IsChamberMasked(Int_t det)
+Bool_t AliTRDcalibDB::IsChamberNoData(Int_t det)
 {
   //
   // Returns status, see name of functions for details ;-)
@@ -910,84 +1755,148 @@ Bool_t AliTRDcalibDB::IsChamberMasked(Int_t det)
     return -1;
   }
 
-  return cal->IsMasked(det);
+  return cal->IsNoData(det);
 
 }
 
 //_____________________________________________________________________________
-const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(AliTRDrecoParam::AliTRDpidMethod method)
+Bool_t AliTRDcalibDB::IsHalfChamberNoData(Int_t det, Int_t side)
 {
   //
-  // Returns the object storing the distributions for PID with likelihood
+  // Returns status, see name of functions for details ;-)
   //
 
-  switch(method) {
-  case AliTRDrecoParam::kLQPID: 
-    return dynamic_cast<const AliTRDCalPID *>(GetCachedCDBObject(kIDPIDLQ));
-  case AliTRDrecoParam::kNNPID: 
-    return dynamic_cast<const AliTRDCalPID *>(GetCachedCDBObject(kIDPIDNN));
+  const AliTRDCalChamberStatus     * cal = dynamic_cast<const AliTRDCalChamberStatus *> 
+                                           (GetCachedCDBObject(kIDChamberStatus));
+  if (!cal) {
+    return -1;
   }
 
-  return 0x0;
+  return side > 0 ? cal->IsNoDataSideB(det) : cal->IsNoDataSideA(det);
 
 }
 
 //_____________________________________________________________________________
-const AliTRDCalMonitoring *AliTRDcalibDB::GetMonitoringObject()
+Bool_t AliTRDcalibDB::IsChamberBadCalibrated(Int_t det)
 {
   //
-  // Returns the object storing the monitoring data
+  // Returns status, see name of functions for details ;-)
   //
 
-  return dynamic_cast<const AliTRDCalMonitoring *> 
-         (GetCachedCDBObject(kIDMonitoringData));
-   
+  const AliTRDCalChamberStatus     * cal = dynamic_cast<const AliTRDCalChamberStatus *> 
+                                           (GetCachedCDBObject(kIDChamberStatus));
+  if (!cal) {
+    return -1;
+  }
+
+  return cal->IsBadCalibrated(det);
+
 }
 
 //_____________________________________________________________________________
-Float_t AliTRDcalibDB::GetOmegaTau(Float_t vdrift, Float_t bz)
+Bool_t AliTRDcalibDB::IsChamberNotCalibrated(Int_t det)
 {
   //
-  // Returns omega*tau (tan(Lorentz-angle)) for a given drift velocity <vdrift> 
-  // and a B-field <bz> for Xe/CO2 (15%).
-  // The values are according to a GARFIELD simulation.
-  //
-  // This function basically does not belong to the calibration class.
-  // It should be moved somewhere else. 
-  // However, currently it is in use by simulation and reconstruction.
+  // Returns status, see name of functions for details ;-)
   //
-  
-  Float_t fieldAbs = TMath::Abs(bz);
-  Float_t fieldSgn = (bz > 0.0) ? 1.0 : -1.0;
 
-  const Int_t kNb = 5;
-  Float_t p0[kNb] = {  0.004810,  0.007412,  0.010252,  0.013409,  0.016888 };
-  Float_t p1[kNb] = {  0.054875,  0.081534,  0.107333,  0.131983,  0.155455 };
-  Float_t p2[kNb] = { -0.008682, -0.012896, -0.016987, -0.020880, -0.024623 };
-  Float_t p3[kNb] = {  0.000155,  0.000238,  0.000330,  0.000428,  0.000541 };
+  const AliTRDCalChamberStatus     * cal = dynamic_cast<const AliTRDCalChamberStatus *> 
+                                           (GetCachedCDBObject(kIDChamberStatus));
+  if (!cal) {
+    return -1;
+  }
+
+  return cal->IsNotCalibrated(det);
 
-  // No ExB if field is too small (or zero)
-  if (fieldAbs < 0.01) {
+}
 
-    return 0.0;
+//_____________________________________________________________________________
+const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(AliTRDpidUtil::ETRDPIDMethod method)
+{
+  //
+  // Returns the object storing the distributions for PID with likelihood
+  //
 
+  switch(method) {
+  case AliTRDpidUtil::kLQ: 
+    return dynamic_cast<const AliTRDCalPID *>(GetCachedCDBObject(kIDPIDLQ));
+  case AliTRDpidUtil::kNN: 
+    return dynamic_cast<const AliTRDCalPID *>(GetCachedCDBObject(kIDPIDNN));
+  case AliTRDpidUtil::kESD:
+    return 0x0; // To avoid compiler warnings 
   }
-  // Calculate ExB from parametrization
-  else {
 
-    Int_t ib = ((Int_t) (10 * (fieldAbs - 0.15)));
-    ib       = TMath::Max(  0,ib);
-    ib       = TMath::Min(kNb,ib);
+  return 0x0;
 
-    Float_t alphaL = p0[ib] 
-                   + p1[ib] * vdrift
-                   + p2[ib] * vdrift*vdrift
-                   + p3[ib] * vdrift*vdrift*vdrift;
+}
 
-    return TMath::Tan(fieldSgn * alphaL);
+//_____________________________________________________________________________
+AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMethod /*method*/)
+{
+  //
+  // Returns the PID response object for 1D-LQ
+  //
+
+  if (!fPIDResponse) {
+    AliDebug(2, "Setting new PID response.");
+
+    fPIDResponse = new AliTRDPIDResponse;
+
+    // Load Reference Histos from OCDB
+//    if(method == AliTRDPIDResponse::kLQ1D){
+    //fPIDResponse->SetPIDmethod(AliTRDPIDResponse::kLQ1D);
+    const TObjArray *references = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDPIDLQ1D));
+
+    TIter refs(references);
+    TObject *obj = NULL;
+    AliTRDPIDReference *ref = NULL;
+    Bool_t hasReference = kFALSE;
+    while ((obj = refs())){
+      if ((ref = dynamic_cast<AliTRDPIDReference *>(obj))){
+        AliDebug(2, "Setting new PID response object.");
+        TDirectory *bkpdir = gDirectory;
+        gROOT->cd();
+        AliTRDPIDResponseObject *ro = new AliTRDPIDResponseObject;
+        ro->SetPIDReference(ref);
+        fPIDResponse->SetPIDResponseObject(ro);
+        hasReference = kTRUE;
+        gDirectory = bkpdir;
+        break;
+      }
+    }
 
+    if (!hasReference) {
+      AliError("Reference histograms not found in the OCDB");
+    }
   }
 
+//  }
+
+  return fPIDResponse;
+
+}
+
+//_____________________________________________________________________________
+const AliTRDCalTrkAttach* AliTRDcalibDB::GetAttachObject()
+{
+  //
+  // Returns the object storing likelihood distributions for cluster to track attachment
+  //
+
+  return dynamic_cast<const AliTRDCalTrkAttach*>(GetCachedCDBObject(kIDAttach));
+
+}
+
+//_____________________________________________________________________________
+const AliTRDCalMonitoring *AliTRDcalibDB::GetMonitoringObject()
+{
+  //
+  // Returns the object storing the monitoring data
+  //
+
+  return dynamic_cast<const AliTRDCalMonitoring *> 
+         (GetCachedCDBObject(kIDMonitoringData));
+   
 }
 
 //_____________________________________________________________________________
@@ -1149,7 +2058,7 @@ void AliTRDcalibDB::SamplePRF()
 
 //_____________________________________________________________________________
 Int_t AliTRDcalibDB::PadResponse(Double_t signal, Double_t dist
-                                , Int_t layer, Double_t *pad) const
+                               , Int_t layer, Double_t *pad) const
 {
   //
   // Applies the pad response
@@ -1157,7 +2066,7 @@ Int_t AliTRDcalibDB::PadResponse(Double_t signal, Double_t dist
   // something dependent on calibration values
   //
 
-  Int_t iBin  = ((Int_t) (( - dist - fPRFlo) / fPRFwid));
+  Int_t iBin  = ((Int_t) ((-dist - fPRFlo) / fPRFwid));
   Int_t iOff  = layer * fPRFbin;
 
   Int_t iBin0 = iBin - fPRFpad + iOff;
@@ -1187,3 +2096,102 @@ Int_t AliTRDcalibDB::PadResponse(Double_t signal, Double_t dist
   }
 
 }
+
+
+AliTRDtrapConfig* AliTRDcalibDB::GetTrapConfig()
+{
+  // return an existing TRAPconfig or load it from the OCDB
+  // in case of failure, a default TRAPconfig is created
+
+  if (fTrapConfig)
+    return fTrapConfig;
+  else {
+    if ((fTrapConfigName.Length() <= 0) || (fTrapConfigVersion.Length() <= 0)) {
+      // query the configuration to be used
+      this->GetGlobalConfiguration(fTrapConfigName);
+      this->GetGlobalConfigurationVersion(fTrapConfigVersion);
+    }
+
+    // try to load the requested configuration
+    this->LoadTrapConfig(fTrapConfigName, fTrapConfigVersion);
+
+    // if we still don't have a valid TRAPconfig, create a default one
+    if (!fTrapConfig) {
+      AliWarning("Falling back to default configuration");
+      static AliTRDtrapConfig trapConfigDefault("default", "default TRAP configuration");
+      fTrapConfig = &trapConfigDefault;
+      AliTRDtrapConfigHandler cfgHandler(fTrapConfig);
+      cfgHandler.Init();
+      cfgHandler.LoadConfig();
+    }
+
+    AliInfo(Form("using TRAPconfig \"%s\"", fTrapConfig->GetTitle()));
+
+    // we still have to load the gain tables
+    // if the gain filter is active
+    if (HasOnlineFilterGain()) {
+      const Int_t nDets = AliTRDgeometry::Ndet();
+      const Int_t nMcms = AliTRDgeometry::MCMmax();
+      const Int_t nChs  = AliTRDgeometry::ADCmax();
+
+      // gain factors are per MCM
+      // allocate the registers accordingly
+      for (Int_t ch = 0; ch < nChs; ++ch) {
+       AliTRDtrapConfig::TrapReg_t regFGAN = (AliTRDtrapConfig::TrapReg_t) (AliTRDtrapConfig::kFGA0 + ch);
+       AliTRDtrapConfig::TrapReg_t regFGFN = (AliTRDtrapConfig::TrapReg_t) (AliTRDtrapConfig::kFGF0 + ch);
+       fTrapConfig->SetTrapRegAlloc(regFGAN, AliTRDtrapConfig::kAllocByMCM);
+       fTrapConfig->SetTrapRegAlloc(regFGFN, AliTRDtrapConfig::kAllocByMCM);
+      }
+
+      for (Int_t iDet = 0; iDet < nDets; ++iDet) {
+       AliTRDCalOnlineGainTableROC *gainTbl = GetOnlineGainTableROC(iDet);
+       if (gainTbl) {
+         const Int_t nRobs = AliTRDgeometry::GetStack(iDet) == 2 ?
+           AliTRDgeometry::ROBmaxC0() : AliTRDgeometry::ROBmaxC1();
+         for (Int_t rob = 0; rob < nRobs; ++rob) {
+           for (Int_t mcm = 0; mcm < nMcms; ++mcm) {
+             AliTRDCalOnlineGainTableMCM *gainTblMCM = gainTbl->GetGainTableMCM(rob, mcm);
+
+             // set ADC reference voltage
+             fTrapConfig->SetTrapReg(AliTRDtrapConfig::kADCDAC, gainTblMCM->GetAdcdac(), iDet, rob, mcm);
+
+             // set constants channel-wise
+             for (Int_t ch = 0; ch < nChs; ++ch) {
+               AliTRDtrapConfig::TrapReg_t regFGAN = (AliTRDtrapConfig::TrapReg_t) (AliTRDtrapConfig::kFGA0 + ch);
+               AliTRDtrapConfig::TrapReg_t regFGFN = (AliTRDtrapConfig::TrapReg_t) (AliTRDtrapConfig::kFGF0 + ch);
+               fTrapConfig->SetTrapReg(regFGAN, gainTblMCM->GetFGAN(ch), iDet, rob, mcm);
+               fTrapConfig->SetTrapReg(regFGFN, gainTblMCM->GetFGFN(ch), iDet, rob, mcm);
+             }
+           }
+         }
+       }
+      }
+    }
+
+    return fTrapConfig;
+  }
+}
+
+
+AliTRDtrapConfig* AliTRDcalibDB::LoadTrapConfig(const TString &name, const TString &version)
+{
+  // try to load the specified configuration from the OCDB
+  // if it fails, or it does not exist, return null
+
+  AliInfo(Form("looking for TRAPconfig \"%s.%s\"", name.Data(), version.Data()));
+
+  const AliTRDCalTrapConfig *caltrap = dynamic_cast<const AliTRDCalTrapConfig*> (GetCachedCDBObject(kIDTrapConfig));
+
+  if (caltrap) {
+    TString configName(name);
+    configName.Append(".");
+    configName.Append(version);
+    fTrapConfig = caltrap->Get(configName);
+  }
+  else {
+    fTrapConfig = 0x0;
+    AliError("No TRAPconfig entry found!");
+  }
+
+  return fTrapConfig;
+}