Skip to content

[PWGCF] FemtoUniverse: Add PRF boost to identical pairs, remove qTmin selection #14615

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
prchakra merged 6 commits into from
Jan 24, 2026
+245 −237
Merged

[PWGCF] FemtoUniverse: Add PRF boost to identical pairs, remove qTmin selection #14615

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
f5e6a26
Boost the space part of phase to PRF
alibuild Jan 22, 2026
62533e7
Add q comp. in PRF for identical pairs
alibuild Jan 24, 2026
73e671d
Remove qTmin cut, add q comp. in PRF for identical pairs
alibuild Jan 24, 2026
a22eddd
Add proper name for identical pair in LCMS
alibuild Jan 24, 2026
7bc5d2a
Adapt to updated newpairfunc
alibuild Jan 24, 2026
faa71d0
Fix clang-format issues
alibuild Jan 24, 2026
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
28 changes: 21 additions & 7 deletions PWGCF/FemtoUniverse/Core/FemtoUniverse3DContainer.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,13 +24,13 @@
#include <Framework/Logger.h>

#include "Math/Vector4D.h"
#include "TDatabasePDG.h"

Check failure on line 27 in PWGCF/FemtoUniverse/Core/FemtoUniverse3DContainer.h

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
#include "TMath.h"

#include <string>
#include <vector>

using namespace o2::framework;

Check failure on line 33 in PWGCF/FemtoUniverse/Core/FemtoUniverse3DContainer.h

View workflow job for this annotation

GitHub Actions / O2 linter

[using-directive] 

Do not put using directives at global scope in headers.

namespace o2::analysis::femto_universe
{
Expand Down Expand Up @@ -158,8 +158,8 @@
/// \param pdg2 PDG code of particle two
void setPDGCodes(const int pdg1, const int pdg2)
{
mMassOne = TDatabasePDG::Instance()->GetParticle(pdg1)->Mass();

Check failure on line 161 in PWGCF/FemtoUniverse/Core/FemtoUniverse3DContainer.h

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
mMassTwo = TDatabasePDG::Instance()->GetParticle(pdg2)->Mass();

Check failure on line 162 in PWGCF/FemtoUniverse/Core/FemtoUniverse3DContainer.h

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
mPDGOne = pdg1;
mPDGTwo = pdg2;
}
Expand Down Expand Up @@ -216,22 +216,36 @@
/// \param isiden Choosing identical or non-identical pairs
/// \param islcm Choosing LCMS or PRF
template <bool isMC, typename T>
void setPair(T const& part1, T const& part2, const int mult, bool use3dplots, bool isiden)
void setPair(T const& part1, T const& part2, const int mult, bool use3dplots, bool isIdenLCMS, bool isWeight, bool isIdenPRF)
{
std::vector<double> f3d;
const float kT = FemtoUniverseMath::getkT(part1, mMassOne, part2, mMassTwo);
const float mT = FemtoUniverseMath::getmT(part1, mMassOne, part2, mMassTwo);

f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo, isiden);
f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo, isIdenLCMS, isWeight, isIdenPRF);

float varout = 0.0;
float varside = 0.0;
float varlong = 0.0;

if (isIdenLCMS) {
varout = f3d[1];
varside = f3d[2];
varlong = f3d[3];
} else if (isIdenPRF) {
varout = f3d[6];
varside = f3d[7];
varlong = f3d[8];
}

const float femtoObs1D = f3d[0];
const float femtoObsKout = f3d[1];
const float femtoObsKside = f3d[2];
const float femtoObsKlong = f3d[3];
const float femtoObsKout = varout;
const float femtoObsKside = varside;
const float femtoObsKlong = varlong;

