Move parameters from the source code into the header file

reco/detectors/trd/CbmTrdModuleRecR.cxx

@@ -562,54 +562,17 @@ CbmTrdHit* CbmTrdModuleRecR::MakeHit(Int_t clusterId,
     hit_pos[iDim] = hit_posV[iDim];
   }
 
-  if (EB) {
-    if (errorclass == 0)
-      xVar = 0.0258725;
-    else if (errorclass == 1)
-      xVar = 0.0267693;
-    else if (errorclass == 2)
-      xVar = 0.0344325;
-    else if (errorclass == 3)
-      xVar = 0.0260322;
-    else if (errorclass == 4)
-      xVar = 0.040115;
-
-    if (errorclass == 0)
-      yVar = 0.024549;
-    else if (errorclass == 1)
-      yVar = 0.025957;
-    else if (errorclass == 2)
-      yVar = 0.0250713;
-    else if (errorclass == 3)
-      yVar = 0.0302682;
-    else if (errorclass == 4)
-      yVar = 0.0291146;
-  } else {
-    if (EBP) time -= 46;  // due to the event time of 0 in the EB
-    // mode and the ULong in the the digi time
-    // TODO: move variables to parameter file
-    if (errorclass == 0)
-      xVar = 0.0426313;
-    else if (errorclass == 1)
-      xVar = 0.0426206;
-    else if (errorclass == 2)
-      xVar = 0.0636962;
-    else if (errorclass == 3)
-      xVar = 0.038981;
-    else if (errorclass == 4)
-      xVar = 0.0723851;
-
-    if (errorclass == 0)
-      yVar = 0.0401438;
-    else if (errorclass == 1)
-      yVar = 0.0407502;
-    else if (errorclass == 2)
-      yVar = 0.0397242;
-    else if (errorclass == 3)
-      yVar = 0.0519485;
-    else if (errorclass == 4)
-      yVar = 0.0504586;
-  }
+   if (EB) {
+     xVar = kxVar_Value[0][errorclass];
+     yVar = kyVar_Value[0][errorclass];
+   } else {
+     if (EBP)
+       time -=
+         46;  //due to the event time of 0 in the EB mode and the ULong in the the digi time
+         //TODO: move to parameter file
+     xVar = kxVar_Value[1][errorclass];
+     yVar = kyVar_Value[1][errorclass];
+   }
 
   TVector3 cluster_pad_dposV(xVar, yVar, 0);
 

reco/detectors/trd/CbmTrdModuleRecR.h

@@ -60,6 +60,18 @@ private:
   Int_t fDigiCounter;  // digits over threshold
   Int_t fCheck = 0;
 
+  // different error classes for the position resolution based on the simulation results
+  // the error classes are defined for the different module types
+  // TODO: move to parameter file 
+  static constexpr Double_t kxVar_Value[2][5] = { 
+               {0.0258725, 0.0267693, 0.0344325, 0.0260322, 0.040115}, 
+               {0.0426313, 0.0426206, 0.0636962, 0.038981, 0.0723851} 
+                             };
+  static constexpr Double_t kyVar_Value[2][5] = { 
+               {0.024549, 0.025957, 0.0250713, 0.0302682, 0.0291146},
+               {0.0401438, 0.0407502, 0.0397242, 0.0519485, 0.0504586}
+                              };
+
   std::deque<std::tuple<Int_t, Bool_t, const CbmTrdDigi*>>
     fDigiMap;  //map to sort all digis from the Array into a deque; different module are separated; the tuple contains the digi indice, a bool to flag processed digis and the digi itself
   std::deque<std::deque<std::pair<Int_t, const CbmTrdDigi*>>>