Skip to content

Commit ef06573

Browse files
jinhyunnijinhyunni
authored
Merge pull request #9 from jinhyunni/xic0-update
[xic0omegac0] Re-update selections
2 parents 4a7dc21 + f8dc7c4 commit ef06573

File tree

2 files changed

+22
-23
lines changed

2 files changed

+22
-23
lines changed

‎PWGHF/Core/SelectorCuts.h‎

Lines changed: 14 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -765,7 +765,7 @@ static const std::vector<std::string> labelsCutVar = {"pT Ka from Omegac"};
765765
namespace hf_cuts_xic_to_xi_pi
766766
{
767767
static constexpr int NBinsPt = 12;
768-
static constexpr int NCutVars = 28;
768+
static constexpr int NCutVars = 27;
769769
// default values for the pT bin edges (can be used to configure histogram axis)
770770
// offset by 1 from the bin numbers in cuts array
771771
constexpr double BinsPt[NBinsPt + 1] = {
@@ -786,18 +786,18 @@ constexpr double BinsPt[NBinsPt + 1] = {
786786
const auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
787787

788788
// default values for the cuts
789-
constexpr double Cuts[NBinsPt][NCutVars] = {{0.2, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 0 < pt < 1 */
790-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 1 < pt < 2 */
791-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 2 < pt < 3 */
792-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 3 < pt < 4 */
793-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 4 < pt < 5 */
794-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 5 < pt < 6 */
795-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 6 < pt < 8 */
796-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 8 < pt < 10 */
797-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 10 < pt < 12 */
798-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 12 < pt < 16 */
799-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 16 < pt < 24 */
800-
{0.5, 0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}}; /* 24 < pt < 50 -> Implemented just for test*/
789+
constexpr double Cuts[NBinsPt][NCutVars] = {{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 0 < pt < 1 */
790+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 1 < pt < 2 */
791+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 2 < pt < 3 */
792+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 3 < pt < 4 */
793+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 4 < pt < 5 */
794+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 5 < pt < 6 */
795+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 6 < pt < 8 */
796+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 8 < pt < 10 */
797+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 10 < pt < 12 */
798+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 12 < pt < 16 */
799+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}, /* 16 < pt < 24 */
800+
{0.97, 0.97, 0.99, 0.99, 2.0, 1.0, 1.0, 0.04, 0.06, 0.06, 0.05, 0.3, 70, 60, 100, 120, 250, 250, 0.4, 100, 300, 0., 0., 1.5, 0., 0., 0.4}}; /* 24 < pt < 50 -> Implemented just for test*/
801801

802802
// row labels
803803
static const std::vector<std::string> labelsPt = {
@@ -815,8 +815,7 @@ static const std::vector<std::string> labelsPt = {
815815
"pT bin 11"};
816816

817817
// column labels
818-
static const std::vector<std::string> labelsCutVar = {"ptPiFromCharmBaryon",
819-
"cosPACasc", "cosPAV0",
818+
static const std::vector<std::string> labelsCutVar = {"cosPACasc", "cosPAV0",
820819
"cosPaCascToXic", "cosPaV0ToCasc",
821820
"dcaCharmBaryonDau", "dcaCascDau", "dcaV0Dau",
822821
"dcaXYToPvCascDau", "dcaXYToPvV0Dau0", "dcaXYToPvV0Dau1",

‎PWGHF/TableProducer/candidateSelectorToXiPiQa.cxx‎

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -359,7 +359,7 @@ struct HfCandidateSelectorToXiPiQa {
359359

360360
// pT: Bachelor
361361
double ptPiFromCasc = RecoDecay::sqrtSumOfSquares(candidate.pxBachFromCasc(), candidate.pyBachFromCasc());
362-
if (ptPiFromCasc < ptPiFromCascMin) {
362+
if (std::abs(ptPiFromCasc) < ptPiFromCascMin) {
363363
return false;
364364
}
365365
registry.fill(HIST("hSelStatusLf"), 11.0);
@@ -446,7 +446,7 @@ struct HfCandidateSelectorToXiPiQa {
446446
registry.fill(HIST("hSelStatusHf"), 0.0);
447447

448448
// eta selection on charm bayron bachelor
449-
if (candidate.etaBachFromCharmBaryon() > etaTrackCharmBachMax) {
449+
if (std::abs(candidate.etaBachFromCharmBaryon()) > etaTrackCharmBachMax) {
450450
return false;
451451
}
452452
registry.fill(HIST("hSelStatusHf"), 1.0);
@@ -458,7 +458,7 @@ struct HfCandidateSelectorToXiPiQa {
458458

459459
// pT: Charm Bachelor
460460
double ptPiFromCharmBaryon = RecoDecay::pt(candidate.pxBachFromCharmBaryon(), candidate.pyBachFromCharmBaryon());
461-
if (ptPiFromCharmBaryon < cuts->get("ptPiFromCharmBaryon")) {
461+
if (std::abs(ptPiFromCharmBaryon) < ptPiFromCharmBaryonMin) {
462462
return false;
463463
}
464464
registry.fill(HIST("hSelStatusHf"), 3.0);
@@ -504,11 +504,11 @@ struct HfCandidateSelectorToXiPiQa {
504504
registry.fill(HIST("hSelStatusHf"), 10.0);
505505
} else {
506506
// Impact parameter(DCA?)
507-
if (std::abs(candidate.impactParBachFromCharmBaryonXY()) < impactParXYPiFromCharmBaryonMin || std::abs(candidate.impactParBachFromCharmBaryonXY()) > impactParXYPiFromCharmBaryonMax) {
507+
if ((std::abs(candidate.impactParBachFromCharmBaryonXY()) < impactParXYPiFromCharmBaryonMin) || (std::abs(candidate.impactParBachFromCharmBaryonXY()) > impactParXYPiFromCharmBaryonMax)) {
508508
return false;
509509
}
510510
registry.fill(HIST("hSelStatusHf"), 4.0);
511-
if (std::abs(candidate.impactParBachFromCharmBaryonZ()) < impactParZPiFromCharmBaryonMin || std::abs(candidate.impactParBachFromCharmBaryonZ()) > impactParZPiFromCharmBaryonMax) {
511+
if ((std::abs(candidate.impactParBachFromCharmBaryonZ()) < impactParZPiFromCharmBaryonMin) || (std::abs(candidate.impactParBachFromCharmBaryonZ()) > impactParZPiFromCharmBaryonMax)) {
512512
return false;
513513
}
514514
registry.fill(HIST("hSelStatusHf"), 5.0);
@@ -683,13 +683,13 @@ struct HfCandidateSelectorToXiPiQa {
683683
double invMassCascade = candidate.invMassCascade();
684684
double invMassCharmBaryon = candidate.invMassCharmBaryon();
685685

686-
if ((invMassLambda - o2::constants::physics::MassLambda0) < v0MassWindow) {
686+
if (std::abs(invMassLambda - o2::constants::physics::MassLambda0) < v0MassWindow) {
687687
statusInvMassLambda = true;
688688
}
689-
if ((invMassCascade - o2::constants::physics::MassXiMinus) < cascMassWindow) {
689+
if (std::abs(invMassCascade - o2::constants::physics::MassXiMinus) < cascMassWindow) {
690690
statusInvMassCascade = true;
691691
}
692-
if ((invMassCharmBaryon > invMassCharmBaryonMin) && (invMassCharmBaryon < invMassCharmBaryonMax)) {
692+
if ((invMassCharmBaryon >= invMassCharmBaryonMin) && (invMassCharmBaryon <= invMassCharmBaryonMax)) {
693693
statusInvMassCharmBaryon = true;
694694
}
695695

0 commit comments

Comments
 (0)