Skip to content

[PWGLF] generalized k892hadronphotonBkg to work for LambdaStar & bug fixes #17480

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

Open
Oussama3030 wants to merge 2 commits into
base: master
Choose a base branch
from
+354 −151
Open

[PWGLF] generalized k892hadronphotonBkg to work for LambdaStar & bug fixes #17480

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a0b1ba3
generalized k892hadronphotonBkg to work for LambdaStar & bug fixes
Oussama3030 Aug 17, 2026
bf1bf98
minor fixes
Oussama3030 Aug 18, 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
2 changes: 1 addition & 1 deletion PWGLF/DataModel/LFSigmaTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@
DECLARE_SOA_COLUMN(X, x, float);
DECLARE_SOA_COLUMN(Y, y, float);
DECLARE_SOA_COLUMN(Z, z, float);
DECLARE_SOA_COLUMN(DCADaughters, dcadaughters, float);

Check failure on line 41 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.

DECLARE_SOA_COLUMN(PhotonPx, photonPx, float);
DECLARE_SOA_COLUMN(PhotonPy, photonPy, float);
Expand All @@ -49,7 +49,7 @@
DECLARE_SOA_COLUMN(LambdaPy, lambdaPy, float);
DECLARE_SOA_COLUMN(LambdaPz, lambdaPz, float);
DECLARE_SOA_COLUMN(LambdaMass, lambdaMass, float);
DECLARE_SOA_COLUMN(AntiLambdaMass, antilambdaMass, float);

Check failure on line 52 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.

//______________________________________________________
// DYNAMIC COLUMNS
Expand Down Expand Up @@ -100,7 +100,7 @@
DECLARE_SOA_DYNAMIC_COLUMN(Radius, radius, //! Sigma0 decay radius (2D, centered at zero)
[](float x, float y) -> float { return RecoDecay::sqrtSumOfSquares(x, y); });

DECLARE_SOA_DYNAMIC_COLUMN(OPAngle, opAngle,

Check failure on line 103 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
[](float photonPx, float photonPy, float photonPz, float lambdaPx, float lambdaPy, float lambdaPz) {
TVector3 v1(photonPx, photonPy, photonPz);
TVector3 v2(lambdaPx, lambdaPy, lambdaPz);
Expand All @@ -113,7 +113,7 @@
return RecoDecay::sqrtSumOfSquares(photonPx, photonPy);
});

DECLARE_SOA_DYNAMIC_COLUMN(PhotonP, photonp, //! Total momentum in GeV/c

Check failure on line 116 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
[](float photonPx, float photonPy, float photonPz) -> float {
return RecoDecay::sqrtSumOfSquares(photonPx, photonPy, photonPz);
});
Expand All @@ -137,7 +137,7 @@
return RecoDecay::sqrtSumOfSquares(lambdaPx, lambdaPy);
});

