Skip to content

Commit f4047c3

Browse files
lmattei01lmattei01
committed
Generalize Jpis-D0 associated table maker to also include Dplus
1 parent 1fbe0d1 commit f4047c3

3 files changed

Lines changed: 248 additions & 57 deletions

File tree

‎PWGDQ/Core/VarManager.h‎

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1214,7 +1214,8 @@ class VarManager : public TObject
12141214
enum DileptonCharmHadronTypes {
12151215
kJPsi = 0,
12161216
kD0ToPiK,
1217-
kD0barToKPi
1217+
kD0barToKPi,
1218+
kDplusToPiKPi
12181219
};
12191220

12201221
enum EventFilters {
@@ -6608,6 +6609,13 @@ void VarManager::FillSingleDileptonCharmHadron(Cand const& candidate, H hfHelper
66086609
values[kRapCharmHadron] = hfHelper.yD0(candidate);
66096610
values[kBdtCharmHadron] = static_cast<float>(bdtScoreCharmHad);
66106611
}
6612+
if constexpr (partType == kDplusToPiKPi) {
6613+
values[kMassCharmHadron] = hfHelper.invMassDplusToPiKPi(candidate);
6614+
values[kPtCharmHadron] = candidate.pt();
6615+
values[kPhiCharmHadron] = candidate.phi();
6616+
values[kRapCharmHadron] = hfHelper.yDplus(candidate);
6617+
values[kBdtCharmHadron] = static_cast<float>(bdtScoreCharmHad);
6618+
}
66116619
}
66126620

66136621
template <int partTypeCharmHad, typename DQ, typename HF, typename H, typename T>

‎PWGDQ/DataModel/ReducedInfoTables.h‎

Lines changed: 34 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1338,32 +1338,44 @@ DECLARE_SOA_COLUMN(DeltaY, deltaY, float);
13381338
DECLARE_SOA_COLUMN(DeltaPhi, deltaPhi, float); //!
13391339
DECLARE_SOA_COLUMN(NumItsClsDmesProng0, numItsClsDmesProng0, int); //!
13401340
DECLARE_SOA_COLUMN(NumItsClsDmesProng1, numItsClsDmesProng1, int); //!
1341+
DECLARE_SOA_COLUMN(NumItsClsDmesProng2, numItsClsDmesProng2, int); //!
13411342
DECLARE_SOA_COLUMN(NumTpcCrossedRowsDmesProng0, numTpcCrossedRowsDmesProng0, int); //!
13421343
DECLARE_SOA_COLUMN(NumTpcCrossedRowsDmesProng1, numTpcCrossedRowsDmesProng1, int); //!
1344+
DECLARE_SOA_COLUMN(NumTpcCrossedRowsDmesProng2, numTpcCrossedRowsDmesProng2, int); //!
13431345
DECLARE_SOA_COLUMN(EtaDmesProng0, etaDmesProng0, float); //!
13441346
DECLARE_SOA_COLUMN(EtaDmesProng1, etaDmesProng1, float); //!
1347+
DECLARE_SOA_COLUMN(EtaDmesProng2, etaDmesProng2, float); //!
13451348
DECLARE_SOA_COLUMN(PtDmesProng0, ptDmesProng0, float); //!
13461349
DECLARE_SOA_COLUMN(PtDmesProng1, ptDmesProng1, float); //!
1350+
DECLARE_SOA_COLUMN(PtDmesProng2, ptDmesProng2, float); //!
13471351
DECLARE_SOA_COLUMN(MinNumItsClsDmesProng, minNumItsClsDmesProng, int); //!
13481352
DECLARE_SOA_COLUMN(MinNumTpcCrossedRowsDmesProng, minNumTpcCrossedRowsDmesProng, int); //!
13491353
DECLARE_SOA_COLUMN(MinAbsEtaDmesProng, minAbsEtaDmesProng, float); //!
13501354
DECLARE_SOA_COLUMN(MinPtDmesProng, minPtDmesProng, float); //!
13511355
DECLARE_SOA_COLUMN(NumSigmaTpcPiProng0, numSigmaTpcPiProng0, float); //!
13521356
DECLARE_SOA_COLUMN(NumSigmaTpcPiProng1, numSigmaTpcPiProng1, float); //!
1357+
DECLARE_SOA_COLUMN(NumSigmaTpcPiProng2, numSigmaTpcPiProng2, float); //!
13531358
DECLARE_SOA_COLUMN(NumSigmaTofPiProng0, numSigmaTofPiProng0, float); //!
13541359
DECLARE_SOA_COLUMN(NumSigmaTofPiProng1, numSigmaTofPiProng1, float); //!
1360+
DECLARE_SOA_COLUMN(NumSigmaTofPiProng2, numSigmaTofPiProng2, float); //!
13551361
DECLARE_SOA_COLUMN(NumSigmaTpcKaProng0, numSigmaTpcKaProng0, float); //!
13561362
DECLARE_SOA_COLUMN(NumSigmaTpcKaProng1, numSigmaTpcKaProng1, float); //!
1363+
DECLARE_SOA_COLUMN(NumSigmaTpcKaProng2, numSigmaTpcKaProng2, float); //!
13571364
DECLARE_SOA_COLUMN(NumSigmaTofKaProng0, numSigmaTofKaProng0, float); //!
13581365
DECLARE_SOA_COLUMN(NumSigmaTofKaProng1, numSigmaTofKaProng1, float); //!
1366+
DECLARE_SOA_COLUMN(NumSigmaTofKaProng2, numSigmaTofKaProng2, float); //!
13591367
DECLARE_SOA_DYNAMIC_COLUMN(NumSigmaTpcTofPiProng0, numSigmaTpcTofPiProng0, //!
13601368
[](float tpcNSigmaPi0, float tofNSigmaPi0) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaPi0, tofNSigmaPi0); });
13611369
DECLARE_SOA_DYNAMIC_COLUMN(NumSigmaTpcTofPiProng1, numSigmaTpcTofPiProng1, //!
13621370
[](float tpcNSigmaPi1, float tofNSigmaPi1) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaPi1, tofNSigmaPi1); });
1371+
DECLARE_SOA_DYNAMIC_COLUMN(NumSigmaTpcTofPiProng2, numSigmaTpcTofPiProng2, //!
1372+
[](float tpcNSigmaPi2, float tofNSigmaPi2) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaPi2, tofNSigmaPi2); });
13631373
DECLARE_SOA_DYNAMIC_COLUMN(NumSigmaTpcTofKaProng0, numSigmaTpcTofKaProng0, //!
13641374
[](float tpcNSigmaKa0, float tofNSigmaKa0) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaKa0, tofNSigmaKa0); });
13651375
DECLARE_SOA_DYNAMIC_COLUMN(NumSigmaTpcTofKaProng1, numSigmaTpcTofKaProng1, //!
13661376
[](float tpcNSigmaKa1, float tofNSigmaKa1) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaKa1, tofNSigmaKa1); });
1377+
DECLARE_SOA_DYNAMIC_COLUMN(NumSigmaTpcTofKaProng2, numSigmaTpcTofKaProng2, //!
1378+
[](float tpcNSigmaKa2, float tofNSigmaKa2) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaKa2, tofNSigmaKa2); });
13671379
} // namespace jpsidmescorr
13681380

