added histogram to monitor run dependent variations of vxy and commented out few lines of code and a histogram

Author: Prottay Das

PWGLF/TableProducer/Common/spvector.cxx

@@ -139,12 +139,14 @@ struct spvector {
   Configurable<bool> followpub{"followpub", true, "flag to use alphaZDC"};
   Configurable<bool> useGainCallib{"useGainCallib", false, "use gain calibration"};
   Configurable<bool> useRecentereSp{"useRecentereSp", false, "use Recentering with Sparse or THn"};
+  Configurable<bool> useRecentereVxy{"useRecentereVxy", false, "use Recentering with Sparse or THn for Vxy"};
   Configurable<bool> useRecenterefineSp{"useRecenterefineSp", false, "use fine Recentering with Sparse or THn"};
   Configurable<bool> useRecenteresqSp{"useRecenteresqSp", false, "use Recenteringsq with Sparse or THn"};
   Configurable<bool> recwitherror{"recwitherror", false, "use Recentering with error"};
   Configurable<bool> recfinewitherror{"recfinewitherror", false, "use Recentering fine with error"};
   Configurable<std::string> ConfGainPath{"ConfGainPath", "Users/p/prottay/My/Object/NewPbPbpass4_10092024/gaincallib", "Path to gain calibration"};
   Configurable<std::string> ConfRecentereSp{"ConfRecentereSp", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for recentere"};
+  Configurable<std::string> ConfRecentereVxy{"ConfRecentereVxy", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for recentereVxy"};
   Configurable<std::string> ConfRecenterecentSp{"ConfRecenterecentSp", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for cent recentere"};
   Configurable<std::string> ConfRecenterevxSp{"ConfRecenterevxSp", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for vx recentere"};
   Configurable<std::string> ConfRecenterevySp{"ConfRecenterevySp", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn Path for vy recentere"};
@@ -187,6 +189,7 @@ struct spvector {
     //  AxisSpec amplitudeZDC = {configAxisZDCgain, "ZDC amplitude"};
     // AxisSpec amplitudeZDC = {ZDCgainNbins, lbinZDCgain, hbinZDCgain, "ZDC amplitude"};
     AxisSpec channelZDCAxis = {8, 0.0, 8.0, "ZDC tower"};
+    AxisSpec binAxis = {2, 0.0, 2.0, "Vx and vy bin positions"};
     AxisSpec qxZDCAxis = {QxyNbins, lbinQxy, hbinQxy, "Qx"};
     AxisSpec phiAxis = {PhiNbins, -6.28, 6.28, "phi"};
     AxisSpec vzAxis = {VzNbins, lbinVz, hbinVz, "vz"};
@@ -248,7 +251,8 @@ struct spvector {
     // histos.add("ZDCAmpCommon", "ZDCAmpCommon", kTProfile3D, {{2,0.0,2.0}, vzfineAxis, centfineAxis});
 
     if (QA) {
-      histos.add("Vz", "Vz", kTH1F, {vzfineAxis});
+      histos.add("hVz", "hVz", kTH1F, {vzfineAxis});
+      histos.add("hpVxVy", "hpVxVy", kTProfile, {binAxis});
       histos.add("hpCosPsiAPsiC", "hpCosPsiAPsiC", kTProfile, {centfineAxis});
       histos.add("hpSinPsiAPsiC", "hpSinPsiAPsiC", kTProfile, {centfineAxis});
     }
@@ -275,8 +279,8 @@ struct spvector {
 
   int currentRunNumber = -999;
   int lastRunNumber = -999;
-  // TH3D* gainprofile;
   TH2D* gainprofile;
+  TProfile* hrecentereVxy;
   THnF* hrecentereSp;
   TH2F* hrecenterecentSp;
   TH2F* hrecenterevxSp;
@@ -439,7 +443,16 @@ struct spvector {
       }
 
       histos.fill(HIST("hCentrality"), centrality);
-      histos.fill(HIST("Vz"), vz);
+      histos.fill(HIST("hVz"), vz);
+
+      histos.fill(HIST("hpVxVy"), 0.5, vx);
+      histos.fill(HIST("hpVxVy"), 1.5, vy);
+
+      if (useRecentereVxy && (currentRunNumber != lastRunNumber)) {
+        hrecentereVxy = ccdb->getForTimeStamp<TProfile>(ConfRecentereVxy.value, bc.timestamp());
+        vx = vx - hrecentereVxy->GetBinContent(1);
+        vy = vy - hrecentereVxy->GetBinContent(2);
+      }
 
       if (useRecentereSp && (currentRunNumber != lastRunNumber)) {
         hrecentereSp = ccdb->getForTimeStamp<THnF>(ConfRecentereSp.value, bc.timestamp());
@@ -521,7 +534,7 @@ struct spvector {
           }
         }
       }
-
+      /*
       if (useRecenterefineSp && hrecenterecentSp) {
 
         double meanxAcent = hrecenterecentSp->GetBinContent(hrecenterecentSp->FindBin(centrality, 0.5));
@@ -653,7 +666,7 @@ struct spvector {
           qyZDCC = qyZDCC / meanyCvzerror;
         }
       }
-
+      */
       // LOG(info) << "**********qxa values in spvector************" << qxZDCA<<" "<<centrality<<" "<<vx<<" "<<vy<<" "<<vz;
 
       psiZDCC = 1.0 * TMath::ATan2(qyZDCC, qxZDCC);

PWGLF/Tasks/Strangeness/lambdapolsp.cxx

@@ -395,8 +395,7 @@ struct lambdapolsp {
       return;
     }
     auto centrality = collision.centFT0C();
-
-    histos.fill(HIST("hCentrality0"), centrality);
+    // histos.fill(HIST("hCentrality0"), centrality);
     if (!collision.triggereventsp()) {
       return;
     }