Skip to content

Commit f8dc7c4

Browse files
jinhyunnijinhyunni
committed
[xic0omegac0] Updated some selections
1 parent 2ea4baa commit f8dc7c4

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
@@ -756,7 +756,7 @@ static const std::vector<std::string> labelsCutVar = {"pT Ka from Omegac"};
756756
namespace hf_cuts_xic_to_xi_pi
757757
{
758758
static constexpr int NBinsPt = 12;
759-
static constexpr int NCutVars = 28;
759+
static constexpr int NCutVars = 27;
760760
// default values for the pT bin edges (can be used to configure histogram axis)
761761
// offset by 1 from the bin numbers in cuts array
762762
constexpr double BinsPt[NBinsPt + 1] = {
@@ -777,18 +777,18 @@ constexpr double BinsPt[NBinsPt + 1] = {
777777
const auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
778778

779779
// default values for the cuts
780-
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 */
781-
{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 */
782-
{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 */
783-
{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 */
784-
{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 */
785-
{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 */
786-
{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 */
787-
{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 */
788-
{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 */
789-
{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 */
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}, /* 16 < pt < 24 */
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}}; /* 24 < pt < 50 -> Implemented just for test*/
780+
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 */
781+
{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 */
782+
{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 */
783+
{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 */
784+
{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 */
785+
{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 */
786+
{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 */
787+
{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 */
788+
{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 */
789+
{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 */
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}, /* 16 < pt < 24 */
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}}; /* 24 < pt < 50 -> Implemented just for test*/
792792

793793
// row labels
794794
static const std::vector<std::string> labelsPt = {
@@ -806,8 +806,7 @@ static const std::vector<std::string> labelsPt = {
806806
"pT bin 11"};
807807

808808
// column labels
809-
static const std::vector<std::string> labelsCutVar = {"ptPiFromCharmBaryon",
810-
"cosPACasc", "cosPAV0",
809+
static const std::vector<std::string> labelsCutVar = {"cosPACasc", "cosPAV0",
811810
"cosPaCascToXic", "cosPaV0ToCasc",
812811
"dcaCharmBaryonDau", "dcaCascDau", "dcaV0Dau",
813812
"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)