if (mHistogramRegistry) {
setPairBase<o2::aod::femtouniverse_mc_particle::MCType::kRecon>(femtoObsKout, femtoObsKside, femtoObsKlong, femtoObs1D, kT, mT, part1, part2, mult, use3dplots, isiden);
if (!isiden) {
setPairBase<o2::aod::femtouniverse_mc_particle::MCType::kRecon>(femtoObsKout, femtoObsKside, femtoObsKlong, femtoObs1D, kT, mT, part1, part2, mult, use3dplots, isIdenLCMS);
if (!isIdenLCMS) {
if (femtoObsKout > 0.0) {
mHistogramRegistry->fill(HIST(FolderSuffix[EventType]) + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[o2::aod::femtouniverse_mc_particle::MCType::kRecon]) + HIST("/KStarOutP"), femtoObs1D);
} else {
Expand Down
14 changes: 11 additions & 3 deletions PWGCF/FemtoUniverse/Core/FemtoUniverseMath.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@

#include "Math/Boost.h"
#include "Math/Vector4D.h"
#include "TLorentzVector.h"

Check failure on line 24 in PWGCF/FemtoUniverse/Core/FemtoUniverseMath.h

View workflow job for this annotation

GitHub Actions / O2 linter

[root/lorentz-vector] 

Do not use the TLorentzVector legacy class. Use std::array with RecoDecay methods or the ROOT::Math::LorentzVector template instead.
#include "TMath.h"

#include <algorithm>
Expand Down Expand Up @@ -177,7 +177,7 @@
/// \param mass2 Mass of particle 2
/// \param isiden Identical or non-identical particle pair
template <typename T>
static std::vector<double> newpairfunc(const T& part1, const float mass1, const T& part2, const float mass2, bool isiden, bool isWeight = 0)
static std::vector<double> newpairfunc(const T& part1, const float mass1, const T& part2, const float mass2, bool isIdenLCMS, bool isWeight, bool isIdenPRF)
{
const double e1 = std::sqrt(std::pow(part1.px(), 2) + std::pow(part1.py(), 2) + std::pow(part1.pz(), 2) + std::pow(mass1, 2));
const double e2 = std::sqrt(std::pow(part2.px(), 2) + std::pow(part2.py(), 2) + std::pow(part2.pz(), 2) + std::pow(mass2, 2));
Expand Down Expand Up @@ -240,7 +240,7 @@
const double qinv = std::sqrt(fDKOutPRF * fDKOutPRF + fDKSidePRF * fDKSidePRF + fDKLongPRF * fDKLongPRF);
const double kstar = std::sqrt(fKOut * fKOut + fDKSide * fDKSide + fDKLong * fDKLong);

if (isiden) {
if (isIdenLCMS) {
vect.push_back(qinv);
vect.push_back(fDKOutLCMS);
vect.push_back(fDKSideLCMS);
Expand All @@ -251,9 +251,10 @@
vect.push_back(fDKOut);
vect.push_back(fDKSide);
vect.push_back(fDKLong);
vect.push_back(1.0);
}

if (isiden && isWeight) {
if (isIdenLCMS && isWeight) {
const double x1_lcms = (vecspace_part1.x() * tPx + vecspace_part1.y() * tPy) / tPt;
const double y1_lcms = (-vecspace_part1.x() * tPy + vecspace_part1.y() * tPx) / tPt;
const double x2_lcms = (vecspace_part2.x() * tPx + vecspace_part2.y() * tPy) / tPt;
Expand All @@ -277,6 +278,13 @@
} else {
vect.push_back(1.0);
}

if (isIdenPRF) {
vect.push_back(fDKOutPRF);
vect.push_back(fDKSidePRF);
vect.push_back(fDKLongPRF);
}

return vect;
}
};
Expand Down
33 changes: 23 additions & 10 deletions PWGCF/FemtoUniverse/Core/FemtoUniversePairSHCentMultKt.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -180,7 +180,7 @@ class PairSHCentMultKt
/// \param ktval kT value
template <typename T>
void fillMultNumDen(T const& part1, T const& part2, uint8_t ChosenEventType,
int maxl, int multval, float ktval, bool isiden, bool isqinvfill)
int maxl, int multval, float ktval, bool isIdenLCMS, bool isqinvfill, bool isWeight, bool isIdenPRF)
{
int multbinval;
int absmultval = multval;
Expand All @@ -197,7 +197,7 @@ class PairSHCentMultKt
return;
}
// std::cout<<"multbinval "<<multbinval<<std::endl;
fillkTNumDen(part1, part2, ChosenEventType, maxl, multbinval, ktval, isiden, isqinvfill);
fillkTNumDen(part1, part2, ChosenEventType, maxl, multbinval, ktval, isIdenLCMS, isqinvfill, isWeight, isIdenPRF);
}

/// Templated function to access different kT directory and call addEventPair
Expand All @@ -209,7 +209,7 @@ class PairSHCentMultKt
/// \param ktval kT value
template <typename T>
void fillkTNumDen(T const& part1, T const& part2, uint8_t ChosenEventType,
int maxl, int multval, float ktval, bool isiden, bool isqinvfill)
int maxl, int multval, float ktval, bool isIdenLCMS, bool isqinvfill, bool isWeight, bool isIdenPRF)
{
int ktbinval = -1;
if (ktval >= ktBins[0] && ktval < ktBins[1]) {
Expand All @@ -229,7 +229,7 @@ class PairSHCentMultKt
} else {
return;
}
addEventPair(part1, part2, ChosenEventType, maxl, multval, ktbinval, isiden, isqinvfill);
addEventPair(part1, part2, ChosenEventType, maxl, multval, ktbinval, isIdenLCMS, isqinvfill, isWeight, isIdenPRF);
}

/// Set the PDG codes of the two particles involved
Expand Down Expand Up @@ -262,20 +262,33 @@ class PairSHCentMultKt
/// \param ktval kT value
template <typename T>
void addEventPair(T const& part1, T const& part2, uint8_t ChosenEventType,
int /*maxl*/, int multval, int ktval, bool isiden, bool isqinvfill)
int /*maxl*/, int multval, int ktval, bool isIdenLCMS, bool isqinvfill, bool isWeight, bool isIdenPRF)
{
int fMultBin = multval;
int fKtBin = ktval;
std::vector<std::complex<double>> fYlmBuffer(kMaxJM);
std::vector<double> f3d;
setPionPairMass();
f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo,
isiden);

const float qout = f3d[1];
const float qside = f3d[2];
const float qlong = f3d[3];
f3d = FemtoUniverseMath::newpairfunc(part1, mMassOne, part2, mMassTwo,
isIdenLCMS, isWeight, isIdenPRF);
double varout = 0.0;
double varside = 0.0;
double varlong = 0.0;

if (isIdenLCMS) {
varout = f3d[1];
varside = f3d[2];
varlong = f3d[3];
} else if (isIdenPRF) {
varout = f3d[6];
varside = f3d[7];
varlong = f3d[8];
}

const double qout = varout;
const double qside = varside;
const double qlong = varlong;
double kv = std::sqrt(qout * qout + qside * qside + qlong * qlong);

// int nqbin = fbinctn[0][0]->GetXaxis()->FindFixBin(kv);
Expand Down
Loading
Loading