DECLARE_SOA_DYNAMIC_COLUMN(LambdaP, lambdap, //! Total momentum in GeV/c

Check failure on line 140 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
[](float lambdaPx, float lambdaPy, float lambdaPz) -> float {
return RecoDecay::sqrtSumOfSquares(lambdaPx, lambdaPy, lambdaPz);
});
Expand Down Expand Up @@ -196,21 +196,21 @@
namespace kstarCore
{
DECLARE_SOA_COLUMN(PhotonV0ID, photonV0ID, int);
DECLARE_SOA_COLUMN(KShortV0ID, kshortV0ID, int);

Check failure on line 199 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.

DECLARE_SOA_COLUMN(X, x, float);
DECLARE_SOA_COLUMN(Y, y, float);
DECLARE_SOA_COLUMN(Z, z, float);
DECLARE_SOA_COLUMN(DCADaughters, dcadaughters, float);

Check failure on line 204 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.

DECLARE_SOA_COLUMN(PhotonPx, photonPx, float);
DECLARE_SOA_COLUMN(PhotonPy, photonPy, float);
DECLARE_SOA_COLUMN(PhotonPz, photonPz, float);
DECLARE_SOA_COLUMN(PhotonMass, photonMass, float);

DECLARE_SOA_COLUMN(KShortPx, kshortPx, float);

Check failure on line 211 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
DECLARE_SOA_COLUMN(KShortPy, kshortPy, float);

Check failure on line 212 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
DECLARE_SOA_COLUMN(KShortPz, kshortPz, float);

Check failure on line 213 in PWGLF/DataModel/LFSigmaTables.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
DECLARE_SOA_COLUMN(KShortMass, kshortMass, float);

// Dynamic Columns
Expand Down Expand Up @@ -850,7 +850,7 @@

DECLARE_SOA_DYNAMIC_COLUMN(KStarMCY, kstarMCY,
[](float photonMCPx, float photonMCPy, float photonMCPz, float kshortMCPx, float kshortMCPy, float kshortMCPz) -> float {
return RecoDecay::y(std::array{photonMCPx + kshortMCPx, photonMCPy + kshortMCPy, photonMCPz + kshortMCPz}, o2::constants::physics::MassKaonNeutral);
return RecoDecay::y(std::array{photonMCPx + kshortMCPx, photonMCPy + kshortMCPy, photonMCPz + kshortMCPz}, o2::constants::physics::MassK0Star892);
});

DECLARE_SOA_DYNAMIC_COLUMN(MCPhi, mcphi, //! Phi in the range [0, 2pi)
Expand Down
14 changes: 7 additions & 7 deletions PWGLF/TableProducer/Strangeness/sigma0builder.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1620,13 +1620,13 @@ struct sigma0builder {

// Fill with properties
GenInfo.IsPrimary = mcParticle.isPhysicalPrimary();
GenInfo.IsV0Lambda = mcParticle.pdgCode() == PDG_t::kLambda0; // 3122
GenInfo.IsV0AntiLambda = mcParticle.pdgCode() == PDG_t::kLambda0Bar; //-3122
GenInfo.IsV0KShort = mcParticle.pdgCode() == PDG_t::kK0Short; // 310
GenInfo.IsPi0 = mcParticle.pdgCode() == PDG_t::kPi0; // 111;
GenInfo.IsSigma0 = mcParticle.pdgCode() == PDG_t::kSigma0; // PDG_t::kSigma0
GenInfo.IsAntiSigma0 = mcParticle.pdgCode() == PDG_t::kSigma0Bar; //-3212
GenInfo.IsKStar = mcParticle.pdgCode() == o2::constants::physics::Pdg::kK0Star892; // 313;
GenInfo.IsV0Lambda = mcParticle.pdgCode() == PDG_t::kLambda0; // 3122
GenInfo.IsV0AntiLambda = mcParticle.pdgCode() == PDG_t::kLambda0Bar; //-3122
GenInfo.IsV0KShort = std::abs(mcParticle.pdgCode()) == PDG_t::kK0Short; // 310
GenInfo.IsPi0 = mcParticle.pdgCode() == PDG_t::kPi0; // 111;
GenInfo.IsSigma0 = mcParticle.pdgCode() == PDG_t::kSigma0; // PDG_t::kSigma0
GenInfo.IsAntiSigma0 = mcParticle.pdgCode() == PDG_t::kSigma0Bar; //-3212
GenInfo.IsKStar = std::abs(mcParticle.pdgCode()) == o2::constants::physics::Pdg::kK0Star892; // 313;
GenInfo.IsProducedByGenerator = mcParticle.producedByGenerator();
GenInfo.MCProcess = mcParticle.getProcess();
GenInfo.MCPt = mcParticle.pt();
Expand Down
12 changes: 7 additions & 5 deletions PWGLF/Tasks/Resonances/k892hadronphoton.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -809,6 +809,8 @@ struct k892hadronphoton {
int trkCode = 10; // 1: TPC-only, 2: TPC+Something, 3: ITS-Only, 4: ITS+TPC + Something, 10: anything else
auto photonPosTrackCode = kstar.photonPosTrackCode();
auto photonNegTrackCode = kstar.photonNegTrackCode();
auto kshortPosTrackCode = kstar.kshortPosTrackCode();
auto kshortNegTrackCode = kstar.kshortNegTrackCode();

if (isGamma) {
if (photonPosTrackCode == DauTPCOnly && photonNegTrackCode == DauTPCOnly)
Expand All @@ -820,13 +822,13 @@ struct k892hadronphoton {
if (photonPosTrackCode == DauITSTrackerOnly || photonNegTrackCode == DauITSTrackerOnly)
trkCode = 4;
} else {
if (photonPosTrackCode == DauTPCOnly && photonNegTrackCode == DauTPCOnly)
if (kshortPosTrackCode == DauTPCOnly && kshortNegTrackCode == DauTPCOnly)
trkCode = 1;
if ((photonPosTrackCode != DauTPCOnly && photonNegTrackCode == DauTPCOnly) || (photonPosTrackCode == DauTPCOnly && photonNegTrackCode != DauTPCOnly))
if ((kshortPosTrackCode != DauTPCOnly && kshortNegTrackCode == DauTPCOnly) || (kshortPosTrackCode == DauTPCOnly && kshortNegTrackCode != DauTPCOnly))
trkCode = 2;
if (photonPosTrackCode == DauITSTPC && photonNegTrackCode == DauITSTPC)
if (kshortPosTrackCode == DauITSTPC && kshortNegTrackCode == DauITSTPC)
trkCode = 3;
if (photonPosTrackCode == DauITSTrackerOnly || photonNegTrackCode == DauITSTrackerOnly)
if (kshortPosTrackCode == DauITSTrackerOnly || kshortNegTrackCode == DauITSTrackerOnly)
trkCode = 4;
}

Expand Down Expand Up @@ -899,7 +901,7 @@ struct k892hadronphoton {
if (std::abs(kshortPDGCodeGrandMother) == o2::constants::physics::Pdg::kK0Star892 && std::abs(photonPDGCodeGrandMother) == o2::constants::physics::Pdg::kK0Star892 && (photonGlobalIndexGrandMother == kshortGlobalIndexGrandMother)) // K*(892)0
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_KStarPionKaon"), kstarpT, kstarMass);

if (photonGlobalIndexGrandMother == kshortGlobalIndexGrandMother) {
if ((photonGlobalIndexGrandMother == kshortGlobalIndexGrandMother) && (photonPDGCodeGrandMother != o2::constants::physics::Pdg::kK0Star892) && (photonPDGCodeGrandMother != PDG_t::kProton)) {
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_PionKaon"), kstarpT, kstarMass);
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPionKaonMother"), kshortPDGCodeGrandMother, photonPDGCodeGrandMother);
}
Expand Down
Loading
Loading