From f90cd618d26326fe8edf07ae16eebf554ded650b Mon Sep 17 00:00:00 2001 From: zjxiongOvO Date: Mon, 4 May 2026 11:10:41 +0200 Subject: [PATCH] add acceptance cut for pbpb --- PWGDQ/Core/CutsLibrary.cxx | 9 +++++++++ PWGDQ/Core/VarManager.cxx | 1 + PWGDQ/Core/VarManager.h | 4 ++++ 3 files changed, 14 insertions(+) diff --git a/PWGDQ/Core/CutsLibrary.cxx b/PWGDQ/Core/CutsLibrary.cxx index 5c0a0a67de9..a5adf320df5 100644 --- a/PWGDQ/Core/CutsLibrary.cxx +++ b/PWGDQ/Core/CutsLibrary.cxx @@ -4622,6 +4622,15 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName) return cut; } + if (!nameStr.compare("acceptance_PbPb5360")) { + cut->AddCut(VarManager::kMCY, -0.9, 0.9); + cut->AddCut(VarManager::kMCP1, 1.0, 1000.0); + cut->AddCut(VarManager::kMCP2, 1.0, 1000.0); + cut->AddCut(VarManager::kMCEta1, -0.9, 0.9); + cut->AddCut(VarManager::kMCEta2, -0.9, 0.9); + return cut; + } + // --------------------------------------------------- // MC generated particle acceptance cuts diff --git a/PWGDQ/Core/VarManager.cxx b/PWGDQ/Core/VarManager.cxx index b236b849b1a..fa081babced 100644 --- a/PWGDQ/Core/VarManager.cxx +++ b/PWGDQ/Core/VarManager.cxx @@ -1953,6 +1953,7 @@ void VarManager::SetDefaultVarNames() fgVarNamesMap["kMCEventWeight"] = kMCEventWeight; fgVarNamesMap["kMCEventImpParam"] = kMCEventImpParam; fgVarNamesMap["kMCEventPlaneAngle"] = kMCEventPlaneAngle; + fgVarNamesMap["kMCEventCentrFT0C"] = kMCEventCentrFT0C; fgVarNamesMap["kQ1ZNAX"] = kQ1ZNAX; fgVarNamesMap["kQ1ZNAY"] = kQ1ZNAY; fgVarNamesMap["kQ1ZNCX"] = kQ1ZNCX; diff --git a/PWGDQ/Core/VarManager.h b/PWGDQ/Core/VarManager.h index 6d45eb620cf..37e0523cd91 100644 --- a/PWGDQ/Core/VarManager.h +++ b/PWGDQ/Core/VarManager.h @@ -715,8 +715,10 @@ class VarManager : public TObject // MC pair variables kMCPt1, kMCEta1, + kMCP1, kMCPt2, kMCEta2, + kMCP2, kMCCosThetaHE, kMCPhiHE, kMCPhiTildeHE, @@ -4201,6 +4203,8 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values) values[kMCPt2] = t2.pt(); values[kMCEta1] = t1.eta(); values[kMCEta2] = t2.eta(); + values[kMCP1] = t1.p(); + values[kMCP2] = t2.p(); // polarization parameters bool useHE = fgUsedVars[kMCCosThetaHE] || fgUsedVars[kMCPhiHE]; // helicity frame