Skip to content

Commit ac62b6b

Browse files
FDUEnrichFDUEnrich
committed
cut pt(final)
1 parent df9a8dc commit ac62b6b

File tree

1 file changed

+45
-37
lines changed

1 file changed

+45
-37
lines changed

‎PWGHF/TableProducer/trackIndexSkimCreator.cxx‎

Lines changed: 45 additions & 37 deletions
Original file line numberDiff line numberDiff line change
@@ -3236,10 +3236,10 @@ struct HfTrackIndexSkimCreatorLfCascades {
32363236
// selections have been set to run2 lambda dedicated cuts
32373237
// selections for cascade have been set to the loosest value between xi and omega
32383238
// a tolerance has been added to be more conservative ptMinOmegaczeroToOmegaKaLfCasc ptMinXicZeroOmegacZeroToXiPiLfCasc
3239-
Configurable<float> ptMinOmegacZeroToOmegaPiLfCasc{"ptMinOmegacZeroToOmegaPiLfCasc", 2.0, "min. pT for Omegaczero in Omega + Pi decays"};
3240-
Configurable<float> ptMinOmegaczeroToOmegaKaLfCasc{"ptMinOmegaczeroToOmegaKaLfCasc", 2.0, "min. pT for Omegaczero in Omega + Ka decays"};
3241-
Configurable<float> ptMinXicZeroOmegacZeroToXiPiLfCasc{"ptMinXicZeroOmegacZeroToXiPiLfCasc", 2.0, "min. pT for XicZeroOmegacZero in Xi + Pi decays"};
3242-
Configurable<float> ptMinXicplusLfCasc{"ptMinXicplusLfCasc", 2.0, "min. pT for Xicplus in Xi + Pi + Pi decays"};
3239+
Configurable<float> ptMinOmegacZeroToOmegaPiLfCasc{"ptMinOmegacZeroToOmegaPiLfCasc", 0.f, "min. pT for Omegaczero in Omega + Pi decays"};
3240+
Configurable<float> ptMinOmegaczeroToOmegaKaLfCasc{"ptMinOmegaczeroToOmegaKaLfCasc", 0.f, "min. pT for Omegaczero in Omega + Ka decays"};
3241+
Configurable<float> ptMinXicZeroOmegacZeroToXiPiLfCasc{"ptMinXicZeroOmegacZeroToXiPiLfCasc", 0.f, "min. pT for XicZeroOmegacZero in Xi + Pi decays"};
3242+
Configurable<float> ptMinXicplusLfCasc{"ptMinXicplusLfCasc", 0.f, "min. pT for Xicplus in Xi + Pi + Pi decays"};
32433243
Configurable<float> v0TransvRadius{"v0TransvRadius", 1.0, "V0 radius in xy plane"}; // 1.2 (xi) and 1.1 (omega) in run2
32443244
Configurable<float> cascTransvRadius{"cascTransvRadius", 0.4, "Cascade radius in xy plane"}; // 0.5 cm (xi) and 0.6 (omega) in run2
32453245
Configurable<float> dcaBachToPv{"dcaBachToPv", 0.03, "DCA Bach To PV"}; // 0.04 in run2
@@ -3363,15 +3363,15 @@ struct HfTrackIndexSkimCreatorLfCascades {
33633363
registry.add("hDCACascDau", "hDCACascDau", {HistType::kTH1D, {{500, 0.0f, 5.0f, "cm^{2}"}}});
33643364
registry.add("hLambdaMass", "hLambdaMass", {HistType::kTH1D, {{400, 0.916f, 1.316f, "Inv. Mass (GeV/c^{2})"}}});
33653365

3366-
// pT rej
3367-
registry.add("hPtXicZeroOmegacZeroToXiPi", "Omegac/Xic to Xi Pi tracks selected by pT;#it{p}_{T}^{track} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3368-
registry.add("hPtCutsOmegacZeroToOmegaPi", "Omegac to Omega Pi tracks selected by pT;#it{p}_{T}^{track} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3369-
registry.add("hPtCutsOmegacZeroToOmegaKa", "Omegac to Omega Ka tracks selected by pT;#it{p}_{T}^{track} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3370-
registry.add("hPtCutsXicPlusToXiPiPi", "Xicplus to Xi Pi Pi tracks selected by pT;#it{p}_{T}^{track} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3371-
registry.add("hRejpTStatusXicZeroOmegacZeroToXiPi", "XicZeroOmegacZeroToXiPi rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}}); // pass dcafitter --> 0, pT>pTmin --> 1
3372-
registry.add("hRejpTStatusOmegacZeroToOmegaPi", "OmegacZeroToOmegaPi rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}});
3373-
registry.add("hRejpTStatusOmegacZeroToOmegaKa", "OmegacZeroToOmegaKa rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}});
3374-
registry.add("hRejpTStatusXicPlusToXiPiPi", "XicPlusToXiPiPi rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}});
3366+
//pT rej
3367+
registry.add("hPtCutsXicZeroOmegacZeroToXiPi", "Omegac/Xic to Xi Pi tracks selected by pT;#it{p}_{T} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3368+
registry.add("hPtCutsOmegacZeroToOmegaPi", "Omegac to Omega Pi tracks selected by pT;#it{p}_{T} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3369+
registry.add("hPtCutsOmegacZeroToOmegaKa", "Omegac to Omega Ka tracks selected by pT;#it{p}_{T} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3370+
registry.add("hPtCutsXicPlusToXiPiPi", "Xicplus to Xi Pi Pi tracks selected by pT;#it{p}_{T} (GeV/#it{c});entries", {HistType::kTH1D, {{500, 0., 50.}}});
3371+
registry.add("hRejpTStatusXicZeroOmegacZeroToXiPi", "XicZeroOmegacZeroToXiPi rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}}); //pass dcafitter --> 0, pT>pTmin --> 1
3372+
registry.add("hRejpTStatusOmegacZeroToOmegaPi", "OmegacZeroToOmegaPi rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}});
3373+
registry.add("hRejpTStatusOmegacZeroToOmegaKa", "OmegacZeroToOmegaKa rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}});
3374+
registry.add("hRejpTStatusXicPlusToXiPiPi", "XicPlusToXiPiPi rejected by pT status;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}});
33753375

33763376
// mass spectra
33773377
registry.add("hMassXicZeroOmegacZeroToXiPi", "2-prong candidates;inv. mass (#Xi #pi) (GeV/#it{c}^{2});entries", {HistType::kTH1D, {{500, 2., 3.}}});
@@ -3570,23 +3570,25 @@ struct HfTrackIndexSkimCreatorLfCascades {
35703570
std::array<float, 3> pVecPion1XiHyp = {0.};
35713571
df2.getTrack(0).getPxPyPzGlo(pVecXi);
35723572
df2.getTrack(1).getPxPyPzGlo(pVecPion1XiHyp);
3573-
float p_T_Xic = std::hypot(pVecXi[0] + pVecPion1XiHyp[0], pVecXi[1] + pVecPion1XiHyp[1]);
3573+
float ptXic = RecoDecay::pt(pVecXi, pVecPion1XiHyp);
35743574

35753575
std::array<std::array<float, 3>, 2> arrMomToXi = {pVecXi, pVecPion1XiHyp};
35763576
auto mass2ProngXiHyp = RecoDecay::m(arrMomToXi, arrMass2Prong[hf_cand_casc_lf::DecayType2Prong::XiczeroOmegaczeroToXiPi]);
35773577

3578-
if ((std::abs(casc.mXi() - massXi) < config.cascadeMassWindow) && (mass2ProngXiHyp >= config.massXiPiMin) && (mass2ProngXiHyp <= config.massXiPiMax)) {
3579-
registry.fill(HIST("hRejpTStatusXicZeroOmegacZeroToXiPi"), 0);
3580-
if (p_T_Xic >= config.ptMinXicZeroOmegacZeroToXiPiLfCasc) {
3578+
if ((std::abs(casc.mXi() - massXi) < config.cascadeMassWindow) && (mass2ProngXiHyp >= config.massXiPiMin) && (mass2ProngXiHyp <= config.massXiPiMax))
3579+
{
3580+
registry.fill(HIST("hRejpTStatusXicZeroOmegacZeroToXiPi"),0);
3581+
if (ptXic >= config.ptMinXicZeroOmegacZeroToXiPiLfCasc)
3582+
{
35813583
SETBIT(hfFlag, aod::hf_cand_casc_lf::DecayType2Prong::XiczeroOmegaczeroToXiPi);
3582-
registry.fill(HIST("hRejpTStatusXicZeroOmegacZeroToXiPi"), 1);
3584+
registry.fill(HIST("hRejpTStatusXicZeroOmegacZeroToXiPi"),1);
35833585
}
35843586
}
35853587

35863588
// fill histograms
35873589
if (config.fillHistograms && (TESTBIT(hfFlag, aod::hf_cand_casc_lf::DecayType2Prong::XiczeroOmegaczeroToXiPi))) {
35883590
registry.fill(HIST("hMassXicZeroOmegacZeroToXiPi"), mass2ProngXiHyp);
3589-
registry.fill(HIST("hPtXicZeroOmegacZeroToXiPi"), p_T_Xic);
3591+
registry.fill(HIST("hPtCutsXicZeroOmegacZeroToXiPi"), ptXic);
35903592
}
35913593
} else if (df2.isPropagationFailure()) {
35923594
LOGF(info, "Exception caught: failed to propagate tracks (2prong - xi) to charm baryon decay vtx");
@@ -3617,25 +3619,30 @@ struct HfTrackIndexSkimCreatorLfCascades {
36173619
std::array<float, 3> pVecCharmBachelor1OmegaHyp = {0.};
36183620
df2.getTrack(0).getPxPyPzGlo(pVecOmega);
36193621
df2.getTrack(1).getPxPyPzGlo(pVecCharmBachelor1OmegaHyp);
3620-
float p_T_Omegac = std::hypot(pVecOmega[0] + pVecCharmBachelor1OmegaHyp[0], pVecOmega[1] + pVecCharmBachelor1OmegaHyp[1]);
3622+
float ptOmegac = RecoDecay::pt(pVecOmega, pVecCharmBachelor1OmegaHyp);
36213623

36223624
std::array<std::array<float, 3>, 2> arrMomToOmega = {pVecOmega, pVecCharmBachelor1OmegaHyp};
36233625
auto mass2ProngOmegaPiHyp = RecoDecay::m(arrMomToOmega, arrMass2Prong[hf_cand_casc_lf::DecayType2Prong::OmegaczeroToOmegaPi]);
36243626
auto mass2ProngOmegaKHyp = RecoDecay::m(arrMomToOmega, arrMass2Prong[hf_cand_casc_lf::DecayType2Prong::OmegaczeroToOmegaK]);
36253627

3626-
if (std::abs(casc.mOmega() - massOmega) < config.cascadeMassWindow) {
3627-
if ((mass2ProngOmegaPiHyp >= config.massOmegaCharmBachelorMin) && (mass2ProngOmegaPiHyp <= config.massOmegaCharmBachelorMax)) {
3628-
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaPi"), 0);
3629-
if (p_T_Omegac >= config.ptMinOmegacZeroToOmegaPiLfCasc) {
3628+
if (std::abs(casc.mOmega() - massOmega) < config.cascadeMassWindow)
3629+
{
3630+
if ((mass2ProngOmegaPiHyp >= config.massOmegaCharmBachelorMin) && (mass2ProngOmegaPiHyp <= config.massOmegaCharmBachelorMax))
3631+
{
3632+
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaPi"),0);
3633+
if (ptOmegac >= config.ptMinOmegacZeroToOmegaPiLfCasc)
3634+
{
36303635
SETBIT(hfFlag, aod::hf_cand_casc_lf::DecayType2Prong::OmegaczeroToOmegaPi);
3631-
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaPi"), 1);
3636+
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaPi"),1);
36323637
}
36333638
}
3634-
if ((mass2ProngOmegaKHyp >= config.massOmegaCharmBachelorMin) && (mass2ProngOmegaKHyp <= config.massOmegaCharmBachelorMax)) {
3635-
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaKa"), 0);
3636-
if (p_T_Omegac >= config.ptMinOmegaczeroToOmegaKaLfCasc) {
3639+
if ((mass2ProngOmegaKHyp >= config.massOmegaCharmBachelorMin) && (mass2ProngOmegaKHyp <= config.massOmegaCharmBachelorMax))
3640+
{
3641+
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaKa"),0);
3642+
if (ptOmegac >= config.ptMinOmegaczeroToOmegaKaLfCasc)
3643+
{
36373644
SETBIT(hfFlag, aod::hf_cand_casc_lf::DecayType2Prong::OmegaczeroToOmegaK);
3638-
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaKa"), 1);
3645+
registry.fill(HIST("hRejpTStatusOmegacZeroToOmegaKa"),1);
36393646
}
36403647
}
36413648
}
@@ -3644,11 +3651,11 @@ struct HfTrackIndexSkimCreatorLfCascades {
36443651
if (config.fillHistograms) {
36453652
if (TESTBIT(hfFlag, aod::hf_cand_casc_lf::DecayType2Prong::OmegaczeroToOmegaPi)) {
36463653
registry.fill(HIST("hMassOmegacZeroToOmegaPi"), mass2ProngOmegaPiHyp);
3647-
registry.fill(HIST("hPtCutsOmegacZeroToOmegaPi"), p_T_Omegac);
3654+
registry.fill(HIST("hPtCutsOmegacZeroToOmegaPi"), ptOmegac);
36483655
}
36493656
if (TESTBIT(hfFlag, aod::hf_cand_casc_lf::DecayType2Prong::OmegaczeroToOmegaK)) {
36503657
registry.fill(HIST("hMassOmegacZeroToOmegaK"), mass2ProngOmegaKHyp);
3651-
registry.fill(HIST("hPtCutsOmegacZeroToOmegaKa"), p_T_Omegac);
3658+
registry.fill(HIST("hPtCutsOmegacZeroToOmegaKa"), ptOmegac);
36523659
}
36533660
}
36543661
} else if (df2.isPropagationFailure()) {
@@ -3721,23 +3728,24 @@ struct HfTrackIndexSkimCreatorLfCascades {
37213728
df3.getTrack(0).getPxPyPzGlo(pVec1); // take the momentum at the Xic vertex
37223729
df3.getTrack(1).getPxPyPzGlo(pVec2);
37233730
df3.getTrack(2).getPxPyPzGlo(pVec3);
3724-
float p_T_Xic3Prong = std::hypot(pVec1[0] + pVec2[0] + pVec3[0], pVec1[1] + pVec2[1] + pVec3[1]);
3731+
float ptXic3Prong = RecoDecay::pt(pVec1, pVec2, pVec3);
37253732

37263733
std::array<std::array<float, 3>, 3> arr3Mom = {pVec1, pVec2, pVec3};
37273734
auto mass3Prong = RecoDecay::m(arr3Mom, arrMass3Prong[hf_cand_casc_lf::DecayType3Prong::XicplusToXiPiPi]);
37283735

3729-
if ((std::abs(casc.mXi() - massXi) < config.cascadeMassWindow) && (mass3Prong >= config.massXiPiPiMin) && (mass3Prong <= config.massXiPiPiMax)) {
3730-
registry.fill(HIST("hRejpTStatusXicPlusToXiPiPi"), 0);
3731-
if (p_T_Xic3Prong >= config.ptMinXicplusLfCasc) {
3736+
if ((std::abs(casc.mXi() - massXi) < config.cascadeMassWindow) && (mass3Prong >= config.massXiPiPiMin) && (mass3Prong <= config.massXiPiPiMax))
3737+
{ registry.fill(HIST("hRejpTStatusXicPlusToXiPiPi"),0);
3738+
if (ptXic3Prong >= config.ptMinXicplusLfCasc)
3739+
{
37323740
SETBIT(hfFlag, aod::hf_cand_casc_lf::DecayType3Prong::XicplusToXiPiPi);
3733-
registry.fill(HIST("hRejpTStatusXicPlusToXiPiPi"), 1);
3741+
registry.fill(HIST("hRejpTStatusXicPlusToXiPiPi"),1);
37343742
}
37353743
}
37363744

37373745
// fill histograms
37383746
if (config.fillHistograms && (TESTBIT(hfFlag, aod::hf_cand_casc_lf::DecayType3Prong::XicplusToXiPiPi))) {
37393747
registry.fill(HIST("hMassXicPlusToXiPiPi"), mass3Prong);
3740-
registry.fill(HIST("hPtCutsXicPlusToXiPiPi"), p_T_Xic3Prong);
3748+
registry.fill(HIST("hPtCutsXicPlusToXiPiPi"), ptXic3Prong);
37413749
}
37423750
} else if (df3.isPropagationFailure()) {
37433751
LOGF(info, "Exception caught: failed to propagate tracks (3prong) to charm baryon decay vtx");

0 commit comments

Comments
 (0)