Skip to content
Merged
Changes from all commits
Commits
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
96 changes: 48 additions & 48 deletions PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -651,7 +651,7 @@ struct correlateStrangeness {
if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) {
return;
}
if (!collision.isInelGt0()&&selectINELgtZERO){
if (!collision.isInelGt0() && selectINELgtZERO) {
return;
}
// ________________________________________________
Expand Down Expand Up @@ -713,7 +713,7 @@ struct correlateStrangeness {
if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) {
return;
}
if (!collision.isInelGt0()&&selectINELgtZERO){
if (!collision.isInelGt0() && selectINELgtZERO) {
return;
}
// ________________________________________________
Expand Down Expand Up @@ -771,7 +771,7 @@ struct correlateStrangeness {
if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) {
return;
}
if (!collision.isInelGt0()&&selectINELgtZERO){
if (!collision.isInelGt0() && selectINELgtZERO) {
return;
}
// ________________________________________________
Expand Down Expand Up @@ -813,25 +813,25 @@ struct correlateStrangeness {
continue;
if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0])
continue;
if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){
if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) {
continue;

if (!doprocessMixedEventHCascades) {
if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
if (!doprocessMixedEventHCascades) {
if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
}
histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ());
histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ());
histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()}));
}
histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ());
histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ());
histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()}));
// ________________________________________________
// Do slicing
auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex());
auto slicedAssocV0s = associatedV0s.sliceBy(collisionSliceV0s, collision2.globalIndex());
// ________________________________________________
// Do hadron - V0 correlations
fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M());
}
// ________________________________________________
// Do slicing
auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex());
auto slicedAssocV0s = associatedV0s.sliceBy(collisionSliceV0s, collision2.globalIndex());
// ________________________________________________
// Do hadron - V0 correlations
fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M());
}
}
void processMixedEventHCascades(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks,
Expand All @@ -850,24 +850,24 @@ struct correlateStrangeness {
continue;
if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0])
continue;
if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){
if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) {
continue;

if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
}
if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
}

histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ());
histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ());
histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()}));
// ________________________________________________
// Do slicing
auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex());
auto slicedAssocCascades = associatedCascades.sliceBy(collisionSliceCascades, collision2.globalIndex());
// ________________________________________________
// Do hadron - cascade correlations
fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M());
}
histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ());
histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ());
histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()}));
// ________________________________________________
// Do slicing
auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex());
auto slicedAssocCascades = associatedCascades.sliceBy(collisionSliceCascades, collision2.globalIndex());
// ________________________________________________
// Do hadron - cascade correlations
fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M());
}
}
void processMixedEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks,
Expand All @@ -884,24 +884,24 @@ struct correlateStrangeness {
continue;
if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0])
continue;
if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){
if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) {
continue;

if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
}
if (collision1.globalIndex() == collision2.globalIndex()) {
histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting
}

histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ());
histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ());
histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()}));
// ________________________________________________
// Do slicing
auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex());
auto slicedAssocPions = assocPions.sliceBy(collisionSlicePions, collision2.globalIndex());
// ________________________________________________
// Do hadron - cascade correlations
fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M());
}
histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ());
histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ());
histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()}));
// ________________________________________________
// Do slicing
auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex());
auto slicedAssocPions = assocPions.sliceBy(collisionSlicePions, collision2.globalIndex());
// ________________________________________________
// Do hadron - cascade correlations
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)
{
Expand Down