Skip to content

[PWGLF] Add SigmaPlus PCM tables #17367

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
romainschotter merged 1 commit into from
Aug 6, 2026
+67 −14
Merged

[PWGLF] Add SigmaPlus PCM tables #17367

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
34 changes: 32 additions & 2 deletions PWGLF/DataModel/LFKinkDecayTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -166,7 +166,9 @@ DECLARE_SOA_COLUMN(XDecVtx, xDecVtx, float); //! Decay vertex of t
DECLARE_SOA_COLUMN(YDecVtx, yDecVtx, float); //! Decay vertex of the candidate (y direction)
DECLARE_SOA_COLUMN(ZDecVtx, zDecVtx, float); //! Decay vertex of the candidate (z direction)
DECLARE_SOA_COLUMN(Radius, radius, float); //! Decay radius of the candidate (cm)
DECLARE_SOA_COLUMN(FlightDistance, flightDistance, float); //! Flight distance of the candidate (PV to decay vertex, cm)
DECLARE_SOA_COLUMN(DcaProtonGamma, dcaProtonGamma, float); //! DCA between proton and photon at the fitted vertex (cm)
DECLARE_SOA_COLUMN(Chi2, chi2, float); //! chi2 of the proton-photon vertex fit

DECLARE_SOA_COLUMN(PxProton, pxProton, float); //! Px of the proton
DECLARE_SOA_COLUMN(PyProton, pyProton, float); //! Py of the proton
Expand All @@ -183,6 +185,24 @@ DECLARE_SOA_COLUMN(NSigmaTOFProton, nSigmaTOFProton, float); //! TOF nSigma of t
DECLARE_SOA_COLUMN(NSigmaTPCElPos, nSigmaTPCElPos, float); //! TPC nSigma_el of the photon's positive daughter
DECLARE_SOA_COLUMN(NSigmaTPCElNeg, nSigmaTPCElNeg, float); //! TPC nSigma_el of the photon's negative daughter

DECLARE_SOA_COLUMN(PhotonMass, photonMass, float); //! Invariant mass of the measured photon (V0) candidate (GeV/c^2)
DECLARE_SOA_COLUMN(PhotonAlpha, photonAlpha, float); //! Armenteros-Podolanski alpha of the measured photon
DECLARE_SOA_COLUMN(PhotonQt, photonQt, float); //! Armenteros-Podolanski qT of the measured photon (GeV/c)
DECLARE_SOA_COLUMN(PhotonConvRadius, photonConvRadius, float); //! Conversion radius of the measured photon (cm)
DECLARE_SOA_COLUMN(PhotonOpeningAngle, photonOpeningAngle, float); //! Opening angle between the photon's e+e- daughters (rad)
DECLARE_SOA_COLUMN(PhotonPointingAngle, photonPointingAngle, float); //! Angle between the photon momentum and the line from its conversion point to the candidate decay vertex (rad)
DECLARE_SOA_COLUMN(PhotonDcaToPV, photonDcaToPV, float); //! DCA of the photon's flight line to the primary vertex (cm)

DECLARE_SOA_COLUMN(ProtonItsNCls, protonItsNCls, uint8_t); //! Number of ITS clusters of the proton track
DECLARE_SOA_COLUMN(ProtonTpcNCls, protonTpcNCls, int16_t); //! Number of found TPC clusters of the proton track
DECLARE_SOA_COLUMN(ProtonDcaXY, protonDcaXY, float); //! DCA of the proton track to the primary vertex, xy (cm)
DECLARE_SOA_COLUMN(ProtonDcaZ, protonDcaZ, float); //! DCA of the proton track to the primary vertex, z (cm)

DECLARE_SOA_COLUMN(PhotonPosItsNCls, photonPosItsNCls, uint8_t); //! Number of ITS clusters of the photon's positive daughter
DECLARE_SOA_COLUMN(PhotonPosTpcNCls, photonPosTpcNCls, int16_t); //! Number of found TPC clusters of the photon's positive daughter
DECLARE_SOA_COLUMN(PhotonNegItsNCls, photonNegItsNCls, uint8_t); //! Number of ITS clusters of the photon's negative daughter
DECLARE_SOA_COLUMN(PhotonNegTpcNCls, photonNegTpcNCls, int16_t); //! Number of found TPC clusters of the photon's negative daughter

// MC columns
DECLARE_SOA_COLUMN(IsSignal, isSignal, bool); //! True if proton and photon are MC-truth matched to the same Sigma+

