]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDcalibDB.cxx
Added some plots
[u/mrichter/AliRoot.git] / TRD / AliTRDcalibDB.cxx
index 2559f83c24cee93371882c095d49981c3ef5ae16..927b4f855a1507740359010e278728e3a4f73e23 100644 (file)
@@ -28,6 +28,7 @@
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
+#include <TROOT.h>
 #include <TClonesArray.h>
 #include <TObjArray.h>
 
 #include "AliLog.h"
 
 #include "AliTRDPIDReference.h"
+#include "AliTRDPIDResponseObject.h"
 #include "AliTRDcalibDB.h"
 #include "AliTRDtrapConfig.h"
 #include "AliTRDtrapConfigHandler.h"
 #include "AliTRDCommonParam.h"
+#include "AliTRDgeometry.h"
 
 #include "Cal/AliTRDCalROC.h"
 #include "Cal/AliTRDCalPad.h"
@@ -279,8 +282,14 @@ 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;
+      default:
+       AliError(Form("unknown gaintable requested with ID"));
       }
       break;
 
@@ -312,7 +321,7 @@ const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
       return CacheCDBEntry(kIDFEE               ,"TRD/Calib/FEE"); 
       break;
     case kIDTrapConfig :
-      return CacheCDBEntry(kIDFEE               ,"TRD/Calib/TrapConfig"); 
+      return CacheCDBEntry(kIDTrapConfig        ,"TRD/Calib/TrapConfig"); 
       break;
     case kIDDCS :
       return CacheCDBEntry(kIDDCS               ,"TRD/Calib/DCS");
@@ -369,6 +378,8 @@ const TObject *AliTRDcalibDB::CacheCDBEntry(Int_t id, const char *cdbPath)
     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()));
     }
   } 
   
@@ -754,7 +765,7 @@ Float_t AliTRDcalibDB::GetOnlineGainFactor(Int_t det, Int_t col, Int_t row)
   //
 
   if (!HasOnlineFilterGain()) {
-    return 0x0;
+    return -1;
   }
   
   const AliTRDCalOnlineGainTable *calOnline 
@@ -902,23 +913,12 @@ 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;
-
-  // Get the corresponding parameter
-  TString tbstr = "";
-  GetDCSConfigParOption(kTimebin, 0, tbstr);
-
   // Check if there is any content in the string first
   if (tbstr.Length() == 0) {
     AliError("Parameter for number of timebins is empty!");
-    return nUndef;
+    return -1;
   }
 
   // Check if we have the correct config parameter
@@ -926,20 +926,36 @@ Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
   TString tbsubstr = tbstr(0,2);
   if (!tbsubstr.EqualTo(tbident)) {
     AliError(Form("Parameter for number of timebins is corrupted (%s)!", tbstr.Data()));
-    return nUndef;
+    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 nUndef;
+    return -1;
   }
 
-  Int_t ntb = tbstr.Atoi();
-  AliInfo(Form("Number of timebins from CDB: %d", ntb));
+  return tbstr.Atoi();
 
-  return ntb;
+}
+
+//_____________________________________________________________________________
+Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
+{
+  //
+  // Returns number of time bins from the DCS
+  //
+
+  // 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);
 
 }
 
@@ -950,7 +966,9 @@ void AliTRDcalibDB::GetFilterType(TString &filterType)
   // Returns the filter type
   //
 
-  GetDCSConfigParOption(kFltrSet, 0, filterType);
+  TString cfgname = "";
+  GetGlobalConfiguration(cfgname);
+  GetDCSConfigParOption(cfgname, kFltrSet, 0, filterType);
 
 }
 
@@ -997,35 +1015,30 @@ Int_t AliTRDcalibDB::GetOnlineGainTableID()
         break;
       }
     }
-    if (tableName.CompareTo("Krypton_2011-01")               == 0) {
+    if (tableName.CompareTo("Krypton_2011-01")               == 0)
       fOnlineGainTableID = 1;
-      return fOnlineGainTableID;
-    }
-    if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2011-01") == 0) {
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2011-01") == 0)
       fOnlineGainTableID = 2;
-      return fOnlineGainTableID;
-    }
-    if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2011-01") == 0) {
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2011-01") == 0)
       fOnlineGainTableID = 3;
-      return fOnlineGainTableID;
-    }
-    if (tableName.CompareTo("Krypton_2011-02")               == 0) {
+    else if (tableName.CompareTo("Krypton_2011-02")               == 0)
       fOnlineGainTableID = 4;
-      return fOnlineGainTableID;
-    }
-    if (tableName.CompareTo("Krypton_2011-03")               == 0) {
+    else if (tableName.CompareTo("Krypton_2011-03")               == 0)
       fOnlineGainTableID = 5;
-      return fOnlineGainTableID;
-    }
-    if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2012-01") == 0) {
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2012-01") == 0)
       fOnlineGainTableID = 6;
-      return fOnlineGainTableID;
-    }
-    if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2012-01") == 0) {
+    else if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2012-01") == 0)
       fOnlineGainTableID = 7;
-      return fOnlineGainTableID;
-    }
+    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 {
 
@@ -1163,7 +1176,40 @@ void AliTRDcalibDB::GetGlobalConfigurationVersion(TString &version)
 }
 
 //_____________________________________________________________________________
