Skip to content
14 changes: 14 additions & 0 deletions PWGLF/Tasks/Strangeness/kinkAnalysis.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -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<o2::base::MatLayerCylSet>(cfgLUTpath));
ft2.setMaxChi2(5);
Expand Down Expand Up @@ -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);
Expand Down