diff --git a/PWGDQ/Core/VarManager.cxx b/PWGDQ/Core/VarManager.cxx
index df1eaa466bd..9e1590a3d91 100644
--- a/PWGDQ/Core/VarManager.cxx
+++ b/PWGDQ/Core/VarManager.cxx
@@ -1726,6 +1726,8 @@ void VarManager::SetDefaultVarNames()
   fgVariableUnits[kCos2ThetaStarFT0C] = "";
   fgVariableNames[kCosThetaStarRandom] = "cos#it{#theta}^{*}_{Random}";
   fgVariableUnits[kCosThetaStarRandom] = "";
+  fgVariableNames[kAbsCosThetaStarRandom] = "|cos#it{#theta}^{*}_{Random}|";
+  fgVariableUnits[kAbsCosThetaStarRandom] = "";
   fgVariableNames[kCos2ThetaStarRandom] = "cos^{2}#it{#theta}^{*}_{Random}";
   fgVariableUnits[kCos2ThetaStarRandom] = "";
   fgVariableNames[kMCCosThetaStar] = "cos#it{#theta}^{*}_{MC}";
@@ -2587,6 +2589,7 @@ void VarManager::SetDefaultVarNames()
   fgVarNamesMap["kAbsCosThetaStarFT0C"] = kAbsCosThetaStarFT0C;
   fgVarNamesMap["kCos2ThetaStarFT0C"] = kCos2ThetaStarFT0C;
   fgVarNamesMap["kCosThetaStarRandom"] = kCosThetaStarRandom;
+  fgVarNamesMap["kAbsCosThetaStarRandom"] = kAbsCosThetaStarRandom;
   fgVarNamesMap["kCos2ThetaStarRandom"] = kCos2ThetaStarRandom;
   fgVarNamesMap["kMCCosThetaStar"] = kMCCosThetaStar;
   fgVarNamesMap["kPairWeight"] = kPairWeight;
diff --git a/PWGDQ/Core/VarManager.h b/PWGDQ/Core/VarManager.h
index 665d5f6bef9..c9e0aa15531 100644
--- a/PWGDQ/Core/VarManager.h
+++ b/PWGDQ/Core/VarManager.h
@@ -840,6 +840,7 @@ class VarManager : public TObject
     kAbsCosThetaStarFT0C,
     kCos2ThetaStarFT0C,
     kCosThetaStarRandom,
+    kAbsCosThetaStarRandom,
     kCos2ThetaStarRandom,
     kCosPhiVP,
     kPhiVP,
@@ -3919,6 +3920,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
     // Randomize the event plane angle to check the unpolarized contribution
     ROOT::Math::XYZVector zaxisRandom = ROOT::Math::XYZVector(TMath::Cos(values[kRandomPsi2]), TMath::Sin(values[kRandomPsi2]), 0).Unit();
     values[kCosThetaStarRandom] = v_CM.Dot(zaxisRandom);
+    values[kAbsCosThetaStarRandom] = TMath::Abs(values[kCosThetaStarRandom]);
     values[kCos2ThetaStarRandom] = values[kCosThetaStarRandom] * values[kCosThetaStarRandom];
 
     // if the truth event plane angle is available, calculate the cos(theta*) with respect to the true event plane angle for comparison
@@ -4522,6 +4524,7 @@ void VarManager::FillPairMEAcrossTFs(T const& t1, T const& t2, float* values)
       // Randomize the event plane angle to check the unpolarized contribution
       ROOT::Math::XYZVector zaxisRandom = ROOT::Math::XYZVector(TMath::Cos(values[kRandomPsi2]), TMath::Sin(values[kRandomPsi2]), 0).Unit();
       values[kCosThetaStarRandom] = v_CM.Dot(zaxisRandom);
+      values[kAbsCosThetaStarRandom] = std::abs(values[kCosThetaStarRandom]);
       values[kCos2ThetaStarRandom] = values[kCosThetaStarRandom] * values[kCosThetaStarRandom];
     }
 
@@ -4717,6 +4720,7 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
     // Randomize the event plane angle to check the unpolarized contribution
     ROOT::Math::XYZVector zaxisRandom = ROOT::Math::XYZVector(TMath::Cos(values[kRandomPsi2]), TMath::Sin(values[kRandomPsi2]), 0).Unit();
     values[kCosThetaStarRandom] = v_CM.Dot(zaxisRandom);
+    values[kAbsCosThetaStarRandom] = std::abs(values[kCosThetaStarRandom]);
     values[kCos2ThetaStarRandom] = values[kCosThetaStarRandom] * values[kCosThetaStarRandom];
 
     // truth event plane angle