[PWGLF] updated code with check of fine time rec (#16835)

Co-authored-by: Prottay Das <prottay@alipap1.cern.ch>

Author: prottayCMT

PWGLF/TableProducer/Common/spvector.cxx

@@ -127,10 +127,13 @@ struct spvector {
   Configurable<bool> useCallibvertex{"useCallibvertex", false, "use calibration for vxy"};
   Configurable<bool> coarse1{"coarse1", false, "RE1"};
   Configurable<bool> fine1{"fine1", false, "REfine1"};
+  Configurable<bool> finetime1{"finetime1", false, "REfinetime1"};
   Configurable<bool> coarse2{"coarse2", false, "RE2"};
   Configurable<bool> fine2{"fine2", false, "REfine2"};
+  Configurable<bool> finetime2{"finetime2", false, "REfinetime2"};
   Configurable<bool> coarse3{"coarse3", false, "RE3"};
   Configurable<bool> fine3{"fine3", false, "REfine3"};
+  Configurable<bool> finetime3{"finetime3", false, "REfinetime3"};
   Configurable<bool> coarse4{"coarse4", false, "RE4"};
   Configurable<bool> fine4{"fine4", false, "REfine4"};
   Configurable<bool> coarse5{"coarse5", false, "RE5"};
@@ -139,6 +142,7 @@ struct spvector {
   Configurable<bool> fine6{"fine6", false, "REfine6"};
   Configurable<bool> useRecentereSp{"useRecentereSp", false, "use Recentering with Sparse or THn"};
   Configurable<bool> useRecenterefineSp{"useRecenterefineSp", false, "use fine Recentering with THn"};
+  Configurable<bool> useTimeRecentering{"useTimeRecentering", false, "Use residual time recentering"};
   Configurable<std::string> ConfGainPath{"ConfGainPath", "Users/p/prottay/My/Object/NewPbPbpass4_10092024/gaincallib", "Path to gain calibration"};
   Configurable<std::string> ConfGainPathvxy{"ConfGainPathvxy", "Users/p/prottay/My/Object/swapcoords/PbPbpass4_20112024/recentervert", "Path to gain calibration for vxy"};
   Configurable<std::string> ConfRecentereSp{"ConfRecentereSp", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for recentere"};
@@ -173,8 +177,9 @@ struct spvector {
   Configurable<std::string> ConfRecenterevzSp6{"ConfRecenterevzSp6", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for vz recentere6"};
   Configurable<std::string> ConfShiftC{"ConfShiftC", "Users/p/prottay/My/Object/Testinglocaltree/shiftcallib2", "Path to shift C"};
   Configurable<std::string> ConfShiftA{"ConfShiftA", "Users/p/prottay/My/Object/Testinglocaltree/shiftcallib2", "Path to shift A"};
-  Configurable<bool> useTimeRecentering{"useTimeRecentering", false, "Use residual time recentering"};
-  Configurable<std::string> confRecentereTimeSp{"confRecentereTimeSp", "Users/p/prottay/My/Object/GCwithoutcfactorgoodVztimedep/From676541/TestDDlocal/2024PbPbpass3_23062026/recenterlast2", "Path to time recentering map"};
+  Configurable<std::string> confRecentereTimeSp1{"confRecentereTimeSp1", "Users/p/prottay/My/Object/GCwithoutcfactorgoodVztimedep/From676541/TestDDlocal/2024PbPbpass3_23062026/recenterlast2", "Path to time recentering map 1"};
+  Configurable<std::string> confRecentereTimeSp2{"confRecentereTimeSp2", "Users/p/prottay/My/Object/GCwithoutcfactorgoodVztimedep/From676541/TestDDlocal/2024PbPbpass3_23062026/recenterlast2", "Path to time recentering map 2"};
+  Configurable<std::string> confRecentereTimeSp3{"confRecentereTimeSp3", "Users/p/prottay/My/Object/GCwithoutcfactorgoodVztimedep/From676541/TestDDlocal/2024PbPbpass3_23062026/recenterlast3", "Path to time recentering map 3"};
 
   // Event selection cuts - Alex
   /*
@@ -314,7 +319,9 @@ struct spvector {
   std::array<TH2F*, 6> hrecenterevzSpA;   // Array of 5 histograms
   TProfile3D* shiftprofileA;
   TProfile3D* shiftprofileC;
-  TH2F* hrecentereTimeSp = nullptr;
+  TH2F* hrecentereTimeSp1 = nullptr;
+  TH2F* hrecentereTimeSp2 = nullptr;
+  TH2F* hrecentereTimeSp3 = nullptr;
 
   Bool_t Correctcoarse(const THnF* hrecentereSp, auto centrality, auto vx, auto vy, auto vz, auto& qxZDCA, auto& qyZDCA, auto& qxZDCC, auto& qyZDCC)
   {
@@ -652,6 +659,14 @@ struct spvector {
         resfine = Correctfine(hrecenterecentSpA[0], hrecenterevxSpA[0], hrecenterevySpA[0], hrecenterevzSpA[0], centrality, vx, vy, vz, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
       }
 
+      if (finetime1 && (currentRunNumber != lastRunNumber)) {
+        hrecentereTimeSp1 = ccdb->getForTimeStamp<TH2F>(confRecentereTimeSp1.value, bc.timestamp());
+      }
+      bool restime = false;
+      if (useTimeRecentering) {
+        restime = Correcttime(hrecentereTimeSp1, timeMin, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
+      }
+
       if (coarse2) {
         if (useRecentereSp && (currentRunNumber != lastRunNumber)) {
           hrecentereSpA[1] = ccdb->getForTimeStamp<THnF>(ConfRecentereSp2.value, bc.timestamp());
@@ -669,6 +684,13 @@ struct spvector {
         resfine = Correctfine(hrecenterecentSpA[1], hrecenterevxSpA[1], hrecenterevySpA[1], hrecenterevzSpA[1], centrality, vx, vy, vz, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
       }
 
+      if (finetime2 && (currentRunNumber != lastRunNumber)) {
+        hrecentereTimeSp2 = ccdb->getForTimeStamp<TH2F>(confRecentereTimeSp2.value, bc.timestamp());
+      }
+      if (useTimeRecentering) {
+        restime = Correcttime(hrecentereTimeSp2, timeMin, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
+      }
+
       if (coarse3) {
         if (useRecentereSp && (currentRunNumber != lastRunNumber)) {
           hrecentereSpA[2] = ccdb->getForTimeStamp<THnF>(ConfRecentereSp3.value, bc.timestamp());
@@ -686,6 +708,13 @@ struct spvector {
         resfine = Correctfine(hrecenterecentSpA[2], hrecenterevxSpA[2], hrecenterevySpA[2], hrecenterevzSpA[2], centrality, vx, vy, vz, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
       }
 
+      if (finetime3 && (currentRunNumber != lastRunNumber)) {
+        hrecentereTimeSp3 = ccdb->getForTimeStamp<TH2F>(confRecentereTimeSp3.value, bc.timestamp());
+      }
+      if (useTimeRecentering) {
+        restime = Correcttime(hrecentereTimeSp3, timeMin, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
+      }
+
       if (coarse4) {
         if (useRecentereSp && (currentRunNumber != lastRunNumber)) {
           hrecentereSpA[3] = ccdb->getForTimeStamp<THnF>(ConfRecentereSp4.value, bc.timestamp());
@@ -741,16 +770,7 @@ struct spvector {
         LOG(info) << "Histograms are null";
       }
 
-      if (useTimeRecentering && (currentRunNumber != lastRunNumber)) {
-        hrecentereTimeSp = ccdb->getForTimeStamp<TH2F>(confRecentereTimeSp.value, bc.timestamp());
-      }
-
-      bool restime = false;
-      if (useTimeRecentering) {
-        restime = Correcttime(hrecentereTimeSp, timeMin, qxZDCA, qyZDCA, qxZDCC, qyZDCC);
-      }
-
-      if (restime == 0) {
+      if (useTimeRecentering && restime == 0 && check == 0) {
         LOG(info) << "Histograms are null";
       }