Skip to content

Commit 4c4da05

Browse files
committed
[xic0omegac0] Fix inconsistency in selector
1 parent df3f047 commit 4c4da05

1 file changed

Lines changed: 17 additions & 10 deletions

File tree

PWGHF/TableProducer/candidateSelectorToXiPiQa.cxx

Lines changed: 17 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -251,10 +251,11 @@ struct HfCandidateSelectorToXiPiQa {
251251
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(5, "cosPaCascToXic");
252252
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(6, "kfDcaXYPiFromXic");
253253
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(7, "chi2GeoXic");
254-
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(8, "chi2TopoPiFromXicToPv");
255-
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(9, "chi2TopoCascToXic");
256-
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(10, "decayLenXYXic");
257-
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(11, "cTauXic");
254+
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(8, "chi2TopoXicToPv");
255+
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(9, "chi2TopoPiFromXicToPv");
256+
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(10, "chi2TopoCascToXic");
257+
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(11, "decayLenXYXic");
258+
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(12, "cTauXic");
258259
}
259260

260261
// invarinat mass histograms
@@ -372,7 +373,7 @@ struct HfCandidateSelectorToXiPiQa {
372373
registry.fill(HIST("hSelStatusLf"), 12.0);
373374

374375
// DCA
375-
if (candidate.kfDcaXYCascToPv() > cuts->get(inputPtBin, "kfDcaXYCascToPv")) {
376+
if (std::abs(candidate.kfDcaXYCascToPv()) > cuts->get(inputPtBin, "kfDcaXYCascToPv")) {
376377
return false;
377378
}
378379
registry.fill(HIST("hSelStatusLf"), 13.0);
@@ -481,26 +482,30 @@ struct HfCandidateSelectorToXiPiQa {
481482
return false;
482483
}
483484
registry.fill(HIST("hSelStatusHf"), 6.0);
484-
if (candidate.chi2TopoPiFromXicToPv() < 0 || candidate.chi2TopoPiFromXicToPv() > cuts->get(inputPtBin, "chi2TopoXicToPv")) {
485+
if (candidate.chi2TopoXicToPv() < 0 || candidate.chi2TopoXicToPv() > cuts->get(inputPtBin, "chi2TopoXicToPv")) {
485486
return false;
486487
}
487488
registry.fill(HIST("hSelStatusHf"), 7.0);
488-
if (candidate.chi2TopoCascToXic() < 0 || candidate.chi2TopoCascToXic() > cuts->get(inputPtBin, "chi2TopoCascToXic")) {
489+
if (candidate.chi2TopoPiFromXicToPv() < 0 || candidate.chi2TopoPiFromXicToPv() > cuts->get(inputPtBin, "chi2TopoPiFromXicToPv")) {
489490
return false;
490491
}
491492
registry.fill(HIST("hSelStatusHf"), 8.0);
493+
if (candidate.chi2TopoCascToXic() < 0 || candidate.chi2TopoCascToXic() > cuts->get(inputPtBin, "chi2TopoCascToXic")) {
494+
return false;
495+
}
496+
registry.fill(HIST("hSelStatusHf"), 9.0);
492497

493498
// Decay Length
494499
if (std::abs(candidate.decayLenXYXic()) > cuts->get(inputPtBin, "decayLenXYXic")) {
495500
return false;
496501
}
497-
registry.fill(HIST("hSelStatusHf"), 9.0);
502+
registry.fill(HIST("hSelStatusHf"), 10.0);
498503

499504
// ctau
500505
if (std::abs(candidate.cTauXic()) > cuts->get(inputPtBin, "cTauXic")) {
501506
return false;
502507
}
503-
registry.fill(HIST("hSelStatusHf"), 10.0);
508+
registry.fill(HIST("hSelStatusHf"), 11.0);
504509
} else {
505510
// Impact parameter(DCA?)
506511
if ((std::abs(candidate.impactParBachFromCharmBaryonXY()) < impactParXYPiFromCharmBaryonMin) || (std::abs(candidate.impactParBachFromCharmBaryonXY()) > impactParXYPiFromCharmBaryonMax)) {
@@ -715,7 +720,9 @@ struct HfCandidateSelectorToXiPiQa {
715720
trackPiFromCharm.tpcNSigmaPi(), trackPiFromCasc.tpcNSigmaPi(), trackPiFromLam.tpcNSigmaPi(), trackPrFromLam.tpcNSigmaPr(),
716721
trackPiFromCharm.tofNSigmaPi(), trackPiFromCasc.tofNSigmaPi(), trackPiFromLam.tofNSigmaPi(), trackPrFromLam.tofNSigmaPr());
717722
} else {
718-
if (!statusPidCharmBaryon || !statusInvMassCharmBaryon) {
723+
// Convert resultSelectrions to false if one of pid status and inv status is false
724+
// This is only needed for KF(currently) since output selection table doesn't contain those information
725+
if (!statusPidCharmBaryon || !statusInvMassLambda || !statusInvMassCascade || !statusInvMassCharmBaryon) {
719726
resultSelections = false;
720727
}
721728
hfSelToXiPiKf(resultSelections,

0 commit comments

Comments
 (0)