Skip to content

PWGCF: Change combinations policy for non identical V0s in V0 task #7325

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
saganatt merged 2 commits into from
Aug 16, 2024
Merged

PWGCF: Change combinations policy for non identical V0s in V0 task #7325

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
6be0abe
Change combinations policy for non identical pairs, change congigurab…
alicja-pp Aug 16, 2024
7d26872
Fix one configurable description
alicja-pp Aug 16, 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
2 changes: 1 addition & 1 deletion PWGCF/FemtoUniverse/Core/FemtoUniverseEventHisto.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@ class FemtoUniverseEventHisto
mHistogramRegistry->add("Event/zvtxhist", "; vtx_{z} (cm); Entries", kTH1F, {{300, -12.5, 12.5}});
mHistogramRegistry->add("Event/MultV0M", "; vMultV0M; Entries", kTH1F, {{16384, 0, 32768}});
mHistogramRegistry->add("Event/MultNTr", "; vMultNTr; Entries", kTH1F, {{200, 0, 200}});
mHistogramRegistry->add("Event/MultNTrVSMultV0M", "; vMultNTr; MultV0M", kTH2F, {{200, 0, 200}, {16384, 0, 32768}});
mHistogramRegistry->add("Event/MultNTrVSMultV0M", "; vMultNTr; MultV0M", kTH2F, {{4000, 0, 4000}, {32768, 0, 32768}});
}

