Skip to content

[Trigger,PWGLF] Add sigma+ - proton trigger #14825

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
dsekihat merged 3 commits into from
Feb 6, 2026
+166 −11
Merged

[Trigger,PWGLF] Add sigma+ - proton trigger #14825

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
383c54e
Add sigma-proton trigger
Feb 5, 2026
305ff1a
Please consider the following formatting changes
alibuild Feb 5, 2026
0bb90a6
Merge pull request #40 from alibuild/alibot-cleanup-14825
fmazzasc Feb 5, 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
174 changes: 164 additions & 10 deletions EventFiltering/PWGLF/strangenessFilter.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@

#include "../filterTables.h"

#include "PWGLF/DataModel/LFKinkDecayTables.h"
#include "PWGLF/DataModel/LFParticleIdentification.h"
#include "PWGLF/DataModel/LFStrangenessTables.h"
#include "PWGLF/Utils/strangenessBuilderHelper.h"
Expand All @@ -26,6 +27,7 @@
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Multiplicity.h"
#include "Common/DataModel/PIDResponseTOF.h"
#include "Common/DataModel/TrackSelectionTables.h"

#include "CCDB/BasicCCDBManager.h"
Expand All @@ -42,6 +44,8 @@
#include "ReconstructionDataFormats/TrackParametrization.h"

#include "TVector3.h"
#include <Math/GenVector/Boost.h>
#include <TLorentzVector.h>

#include <cmath>

