]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDcalibDB.cxx
Added some plots
[u/mrichter/AliRoot.git] / TRD / AliTRDcalibDB.cxx
index aba9909ca8621f7d93a2524b520fe72d9566393b..927b4f855a1507740359010e278728e3a4f73e23 100644 (file)
@@ -42,6 +42,7 @@
 #include "AliTRDtrapConfig.h"
 #include "AliTRDtrapConfigHandler.h"
 #include "AliTRDCommonParam.h"
+#include "AliTRDgeometry.h"
 
 #include "Cal/AliTRDCalROC.h"
 #include "Cal/AliTRDCalPad.h"
@@ -287,6 +288,8 @@ const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
          // Online gain table ID 8
           return CacheCDBEntry(kIDOnlineGainFactor  ,"TRD/Calib/Krypton_2012-01"); 
           break;
+      default:
+       AliError(Form("unknown gaintable requested with ID"));
       }
       break;
 
@@ -318,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");
@@ -375,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()));
     }
   } 
   
@@ -1010,39 +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;
-    }
-    if (tableName.CompareTo("Krypton_2011-03")               == 0) {
+    else if (tableName.CompareTo("Krypton_2012-01")               == 0)
       fOnlineGainTableID = 8;
-      return fOnlineGainTableID;
-    }
+    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 {
 
@@ -1180,14 +1176,15 @@ void AliTRDcalibDB::GetGlobalConfigurationVersion(TString &version)
 }
 
 //_____________________________________________________________________________
-Int_t AliTRDcalibDB::GetNumberOfParsDCS(TString cname)
+Int_t AliTRDcalibDB::GetNumberOfParsDCS(TString cname, Char_t delimiter)
 {
   // 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"
-
+  //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;
 }
 
 //_____________________________________________________________________________
@@ -1195,22 +1192,20 @@ 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();
+  Char_t cdelim = '_', // define the delimiters
+         odelim = '-';
+  Int_t nconfig = GetNumberOfParsDCS(cname, cdelim);
 
   // protect
-  if ((nconfig == 0) || ((nconfig-1) < cfgType)) {
+  if ((nconfig == -1) || (nconfig < 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
-
+  TObjArray *carr = cname.Tokenize(cdelim);
+  Int_t nopt = GetNumberOfParsDCS(((TObjString*)carr->At(cfgType))->GetString(), odelim);
+  carr->Delete(); delete carr;
+  return nopt;
 }
 
 //_____________________________________________________________________________
@@ -1222,36 +1217,35 @@ void AliTRDcalibDB::GetDCSConfigParOption(TString cname, Int_t cfgType, Int_t op
   // option >  0 returns the optional parameter Nr. (option) of the configuration (cfgType)
   //
 
-  // define the delimiters
-  TString cdelim = "_";
-  TString odelim = "-";
-
-  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) {
+  if (nconfig == -1) {
     AliError("DCS configuration name empty!");
     cfgo = "";
     return;
-  } else if ((nconfig-1) < cfgType) {
+  } else if (nconfig < cfgType) {
     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();
-
+  TObjArray *carr = cname.Tokenize(cdelim);
+  TString cfgString(((TObjString*)carr->At(cfgType))->GetString());
+  Int_t noptions = GetNumberOfParsDCS(cfgString, odelim);
   // protect
-  if ((noptions-1) < option) {
+  if (noptions < option) {
     AliError("Not enough options in DCS configuration name!");
     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;
 
 }
@@ -1513,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)
 {
@@ -1547,7 +1558,7 @@ AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMetho
 
     // Load Reference Histos from OCDB
 //    if(method == AliTRDPIDResponse::kLQ1D){
-    fPIDResponse->SetPIDmethod(AliTRDPIDResponse::kLQ1D);
+    //fPIDResponse->SetPIDmethod(AliTRDPIDResponse::kLQ1D);
     const TObjArray *references = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDPIDLQ1D));
 
     TIter refs(references);
@@ -1811,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
@@ -1823,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();
@@ -1831,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;
   }
 }