Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions PWGDQ/Core/VarManager.cxx
Original file line numberDiff line numberDiff line change
Expand Up@@ -837,6 +837,8 @@ void VarManager::SetDefaultVarNames()
fgVariableUnits[kCosThetaCS] = "";
fgVariableNames[kPhiCS] = "#varphi_{CS}";
fgVariableUnits[kPhiCS] = "rad.";
fgVariableNames[kDeltaPhiPair2] = "#Delta#phi";
fgVariableUnits[kDeltaPhiPair2] = "rad.";
fgVariableNames[kPsiPair] = "#Psi_{pair}";
fgVariableUnits[kPsiPair] = "rad.";
fgVariableNames[kDeltaPhiPair] = "#Delta#phi";
Expand Down
11 changes: 11 additions & 0 deletions PWGDQ/Core/VarManager.h
Original file line numberDiff line numberDiff line change
Expand Up@@ -569,6 +569,7 @@ class VarManager : public TObject
kCosThetaCS,
kPhiHE,
kPhiCS,
kDeltaPhiPair2,
kPsiPair,
kDeltaPhiPair,
kOpeningAngle,
Expand DownExpand Up@@ -2347,6 +2348,16 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
double Ptot2 = TMath::Sqrt(v2.Px() * v2.Px() + v2.Py() * v2.Py() + v2.Pz() * v2.Pz());
values[kDeltaPtotTracks] = Ptot1 - Ptot2;

if (fgUsedVars[kDeltaPhiPair2]) {
double phipair2 = v1.Phi() - v2.Phi();
if (phipair2 > 3 * TMath::Pi() / 2) {
values[kDeltaPhiPair2] = phipair2 - 2 * TMath::Pi();
}
if (phipair2 < -TMath::Pi() / 2) {
values[kDeltaPhiPair2] = phipair2 + 2 * TMath::Pi();
}
}

if (fgUsedVars[kPsiPair]) {
values[kDeltaPhiPair] = (t1.sign() * fgMagField > 0.) ? (v1.Phi() - v2.Phi()) : (v2.Phi() - v1.Phi());
double xipair = TMath::ACos((v1.Px() * v2.Px() + v1.Py() * v2.Py() + v1.Pz() * v2.Pz()) / v1.P() / v2.P());
Expand Down