Skip to content

Commit e27a654

Browse files
zhangbiao-phyzhangbiao-phy
committed
add more Pid and Lc mass in the derived table
1 parent 3e7c82c commit e27a654

File tree

2 files changed

+15
-16
lines changed

2 files changed

+15
-16
lines changed

‎PWGHF/D2H/Tasks/taskCd.cxx‎

Lines changed: 14 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -63,7 +63,8 @@ namespace o2::aod
6363
namespace full
6464
{
6565
// Candidate kinematics
66-
DECLARE_SOA_COLUMN(M, m, float); //! Invariant mass of candidate (GeV/c^2)
66+
DECLARE_SOA_COLUMN(MCd, mcd, float); //! Invariant mass of cd candidate (GeV/c^2)
67+
DECLARE_SOA_COLUMN(MLc, mlc, float); //! Invariant mass of lc candidate (GeV/c^2)
6768
DECLARE_SOA_COLUMN(Pt, pt, float); //! Transverse momentum of candidate (GeV/c)
6869
DECLARE_SOA_COLUMN(Eta, eta, float); //! eta of candidate (GeV/c)
6970
DECLARE_SOA_COLUMN(Phi, phi, float); //! phi of candidate (GeV/c)
@@ -100,7 +101,8 @@ DECLARE_SOA_COLUMN(TimeStamp, timeStamp, int64_t); //! Timestamp fo
100101

101102
// Full table: include ALL columns declared above
102103
DECLARE_SOA_TABLE(HfCandCd, "AOD", "HFCANDCD",
103-
full::M,
104+
full::MCd,
105+
full::MLc,
104106
full::Pt,
105107
full::Eta,
106108
full::Phi,
@@ -116,13 +118,10 @@ DECLARE_SOA_TABLE(HfCandCd, "AOD", "HFCANDCD",
116118
full::NSigmaTpcDe,
117119
full::NSigmaItsDe,
118120
full::NSigmaTofDe,
119-
full::NItsClusters,
120121
full::NItsNClusterSize,
121-
full::NTpcClusters,
122122
full::NTpcSignalsDe,
123123
full::NTpcSignalsPi,
124124
full::NTpcSignalsKa,
125-
full::NItsSignalsDe,
126125
full::CandidateSelFlag,
127126
full::Cent);
128127
} // namespace o2::aod
@@ -219,6 +218,7 @@ struct HfTaskCd {
219218
registry.add("Data/hImpParErrProng1", "3-prong candidates;prong 1 impact parameter error (cm);entries", {HistType::kTH2F, {{100, -1., 1.}, {binsPt, "#it{p}_{T} (GeV/#it{c})"}}});
220219
registry.add("Data/hImpParErrProng2", "3-prong candidates;prong 2 impact parameter error (cm);entries", {HistType::kTH2F, {{100, -1., 1.}, {binsPt, "#it{p}_{T} (GeV/#it{c})"}}});
221220
registry.add("Data/hNsigmaTPCDeVsP", "deuteron;#it{p} (GeV/#it{c}); n#sigma^{TPC}_{d}", {HistType::kTH2F, {{200, -10.f, 10.f}, {200, -6.f, 6.f}}});
221+
registry.add("Data/hNsigmaTPCPrVsP", "proton;#it{p} (GeV/#it{c}); n#sigma^{TPC}_{p}", {HistType::kTH2F, {{200, -10.f, 10.f}, {200, -6.f, 6.f}}});
222222
registry.add("Data/hNsigmaTOFDeVsP", "deuteron;#it{p} (GeV/#it{c}); n#sigma^{TOF}_{d}", {HistType::kTH2F, {{200, -10.f, 10.f}, {200, -6.f, 6.f}}});
223223
registry.add("Data/hNsigmaITSDeVsP", "deuteron;#it{p} (GeV/#it{c}); n#sigma^{ITS}_{d}", {HistType::kTH2F, {{200, -10.f, 10.f}, {200, -6.f, 6.f}}});
224224
registry.add("Data/hTPCSignalDeVsP", "deuteron;#it{p} (GeV/#it{c}); TPC signals", {HistType::kTH2F, {{200, -10.f, 10.f}, {2000, 0, 2000}}});
@@ -290,11 +290,14 @@ struct HfTaskCd {
290290
const auto cpa = candidate.cpa();
291291
const auto cpaXY = candidate.cpaXY();
292292
float invMassCd = 0.f;
293+
float invMassLc = 0.f;
293294
if (candidate.isSelCdToDeKPi() >= selectionFlagCd) {
294295
invMassCd = HfHelper::invMassCdToDeKPi(candidate);
296+
invMassLc = HfHelper::invMassLcToPKPi(candidate);
295297
}
296298
if (candidate.isSelCdToPiKDe() >= selectionFlagCd) {
297299
invMassCd = HfHelper::invMassCdToPiKDe(candidate);
300+
invMassLc = HfHelper::invMassLcToPiKP(candidate);
298301
}
299302

300303
if (candidate.isSelCdToDeKPi() >= selectionFlagCd) {
@@ -356,13 +359,11 @@ struct HfTaskCd {
356359
if (fillTree) {
357360
int candFlag = -999;
358361

359-
float nSigmaTpcDe = 0.f, nSigmaTpcKa = 0.f, nSigmaTpcPi = 0.f;
362+
float nSigmaTpcDe = 0.f, nSigmaTpcKa = 0.f, nSigmaTpcPi = 0.f, nSigmaTpcPr = 0.f;
360363
float nSigmaItsDe = 0.f;
361364
float nSigmaTofDe = 0.f, nSigmaTofKa = 0.f, nSigmaTofPi = 0.f;
362365

363-
int itsNClusterDe = 0;
364366
int itsNClusterSizeDe = 0;
365-
int tpcNClusterDe = 0;
366367

367368
float tpcSignalsDe = 0.f;
368369
float tpcSignalsPi = 0.f;
@@ -393,13 +394,12 @@ struct HfTaskCd {
393394
pSignedDe = prong0.tpcInnerParam() * prong0.sign();
394395
pSignedPi = prong2.tpcInnerParam() * prong2.sign();
395396
nSigmaTpcDe = candidate.nSigTpcDe0();
397+
nSigmaTpcPr = candidate.nSigTpcPr0();
396398
nSigmaTofDe = candidate.nSigTofDe0();
397399
nSigmaTpcPi = candidate.nSigTpcPi2();
398400
nSigmaTofPi = candidate.nSigTofPi2();
399401
nSigmaItsDe = prong0Its.itsNSigmaDe();
400-
itsNClusterDe = prong0.itsNCls();
401402
itsNClusterSizeDe = prong0.itsClusterSizes();
402-
tpcNClusterDe = prong0.tpcNClsCrossedRows();
403403
tpcSignalsDe = prong0.tpcSignal();
404404
tpcSignalsPi = prong2.tpcSignal();
405405
itsSignalsDe = itsSignal(prong0);
@@ -408,20 +408,20 @@ struct HfTaskCd {
408408
pSignedDe = prong2.tpcInnerParam() * prong2.sign();
409409
pSignedPi = prong0.tpcInnerParam() * prong0.sign();
410410
nSigmaTpcDe = candidate.nSigTpcDe2();
411+
nSigmaTpcPr = candidate.nSigTpcPr2();
411412
nSigmaTofDe = candidate.nSigTofDe2();
412413
nSigmaTpcPi = candidate.nSigTpcPi0();
413414
nSigmaTofPi = candidate.nSigTofPi0();
414415
nSigmaItsDe = prong2Its.itsNSigmaDe();
415-
itsNClusterDe = prong2.itsNCls();
416416
itsNClusterSizeDe = prong2.itsClusterSizes();
417-
tpcNClusterDe = prong2.tpcNClsCrossedRows();
418417
tpcSignalsDe = prong2.tpcSignal();
419418
tpcSignalsPi = prong0.tpcSignal();
420419
itsSignalsDe = itsSignal(prong2);
421420
}
422421

423422
// PID QA
424423
registry.fill(HIST("Data/hNsigmaTPCDeVsP"), pSignedDe, nSigmaTpcDe);
424+
registry.fill(HIST("Data/hNsigmaTPCPrVsP"), pSignedDe, nSigmaTpcPr);
425425
registry.fill(HIST("Data/hNsigmaTOFDeVsP"), pSignedDe, nSigmaTofDe);
426426
registry.fill(HIST("Data/hNsigmaITSDeVsP"), pSignedDe, nSigmaItsDe);
427427
registry.fill(HIST("Data/hTPCSignalDeVsP"), pSignedDe, tpcSignalsDe);
@@ -435,6 +435,7 @@ struct HfTaskCd {
435435

436436
rowCandCd(
437437
invMassCd,
438+
invMassLc,
438439
pt,
439440
eta,
440441
phi,
@@ -448,15 +449,13 @@ struct HfTaskCd {
448449
cpa,
449450
chi2PCA,
450451
nSigmaTpcDe,
452+
nSigmaTpcPr,
451453
nSigmaItsDe,
452454
nSigmaTofDe,
453-
itsNClusterDe,
454455
itsNClusterSizeDe,
455-
tpcNClusterDe,
456456
tpcSignalsDe,
457457
tpcSignalsPi,
458458
tpcSignalsKa,
459-
itsSignalsDe,
460459
candFlag,
461460
cent);
462461
}

‎PWGHF/DataModel/CandidateReconstructionTables.h‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -837,7 +837,7 @@ DECLARE_SOA_TABLE(HfCand3Prong2PidAl, "AOD", "HFCAND3P2PIDAL", //!
837837
using HfCand3Prong = HfCand3ProngExt;
838838
using HfCand3ProngWPidPiKaPr = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidPr, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidPr, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidPr, HfCand3Prong2PidKa>;
839839
using HfCand3ProngWPidPiKa = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidKa>;
840-
using HfCand3ProngWPidPiKaDe = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidDe, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidDe, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidDe, HfCand3Prong2PidKa>;
840+
using HfCand3ProngWPidPiKaDe = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidDe, HfCand3Prong0PidPr, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidDe, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidDe, HfCand3Prong2PidPr, HfCand3Prong2PidKa>;
841841
using HfCand3ProngWPidPiKaTr = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidTr, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidTr, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidTr, HfCand3Prong2PidKa>;
842842
using HfCand3ProngWPidPiKaHe = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidHe, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidHe, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidHe, HfCand3Prong2PidKa>;
843843
using HfCand3ProngWPidPiKaAl = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidAl, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidAl, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidAl, HfCand3Prong2PidKa>;

0 commit comments

Comments
 (0)