Expand Down Expand Up @@ -228,12 +248,17 @@ DECLARE_SOA_DYNAMIC_COLUMN(MassSigmaPlus, massSigmaPlus, //! Invariant mass of t

DECLARE_SOA_TABLE(SigmaPlusCands, "AOD", "SIGMAPLUSCANDS",
sigmapluscand::XDecVtx, sigmapluscand::YDecVtx, sigmapluscand::ZDecVtx,
sigmapluscand::Radius, sigmapluscand::DcaProtonGamma,
sigmapluscand::Radius, sigmapluscand::FlightDistance,
sigmapluscand::DcaProtonGamma, sigmapluscand::Chi2,
sigmapluscand::PxProton, sigmapluscand::PyProton, sigmapluscand::PzProton,
sigmapluscand::PxGamma1, sigmapluscand::PyGamma1, sigmapluscand::PzGamma1,
sigmapluscand::PxGamma2, sigmapluscand::PyGamma2, sigmapluscand::PzGamma2,
sigmapluscand::NSigmaTPCProton, sigmapluscand::NSigmaTOFProton,
sigmapluscand::NSigmaTPCElPos, sigmapluscand::NSigmaTPCElNeg,
sigmapluscand::PhotonMass, sigmapluscand::PhotonAlpha, sigmapluscand::PhotonQt, sigmapluscand::PhotonConvRadius,
sigmapluscand::PhotonOpeningAngle, sigmapluscand::PhotonPointingAngle, sigmapluscand::PhotonDcaToPV,
sigmapluscand::ProtonItsNCls, sigmapluscand::ProtonTpcNCls, sigmapluscand::ProtonDcaXY, sigmapluscand::ProtonDcaZ,
sigmapluscand::PhotonPosItsNCls, sigmapluscand::PhotonPosTpcNCls, sigmapluscand::PhotonNegItsNCls, sigmapluscand::PhotonNegTpcNCls,

// dynamic columns
sigmapluscand::PxSigmaPlus<sigmapluscand::PxProton, sigmapluscand::PxGamma1, sigmapluscand::PxGamma2>,
Expand All @@ -244,12 +269,17 @@ DECLARE_SOA_TABLE(SigmaPlusCands, "AOD", "SIGMAPLUSCANDS",

DECLARE_SOA_TABLE(SigmaPlusCandsMC, "AOD", "SIGMAPLUSMC",
sigmapluscand::XDecVtx, sigmapluscand::YDecVtx, sigmapluscand::ZDecVtx,
sigmapluscand::Radius, sigmapluscand::DcaProtonGamma,
sigmapluscand::Radius, sigmapluscand::FlightDistance,
sigmapluscand::DcaProtonGamma, sigmapluscand::Chi2,
sigmapluscand::PxProton, sigmapluscand::PyProton, sigmapluscand::PzProton,
sigmapluscand::PxGamma1, sigmapluscand::PyGamma1, sigmapluscand::PzGamma1,
sigmapluscand::PxGamma2, sigmapluscand::PyGamma2, sigmapluscand::PzGamma2,
sigmapluscand::NSigmaTPCProton, sigmapluscand::NSigmaTOFProton,
sigmapluscand::NSigmaTPCElPos, sigmapluscand::NSigmaTPCElNeg,
sigmapluscand::PhotonMass, sigmapluscand::PhotonAlpha, sigmapluscand::PhotonQt, sigmapluscand::PhotonConvRadius,
sigmapluscand::PhotonOpeningAngle, sigmapluscand::PhotonPointingAngle, sigmapluscand::PhotonDcaToPV,
sigmapluscand::ProtonItsNCls, sigmapluscand::ProtonTpcNCls, sigmapluscand::ProtonDcaXY, sigmapluscand::ProtonDcaZ,
sigmapluscand::PhotonPosItsNCls, sigmapluscand::PhotonPosTpcNCls, sigmapluscand::PhotonNegItsNCls, sigmapluscand::PhotonNegTpcNCls,
sigmapluscand::IsSignal,
sigmapluscand::ProtonPdgCode, sigmapluscand::ProtonMotherPdgCode,
sigmapluscand::GammaPdgCode, sigmapluscand::GammaMotherPdgCode, sigmapluscand::GammaGMotherPdgCode,
Expand Down
47 changes: 35 additions & 12 deletions PWGLF/TableProducer/Strangeness/sigmaplusbuilder.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,7 @@
#include "Common/Core/trackUtilities.h"
#include "Common/DataModel/PIDResponseTOF.h"
#include "Common/DataModel/PIDResponseTPC.h"
#include "Common/DataModel/TrackSelectionTables.h"

#include <CCDB/BasicCCDBManager.h>
#include <CommonConstants/PhysicsConstants.h>
Expand Down Expand Up @@ -47,7 +48,7 @@
using namespace o2;
using namespace o2::framework;

using TracksFull = soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksCovIU,
using TracksFull = soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::TracksDCA,
aod::pidTPCFullEl, aod::pidTPCFullPr, aod::pidTOFFullPr>;
using TracksFullMC = soa::Join<TracksFull, aod::McTrackLabels>;
using CollisionsFull = aod::Collisions;
Expand All @@ -56,7 +57,7 @@ using CollisionsFullMC = soa::Join<aod::Collisions, aod::McCollisionLabels>;
struct Sigmaplusbuilder {

// photon (PCM) selection
Configurable<float> photonMaxMass{"photonMaxMass", 0.10, "Max photon mass (GeV/c^2)"};
Configurable<float> photonMaxMass{"photonMaxMass", 0.20, "Max photon mass (GeV/c^2)"};
Configurable<float> photonMinRapidity{"photonMinRapidity", -0.8, "Min photon rapidity"};
Configurable<float> photonMaxRapidity{"photonMaxRapidity", 0.8, "Max photon rapidity"};
Configurable<float> photonDauEtaMin{"photonDauEtaMin", -0.8, "Min eta of photon daughter tracks"};
Expand All @@ -65,9 +66,9 @@ struct Sigmaplusbuilder {
Configurable<float> photonMaxRadius{"photonMaxRadius", 115., "Max photon conversion radius (cm)"};
Configurable<float> photonMinV0cospa{"photonMinV0cospa", 0.80, "Min V0 CosPA"};
Configurable<float> photonMaxDCAV0Dau{"photonMaxDCAV0Dau", 3.5, "Max DCA between photon daughters (cm)"};
Configurable<float> photonMaxQt{"photonMaxQt", 0.08, "Max Armenteros qT for photons (GeV/c)"};
Configurable<float> photonMaxQt{"photonMaxQt", 0.15, "Max Armenteros qT for photons (GeV/c)"};
Configurable<float> photonMaxAlpha{"photonMaxAlpha", 1.0, "Max |Armenteros alpha| for photons"};
Configurable<float> photonMaxTPCNSigmaEl{"photonMaxTPCNSigmaEl", 7, "Max |TPC nSigma_el| for photon daughters"};
Configurable<float> photonMaxTPCNSigmaEl{"photonMaxTPCNSigmaEl", 15, "Max |TPC nSigma_el| for photon daughters"};

// proton selection
Configurable<float> protonMinPt{"protonMinPt", 0.3, "Minimum proton pT (GeV/c)"};
Expand Down Expand Up @@ -279,10 +280,6 @@ struct Sigmaplusbuilder {
auto posTrack = v0.template posTrack_as<TTracks>();
auto negTrack = v0.template negTrack_as<TTracks>();

// true if this V0's legs genuinely share a gamma mother whose mother is a
// pi0 whose mother is a Sigma+ -- regardless of whether it passes the
// selection below. Lets us see exactly which cut, if any, throws away
// genuine Sigma+ photons, instead of only knowing the final accept/reject.
bool isSignal = false;
if constexpr (IsMC) {
if (posTrack.has_mcParticle() && negTrack.has_mcParticle()) {
Expand Down Expand Up @@ -557,7 +554,8 @@ struct Sigmaplusbuilder {
}
fillCandStep(1); // Vertex fit

float dcaProtonGamma = std::sqrt(fitter.getChi2AtPCACandidate());
float fitChi2 = fitter.getChi2AtPCACandidate();
float dcaProtonGamma = std::sqrt(fitChi2);
histos.fill(HIST("Candidate/hDcaProtonGamma"), dcaProtonGamma);
if constexpr (IsMC) {
if (isSignal) {
Expand Down Expand Up @@ -587,6 +585,7 @@ struct Sigmaplusbuilder {
// flight direction n and the decay-plane basis n, eIn, eOut
std::array<float, 3> flightVec{secVtx[0] - pv[0], secVtx[1] - pv[1], secVtx[2] - pv[2]};
std::array<float, 3> nHat = normalize3(flightVec);
float flightDistance = std::sqrt(dot3(flightVec, flightVec));

std::array<float, 3> pProton;
std::array<float, 3> pGamma1;
Expand Down Expand Up @@ -725,14 +724,34 @@ struct Sigmaplusbuilder {
}
fillCandStep(6); // filled

// photon (V0) opening angle: angle between the e+/e- daughter momenta at their own reference point
std::array<float, 3> pPosDau{posTrack.px(), posTrack.py(), posTrack.pz()};
std::array<float, 3> pNegDau{negTrack.px(), negTrack.py(), negTrack.pz()};
float photonOpeningAngle = std::acos(std::clamp(dot3(pPosDau, pNegDau) / std::sqrt(dot3(pPosDau, pPosDau) * dot3(pNegDau, pNegDau)), -1.f, 1.f));

// photon pointing angle: angle between the fitted photon momentum and the line from its conversion point to the p-gamma decay vertex
std::array<float, 3> convToDecVtx{secVtx[0] - photon.x(), secVtx[1] - photon.y(), secVtx[2] - photon.z()};
float photonPointingAngle = std::acos(std::clamp(dot3(pGamma1, convToDecVtx) / std::sqrt(dot3(pGamma1, pGamma1) * dot3(convToDecVtx, convToDecVtx)), -1.f, 1.f));

// photon DCA to PV: distance from the PV to the line through the conversion point along the photon momentum direction
std::array<float, 3> convPoint{photon.x(), photon.y(), photon.z()};
std::array<float, 3> photonDir = normalize3({photon.px(), photon.py(), photon.pz()});
std::array<float, 3> pvToConv{pv[0] - convPoint[0], pv[1] - convPoint[1], pv[2] - convPoint[2]};
std::array<float, 3> pvToConvCrossDir = cross3(pvToConv, photonDir);
float photonDcaToPV = std::sqrt(dot3(pvToConvCrossDir, pvToConvCrossDir));

if constexpr (IsMC) {
sigmaPlusCandsMC(secVtx[0], secVtx[1], secVtx[2],
radius, dcaProtonGamma,
radius, flightDistance, dcaProtonGamma, fitChi2,
pProton[0], pProton[1], pProton[2],
pGamma1[0], pGamma1[1], pGamma1[2],
bestMomGamma2[0], bestMomGamma2[1], bestMomGamma2[2],
protonTrack.tpcNSigmaPr(), protonTrack.tofNSigmaPr(),
posTrack.tpcNSigmaEl(), negTrack.tpcNSigmaEl(),
photon.mGamma(), photon.alpha(), photon.qtarm(), photon.v0radius(),
photonOpeningAngle, photonPointingAngle, photonDcaToPV,
protonTrack.itsNCls(), protonTrack.tpcNClsFound(), protonTrack.dcaXY(), protonTrack.dcaZ(),
posTrack.itsNCls(), posTrack.tpcNClsFound(), negTrack.itsNCls(), negTrack.tpcNClsFound(),
isSignal,
protonPdgCode, protonMotherPdgCode,
gammaPdgCode, gammaMotherPdgCode, gammaGMotherPdgCode,
Expand All @@ -741,12 +760,16 @@ struct Sigmaplusbuilder {
mcTrueMomGamma[0], mcTrueMomGamma[1], mcTrueMomGamma[2]);
} else {
sigmaPlusCands(secVtx[0], secVtx[1], secVtx[2],
radius, dcaProtonGamma,
radius, flightDistance, dcaProtonGamma, fitChi2,
pProton[0], pProton[1], pProton[2],
pGamma1[0], pGamma1[1], pGamma1[2],
bestMomGamma2[0], bestMomGamma2[1], bestMomGamma2[2],
protonTrack.tpcNSigmaPr(), protonTrack.tofNSigmaPr(),
posTrack.tpcNSigmaEl(), negTrack.tpcNSigmaEl());
posTrack.tpcNSigmaEl(), negTrack.tpcNSigmaEl(),
photon.mGamma(), photon.alpha(), photon.qtarm(), photon.v0radius(),
photonOpeningAngle, photonPointingAngle, photonDcaToPV,
protonTrack.itsNCls(), protonTrack.tpcNClsFound(), protonTrack.dcaXY(), protonTrack.dcaZ(),
posTrack.itsNCls(), posTrack.tpcNClsFound(), negTrack.itsNCls(), negTrack.tpcNClsFound());
}
}

Expand Down
Loading