/// Some basic QA of the event
Expand Down
5 changes: 2 additions & 3 deletions PWGCF/FemtoUniverse/Tasks/femtoUniverseEfficiencyBase.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +49,7 @@ struct femtoUniverseEfficiencyBase {

/// Particle 1
Configurable<int32_t> ConfPDGCodePartOne{"ConfPDGCodePartOne", 2212, "Particle 1 - PDG code"};
Configurable<uint8_t> ConfParticleTypePartOne{"ConfParticleTypePartOne", aod::femtouniverseparticle::ParticleType::kTrack, "Particle 1 - particle type"};
Configurable<uint8_t> ConfParticleTypePartOne{"ConfParticleTypePartOne", aod::femtouniverseparticle::ParticleType::kTrack, "Particle 1 - particle type: 0 - track, 2 - V0, 6 - phi"};
Configurable<bool> ConfNoPDGPartOne{"ConfNoPDGPartOne", false, "0: selecting part one by PDG, 1: no PID selection"};
Configurable<float> ConfPtLowPart1{"ConfPtLowPart1", 0.2, "Lower limit for Pt for the first particle"};
Configurable<float> ConfPtHighPart1{"ConfPtHighPart1", 2.5, "Higher limit for Pt for the first particle"};
Expand All @@ -70,7 +70,7 @@ struct femtoUniverseEfficiencyBase {
/// Particle 2
Configurable<bool> ConfIsSame{"ConfIsSame", false, "Pairs of the same particle"};
Configurable<int32_t> ConfPDGCodePartTwo{"ConfPDGCodePartTwo", 333, "Particle 2 - PDG code"};
Configurable<uint8_t> ConfParticleTypePartTwo{"ConfParticleTypePartTwo", aod::femtouniverseparticle::ParticleType::kTrack, "Particle 2 - particle type"};
Configurable<uint8_t> ConfParticleTypePartTwo{"ConfParticleTypePartTwo", aod::femtouniverseparticle::ParticleType::kTrack, "Particle 2 - particle type: 0 - track, 2 - V0, 6 - phi"};
Configurable<bool> ConfNoPDGPartTwo{"ConfNoPDGPartTwo", false, "0: selecting part two by PDG, 1: no PID selection"};
Configurable<float> ConfPtLowPart2{"ConfPtLowPart2", 0.2, "Lower limit for Pt for the second particle"};
Configurable<float> ConfPtHighPart2{"ConfPtHighPart2", 2.5, "Higher limit for Pt for the second particle"};
Expand Down Expand Up @@ -333,7 +333,6 @@ struct femtoUniverseEfficiencyBase {
void doMCRecTrackPhi(PartitionType grouppartsOneMCGen, PartitionType grouppartsTwoGen)
{ // part1 is track and part2 is Phi

/// Histogramming same event
for (auto& part : grouppartsOneMCGen) {
if (part.partType() != ConfParticleTypePartOne || part.sign() != ConfChargePart1 || !IsParticleNSigma(ConfPDGCodePartOne, part.p(), trackCuts.getNsigmaTPC(part, o2::track::PID::Proton), trackCuts.getNsigmaTOF(part, o2::track::PID::Proton), trackCuts.getNsigmaTPC(part, o2::track::PID::Pion), trackCuts.getNsigmaTOF(part, o2::track::PID::Pion), trackCuts.getNsigmaTPC(part, o2::track::PID::Kaon), trackCuts.getNsigmaTOF(part, o2::track::PID::Kaon))) {
continue;
Expand Down
84 changes: 59 additions & 25 deletions PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Extended.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -395,36 +395,70 @@ struct femtoUniversePairTaskTrackV0Extended {
}
}

/// Now build the combinations
for (auto& [p1, p2] : combinations(CombinationsStrictlyUpperIndexPolicy(groupPartsTwo, groupPartsTwo))) {
// Lambda invariant mass cut for p1
if (!invMLambda(p1.mLambda(), p1.mAntiLambda()))
continue;
// Lambda invariant mass cut for p2
if (!invMLambda(p2.mLambda(), p2.mAntiLambda()))
continue;
// track cleaning
if (!pairCleanerV0.isCleanPair(p1, p2, parts)) {
continue;
}
if (ConfIsCPR.value) {
if (pairCloseRejectionV0.isClosePair(p1, p2, parts, magFieldTesla, femtoUniverseContainer::EventType::same)) {
if (ConfV0Type1 == ConfV0Type2) {
/// Now build the combinations for identical V0s
for (auto& [p1, p2] : combinations(CombinationsStrictlyUpperIndexPolicy(groupPartsTwo, groupPartsTwo))) {
// Lambda invariant mass cut for p1
if (!invMLambda(p1.mLambda(), p1.mAntiLambda()))
continue;
// Lambda invariant mass cut for p2
if (!invMLambda(p2.mLambda(), p2.mAntiLambda()))
continue;
// track cleaning
if (!pairCleanerV0.isCleanPair(p1, p2, parts)) {
continue;
}
if (ConfIsCPR.value) {
if (pairCloseRejectionV0.isClosePair(p1, p2, parts, magFieldTesla, femtoUniverseContainer::EventType::same)) {
continue;
}
}
const auto& posChild1 = parts.iteratorAt(p1.index() - 2);
const auto& negChild1 = parts.iteratorAt(p1.index() - 1);
/// Daughters that do not pass this condition are not selected
if (!IsParticleTPC(posChild1, V0ChildTable[ConfV0Type1][0]) || !IsParticleTPC(negChild1, V0ChildTable[ConfV0Type1][1]))
continue;

const auto& posChild2 = parts.iteratorAt(p2.index() - 2);
const auto& negChild2 = parts.iteratorAt(p2.index() - 1);
/// Daughters that do not pass this condition are not selected
if (!IsParticleTPC(posChild2, V0ChildTable[ConfV0Type2][0]) || !IsParticleTPC(negChild2, V0ChildTable[ConfV0Type2][1]))
continue;

sameEventCont.setPair<false>(p1, p2, multCol, ConfUse3D);
}
const auto& posChild1 = parts.iteratorAt(p1.index() - 2);
const auto& negChild1 = parts.iteratorAt(p1.index() - 1);
/// Daughters that do not pass this condition are not selected
if (!IsParticleTPC(posChild1, V0ChildTable[ConfV0Type1][0]) || !IsParticleTPC(negChild1, V0ChildTable[ConfV0Type1][1]))
continue;
} else {
/// Now build the combinations for not identical identical V0s
for (auto& [p1, p2] : combinations(CombinationsFullIndexPolicy(groupPartsTwo, groupPartsTwo))) {
// Lambda invariant mass cut for p1
if (!invMLambda(p1.mLambda(), p1.mAntiLambda()))
continue;
// Lambda invariant mass cut for p2
if (!invMLambda(p2.mLambda(), p2.mAntiLambda()))
continue;
// track cleaning
if (!pairCleanerV0.isCleanPair(p1, p2, parts)) {
continue;
}
if (ConfIsCPR.value) {
if (pairCloseRejectionV0.isClosePair(p1, p2, parts, magFieldTesla, femtoUniverseContainer::EventType::same)) {
continue;
}
}
const auto& posChild1 = parts.iteratorAt(p1.index() - 2);
const auto& negChild1 = parts.iteratorAt(p1.index() - 1);
/// Daughters that do not pass this condition are not selected
if (!IsParticleTPC(posChild1, V0ChildTable[ConfV0Type1][0]) || !IsParticleTPC(negChild1, V0ChildTable[ConfV0Type1][1]))
continue;

const auto& posChild2 = parts.iteratorAt(p2.index() - 2);
const auto& negChild2 = parts.iteratorAt(p2.index() - 1);
/// Daughters that do not pass this condition are not selected
if (!IsParticleTPC(posChild2, V0ChildTable[ConfV0Type2][0]) || !IsParticleTPC(negChild2, V0ChildTable[ConfV0Type2][1]))
continue;
const auto& posChild2 = parts.iteratorAt(p2.index() - 2);
const auto& negChild2 = parts.iteratorAt(p2.index() - 1);
/// Daughters that do not pass this condition are not selected
if (!IsParticleTPC(posChild2, V0ChildTable[ConfV0Type2][0]) || !IsParticleTPC(negChild2, V0ChildTable[ConfV0Type2][1]))
continue;

sameEventCont.setPair<false>(p1, p2, multCol, ConfUse3D);
sameEventCont.setPair<false>(p1, p2, multCol, ConfUse3D);
}
}
}

Expand Down