13691381
DECLARE_SOA_TABLE(RedJpDmDileptons, "AOD", "REDJPDMDILEPTON", //!
@@ -1418,6 +1430,23 @@ DECLARE_SOA_TABLE(RedJpDmDmDau1s, "AOD", "REDJPDMDMDAU1", //!
14181430
jpsidmescorr::NumSigmaTpcTofPiProng1<jpsidmescorr::NumSigmaTpcPiProng1, jpsidmescorr::NumSigmaTofPiProng1>,
14191431
jpsidmescorr::NumSigmaTpcTofKaProng1<jpsidmescorr::NumSigmaTpcKaProng1, jpsidmescorr::NumSigmaTofKaProng1>);
14201432

1433+
DECLARE_SOA_TABLE(RedJpDmDmDau2s, "AOD", "REDJPDMDMDAU2", //!
1434+
jpsidmescorr::PtDmesProng2,
1435+
jpsidmescorr::EtaDmesProng2,
1436+
jpsidmescorr::NumItsClsDmesProng2,
1437+
jpsidmescorr::NumTpcCrossedRowsDmesProng2,
1438+
jpsidmescorr::NumSigmaTpcPiProng2,
1439+
jpsidmescorr::NumSigmaTofPiProng2,
1440+
jpsidmescorr::NumSigmaTpcKaProng2,
1441+
jpsidmescorr::NumSigmaTofKaProng2,
1442+
jpsidmescorr::NumSigmaTpcTofPiProng2<jpsidmescorr::NumSigmaTpcPiProng2, jpsidmescorr::NumSigmaTofPiProng2>,
1443+
jpsidmescorr::NumSigmaTpcTofKaProng2<jpsidmescorr::NumSigmaTpcKaProng2, jpsidmescorr::NumSigmaTofKaProng2>);
1444+
1445+
//General table for Dmeson mass
1446+
DECLARE_SOA_TABLE(RedJpDmDmesMasss, "AOD", "REDJPDMDMESMASS", //???
1447+
jpsidmescorr::MassDmes);
1448+
1449+
//Special table for D0 mass: particle and antiparticle separated
14211450
DECLARE_SOA_TABLE(RedJpDmD0Masss, "AOD", "REDJPDMD0MASS", //!
14221451
jpsidmescorr::MassD0,
14231452
jpsidmescorr::MassD0bar);
@@ -1430,6 +1459,11 @@ DECLARE_SOA_TABLE(RedJpDmDmesBdts, "AOD", "REDJPDMDMESBDT", //!
14301459
jpsidmescorr::BdtPromptMassHypo1,
14311460
jpsidmescorr::BdtNonpromptMassHypo1);
14321461

1462+
DECLARE_SOA_TABLE(RedJpDmDplusBdts, "AOD", "REDJPDMDPLUSBDT",
1463+
jpsidmescorr::BdtBkgMassHypo0,
1464+
jpsidmescorr::BdtPromptMassHypo0,
1465+
jpsidmescorr::BdtNonpromptMassHypo0);
1466+
14331467
DECLARE_SOA_TABLE(RedDleptDmesAll, "AOD", "RTDILPTDMESALL", //!
14341468
reducedpair::Mass,
14351469
jpsidmescorr::MassDmes,

0 commit comments

Comments
 (0)