Skip to content
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

[PWGHF,Trigger] Build V0s and cascades in HF triggers #9175

Merged
merged 5 commits into from
Jan 6, 2025
+499 −148
Merged

[PWGHF,Trigger] Build V0s and cascades in HF triggers #9175

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
8b19ab7
Remove dependence on lambdakzerobuilder
fgrosa Dec 31, 2024
f0473df
Remove dependence on cascadebuilder
fgrosa Jan 1, 2025
20f9416
Fix typo
fgrosa Jan 1, 2025
fb04b0f
Please consider the following formatting changes
alibuild Jan 1, 2025
fb44af8
Merge pull request #110 from alibuild/alibot-cleanup-9175
fgrosa Jan 1, 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
126 changes: 60 additions & 66 deletions EventFiltering/PWGHF/HFFilter.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

Check warning on line 1 in EventFiltering/PWGHF/HFFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/workflow-file] 

Name of a workflow file must match the name of the main struct in it (without the PWG prefix). (Class implementation files should be in "Core" directories.)

Check warning on line 1 in EventFiltering/PWGHF/HFFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pwghf/struct-member-order] 

Declare struct members in the conventional order. See the PWGHF coding guidelines.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
Expand Down Expand Up @@ -127,8 +127,8 @@

// CCDB configuration
o2::ccdb::CcdbApi ccdbApi;
Service<o2::ccdb::BasicCCDBManager> ccdb;

Check warning on line 130 in EventFiltering/PWGHF/HFFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pwghf/struct-member-order] 

HfFilter: Service< appears too early (before end of Configurable<).
Configurable<std::string> url{"ccdb-url", "http://alice-ccdb.cern.ch", "url of the ccdb repository"};

Check warning on line 131 in EventFiltering/PWGHF/HFFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/configurable] 

Use lowerCamelCase for names of configurables and use the same name for the struct member as for the JSON string. (Declare the type and names on the same line.)
int currentRun{0}; // needed to detect if the run changed and trigger update of calibrations etc.

// TPC PID calibrations
Expand All @@ -155,8 +155,9 @@
o2::vertexing::DCAFitterN<3> df3; // fitter for Charm Hadron vertex (3-prong vertex fitter)
o2::vertexing::DCAFitterN<2> dfB; // fitter for Beauty Hadron vertex (2-prong vertex fitter)
o2::vertexing::DCAFitterN<3> dfBtoDstar; // fitter for Beauty Hadron to D* vertex (3-prong vertex fitter)
o2::vertexing::DCAFitterN<2> dfStrangeness; // fitter for V0s and cascades (2-prong vertex fitter)

HistogramRegistry registry{"registry"};

Check warning on line 160 in EventFiltering/PWGHF/HFFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pwghf/struct-member-order] 

HfFilter: HistogramRegistry appears too early (before end of Preslice<).
std::shared_ptr<TH1> hProcessedEvents;

