Skip to content

[PWGLF] Local polarization of Lambda induced by jets in Run3 pp collisions at 13.6 TeV #14234

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 3 commits into from
Dec 16, 2025
Merged

[PWGLF] Local polarization of Lambda induced by jets in Run3 pp collisions at 13.6 TeV #14234

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
a8e949a
Add files via upload
15071832337 Dec 16, 2025
c1583c0
Add files via upload
15071832337 Dec 16, 2025
dcf77c7
Add files via upload
15071832337 Dec 16, 2025
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
27 changes: 19 additions & 8 deletions PWGLF/Tasks/Strangeness/lambdaJetpolarization.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,7 +69,7 @@ struct LambdaJetpolarization {
Configurable<float> etaMax{"etaMax", +0.9f, "eta max"};
Configurable<double> deltaEtaEdge{"deltaEtaEdge", 0.00, "eta gap from the edge"};
// track parameters
Configurable<float> minITSnCls{"minITSnCls", 4.0f, "min number of ITS clusters"};
Configurable<float> minITSnCls{"minITSnCls", 2.0f, "min number of ITS clusters"};
Configurable<float> minTPCnClsFound{"minTPCnClsFound", 80.0f, "min number of found TPC clusters"};
Configurable<float> minNCrossedRowsTPC{"minNCrossedRowsTPC", 80.0f, "min number of TPC crossed rows"};
Configurable<float> minTpcNcrossedRowsOverFindable{"minTpcNcrossedRowsOverFindable", 0.8, "crossed rows/findable"};
Expand All @@ -92,17 +92,18 @@ struct LambdaJetpolarization {
Configurable<double> cfgtrkMinPt{"cfgtrkMinPt", 0.10, "set track min pT"};

// v0 parameters
Configurable<float> v0Ptmin{"v0Ptmin", 0.6f, "Minimum V0 pT"};
Configurable<float> v0cospaMin{"v0cospaMin", 0.995f, "Minimum V0 CosPA"};
Configurable<float> v0cospainit{"v0cospainit", 0.97f, "Minimum V0 CosPA"};
Configurable<float> minimumV0Radius{"minimumV0Radius", 0.2f, "Minimum V0 Radius"};
Configurable<float> minimumV0Radius{"minimumV0Radius", 0.5f, "Minimum V0 Radius"};
Configurable<float> maximumV0Radius{"maximumV0Radius", 100000.0f, "Maximum V0 Radius"};
Configurable<float> dcaV0DaughtersMax{"dcaV0DaughtersMax", 1.0f, "Maximum DCA Daughters"};
Configurable<float> dcanegtoPVmin{"dcanegtoPVmin", 0.1f, "Minimum DCA Neg To PV"};
Configurable<float> dcapostoPVmin{"dcapostoPVmin", 0.1f, "Minimum DCA Pos To PV"};
Configurable<float> v0radius{"v0radius", 0.0, "Radius"};
Configurable<float> dcav0dau{"dcav0dau", 10, "DCA V0 Daughters"};
Configurable<float> dcanegtopv{"dcanegtopv", 0.0, "DCA Neg To PV"};
Configurable<float> dcapostopv{"dcapostopv", 0.0, "DCA Pos To PV"};
Configurable<float> dcav0dau{"dcav0dau", 1, "DCA V0 Daughters"};
Configurable<float> dcanegtopv{"dcanegtopv", 0.05, "DCA Neg To PV"};
Configurable<float> dcapostopv{"dcapostopv", 0.05, "DCA Pos To PV"};

// jet selection
Configurable<float> cfgjetPtMin{"cfgjetPtMin", 8.0, "minimum jet pT cut"};
Expand Down Expand Up @@ -130,7 +131,7 @@ struct LambdaJetpolarization {
Configurable<bool> ifinitpasslambda{"ifinitpasslambda", 0, "ifinitpasslambda"};
Configurable<bool> ifpasslambda{"ifpasslambda", 1, "ifpasslambda"};
Configurable<float> paramArmenterosCut{"paramArmenterosCut", 0.2, "parameter Armenteros Cut"};
Configurable<bool> doArmenterosCut{"doArmenterosCut", 0, "do Armenteros Cut"};
Configurable<bool> doArmenterosCut{"doArmenterosCut", 1, "do Armenteros Cut"};
Configurable<bool> noSameBunchPileUp{"noSameBunchPileUp", true, "reject SameBunchPileUp"};
Configurable<int> v0TypeSelection{"v0TypeSelection", 1, "select on a certain V0 type (leave negative if no selection desired)"};
Configurable<bool> notITSAfterburner{"notITSAfterburner", 0, "notITSAfterburner"};
Expand Down Expand Up @@ -303,6 +304,8 @@ struct LambdaJetpolarization {
registryData.add("hprotonThetaInJetV0", "hprotonThetaInJetV0", kTH1F, {axisTheta});

registryData.add("LambdaQA/TH2FLambdaMassPhiInJet", "TH2FLambdaMassPhiInJet", kTH2F, {{200, 0, TMath::Pi()}, {200, 0.9, 1.2}});
registryData.add("LambdaQA/hArmenterosPreAnalyserCuts", "hArmenterosPreAnalyserCuts", kTH2F, {{1000, -1.0f, 1.0f, "#alpha"}, {1000, 0.0f, 0.30f, "#it{Q}_{T}"}});
registryData.add("AntiLambdaQA/hArmenterosPreAnalyserCuts", "hArmenterosPreAnalyserCuts", kTH2F, {{1000, -1.0f, 1.0f, "#alpha"}, {1000, 0.0f, 0.30f, "#it{Q}_{T}"}});

// Lab frame measures
registryData.add("LambdaQA/TH2FprotonCosThetaInLab", "TH2FprotonCosThetaInLab", kTH2F, {{200, 0.9, 1.2}, {200, -1.0, 1.0}});
Expand Down Expand Up @@ -722,6 +725,10 @@ struct LambdaJetpolarization {
evFlag = 1;
}

if (v0.pt() < v0Ptmin) {
return false;
}

registryData.fill(HIST("QA/hv0sSelection"), 0.5);

if (evSel && evFlag < 1)
Expand Down Expand Up @@ -833,7 +840,9 @@ struct LambdaJetpolarization {
if (collision.isInelGt0()) {
evFlag = 1;
}

if (v0.pt() < v0Ptmin) {
return false;
}
if (evSel && evFlag < 1)
return false;

Expand Down Expand Up @@ -1599,7 +1608,7 @@ struct LambdaJetpolarization {
}
}
}
PROCESS_SWITCH(LambdaJetpolarization, processLongitudinalPolarization, "processLongitudinalPolarization", true);
PROCESS_SWITCH(LambdaJetpolarization, processLongitudinalPolarization, "processLongitudinalPolarization", false);

void processLambdaJetPolarization(SelV0Collisions::iterator const& collision, aod::V0Datas const& fullV0s, StrHadronDaughterTracks const& tracks)
{
Expand Down Expand Up @@ -1701,10 +1710,12 @@ struct LambdaJetpolarization {
if (registryDataAcceptV0Lambda(v0, pos, neg, collision)) {
V0Numbers = V0Numbers + 1;
registryData.fill(HIST("LambdaPtMass"), v0.pt(), v0.mLambda());
registryData.fill(HIST("LambdaQA/hArmenterosPreAnalyserCuts"), v0.alpha(), v0.qtarm());
}
if (registryDataAcceptV0AntiLambda(v0, pos, neg, collision)) {
AntiV0Numbers = AntiV0Numbers + 1;
registryData.fill(HIST("AntiLambdaPtMass"), v0.pt(), v0.mAntiLambda());
registryData.fill(HIST("AntiLambdaQA/hArmenterosPreAnalyserCuts"), v0.alpha(), v0.qtarm());
}
}
registryData.fill(HIST("nV0sPerEvent"), V0Numbers);
Expand Down
Loading