Skip to content

PWGLF: add INEL>0 selection to hStr corr analysis #6761

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
njacazio merged 6 commits into from
Jul 10, 2024
Merged

PWGLF: add INEL>0 selection to hStr corr analysis #6761

Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
3ad0004
PWGLF: add INEL>0 selection to hStr corr analysis
ddobrigk Jul 8, 2024
d2461c6
Please consider the following formatting changes (#309)
alibuild Jul 8, 2024
2d13491
Bugfixes
ddobrigk Jul 9, 2024
1349430
Merge branch 'master' into inelgtzero
ddobrigk Jul 9, 2024
3e9ef56
Merge branch 'master' into inelgtzero
ddobrigk Jul 9, 2024
90ce4e1
Please consider the following formatting changes (#310)
alibuild Jul 9, 2024
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
35 changes: 28 additions & 7 deletions PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -59,6 +59,7 @@ struct correlateStrangeness {
Configurable<bool> doCorrelationOmegaMinus{"doCorrelationOmegaMinus", false, "do OmegaMinus correlation"};
Configurable<bool> doCorrelationOmegaPlus{"doCorrelationOmegaPlus", false, "do OmegaPlus correlation"};
Configurable<bool> doCorrelationPion{"doCorrelationPion", false, "do Pion correlation"};
Configurable<bool> selectINELgtZERO{"selectINELgtZERO", true, "select INEL>0 events"};
Configurable<float> zVertexCut{"zVertexCut", 10, "Cut on PV position"};
Configurable<bool> skipUnderOverflowInTHn{"skipUnderOverflowInTHn", false, "skip under/overflow in THns"};
Configurable<int> mixingParameter{"mixingParameter", 10, "how many events are mixed"};
Expand Down Expand Up @@ -663,7 +664,7 @@ struct correlateStrangeness {
}
}

void processSameEventHV0s(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>::iterator const& collision,
void processSameEventHV0s(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>::iterator const& collision,
aod::AssocV0s const& associatedV0s, aod::TriggerTracks const& triggerTracks,
V0DatasWithoutTrackX const&, aod::V0sLinked const&, TracksComplete const&, aod::BCsWithTimestamps const&)
{
Expand All @@ -678,6 +679,9 @@ struct correlateStrangeness {
if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) {
return;
}
if (!collision.isInelGt0() && selectINELgtZERO) {
return;
}
// ________________________________________________
if (!doprocessSameEventHCascades) {
histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()}));
Expand Down Expand Up @@ -727,7 +731,7 @@ struct correlateStrangeness {
fillCorrelationsV0(triggerTracks, associatedV0s, false, collision.posZ(), collision.centFT0M());
}

void processSameEventHCascades(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>::iterator const& collision,
void processSameEventHCascades(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>::iterator const& collision,
aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks,
V0DatasWithoutTrackX const&, aod::V0sLinked const&, aod::CascDatas const&, TracksComplete const&, aod::BCsWithTimestamps const&)
{
Expand All @@ -742,6 +746,9 @@ struct correlateStrangeness {
if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) {
return;
}
if (!collision.isInelGt0() && selectINELgtZERO) {
return;
}
// ________________________________________________
histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()}));
histos.fill(HIST("EventQA/hMult"), collision.centFT0M());
Expand Down Expand Up @@ -787,7 +794,7 @@ struct correlateStrangeness {
// Do hadron - cascade correlations
fillCorrelationsCascade(triggerTracks, associatedCascades, false, collision.posZ(), collision.centFT0M());
}
void processSameEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>::iterator const& collision,
void processSameEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>::iterator const& collision,
aod::AssocPions const& associatedPions, aod::TriggerTracks const& triggerTracks,
TracksComplete const&, aod::BCsWithTimestamps const&)
{
Expand All @@ -802,6 +809,9 @@ struct correlateStrangeness {
if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) {
return;
}
if (!collision.isInelGt0() && selectINELgtZERO) {
return;
}
// ________________________________________________
if (!doprocessSameEventHCascades && !doprocessSameEventHV0s) {
histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()}));
Expand All @@ -824,7 +834,7 @@ struct correlateStrangeness {
// Do hadron - Pion correlations
fillCorrelationsPion(triggerTracks, associatedPions, false, collision.posZ(), collision.centFT0M());
}
void processMixedEventHV0s(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
void processMixedEventHV0s(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults> const& collisions,
aod::AssocV0s const& associatedV0s, aod::TriggerTracks const& triggerTracks,
V0DatasWithoutTrackX const&, aod::V0sLinked const&, TracksComplete const&, aod::BCsWithTimestamps const&)
{
Expand All @@ -843,6 +853,8 @@ struct correlateStrangeness {
continue;
if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0])
continue;
if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO)
continue;

if (!doprocessMixedEventHCascades) {
if (collision1.globalIndex() == collision2.globalIndex()) {
Expand All @@ -861,7 +873,7 @@ struct correlateStrangeness {
fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M());
}
}
void processMixedEventHCascades(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
void processMixedEventHCascades(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults> const& collisions,
aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks,
V0DatasWithoutTrackX const&, aod::V0sLinked const&, aod::CascDatas const&, TracksComplete const&, aod::BCsWithTimestamps const&)
{
Expand All @@ -880,6 +892,8 @@ struct correlateStrangeness {
continue;
if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0])
continue;
if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO)
continue;

if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
Expand All @@ -897,7 +911,7 @@ struct correlateStrangeness {
fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M());
}
}
void processMixedEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
void processMixedEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults> const& collisions,
aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks,
TracksComplete const&)
{
Expand All @@ -912,6 +926,8 @@ struct correlateStrangeness {
continue;
if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0])
continue;
if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO)
continue;

if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
Expand All @@ -929,7 +945,8 @@ struct correlateStrangeness {
fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M());
}
}
void processMCGenerated(aod::McCollision const&, soa::SmallGroups<soa::Join<aod::McCollisionLabels, aod::Collisions, aod::EvSels, aod::CentFT0Ms>> const& collisions, aod::McParticles const& mcParticles)

void processMCGenerated(aod::McCollision const& mccollision, soa::SmallGroups<soa::Join<aod::McCollisionLabels, aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>> const& collisions, aod::McParticles const& mcParticles)
{

for (auto const& mcParticle : mcParticles) {
Expand Down Expand Up @@ -960,12 +977,14 @@ struct correlateStrangeness {
int bestCollisionFT0Mpercentile = -1;
float bestCollisionVtxZ = 0.0f;
bool bestCollisionSel8 = false;
bool bestCollisionINELgtZERO = false;
for (auto& collision : collisions) {
if (biggestNContribs < collision.numContrib()) {
biggestNContribs = collision.numContrib();
bestCollisionFT0Mpercentile = collision.centFT0M();
bestCollisionSel8 = collision.sel8();
bestCollisionVtxZ = collision.posZ();
bestCollisionINELgtZERO = collision.isInelGt0();
}
}

Expand Down Expand Up @@ -1001,6 +1020,8 @@ struct correlateStrangeness {
return;
if (std::abs(bestCollisionVtxZ) > 10.0f)
return;
if (!bestCollisionINELgtZERO)
return;

for (auto const& mcParticle : mcParticles) {
if (!mcParticle.isPhysicalPrimary()) {
Expand Down