diff --git a/PWGLF/Tasks/Strangeness/kinkAnalysis.cxx b/PWGLF/Tasks/Strangeness/kinkAnalysis.cxx index a019f6ff56b..cdc46e561d5 100644 --- a/PWGLF/Tasks/Strangeness/kinkAnalysis.cxx +++ b/PWGLF/Tasks/Strangeness/kinkAnalysis.cxx @@ -342,6 +342,8 @@ struct kinkAnalysis { histos.add("hNSigmaTrVsPt", "hNSigmaTrVsPt", kTH2F, {axisPtHyp, {100, -5, 5, "nSigmaTPC"}}); histos.add("hpRes", "hpRes", kTH2F, {axisPtHyp, {100, -0.5, 0.5, "p_{T} Res"}}); + histos.add("hDCAdaughterMC", "hDCAdaughterMC", kTH1F, {axisDCAdaugh}); + histos.add("tpcClusterTriton", "tpcClusterTriton", kTH1F, {{100, 0, 300, "TPC clusters"}}); lut = o2::base::MatLayerCylSet::rectifyPtrFromFile(ccdb->get(cfgLUTpath)); ft2.setMaxChi2(5); @@ -803,6 +805,18 @@ struct kinkAnalysis { mass = sqrt((neutronE + pionE) * (neutronE + pionE) - sigmaPabsDC * sigmaPabsDC); + if (isDaughter) + histos.fill(HIST("hDCAdaughterMC"), kinkdcaXY); + + if (particleName == Hypertriton) { + if (trackDgh.tpcNSigmaTr() > 4 || trackDgh.tpcNSigmaTr() < -2) + continue; + + histos.fill(HIST("tpcClusterTriton"), trackDgh.tpcNClsFound()); + if (trackDgh.tpcNClsFound() < 90) + continue; + } + if (particleName == Hypertriton) { histos.fill(HIST("hHypMass"), mass); histos.fill(HIST("hHypMassPt"), mass, sigmaPt);