// QA histos
Expand Down Expand Up @@ -215,6 +216,8 @@
helper.setPtRangeSoftPiSigmaC(ptCuts->get(0u, 4u), ptCuts->get(1u, 4u));
helper.setPtDeltaMassRangeSigmaC(cutsPtDeltaMassCharmReso->get(0u, 6u), cutsPtDeltaMassCharmReso->get(1u, 6u), cutsPtDeltaMassCharmReso->get(0u, 7u), cutsPtDeltaMassCharmReso->get(1u, 7u), cutsPtDeltaMassCharmReso->get(0u, 8u), cutsPtDeltaMassCharmReso->get(1u, 8u), cutsPtDeltaMassCharmReso->get(0u, 9u), cutsPtDeltaMassCharmReso->get(1u, 9u), cutsPtDeltaMassCharmReso->get(2u, 6u), cutsPtDeltaMassCharmReso->get(2u, 7u), cutsPtDeltaMassCharmReso->get(2u, 8u), cutsPtDeltaMassCharmReso->get(2u, 9u));
helper.setPtRangeSoftKaonXicResoToSigmaC(ptCuts->get(0u, 5u), ptCuts->get(1u, 5u));
helper.setVtxConfiguration(dfStrangeness, true); // (DCAFitterN, useAbsDCA)
dfStrangeness.setMatCorrType(matCorr);
if (activateSecVtxForB) {
helper.setVtxConfiguration(df2, false); // (DCAFitterN, useAbsDCA)
helper.setVtxConfiguration(df3, false);
Expand Down Expand Up @@ -327,26 +330,28 @@

using BigTracksMCPID = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::pidTPCFullPr, aod::pidTOFFullPr, aod::pidTPCFullDe, aod::pidTOFFullDe, aod::McTrackLabels>;
using BigTracksPID = soa::Join<aod::Tracks, aod::TracksWCovDcaExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::pidTPCFullPr, aod::pidTOFFullPr, aod::pidTPCFullDe, aod::pidTOFFullDe>;
using TracksIUPID = soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::pidTPCFullPr, aod::pidTOFFullPr, aod::pidTPCFullPi, aod::pidTOFFullPi>;
using CollsWithEvSel = soa::Join<aod::Collisions, aod::EvSels>;

using Hf2ProngsWithMl = soa::Join<aod::Hf2Prongs, aod::Hf2ProngMlProbs>;
using Hf3ProngsWithMl = soa::Join<aod::Hf3Prongs, aod::Hf3ProngMlProbs>;

Preslice<aod::TrackAssoc> trackIndicesPerCollision = aod::track_association::collisionId;
Preslice<aod::V0Datas> v0sPerCollision = aod::v0data::collisionId;
Preslice<aod::V0s> v0sPerCollision = aod::v0data::collisionId;
Preslice<Hf2ProngsWithMl> hf2ProngPerCollision = aod::track_association::collisionId;
Preslice<Hf3ProngsWithMl> hf3ProngPerCollision = aod::track_association::collisionId;
Preslice<aod::CascDatas> cascPerCollision = aod::cascdata::collisionId;
Preslice<aod::Cascades> cascPerCollision = aod::cascdata::collisionId;
Preslice<aod::V0PhotonsKF> photonsPerCollision = aod::v0photonkf::collisionId;

void process(CollsWithEvSel const& collisions,
aod::BCsWithTimestamps const&,
aod::V0Datas const& v0s,
aod::CascDatas const& cascades,
aod::V0s const& v0s,
aod::Cascades const& cascades,
Hf2ProngsWithMl const& cand2Prongs,
Hf3ProngsWithMl const& cand3Prongs,
aod::TrackAssoc const& trackIndices,
BigTracksPID const&,
TracksIUPID const& tracksIU,
aod::V0PhotonsKF const& photons,
aod::V0Legs const&)
{
Expand Down Expand Up @@ -380,6 +385,15 @@
helper.setValuesBB(ccdbApi, bc, std::array{ccdbBBPion.value, ccdbBBAntiPion.value, ccdbBBKaon.value, ccdbBBAntiKaon.value, ccdbBBProton.value, ccdbBBAntiProton.value, ccdbBBProton.value, ccdbBBAntiProton.value}); // dummy for deuteron
}

auto bz = o2::base::Propagator::Instance()->getNominalBz();
dfStrangeness.setBz(bz);
if (activateSecVtxForB) {
df2.setBz(bz);
df3.setBz(bz);
dfB.setBz(bz);
dfBtoDstar.setBz(bz);
}

currentRun = bc.runNumber();
}

