]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDcalibDB.cxx
Patch to load 1DLQ PID references
[u/mrichter/AliRoot.git] / TRD / AliTRDcalibDB.cxx
index 32ea602942a73466760bd78a9b1452e091ff1339..8571a2784845ac233b0a5718803b1378679fc3c0 100644 (file)
@@ -36,6 +36,7 @@
 #include "AliLog.h"
 
 #include "AliTRDPIDReference.h"
+#include "AliTRDPIDResponseObject.h"
 #include "AliTRDcalibDB.h"
 #include "AliTRDtrapConfig.h"
 #include "AliTRDtrapConfigHandler.h"
@@ -279,7 +280,11 @@ const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
           break;
         case 7:
          // Online gain table ID 7
-          return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Gaintbl_Uniform_FGAN8_2012-01"); 
+          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;
       }
       break;
@@ -902,81 +907,49 @@ Float_t AliTRDcalibDB::GetPRFWidth(Int_t det, Int_t col, Int_t row)
 }
   
 //_____________________________________________________________________________
-Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
+Int_t AliTRDcalibDB::ExtractTimeBinsFromString(TString tbstr)
 {
-  //
-  // Returns Number of time bins from the DCS
-  //
+  // default value - has not been set
+  Int_t nUndef = -1;
 
-  Int_t nMixed = -2; // not the same number for all chambers
-  Int_t nUndef = -1; // default value - has not been set!
-  Int_t nTbSor = nUndef;
-  Int_t nTbEor = nUndef;
-  Int_t calver = 0; // Check CalDCS version
-
-  const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
-  if (!dcsArr) {
-    AliError("No DCS object found!");
+  // Check if there is any content in the string first
+  if (tbstr.Length() == 0) {
+    AliError("Parameter for number of timebins is empty!");
     return nUndef;
   }
 
-  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 *calDCSsor = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(0));
-    const AliTRDCalDCS *calDCSeor = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(1));
-    if (!calDCSsor) {
-      // the SOR file is mandatory
-      AliError("NO SOR AliTRDCalDCS object found in CDB file!");
-      return nUndef;
-    }
-    if (!calDCSeor) {
-      // this can happen if the run is shorter than a couple of seconds.
-      AliWarning("NO EOR AliTRDCalDCS object found in CDB file.");
-    }
-
-    // get the numbers
-    nTbSor = calDCSsor->GetGlobalNumberOfTimeBins();
-    if (calDCSeor) nTbEor = calDCSeor->GetGlobalNumberOfTimeBins();
-
-  } else if (calver == 2) {
-    // DCSv2 object
-    const AliTRDCalDCSv2 *calDCSsorv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(0));
-    const AliTRDCalDCSv2 *calDCSeorv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(1));
-    if (!calDCSsorv2) {
-      // the SOR file is mandatory
-      AliError("NO SOR AliTRDCalDCSv2 object found in CDB file!");
-      return nUndef;
-    }
-    if (!calDCSeorv2) {
-      // this can happen if the run is shorter than a couple of seconds.
-      AliWarning("NO EOR AliTRDCalDCSv2 object found in CDB file.");
-    }
+  // 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 nUndef;
+  }
 
-    // get the numbers
-    nTbSor = calDCSsorv2->GetGlobalNumberOfTimeBins();
-    if (calDCSeorv2) nTbEor = calDCSeorv2->GetGlobalNumberOfTimeBins();
+  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 nUndef;
+  }
 
-  } else AliError("NO DCS/DCSv2 OCDB entry found!");
+  return tbstr.Atoi();
 
-  // if they're the same return the value
-  // -2 means mixed, -1: no data, >= 0: good number of time bins
-  if (nTbSor == nTbEor) return nTbSor;
+}
 
-  // if they're differing:
-  if (nTbSor == nMixed || nTbEor == nMixed) {
-    AliWarning("Inconsistent number of time bins found!");
-    return nMixed;
-  }
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
+{
+  //
+  // Returns number of time bins from the DCS
+  //
 
-  // one is undefined, the other ok -> return that one
-  if (nTbSor == nUndef) return nTbEor;
-  if (nTbEor == nUndef) return nTbSor;
+  // Get the corresponding parameter
+  TString cfgstr = "", cfgname = "";
+  GetGlobalConfiguration(cfgname);
+  GetDCSConfigParOption(cfgname, kTimebin, 0, cfgstr);
 
-  // only remains: two different numbers >= 0
-  return nMixed;
+  return ExtractTimeBinsFromString(cfgstr);
 
 }
 
@@ -987,7 +960,9 @@ void AliTRDcalibDB::GetFilterType(TString &filterType)
   // Returns the filter type
   //
 
-  GetDCSConfigParOption(kFltrSet, 0, filterType);
+  TString cfgname = "";
+  GetGlobalConfiguration(cfgname);
+  GetDCSConfigParOption(cfgname, kFltrSet, 0, filterType);
 
 }
 
@@ -1062,6 +1037,10 @@ Int_t AliTRDcalibDB::GetOnlineGainTableID()
       fOnlineGainTableID = 7;
       return fOnlineGainTableID;
     }