-void AliTRDcalibDB::GetDCSConfigParOption(Int_t cfgType, Int_t option, TString &cfgo)
+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
@@ -1171,56 +1217,36 @@ void AliTRDcalibDB::GetDCSConfigParOption(Int_t cfgType, Int_t option, TString &
   // option >  0 returns the optional parameter Nr. (option) of the configuration (cfgType)
   //
 
-  // define the delimiters
-  TString cdelim = "_";
-  TString odelim = "-";
-
-  // get the full configuration name
-  TString cname;
-  GetGlobalConfiguration(cname);
-  TObjArray *carr = cname.Tokenize(cdelim);
-  Int_t nconfig = carr->GetEntries();
+  Char_t cdelim = '_', // define the delimiters
+         odelim = '-';
 
+  Int_t nconfig = GetNumberOfParsDCS(cname, cdelim);
   // protect
-  if (nconfig == 0) {
-    AliError("Bad DCS configuration name!");
+  if (nconfig == -1) {
+    AliError("DCS configuration name empty!");
     cfgo = "";
     return;
-  } else if ((nconfig-1) < cfgType) {
-    AliError("Not enough DCS configuration parameters!");
+  } else if (nconfig < cfgType) {
+    AliError("Not enough parameters in DCS configuration name!");
     cfgo = "";
     return;
   }
 
-  TString fullcfg = ((TObjString*)carr->At(cfgType))->GetString();
-
-  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();
+  TObjArray *carr = cname.Tokenize(cdelim);
+  TString cfgString(((TObjString*)carr->At(cfgType))->GetString());
+  Int_t noptions = GetNumberOfParsDCS(cfgString, odelim);
+  // protect
+  if (noptions < option) {
+    AliError("Not enough options in DCS configuration name!");
+    cfgo = "";
+    carr->Delete(); delete carr;
     return;
-
-  }
-  else {
-
-    if (option != 0) {
-      AliError("Not enough DCS configuration options defined!");
-      cfgo = "";
-      return;
-    }
-    cfgo = fullcfg;
-    return;
-
   }
+  TObjArray *oarr = cfgString.Tokenize(odelim);
+  cfgo = ((TObjString*)oarr->At(option))->GetString();
+  carr->Delete(); delete carr;
+  oarr->Delete(); delete oarr;
+  return;
 
 }
 
@@ -1481,6 +1507,23 @@ Bool_t AliTRDcalibDB::IsChamberBadCalibrated(Int_t det)
 
 }
 
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::IsChamberNotCalibrated(Int_t det)
+{
+  //
+  // Returns status, see name of functions for details ;-)
+  //
+
+  const AliTRDCalChamberStatus     * cal = dynamic_cast<const AliTRDCalChamberStatus *> 
+                                           (GetCachedCDBObject(kIDChamberStatus));
+  if (!cal) {
+    return -1;
+  }
+
+  return cal->IsNotCalibrated(det);
+
+}
+
 //_____________________________________________________________________________
 const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(AliTRDpidUtil::ETRDPIDMethod method)
 {
@@ -1502,18 +1545,20 @@ const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(AliTRDpidUtil::ETRDPIDMethod met
 }
 
 //_____________________________________________________________________________
-AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMethod method)
+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
-    fPIDResponse->SetPIDmethod(method);
+//    if(method == AliTRDPIDResponse::kLQ1D){
+    //fPIDResponse->SetPIDmethod(AliTRDPIDResponse::kLQ1D);
     const TObjArray *references = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDPIDLQ1D));
 
     TIter refs(references);
@@ -1522,8 +1567,14 @@ AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMetho
     Bool_t hasReference = kFALSE;
     while ((obj = refs())){
       if ((ref = dynamic_cast<AliTRDPIDReference *>(obj))){
-        fPIDResponse->Load(ref);
+        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;
       }
     }
@@ -1531,9 +1582,10 @@ AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMetho
     if (!hasReference) {
       AliError("Reference histograms not found in the OCDB");
     }
-
   }
 
+//  }
+
   return fPIDResponse;
 
 }
@@ -1770,10 +1822,8 @@ AliTRDtrapConfig* AliTRDcalibDB::GetTrapConfig()
   else {
     if ((fTrapConfigName.Length() <= 0) || (fTrapConfigVersion.Length() <= 0)) {
       // query the configuration to be used
-      TString configName;
-      this->GetGlobalConfiguration(configName);
-      TString configVersion;
-      this->GetGlobalConfigurationVersion(configVersion);
+      this->GetGlobalConfiguration(fTrapConfigName);
+      this->GetGlobalConfigurationVersion(fTrapConfigVersion);
     }
 
     // try to load the requested configuration
@@ -1782,7 +1832,8 @@ AliTRDtrapConfig* AliTRDcalibDB::GetTrapConfig()
     // if we still don't have a valid TRAPconfig, create a default one
     if (!fTrapConfig) {
       AliWarning("Falling back to default configuration");
-      fTrapConfig = new AliTRDtrapConfig("default", "default TRAP configuration");
+      static AliTRDtrapConfig trapConfigDefault("default", "default TRAP configuration");
+      fTrapConfig = &trapConfigDefault;
       AliTRDtrapConfigHandler cfgHandler(fTrapConfig);
       cfgHandler.Init();
       cfgHandler.LoadConfig();
@@ -1790,6 +1841,47 @@ AliTRDtrapConfig* AliTRDcalibDB::GetTrapConfig()
 
     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;
   }
 }