]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDcalibDB.cxx
- add setter for error in tracklet y
[u/mrichter/AliRoot.git] / TRD / AliTRDcalibDB.cxx
index 4255765c3908fd90d0a5dea8bffb70a5947e5a62..2c412156a67c4949dbd1cd1fa82e0f179faea9f8 100644 (file)
 #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/AliTRDCalDCS.h"
+#include "Cal/AliTRDCalDCSFEE.h"
 #include "Cal/AliTRDCalDCSv2.h"
 #include "Cal/AliTRDCalDCSFEEv2.h"
 #include "Cal/AliTRDCalPID.h"
@@ -287,6 +289,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 +322,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 +379,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()));
     }
   } 
   
@@ -935,6 +941,54 @@ Int_t AliTRDcalibDB::ExtractTimeBinsFromString(TString tbstr)
 
 }
 
+//_____________________________________________________________________________
+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()
 {
@@ -942,15 +996,51 @@ 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);
+  // 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);
 
 }
 
@@ -1010,39 +1100,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 {
 
@@ -1130,7 +1211,206 @@ void AliTRDcalibDB::GetGlobalConfiguration(TString &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)
 {
@@ -1180,14 +1460,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 +1476,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 +1501,37 @@ 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) {
-    AliError("Not enough parameters in DCS configuration name!");
+  } else if (nconfig < cfgType) {
+    AliError(Form("Not enough parameters in DCS configuration name!"
+                 " Name %s",cname.Data()));
     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) {
-    AliError("Not enough options in DCS configuration name!");
+  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;
 
 }
@@ -1828,10 +2108,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
@@ -1840,7 +2118,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();
@@ -1848,6 +2127,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;
   }
 }