Expand Down Expand Up @@ -702,30 +716,21 @@
if (!keepEvent[kV0Charm2P] && isSignalTagged && (TESTBIT(selD0, 0) || TESTBIT(selD0, 1))) {
auto v0sThisCollision = v0s.sliceBy(v0sPerCollision, thisCollId);
for (const auto& v0 : v0sThisCollision) {
auto posTrack = v0.posTrack_as<BigTracksPID>();
auto negTrack = v0.negTrack_as<BigTracksPID>();
auto selV0 = helper.isSelectedV0(v0, std::array{posTrack, negTrack}, collision, activateQA, hV0Selected, hArmPod);
V0Cand v0Cand;
if (!helper.buildV0(v0, tracksIU, collision, dfStrangeness, std::vector{cand2Prong.prong0Id(), cand2Prong.prong1Id()}, v0Cand)) {
continue;
}
auto selV0 = helper.isSelectedV0(v0Cand, activateQA, hV0Selected, hArmPod);
if (!selV0) {
continue;
}

// propagate to PV
gpu::gpustd::array<float, 2> dcaInfo;
std::array<float, 3> pVecV0 = {v0.px(), v0.py(), v0.pz()};
std::array<float, 3> pVecV0Orig = {v0.px(), v0.py(), v0.pz()};
std::array<float, 3> posVecV0 = {v0.x(), v0.y(), v0.z()};
if (!keepEvent[kV0Charm2P] && TESTBIT(selV0, kK0S)) {

auto trackParK0 = o2::track::TrackPar(posVecV0, pVecV0Orig, 0, true);
trackParK0.setPID(o2::track::PID::K0);
trackParK0.setAbsCharge(0);
o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackParK0, 2.f, matCorr, &dcaInfo);
getPxPyPz(trackParK0, pVecV0);

// we first look for a D*+
for (const auto& trackBachelorId : trackIdsThisCollision) { // start loop over tracks
auto trackBachelor = trackBachelorId.track_as<BigTracksPID>();
if (trackBachelor.globalIndex() == trackPos.globalIndex() || trackBachelor.globalIndex() == trackNeg.globalIndex()) {
if (trackBachelor.globalIndex() == trackPos.globalIndex() || trackBachelor.globalIndex() == trackNeg.globalIndex() || trackBachelor.globalIndex() == v0.posTrackId() || trackBachelor.globalIndex() == v0.negTrackId()) {
continue;
}

Expand Down Expand Up @@ -756,10 +761,10 @@
if (activateQA) {
hMassVsPtC[kNCharmParticles]->Fill(ptDStarCand, massDiffDstar);
}
auto pVecReso2Prong = RecoDecay::pVec(pVecDStarCand, pVecV0);
auto pVecReso2Prong = RecoDecay::pVec(pVecDStarCand, v0Cand.mom);
auto ptCand = RecoDecay::pt(pVecReso2Prong);
if (ptCand > cutsPtDeltaMassCharmReso->get(2u, 3u)) {
auto massDStarK0S = RecoDecay::m(std::array{pVecPos, pVecNeg, pVecBachelor, pVecV0}, std::array{massDausD0[0], massDausD0[1], massPi, massK0S});
auto massDStarK0S = RecoDecay::m(std::array{pVecPos, pVecNeg, pVecBachelor, v0Cand.mom}, std::array{massDausD0[0], massDausD0[1], massPi, massK0S});
auto massDiffDsReso = massDStarK0S - massDStarCand;
if (cutsPtDeltaMassCharmReso->get(0u, 3u) < massDiffDsReso && massDiffDsReso < cutsPtDeltaMassCharmReso->get(1u, 3u)) {
if (activateQA) {
Expand All @@ -775,24 +780,19 @@
}
}
if (!keepEvent[kV0Charm2P] && (TESTBIT(selV0, kLambda) || TESTBIT(selV0, kAntiLambda))) { // Xic(3055) and Xic(3080) --> since it occupies only a small bandwidth, we might want to keep also wrong sign pairs
auto trackParLambda = o2::track::TrackPar(posVecV0, pVecV0Orig, 0, true);
trackParLambda.setAbsCharge(0);
trackParLambda.setPID(o2::track::PID::Lambda);
o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackParLambda, 2.f, matCorr, &dcaInfo);
getPxPyPz(trackParLambda, pVecV0);
float massXicStarCand{-999.}, massXicStarBarCand{-999.};
float massDiffXicStarCand{-999.}, massDiffXicStarBarCand{-999.};
bool isRightSignXicStar{false}, isRightSignXicStarBar{false};
auto pVecReso2Prong = RecoDecay::pVec(pVec2Prong, pVecV0);
auto pVecReso2Prong = RecoDecay::pVec(pVec2Prong, v0Cand.mom);
auto ptCand = RecoDecay::pt(pVecReso2Prong);
if (ptCand > cutsPtDeltaMassCharmReso->get(2u, 5u)) {
if (TESTBIT(selD0, 0)) {
massXicStarCand = RecoDecay::m(std::array{pVecPos, pVecNeg, pVecV0}, std::array{massPi, massKa, massLambda});
massXicStarCand = RecoDecay::m(std::array{pVecPos, pVecNeg, v0Cand.mom}, std::array{massPi, massKa, massLambda});
massDiffXicStarCand = massXicStarCand - massD0Cand;
isRightSignXicStar = TESTBIT(selV0, kLambda); // right sign if Lambda
}
if (TESTBIT(selD0, 1)) {
massXicStarBarCand = RecoDecay::m(std::array{pVecPos, pVecNeg, pVecV0}, std::array{massKa, massPi, massLambda});
massXicStarBarCand = RecoDecay::m(std::array{pVecPos, pVecNeg, v0Cand.mom}, std::array{massKa, massPi, massLambda});
massDiffXicStarBarCand = massXicStarBarCand - massD0BarCand;
isRightSignXicStarBar = TESTBIT(selV0, kAntiLambda); // right sign if AntiLambda
}
Expand Down Expand Up @@ -1249,28 +1249,21 @@

if ((!keepEvent[kV0Charm3P] && isGoodDPlus) || (!keepEvent[kSigmaC0K0] && (isGoodLcToPKPi || isGoodLcToPiKP))) {
for (const auto& v0 : v0sThisCollision) {
auto posTrack = v0.posTrack_as<BigTracksPID>();
auto negTrack = v0.negTrack_as<BigTracksPID>();
auto selV0 = helper.isSelectedV0(v0, std::array{posTrack, negTrack}, collision, activateQA, hV0Selected, hArmPod);
V0Cand v0Cand;
if (!helper.buildV0(v0, tracksIU, collision, dfStrangeness, std::vector{cand3Prong.prong0Id(), cand3Prong.prong1Id(), cand3Prong.prong2Id()}, v0Cand)) {
continue;
}
auto selV0 = helper.isSelectedV0(v0Cand, activateQA, hV0Selected, hArmPod);
if (!selV0) {
continue;
}
gpu::gpustd::array<float, 2> dcaInfo;
std::array<float, 3> pVecV0Orig = {v0.px(), v0.py(), v0.pz()};
std::array<float, 3> pVecV0 = {v0.px(), v0.py(), v0.pz()};
std::array<float, 3> posVecV0 = {v0.x(), v0.y(), v0.z()};

// we pair D+ with V0
if (!keepEvent[kV0Charm3P] && isGoodDPlus) {
if (!keepEvent[kV0Charm3P] && TESTBIT(selV0, kK0S)) { // Ds2*
auto trackParK0S = o2::track::TrackPar(posVecV0, pVecV0Orig, 0, true);
trackParK0S.setAbsCharge(0);
trackParK0S.setPID(o2::track::PID::K0);
o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackParK0S, 2.f, matCorr, &dcaInfo);
getPxPyPz(trackParK0S, pVecV0);
auto massDsStarCand = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecV0}, std::array{massPi, massKa, massPi, massK0S});
auto massDsStarCand = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, v0Cand.mom}, std::array{massPi, massKa, massPi, massK0S});
auto massDiffDsStar = massDsStarCand - massDPlusCand;
auto pVecReso3Prong = RecoDecay::pVec(pVec3Prong, pVecV0);
auto pVecReso3Prong = RecoDecay::pVec(pVec3Prong, v0Cand.mom);
auto ptCand = RecoDecay::pt(pVecReso3Prong);
if (ptCand > cutsPtDeltaMassCharmReso->get(2u, 4u)) {
if (cutsPtDeltaMassCharmReso->get(0u, 4u) < massDiffDsStar && massDiffDsStar < cutsPtDeltaMassCharmReso->get(1u, 4u)) {
Expand All @@ -1282,15 +1275,10 @@
}
}
if (!keepEvent[kV0Charm3P] && (TESTBIT(selV0, kLambda) || TESTBIT(selV0, kAntiLambda))) { // Xic(3055) and Xic(3080) --> since it occupies only a small bandwidth, we might want to keep also wrong sign pairs
auto trackParLambda = o2::track::TrackPar(posVecV0, pVecV0Orig, 0, true);
trackParLambda.setAbsCharge(0);
trackParLambda.setPID(o2::track::PID::Lambda);
o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackParLambda, 2.f, matCorr, &dcaInfo);
getPxPyPz(trackParLambda, pVecV0);
auto pVecReso3Prong = RecoDecay::pVec(pVec3Prong, pVecV0);
auto pVecReso3Prong = RecoDecay::pVec(pVec3Prong, v0Cand.mom);
auto ptCand = RecoDecay::pt(pVecReso3Prong);
if (ptCand > cutsPtDeltaMassCharmReso->get(2u, 5u)) {
auto massXicStarCand = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecV0}, std::array{massPi, massKa, massPi, massLambda});
auto massXicStarCand = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, v0Cand.mom}, std::array{massPi, massKa, massPi, massLambda});
auto massDiffXicStar = massXicStarCand - massDPlusCand;
bool isRightSign = ((TESTBIT(selV0, kLambda) && sign3Prong > 0) || (TESTBIT(selV0, kAntiLambda) && sign3Prong < 0));
if (cutsPtDeltaMassCharmReso->get(0u, 5u) < massDiffXicStar && massDiffXicStar < cutsPtDeltaMassCharmReso->get(1u, 5u)) {
Expand Down Expand Up @@ -1319,8 +1307,8 @@
auto globalIndexSoftPi = trackSoftPi.globalIndex();

// exclude tracks already used to build the 3-prong candidate
if (globalIndexSoftPi == trackFirst.globalIndex() || globalIndexSoftPi == trackSecond.globalIndex() || globalIndexSoftPi == trackThird.globalIndex()) {
// do not consider as candidate soft pion a track already used to build the current 3-prong candidate
if (globalIndexSoftPi == trackFirst.globalIndex() || globalIndexSoftPi == trackSecond.globalIndex() || globalIndexSoftPi == trackThird.globalIndex() || globalIndexSoftPi == v0.posTrackId() || globalIndexSoftPi == v0.negTrackId()) {
// do not consider as candidate soft pion a track already used to build the current 3-prong candidate / V0 candidate
continue;
}

Expand Down Expand Up @@ -1353,17 +1341,15 @@
/// and keep it only if it is in the correct mass range

float massSigmaCPKPi{-999.}, massSigmaCPiKP{-999.}, deltaMassXicResoPKPi{-999.}, deltaMassXicResoPiKP{-999.};
std::array<float, 3> pVecPiPosK0s = posTrack.pVector();
std::array<float, 3> pVecPiNegK0s = negTrack.pVector();
float ptSigmaCKaon = RecoDecay::pt(pVecSigmaC, pVecPiPosK0s, pVecPiNegK0s);
float ptSigmaCKaon = RecoDecay::pt(pVecSigmaC, v0Cand.mom);
if (ptSigmaCKaon > cutsPtDeltaMassCharmReso->get(2u, 10u)) {
if (TESTBIT(whichSigmaC, 0)) {
massSigmaCPKPi = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecSoftPi}, std::array{massProton, massKa, massPi, massPi});
deltaMassXicResoPKPi = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecSoftPi, pVecPiPosK0s, pVecPiNegK0s}, std::array{massProton, massKa, massPi, massPi, massPi, massPi}) - massSigmaCPKPi;
deltaMassXicResoPKPi = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecSoftPi, v0Cand.mom}, std::array{massProton, massKa, massPi, massPi, massK0S}) - massSigmaCPKPi;
}
if (TESTBIT(whichSigmaC, 1)) {
massSigmaCPiKP = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecSoftPi}, std::array{massPi, massKa, massProton, massPi});
deltaMassXicResoPiKP = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecSoftPi, pVecPiPosK0s, pVecPiNegK0s}, std::array{massPi, massKa, massProton, massPi, massPi, massPi}) - massSigmaCPiKP;
deltaMassXicResoPiKP = RecoDecay::m(std::array{pVecFirst, pVecSecond, pVecThird, pVecSoftPi, v0Cand.mom}, std::array{massPi, massKa, massProton, massPi, massK0S}) - massSigmaCPiKP;
}