Expand Down Expand Up @@ -78,8 +82,10 @@ struct strangenessFilter {
HistogramRegistry QAHistosTopologicalVariables{"QAHistosTopologicalVariables", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
HistogramRegistry QAHistosTriggerParticles{"QAHistosTriggerParticles", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
HistogramRegistry QAHistosStrangenessTracking{"QAHistosStrangenessTracking", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
HistogramRegistry QAHistosSigma{"QAHistosSigma", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};

HistogramRegistry EventsvsMultiplicity{"EventsvsMultiplicity", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
OutputObj<TH1D> hProcessedEvents{TH1D("hProcessedEvents", "Strangeness - event filtered;; Number of events", 20, -1., 19.)};
OutputObj<TH1D> hProcessedEvents{TH1D("hProcessedEvents", "Strangeness - event filtered;; Number of events", 21, -1., 20.)};
OutputObj<TH1F> hCandidate{TH1F("hCandidate", "; Candidate pass selection; Number of events", 30, 0., 30.)};
OutputObj<TH1F> hEvtvshMinPt{TH1F("hEvtvshMinPt", " Number of h-Omega events with pT_h higher than thrd; min p_{T, trigg} (GeV/c); Number of events", 11, 0., 11.)};

Expand Down Expand Up @@ -164,9 +170,11 @@ struct strangenessFilter {
Configurable<float> massWindowXiNsigma{"massWindowXiNsigma", 6, "Inv. mass window for tracked Xi"};

// Selections criteria for tracks
Configurable<float> hEta{"hEta", 0.9f, "Eta range for trigger particles"};
Configurable<float> hMinPt{"hMinPt", 1.0f, "Min pt for trigger particles"};
Configurable<bool> isTrackFilter{"isTrackFilter", true, "Apply track myTrackSelections"};
struct : ConfigurableGroup {
Configurable<float> hEta{"hEta", 0.9f, "Eta range for trigger particles"};
Configurable<float> hMinPt{"hMinPt", 1.0f, "Min pt for trigger particles"};
Configurable<bool> isTrackFilter{"isTrackFilter", true, "Apply track myTrackSelections"};
} cfgTrackCuts;

// Settings for strangeness tracking filter
Configurable<std::string> ccdbUrl{"ccdbUrl", "http://alice-ccdb.cern.ch", "url of the ccdb repository"};
Expand Down Expand Up @@ -196,6 +204,23 @@ struct strangenessFilter {
Configurable<float> maxCpaTrackedXi{"maxCpaTrackedXi", 1., "Maximum CPA for tracked cascades"};
Configurable<float> minDcaTrackedOmega{"minDcaTrackedOmega", -1., "Minimum DCA for tracked cascades (ST)"};
Configurable<float> maxCpaTrackedOmega{"maxCpaTrackedOmega", 1., "Maximum CPA for tracked cascades (ST)"};

// Settings for sigmaplus filter
struct : ConfigurableGroup {
Configurable<float> nsigmatpcSigma{"nsigmatpcSigma", 3.f, "N Sigmas TPC Sigma"};
Configurable<float> nsigmatofSigma{"nsigmatofSigma", 3.f, "N Sigmas TOF Sigma"};
Configurable<float> minMassSigma{"minMassSigma", 1.15f, "min mass for Sigma"};
Configurable<float> maxMassSigma{"maxMassSigma", 1.25f, "max mass for Sigma"};
Configurable<float> minPtSigma{"minPtSigma", 1.2f, "min Pt for Sigma"};
Configurable<float> minQtAPSigma{"minQtAPSigma", 0.15f, "min QtAP for Sigma"};
Configurable<float> maxQtAPSigma{"maxQtAPSigma", 0.20f, "max QtAP for Sigma"};
Configurable<float> maxDCAtoPVSigma{"maxDCAtoPVSigma", 0.1f, "Max DCA to primary vertex for Sigma candidates (cm)"};
Configurable<float> minRadiusSigma{"minRadiusSigma", 19.f, "Min radius for Sigma+ decay vertex (cm)"};
Configurable<float> minCosPASigma{"minCosPASigma", 0.995f, "Min Cosine of pointing angle for Sigma candidates"};
Configurable<float> minPtProtonTOF{"minPtProtonTOF", 0.75f, "Min Pt for proton to have TOF signal (GeV/c)"};
Configurable<float> maxKStarSigmaProton{"maxKStarSigmaProton", 0.8f, "Max k* for Sigma-Proton pairs (GeV/c)"};
} cfgSigma;

Configurable<LabeledArray<double>> parSigmaMass{
"parSigmaMass",
{stfilter::massSigmaParameters[0], 4, 2,
Expand All @@ -215,7 +240,7 @@ struct strangenessFilter {

bool selectTrack(const auto& track)
{
return track.pt() > hMinPt && std::abs(track.eta()) < hEta && track.tpcNClsCrossedRows() >= tpcmincrossedrows && track.tpcChi2NCl() <= 4.f && track.itsChi2NCl() <= 36.f && (track.itsClusterMap() & 0x7) != 0;
return track.pt() > cfgTrackCuts.hMinPt && std::abs(track.eta()) < cfgTrackCuts.hEta && track.tpcNClsCrossedRows() >= tpcmincrossedrows && track.tpcChi2NCl() <= 4.f && track.itsChi2NCl() <= 36.f && (track.itsClusterMap() & 0x7) != 0;
}
bool selectTrackOHM(const auto& track)
{
Expand Down Expand Up @@ -296,6 +321,48 @@ struct strangenessFilter {
return true;
}

float getAlphaAP(const std::array<float, 3>& momMother, const std::array<float, 3>& momKink)
{
std::array<float, 3> momMissing = {momMother[0] - momKink[0], momMother[1] - momKink[1], momMother[2] - momKink[2]};
float lQlP = std::inner_product(momMother.begin(), momMother.end(), momKink.begin(), 0.f);
float lQlN = std::inner_product(momMother.begin(), momMother.end(), momMissing.begin(), 0.f);
return (lQlP - lQlN) / (lQlP + lQlN);
}

float getQtAP(const std::array<float, 3>& momMother, const std::array<float, 3>& momKink)
{
float dp = std::inner_product(momMother.begin(), momMother.end(), momKink.begin(), 0.f);
float p2V0 = std::inner_product(momMother.begin(), momMother.end(), momMother.begin(), 0.f);
float p2A = std::inner_product(momKink.begin(), momKink.end(), momKink.begin(), 0.f);
return std::sqrt(p2A - dp * dp / p2V0);
}

float getCosPA(const std::array<float, 3>& momMother, const std::array<float, 3>& decayVertex, const std::array<float, 3>& primaryVertex)
{
std::array<float, 3> decayVec = {decayVertex[0] - primaryVertex[0], decayVertex[1] - primaryVertex[1], decayVertex[2] - primaryVertex[2]};
float dotProduct = std::inner_product(momMother.begin(), momMother.end(), decayVec.begin(), 0.f);
float momMotherMag = std::sqrt(std::inner_product(momMother.begin(), momMother.end(), momMother.begin(), 0.f));
float decayVecMag = std::sqrt(std::inner_product(decayVec.begin(), decayVec.end(), decayVec.begin(), 0.f));
return dotProduct / (momMotherMag * decayVecMag);
}

float getKStar(TLorentzVector const& part1, TLorentzVector const& part2)
{
TLorentzVector trackSum, PartOneCMS, PartTwoCMS, trackRelK;
trackSum = part1 + part2;
const float beta = trackSum.Beta();
const float betax = beta * std::cos(trackSum.Phi()) * std::sin(trackSum.Theta());
const float betay = beta * std::sin(trackSum.Phi()) * std::sin(trackSum.Theta());
const float betaz = beta * std::cos(trackSum.Theta());
PartOneCMS.SetXYZM(part1.Px(), part1.Py(), part1.Pz(), part1.M());
PartTwoCMS.SetXYZM(part2.Px(), part2.Py(), part2.Pz(), part2.M());
const ROOT::Math::Boost boostPRF = ROOT::Math::Boost(-betax, -betay, -betaz);
PartOneCMS = boostPRF(PartOneCMS);
PartTwoCMS = boostPRF(PartTwoCMS);
trackRelK = PartOneCMS - PartTwoCMS;
return 0.5 * trackRelK.P();
}

void init(o2::framework::InitContext&)
{
// set V0 parameters in the helper
Expand Down Expand Up @@ -345,6 +412,7 @@ struct strangenessFilter {
hProcessedEvents->GetXaxis()->SetBinLabel(18, aod::filtering::OmegaHighMultTrk::columnLabel());
hProcessedEvents->GetXaxis()->SetBinLabel(19, aod::filtering::HighMultFT0M::columnLabel());
hProcessedEvents->GetXaxis()->SetBinLabel(20, aod::filtering::HighMultTrk::columnLabel());
hProcessedEvents->GetXaxis()->SetBinLabel(21, aod::filtering::SigmaProton::columnLabel());

hCandidate->GetXaxis()->SetBinLabel(1, "All");
hCandidate->GetXaxis()->SetBinLabel(2, "PassBuilderSel");
Expand Down Expand Up @@ -530,6 +598,13 @@ struct strangenessFilter {
QAHistosStrangenessTracking.add("hPtVsMassTrkOmega", "cascades;p_{T} (GeV/#it{c});m (GeV/#it{c}^2)", HistType::kTH2D, {{200, 0., 10.}, {1000, 1.6, 2.1}});
QAHistosStrangenessTracking.add("hPtVsMassTrkXiSelected", "cascades;p_{T} (GeV/#it{c});m (GeV/#it{c}^2)", HistType::kTH2D, {{200, 0., 10.}, {1000, 1.2, 1.7}});
QAHistosStrangenessTracking.add("hPtVsMassTrkOmegaSelected", "cascades;p_{T} (GeV/#it{c});m (GeV/#it{c}^2)", HistType::kTH2D, {{200, 0., 10.}, {1000, 1.6, 2.1}});

// Sigma QA histograms
QAHistosSigma.add("hPtVsMassSigmaPlus", ";p_{T} (GeV/#it{c});m (GeV/#it{c}^2)", HistType::kTH2D, {{20, -5, 5.}, {50, 1.1, 1.3}});
QAHistosSigma.add("hDecayRadiusSigma", ";R (cm)", HistType::kTH1D, {{100, 0., 40.}});
QAHistosSigma.add("hPtNSigmaTPCPrPair", ";p_{T} (GeV/#it{c});#sigma_{TPC}", HistType::kTH2D, {{200, -5., 5.}, {200, -5., 5.}});
QAHistosSigma.add("hPtNSigmaTOFPrPair", ";p_{T} (GeV/#it{c});#sigma_{TOF}", HistType::kTH2D, {{200, -5., 5.}, {200, -5., 5.}});
QAHistosSigma.add("hKStarSigmaPr", ";k*", HistType::kTH1D, {{200, 0., 2}});
}
}

Expand Down Expand Up @@ -566,7 +641,7 @@ struct strangenessFilter {

// Tables
using CollisionCandidates = soa::Join<aod::Collisions, aod::EvSels, aod::MultZeqs, aod::FT0Mults, aod::MultsGlobal>::iterator;
using TrackCandidates = soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::pidTPCLfFullPi, aod::pidTPCLfFullPr, aod::pidTPCLfFullKa>;
using TrackCandidates = soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::pidTPCLfFullPi, aod::pidTPCLfFullPr, aod::pidTPCLfFullKa, aod::pidTOFFullPr>;

float getMassWindow(const stfilter::species s, const float pt, const float nsigma = 6)
{
Expand All @@ -580,14 +655,14 @@ struct strangenessFilter {

void fillTriggerTable(bool keepEvent[])
{
strgtable(keepEvent[0], keepEvent[1], keepEvent[2], keepEvent[3], keepEvent[4], keepEvent[5], keepEvent[6], keepEvent[7], keepEvent[8], keepEvent[9], keepEvent[10], keepEvent[11], keepEvent[12], keepEvent[13], keepEvent[14], keepEvent[15]);
strgtable(keepEvent[0], keepEvent[1], keepEvent[2], keepEvent[3], keepEvent[4], keepEvent[5], keepEvent[6], keepEvent[7], keepEvent[8], keepEvent[9], keepEvent[10], keepEvent[11], keepEvent[12], keepEvent[13], keepEvent[14], keepEvent[15], keepEvent[16]);
}

void process(CollisionCandidates const& collision, TrackCandidates const& tracks, aod::Cascades const& cascadesBase, aod::AssignedTrackedCascades const& trackedCascades, aod::AssignedTrackedV0s const& /*trackedV0s*/, aod::AssignedTracked3Bodys const& /*tracked3Bodys*/, aod::V0s const& v0Base, aod::BCs const&, aod::FT0s const& /*ft0s*/)
void process(CollisionCandidates const& collision, TrackCandidates const& tracks, aod::Cascades const& cascadesBase, aod::AssignedTrackedCascades const& trackedCascades, aod::KinkCands const& kinkCands, aod::AssignedTrackedV0s const& /*trackedV0s*/, aod::AssignedTracked3Bodys const& /*tracked3Bodys*/, aod::V0s const& v0Base, aod::BCs const&, aod::FT0s const& /*ft0s*/)
{
// Is event good? [0] = Omega, [1] = high-pT hadron + Omega, [2] = 2Xi, [3] = 3Xi, [4] = 4Xi, [5] single-Xi, [6] Omega with high radius
// [7] tracked Xi, [8] tracked Omega, [9] Omega + high mult event
bool keepEvent[16]{}; // explicitly zero-initialised
bool keepEvent[17]{}; // explicitly zero-initialised
std::vector<std::array<int64_t, 2>> v0sFromOmegaID;
std::vector<std::array<int64_t, 2>> v0sFromXiID;

Expand Down Expand Up @@ -1109,7 +1184,7 @@ struct strangenessFilter {
// QA tracks
int triggcounterAllEv = 0;
for (auto track : tracks) { // start loop over tracks
if (isTrackFilter && !selectTrack(track)) {
if (cfgTrackCuts.isTrackFilter && !selectTrack(track)) {
continue;
}
triggcounterAllEv++;
Expand Down Expand Up @@ -1373,6 +1448,82 @@ struct strangenessFilter {
}
}
}
// // Sigma - proton trigger definition
for (const auto& kinkCand : kinkCands) {
auto dauTrack = kinkCand.trackDaug_as<TrackCandidates>();
if (!dauTrack.hasTPC() || !dauTrack.hasTOF()) {
continue;
}
if (std::abs(dauTrack.tpcNSigmaPr()) > cfgSigma.nsigmatpcSigma || std::abs(dauTrack.tofNSigmaPr()) > cfgSigma.nsigmatofSigma) {
continue;
}
if (kinkCand.ptMoth() < cfgSigma.minPtSigma) {
continue;
}
if (kinkCand.mSigmaPlus() < cfgSigma.minMassSigma || kinkCand.mSigmaPlus() > cfgSigma.maxMassSigma) {
continue;
}
QAHistosSigma.fill(HIST("hPtVsMassSigmaPlus"), kinkCand.mothSign() * kinkCand.ptMoth(), kinkCand.mSigmaPlus());
std::array<float, 3> momMoth = {kinkCand.pxMoth(), kinkCand.pyMoth(), kinkCand.pzMoth()};
std::array<float, 3> momDaug = {kinkCand.pxDaug(), kinkCand.pyDaug(), kinkCand.pzDaug()};
std::array<float, 3> primaryVtx = {collision.posX(), collision.posY(), collision.posZ()};
std::array<float, 3> decayVtx = {kinkCand.xDecVtx(), kinkCand.yDecVtx(), kinkCand.zDecVtx()};
float qtAP = getQtAP(momMoth, momDaug);
float alphaAP = getAlphaAP(momMoth, momDaug);
float cosPA = getCosPA(momMoth, decayVtx, primaryVtx);
if (alphaAP < 0) {
continue;
}
if (qtAP < cfgSigma.minQtAPSigma || qtAP > cfgSigma.maxQtAPSigma) {
continue;
}
if (cosPA < cfgSigma.minCosPASigma) {
continue;
}
if (std::abs(kinkCand.dcaMothPv()) > cfgSigma.maxDCAtoPVSigma) {
continue;
}
float decRad = std::hypot(kinkCand.xDecVtx(), kinkCand.yDecVtx());
if (decRad < cfgSigma.minRadiusSigma) {
continue;
}
QAHistosSigma.fill(HIST("hDecayRadiusSigma"), decRad);
// pair a proton
bool isProtonPaired = false;
for (auto track : tracks) {
if (track.globalIndex() == dauTrack.globalIndex()) {
continue;
}
if (std::abs(track.tpcNSigmaPr()) > cfgSigma.nsigmatpcSigma) {
continue;
}
QAHistosSigma.fill(HIST("hPtNSigmaTPCPrPair"), track.sign() * track.pt(), track.tpcNSigmaPr());
TLorentzVector sigmaVec, protonVec;
sigmaVec.SetXYZM(kinkCand.pxMoth(), kinkCand.pyMoth(), kinkCand.pzMoth(), o2::constants::physics::MassSigmaPlus);
protonVec.SetXYZM(track.px(), track.py(), track.pz(), o2::constants::physics::MassProton);
float kstar = getKStar(sigmaVec, protonVec);
if (kstar > cfgSigma.maxKStarSigmaProton) {
continue;
}
QAHistosSigma.fill(HIST("hKStarSigmaPr"), kstar);
if (track.pt() < cfgSigma.minPtProtonTOF) {
isProtonPaired = true;
break;
}
if (!track.hasTOF()) {
continue;
}
QAHistosSigma.fill(HIST("hPtNSigmaTOFPrPair"), track.sign() * track.pt(), track.tofNSigmaPr());
if (std::abs(track.tofNSigmaPr()) > cfgSigma.nsigmatofSigma) {
continue;
}
isProtonPaired = true;
break;
}
if (isProtonPaired) {
keepEvent[16] = true;
}
}

// Fill centrality dependent histos
if (keepEvent[0]) {
Expand Down Expand Up @@ -1426,6 +1577,9 @@ struct strangenessFilter {
if (keepEvent[15]) {
hProcessedEvents->Fill(18.5);
}
if (keepEvent[16]) {
hProcessedEvents->Fill(19.5);
}
// Filling the table
fillTriggerTable(keepEvent);
}
Expand Down
3 changes: 2 additions & 1 deletion EventFiltering/filterTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -174,6 +174,7 @@ DECLARE_SOA_COLUMN(LambdaLambda, lambdaLambda, bool); //! at least 2
DECLARE_SOA_COLUMN(OmegaHighMultTrk, hasOmegaHighMultTrk, bool); //! at least 1 Omega + high-mult track event
DECLARE_SOA_COLUMN(HighMultFT0M, hasHighMultFT0M, bool); //! at least 1 Omega + high-mult track event
DECLARE_SOA_COLUMN(HighMultTrk, hasHighMultTrk, bool); //! at least 1 Omega + high-mult track event
DECLARE_SOA_COLUMN(SigmaProton, hasSigmaProton, bool); //! at least 1 Sigma - proton candidate
// F1-proton
DECLARE_SOA_COLUMN(TriggerEventF1Proton, triggereventf1proton, bool); //! F1 - proton femto trigger event

Expand Down Expand Up @@ -334,7 +335,7 @@ using FullJetFilter = FullJetFilters::iterator;

// strangeness (lf)
DECLARE_SOA_TABLE(StrangenessFilters, "AOD", "LFStrgFilters", //!
filtering::Omega, filtering::hadronOmega, filtering::DoubleXi, filtering::TripleXi, filtering::QuadrupleXi, filtering::SingleXiYN, filtering::OmegaLargeRadius, filtering::TrackedXi, filtering::TrackedOmega, filtering::OmegaHighMult, filtering::DoubleOmega, filtering::OmegaXi, filtering::LambdaLambda, filtering::OmegaHighMultTrk, filtering::HighMultFT0M, filtering::HighMultTrk);
filtering::Omega, filtering::hadronOmega, filtering::DoubleXi, filtering::TripleXi, filtering::QuadrupleXi, filtering::SingleXiYN, filtering::OmegaLargeRadius, filtering::TrackedXi, filtering::TrackedOmega, filtering::OmegaHighMult, filtering::DoubleOmega, filtering::OmegaXi, filtering::LambdaLambda, filtering::OmegaHighMultTrk, filtering::HighMultFT0M, filtering::HighMultTrk, filtering::SigmaProton);
using StrangenessFilter = StrangenessFilters::iterator;

// F1 proton
Expand Down
Loading