+    if (tableName.CompareTo("Krypton_2011-03")               == 0) {
+      fOnlineGainTableID = 8;
+      return fOnlineGainTableID;
+    }
 
   } 
   else {
@@ -1084,43 +1063,71 @@ void AliTRDcalibDB::GetGlobalConfiguration(TString &config)
 
   const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
   if(!dcsArr){
+    AliError("No DCS CDB Object available!");
     config = "";
     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;
+  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
 
-  if      (calver == 1) {
+  // 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!");
 
-    // DCS object
-    const AliTRDCalDCS   *calDCS   = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(esor));
-    if(!calDCS){
-      config = "";
-      return;
-    } 
-    config = calDCS->GetGlobalConfigName();
+  // 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 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
-    if(!calDCSv2){
-      config = "";
-      return;
-    } 
-    config = calDCSv2->GetGlobalConfigName();
-
+    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;
+
 }
 
 //_____________________________________________________________________________
@@ -1172,7 +1179,41 @@ void AliTRDcalibDB::GetGlobalConfigurationVersion(TString &version)
 }
 
 //_____________________________________________________________________________
-void AliTRDcalibDB::GetDCSConfigParOption(Int_t cfgType, Int_t option, TString &cfgo)
+Int_t AliTRDcalibDB::GetNumberOfParsDCS(TString cname)
+{
+  // Get the number of configuration parameters from the DCS config
+
+  TString cdelim = "_";
+  TObjArray *carr = cname.Tokenize(cdelim);
+  return carr->GetEntries() - 1; // -1 for the "cf"
+
+}
+
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfOptsDCS(TString cname, Int_t cfgType)
+{
+  // Get the number of options of a given configuration parameter from DCS
+
+  TString cdelim = "_";
+  TString odelim = "-";
+
+  TObjArray *carr = cname.Tokenize(cdelim);
+  Int_t nconfig = carr->GetEntries();
+
+  // protect
+  if ((nconfig == 0) || ((nconfig-1) < cfgType)) {
+    AliError("Not enough parameters in DCS configuration name!");
+    return 0;
+  }
+
+  TString fullcfg = ((TObjString*)carr->At(cfgType))->GetString();
+  TObjArray *oarr = fullcfg.Tokenize(odelim);
+  return oarr->GetEntries() -1; // -1 for the parameter name
+
+}
+
+//_____________________________________________________________________________
+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
@@ -1184,52 +1225,33 @@ void AliTRDcalibDB::GetDCSConfigParOption(Int_t cfgType, Int_t option, TString &
   TString cdelim = "_";
   TString odelim = "-";
 
-  // get the full configuration name
-  TString cname;
-  GetGlobalConfiguration(cname);
   TObjArray *carr = cname.Tokenize(cdelim);
   Int_t nconfig = carr->GetEntries();
 
   // protect
   if (nconfig == 0) {
-    AliError("Bad DCS configuration name!");
+    AliError("DCS configuration name empty!");
     cfgo = "";
     return;
   } else if ((nconfig-1) < cfgType) {
-    AliError("Not enough DCS configuration parameters!");
+    AliError("Not enough parameters in DCS configuration name!");
     cfgo = "";
     return;
   }
 
   TString fullcfg = ((TObjString*)carr->At(cfgType))->GetString();
+  TObjArray *oarr = fullcfg.Tokenize(odelim);
+  Int_t noptions = oarr->GetEntries();
 
-  if (fullcfg.Contains(odelim)) {
-
-    TObjArray *oarr = fullcfg.Tokenize(odelim);
-    Int_t noptions = oarr->GetEntries();
-
-    // protect
-    if ((noptions-1) < option) {
-      AliError("Not enough DCS configuration options defined!");
-      cfgo = "";
-      return;
-    }
-
-    cfgo = ((TObjString*)oarr->At(option))->GetString();
+  // protect
+  if ((noptions-1) < option) {
+    AliError("Not enough options in DCS configuration name!");
+    cfgo = "";
     return;
-
   }
-  else {
 
-    if (option != 0) {
-      AliError("Not enough DCS configuration options defined!");
-      cfgo = "";
-      return;
-    }
-    cfgo = fullcfg;
-    return;
-
-  }
+  cfgo = ((TObjString*)oarr->At(option))->GetString();
+  return;
 
 }
 
@@ -1522,23 +1544,27 @@ AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMetho
     fPIDResponse = new AliTRDPIDResponse;
 
     // Load Reference Histos from OCDB
-    fPIDResponse->SetPIDmethod(method);
-    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))){
-        fPIDResponse->Load(ref);
-        hasReference = kTRUE;
-        break;
+    if(method == AliTRDPIDResponse::kLQ1D){
+      fPIDResponse->SetPIDmethod(method);
+      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))){
+          AliTRDPIDResponseObject *ro = new AliTRDPIDResponseObject;
+          ro->SetPIDReference(ref);
+          fPIDResponse->SetPIDResponseObject(ro);
+          hasReference = kTRUE;
+          break;
+        }
       }
-    }
 
-    if (!hasReference) {
-      AliError("Reference histograms not found in the OCDB");
+      if (!hasReference) {
+        AliError("Reference histograms not found in the OCDB");
+      }
     }
 
   }