bool isPKPiOk = (cutsPtDeltaMassCharmReso->get(0u, 10u) < deltaMassXicResoPKPi && deltaMassXicResoPKPi < cutsPtDeltaMassCharmReso->get(1u, 10u));
Expand Down Expand Up @@ -1404,35 +1390,43 @@
if (!keepEvent[kCharmBarToXiBach]) {
auto cascThisColl = cascades.sliceBy(cascPerCollision, thisCollId);
for (const auto& casc : cascThisColl) {
auto bachelorCasc = casc.bachelor_as<BigTracksPID>();
auto v0DauPos = casc.posTrack_as<BigTracksPID>();
auto v0DauNeg = casc.negTrack_as<BigTracksPID>();

if (!helper.isSelectedCascade(casc, std::array{bachelorCasc, v0DauPos, v0DauNeg}, collision)) {
CascCand cascCand;
if (!helper.buildCascade(casc, v0s, tracksIU, collision, dfStrangeness, {}, cascCand)) {
continue;
}

if (!helper.isSelectedCascade(cascCand)) {
continue;
}

if (activateQA) {
hMassXi->Fill(casc.mXi());
hMassXi->Fill(cascCand.mXi);
}

auto bachelorCascId = casc.bachelorId();
auto v0 = v0s.rawIteratorAt(casc.v0Id());
auto v0DauPosId = v0.posTrackId();
auto v0DauNegId = v0.negTrackId();

auto trackIdsThisCollision = trackIndices.sliceBy(trackIndicesPerCollision, thisCollId);
for (const auto& trackId : trackIdsThisCollision) { // start loop over tracks
auto track = trackId.track_as<BigTracksPID>();

// ask for opposite sign daughters (omegac daughters)
if (track.sign() * bachelorCasc.sign() >= 0) {
if (track.sign() * cascCand.sign > 0) {
continue;
}

// check if track is one of the Xi daughters
if (track.globalIndex() == bachelorCasc.globalIndex() || track.globalIndex() == v0DauPos.globalIndex() || track.globalIndex() == v0DauNeg.globalIndex()) {
if (track.globalIndex() == bachelorCascId || track.globalIndex() == v0DauPosId || track.globalIndex() == v0DauNegId) {
continue;
}

// propagate to PV
gpu::gpustd::array<float, 2> dcaInfo;
std::array<float, 3> pVecCascade = {casc.px(), casc.py(), casc.pz()};
auto trackParCasc = o2::track::TrackPar(std::array{casc.x(), casc.y(), casc.z()}, pVecCascade, bachelorCasc.sign(), true);
std::array<float, 3> pVecCascade = cascCand.mom;
auto trackParCasc = o2::track::TrackPar(cascCand.vtx, cascCand.mom, cascCand.sign, true);
trackParCasc.setPID(o2::track::PID::XiMinus);
o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackParCasc, 2.f, matCorr, &dcaInfo);
getPxPyPz(trackParCasc, pVecCascade);
Expand Down
Loading
Loading