Skip to content

Commit 2abc9d2

Browse files
creetz16alibuild
creetz16
and
alibuild
authored
[PWGLF] Update KF reconstruction of hypertriton 3-body analysis (#8368)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 7a6bcb5 commit 2abc9d2

File tree

3 files changed

+129
-71
lines changed

3 files changed

+129
-71
lines changed

‎PWGLF/DataModel/Vtx3BodyTables.h‎

Lines changed: 34 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -382,14 +382,19 @@ DECLARE_SOA_COLUMN(DCAVtxDaughtersKF, dcavtxdaughterskf, float); //! sum
382382
DECLARE_SOA_COLUMN(Track0Sign, track0sign, float); //! sign of proton daughter track
383383
DECLARE_SOA_COLUMN(Track1Sign, track1sign, float); //! sign of pion daughter track
384384
DECLARE_SOA_COLUMN(Track2Sign, track2sign, float); //! sign of deuteron daughter track
385+
DECLARE_SOA_COLUMN(TPCInnerParamTrack0, tpcinnerparamtrack0, float); //! momentum at inner wall of TPC of proton daughter
386+
DECLARE_SOA_COLUMN(TPCInnerParamTrack1, tpcinnerparamtrack1, float); //! momentum at inner wall of TPC of pion daughter
387+
DECLARE_SOA_COLUMN(TPCInnerParamTrack2, tpcinnerparamtrack2, float); //! momentum at inner wall of TPC of deuteron daughter
385388
// PID
386-
DECLARE_SOA_COLUMN(TPCNSigmaProton, tpcnsigmaproton, float); //! nsigma of TPC PID of the proton daughter
387-
DECLARE_SOA_COLUMN(TPCNSigmaPion, tpcnsigmapion, float); //! nsigma of TPC PID of the pion daughter
388-
DECLARE_SOA_COLUMN(TPCNSigmaDeuteron, tpcnsigmadeuteron, float); //! nsigma of TPC PID of the bachelor daughter
389-
DECLARE_SOA_COLUMN(TPCdEdxProton, tpcdedxproton, float); //! TPC dEdx of the proton daughter
390-
DECLARE_SOA_COLUMN(TPCdEdxPion, tpcdedxpion, float); //! TPC dEdx of the pion daughter
391-
DECLARE_SOA_COLUMN(TPCdEdxDeuteron, tpcdedxdeuteron, float); //! TPC dEdx of the bachelor daughter
392-
DECLARE_SOA_COLUMN(TOFNSigmaDeuteron, tofnsigmadeuteron, float); //! nsigma of TOF PID of the bachelor daughter
389+
DECLARE_SOA_COLUMN(TPCNSigmaProton, tpcnsigmaproton, float); //! nsigma of TPC PID of the proton daughter
390+
DECLARE_SOA_COLUMN(TPCNSigmaPion, tpcnsigmapion, float); //! nsigma of TPC PID of the pion daughter
391+
DECLARE_SOA_COLUMN(TPCNSigmaDeuteron, tpcnsigmadeuteron, float); //! nsigma of TPC PID of the bachelor daughter
392+
DECLARE_SOA_COLUMN(TPCdEdxProton, tpcdedxproton, float); //! TPC dEdx of the proton daughter
393+
DECLARE_SOA_COLUMN(TPCdEdxPion, tpcdedxpion, float); //! TPC dEdx of the pion daughter
394+
DECLARE_SOA_COLUMN(TPCdEdxDeuteron, tpcdedxdeuteron, float); //! TPC dEdx of the bachelor daughter
395+
DECLARE_SOA_COLUMN(TOFNSigmaDeuteron, tofnsigmadeuteron, float); //! nsigma of TOF PID of the bachelor daughter
396+
DECLARE_SOA_COLUMN(ITSClusSizeDeuteron, itsclussizedeuteron, double); //! average ITS cluster size of bachelor daughter
397+
DECLARE_SOA_COLUMN(PIDTrackingDeuteron, pidtrackingdeuteron, uint32_t); //! PID during tracking of bachelor daughter
393398

394399
// Monte Carlo
395400
DECLARE_SOA_COLUMN(GenP, genp, float); //! generated momentum
@@ -410,6 +415,7 @@ DECLARE_SOA_COLUMN(IsReco, isreco, int); //! flag if candidate wa
410415
// V0
411416
DECLARE_SOA_COLUMN(MassV0, massv0, float); //! proton, pion vertex mass
412417
DECLARE_SOA_COLUMN(Chi2MassV0, chi2massv0, float); //! chi2 of proton, pion mass constraint to Lambda mass
418+
DECLARE_SOA_COLUMN(CosPAV0, cospav0, float); //! proton, pion vertex mass
413419

414420
} // namespace kfvtx3body
415421

@@ -432,11 +438,13 @@ DECLARE_SOA_TABLE(KFVtx3BodyDatas, "AOD", "KFVTX3BODYDATA",
432438

433439
// V0
434440
kfvtx3body::MassV0, kfvtx3body::Chi2MassV0,
441+
kfvtx3body::CosPAV0,
435442

436443
// daughters
437-
vtx3body::PxTrack0, vtx3body::PyTrack0, vtx3body::PzTrack0, // proton
438-
vtx3body::PxTrack1, vtx3body::PyTrack1, vtx3body::PzTrack1, // pion
439-
vtx3body::PxTrack2, vtx3body::PyTrack2, vtx3body::PzTrack2, // deuteron
444+
vtx3body::PxTrack0, vtx3body::PyTrack0, vtx3body::PzTrack0, // proton
445+
vtx3body::PxTrack1, vtx3body::PyTrack1, vtx3body::PzTrack1, // pion
446+
vtx3body::PxTrack2, vtx3body::PyTrack2, vtx3body::PzTrack2, // deuteron
447+
kfvtx3body::TPCInnerParamTrack0, kfvtx3body::TPCInnerParamTrack1, kfvtx3body::TPCInnerParamTrack2, // proton, pion, deuteron
440448
kfvtx3body::DCATrack0ToPVKF, kfvtx3body::DCATrack1ToPVKF, kfvtx3body::DCATrack2ToPVKF, kfvtx3body::DCAxyTrack0ToPVKF, kfvtx3body::DCAxyTrack1ToPVKF, kfvtx3body::DCAxyTrack2ToPVKF,
441449
kfvtx3body::DCAxyTrack0ToSVKF, kfvtx3body::DCAxyTrack1ToSVKF, kfvtx3body::DCAxyTrack2ToSVKF,
442450
kfvtx3body::DCAxyTrack0ToTrack1KF, kfvtx3body::DCAxyTrack0ToTrack2KF, kfvtx3body::DCAxyTrack1ToTrack2KF,
@@ -447,6 +455,8 @@ DECLARE_SOA_TABLE(KFVtx3BodyDatas, "AOD", "KFVTX3BODYDATA",
447455
kfvtx3body::TPCNSigmaProton, kfvtx3body::TPCNSigmaPion, kfvtx3body::TPCNSigmaDeuteron,
448456
kfvtx3body::TPCdEdxProton, kfvtx3body::TPCdEdxPion, kfvtx3body::TPCdEdxDeuteron,
449457
kfvtx3body::TOFNSigmaDeuteron,
458+
kfvtx3body::ITSClusSizeDeuteron,
459+
kfvtx3body::PIDTrackingDeuteron,
450460

451461
// dynamic columns
452462
vtx3body::VtxRadius<vtx3body::X, vtx3body::Y>,
@@ -486,33 +496,33 @@ DECLARE_SOA_TABLE(KFVtx3BodyDatasLite, "AOD", "KF3BODYLITE",
486496
// hypertriton candidate
487497
kfvtx3body::Mass,
488498
vtx3body::X, vtx3body::Y, vtx3body::Z,
489-
kfvtx3body::XErr, kfvtx3body::YErr, kfvtx3body::ZErr,
490499
kfvtx3body::Px, kfvtx3body::Py, kfvtx3body::Pz, kfvtx3body::Pt,
491500
kfvtx3body::Sign,
492501
kfvtx3body::DCAVtxToPVKF, kfvtx3body::DCAXYVtxToPVKF,
493502
kfvtx3body::VtxCosPAKF, kfvtx3body::VtxCosXYPAKF,
494-
kfvtx3body::VtxCosPAKFtopo, kfvtx3body::VtxCosXYPAKFtopo,
495503
kfvtx3body::DecayLKF, kfvtx3body::DecayLXYKF, kfvtx3body::DecayLDeltaL,
496504
kfvtx3body::Chi2geoNDF, kfvtx3body::Chi2topoNDF,
497505
kfvtx3body::CTauKFtopo,
498506

499507
// V0
500508
kfvtx3body::MassV0, kfvtx3body::Chi2MassV0,
509+
kfvtx3body::CosPAV0,
501510

502511
// daughters
503-
vtx3body::PxTrack0, vtx3body::PyTrack0, vtx3body::PzTrack0, // proton
504-
vtx3body::PxTrack1, vtx3body::PyTrack1, vtx3body::PzTrack1, // pion
505-
vtx3body::PxTrack2, vtx3body::PyTrack2, vtx3body::PzTrack2, // deuteron
512+
vtx3body::PxTrack0, vtx3body::PyTrack0, vtx3body::PzTrack0, // proton
513+
vtx3body::PxTrack1, vtx3body::PyTrack1, vtx3body::PzTrack1, // pion
514+
vtx3body::PxTrack2, vtx3body::PyTrack2, vtx3body::PzTrack2, // deuteron
515+
kfvtx3body::TPCInnerParamTrack0, kfvtx3body::TPCInnerParamTrack1, kfvtx3body::TPCInnerParamTrack2, // proton, pion, deuteron
506516
kfvtx3body::DCATrack0ToPVKF, kfvtx3body::DCATrack1ToPVKF, kfvtx3body::DCATrack2ToPVKF, kfvtx3body::DCAxyTrack0ToPVKF, kfvtx3body::DCAxyTrack1ToPVKF, kfvtx3body::DCAxyTrack2ToPVKF,
507517
kfvtx3body::DCAxyTrack0ToSVKF, kfvtx3body::DCAxyTrack1ToSVKF, kfvtx3body::DCAxyTrack2ToSVKF,
508518
kfvtx3body::DCAxyTrack0ToTrack1KF, kfvtx3body::DCAxyTrack0ToTrack2KF, kfvtx3body::DCAxyTrack1ToTrack2KF,
509519
kfvtx3body::DCAVtxDaughtersKF,
510-
kfvtx3body::DCAxyTrackPosToPV, kfvtx3body::DCAxyTrackNegToPV, kfvtx3body::DCAxyTrackBachToPV,
511-
kfvtx3body::DCATrackPosToPV, kfvtx3body::DCATrackNegToPV, kfvtx3body::DCATrackBachToPV,
512520
kfvtx3body::Track0Sign, kfvtx3body::Track1Sign, kfvtx3body::Track2Sign, // track sing: proton, pion, deuteron
513521
kfvtx3body::TPCNSigmaProton, kfvtx3body::TPCNSigmaPion, kfvtx3body::TPCNSigmaDeuteron,
514522
kfvtx3body::TPCdEdxProton, kfvtx3body::TPCdEdxPion, kfvtx3body::TPCdEdxDeuteron,
515523
kfvtx3body::TOFNSigmaDeuteron,
524+
kfvtx3body::ITSClusSizeDeuteron,
525+
kfvtx3body::PIDTrackingDeuteron,
516526

517527
// dynamic columns
518528
vtx3body::VtxRadius<vtx3body::X, vtx3body::Y>,
@@ -555,11 +565,13 @@ DECLARE_SOA_TABLE(McKFVtx3BodyDatas, "AOD", "MCKF3BODYDATAS",
555565

556566
// V0
557567
kfvtx3body::MassV0, kfvtx3body::Chi2MassV0,
568+
kfvtx3body::CosPAV0,
558569

559570
// daughters
560-
vtx3body::PxTrack0, vtx3body::PyTrack0, vtx3body::PzTrack0, // proton
561-
vtx3body::PxTrack1, vtx3body::PyTrack1, vtx3body::PzTrack1, // pion
562-
vtx3body::PxTrack2, vtx3body::PyTrack2, vtx3body::PzTrack2, // deuteron
571+
vtx3body::PxTrack0, vtx3body::PyTrack0, vtx3body::PzTrack0, // proton
572+
vtx3body::PxTrack1, vtx3body::PyTrack1, vtx3body::PzTrack1, // pion
573+
vtx3body::PxTrack2, vtx3body::PyTrack2, vtx3body::PzTrack2, // deuteron
574+
kfvtx3body::TPCInnerParamTrack0, kfvtx3body::TPCInnerParamTrack1, kfvtx3body::TPCInnerParamTrack2, // proton, pion, deuteron
563575
kfvtx3body::DCATrack0ToPVKF, kfvtx3body::DCATrack1ToPVKF, kfvtx3body::DCATrack2ToPVKF, kfvtx3body::DCAxyTrack0ToPVKF, kfvtx3body::DCAxyTrack1ToPVKF, kfvtx3body::DCAxyTrack2ToPVKF,
564576
kfvtx3body::DCAxyTrack0ToSVKF, kfvtx3body::DCAxyTrack1ToSVKF, kfvtx3body::DCAxyTrack2ToSVKF,
565577
kfvtx3body::DCAxyTrack0ToTrack1KF, kfvtx3body::DCAxyTrack0ToTrack2KF, kfvtx3body::DCAxyTrack1ToTrack2KF,
@@ -570,6 +582,8 @@ DECLARE_SOA_TABLE(McKFVtx3BodyDatas, "AOD", "MCKF3BODYDATAS",
570582
kfvtx3body::TPCNSigmaProton, kfvtx3body::TPCNSigmaPion, kfvtx3body::TPCNSigmaDeuteron,
571583
kfvtx3body::TPCdEdxProton, kfvtx3body::TPCdEdxPion, kfvtx3body::TPCdEdxDeuteron,
572584
kfvtx3body::TOFNSigmaDeuteron,
585+
kfvtx3body::ITSClusSizeDeuteron,
586+
kfvtx3body::PIDTrackingDeuteron,
573587

574588
// MC information
575589
kfvtx3body::GenP,

‎PWGLF/TableProducer/Nuspex/decay3bodybuilder.cxx‎

Lines changed: 58 additions & 26 deletions
Original file line numberDiff line numberDiff line change
@@ -16,6 +16,8 @@
1616
#include <cmath>
1717
#include <array>
1818
#include <cstdlib>
19+
#include <string>
20+
#include <vector>
1921

2022
#include "Framework/runDataProcessing.h"
2123
#include "Framework/AnalysisTask.h"
@@ -765,6 +767,17 @@ struct decay3bodyBuilder {
765767
tofNSigmaDeuteron = bachelorTOFPID.GetTOFNSigma(trackBach, originalcol, collision);
766768
}
767769

770+
// Average ITS cluster size of deuteron track
771+
double averageClusterSizeDeuteron(0);
772+
int nCls(0);
773+
for (int i = 0; i < 7; i++) {
774+
int clusterSize = trackBach.itsClsSizeInLayer(i);
775+
averageClusterSizeDeuteron += static_cast<double>(clusterSize);
776+
if (clusterSize > 0)
777+
nCls++;
778+
}
779+
averageClusterSizeDeuteron = averageClusterSizeDeuteron / static_cast<double>(nCls);
780+
768781
// track DCAxy and DCAz to PV associated with decay3body
769782
o2::dataformats::VertexBase mPV;
770783
o2::dataformats::DCA mDcaInfoCovPos;
@@ -801,6 +814,18 @@ struct decay3bodyBuilder {
801814
auto TrackNegDca = std::sqrt(TrackNegDcaXY * TrackNegDcaXY + TrackNegDcaZ * TrackNegDcaZ);
802815
auto TrackBachDca = std::sqrt(TrackBachDcaXY * TrackBachDcaXY + TrackBachDcaZ * TrackBachDcaZ);
803816

817+
// daughter track momentum at inner wall of TPC
818+
float tpcInnerParamProton;
819+
float tpcInnerParamPion;
820+
float tpcInnerParamDeuteron = trackBach.tpcInnerParam();
821+
if (isMatter) { // hypertriton (proton, pi-, deuteron)
822+
tpcInnerParamProton = trackPos.tpcInnerParam();
823+
tpcInnerParamPion = trackNeg.tpcInnerParam();
824+
} else if (!isMatter) { // anti-hypertriton (anti-proton, pi+, deuteron)
825+
tpcInnerParamProton = trackNeg.tpcInnerParam();
826+
tpcInnerParamPion = trackPos.tpcInnerParam();
827+
}
828+
804829
// -------- STEP 2: fit vertex with proton and pion --------
805830
// Fit vertex with DCA fitter to find minimization point --> uses material corrections implicitly
806831
if (kfparticleConfigurations.doDCAFitterPreMinimum) {
@@ -829,13 +854,13 @@ struct decay3bodyBuilder {
829854
}
830855
LOG(debug) << "KFParticle objects created from daughter tracks.";
831856

832-
// Construct V0
857+
// Construct V0 as intermediate step
833858
KFParticle KFV0;
834-
int nDaughters = 2;
835-
const KFParticle* Daughters[2] = {&kfpProton, &kfpPion};
859+
int nDaughtersV0 = 2;
860+
const KFParticle* DaughtersV0[2] = {&kfpProton, &kfpPion};
836861
KFV0.SetConstructMethod(2);
837862
try {
838-
KFV0.Construct(Daughters, nDaughters);
863+
KFV0.Construct(DaughtersV0, nDaughtersV0);
839864
} catch (std::runtime_error& e) {
840865
LOG(debug) << "Failed to create V0 vertex from daughter tracks." << e.what();
841866
continue;
@@ -855,31 +880,35 @@ struct decay3bodyBuilder {
855880
}
856881
registry.fill(HIST("hVtx3BodyCounterKFParticle"), kKfVtxV0MassConst);
857882

858-
// -------- STEP 3: fit vertex with V0 and deuteron --------
883+
// -------- STEP 3: fit three body vertex --------
859884
// Create KFParticle object from deuteron track
860885
KFParticle kfpDeuteron;
861886
kfpDeuteron = createKFParticleFromTrackParCov(trackParCovBach, trackBach.sign() * bachelorcharge, constants::physics::MassDeuteron);
862887
LOG(debug) << "KFParticle created from deuteron track.";
863-
// Add deuteron to V0 vertex
888+
// Construct 3body vertex
889+
int nDaughters3body = 3;
890+
const KFParticle* Daughters3body[3] = {&kfpProton, &kfpPion, &kfpDeuteron};
864891
KFParticle KFHt;
865-
KFHt = KFV0;
866892
KFHt.SetConstructMethod(2);
867893
try {
868-
KFHt.AddDaughter(kfpDeuteron);
894+
KFHt.Construct(Daughters3body, nDaughters3body);
869895
} catch (std::runtime_error& e) {
870-
LOG(debug) << "Failed to create Hyper triton from V0 and deuteron." << e.what();
896+
LOG(debug) << "Failed to create Hyper triton 3-body vertex." << e.what();
871897
continue;
872898
}
899+
// transport all daughter tracks to hypertriton vertex
900+
float HtVtx[3] = {0.};
901+
HtVtx[0] = KFHt.GetX();
902+
HtVtx[1] = KFHt.GetY();
903+
HtVtx[2] = KFHt.GetZ();
904+
kfpProton.TransportToPoint(HtVtx);
905+
kfpPion.TransportToPoint(HtVtx);
906+
kfpDeuteron.TransportToPoint(HtVtx);
873907
registry.fill(HIST("hVtx3BodyCounterKFParticle"), kKfVtxhasSV);
874908
LOG(debug) << "Hypertriton vertex constructed.";
875909

876910
// -------- STEP 4: selections after geometrical vertex fit --------
877-
// Get updated daughter tracks
878-
kfpProton.SetProductionVertex(KFHt);
879-
kfpPion.SetProductionVertex(KFHt);
880-
kfpDeuteron.SetProductionVertex(KFHt);
881-
LOG(debug) << "Topo constraint applied to daughters.";
882-
// daughter DCAs
911+
// daughter DCAs with KF
883912
if ((kfpProton.GetDistanceFromParticle(kfpPion) >= kfparticleConfigurations.maxDcaProPi) || (kfpProton.GetDistanceFromParticle(kfpDeuteron) >= kfparticleConfigurations.maxDcaProDeu) || (kfpPion.GetDistanceFromParticle(kfpDeuteron) >= kfparticleConfigurations.maxDcaPiDe)) {
884913
continue;
885914
}
@@ -988,10 +1017,13 @@ struct decay3bodyBuilder {
9881017
KFHtPV.GetLifeTime(),
9891018
// V0
9901019
massV0, chi2massV0,
991-
// daughter momenta
1020+
cpaFromKF(KFV0, kfpv),
1021+
// daughter momenta at vertex
9921022
kfpProton.GetPx(), kfpProton.GetPy(), kfpProton.GetPz(),
9931023
kfpPion.GetPx(), kfpPion.GetPy(), kfpPion.GetPz(),
9941024
kfpDeuteron.GetPx(), kfpDeuteron.GetPy(), kfpDeuteron.GetPz(),
1025+
// daughter momenta at inner wall of TPC
1026+
tpcInnerParamProton, tpcInnerParamPion, tpcInnerParamDeuteron,
9951027
// daughter DCAs KF
9961028
kfpProton.GetDistanceFromVertex(kfpv),
9971029
kfpPion.GetDistanceFromVertex(kfpv),
@@ -1021,31 +1053,33 @@ struct decay3bodyBuilder {
10211053
dEdxProton,
10221054
dEdxPion,
10231055
dEdxDeuteron,
1024-
tofNSigmaDeuteron);
1056+
tofNSigmaDeuteron,
1057+
averageClusterSizeDeuteron,
1058+
trackBach.pidForTracking());
10251059

10261060
if (kfparticleConfigurations.fillCandidateLiteTable) {
10271061
kfvtx3bodydatalite(
10281062
// hypertriton
10291063
massHt,
10301064
KFHt.GetX(), KFHt.GetY(), KFHt.GetZ(),
1031-
KFHt.GetErrX(), KFHt.GetErrY(), KFHt.GetErrZ(),
10321065
KFHt.GetPx(), KFHt.GetPy(), KFHt.GetPz(), KFHt.GetPt(),
10331066
KFHt.GetQ(),
10341067
KFHt.GetDistanceFromVertex(kfpv), KFHt.GetDistanceFromVertexXY(kfpv),
10351068
cpaFromKF(KFHt, kfpv), // before topo constraint
10361069
cpaXYFromKF(KFHt, kfpv),
1037-
cpaFromKF(KFHtPV, kfpv), // after topo constraint
1038-
cpaXYFromKF(KFHtPV, kfpv),
10391070
KFHtPV.GetDecayLength(), KFHtPV.GetDecayLengthXY(), // decay length defined after topological constraint
10401071
KFHtPV.GetDecayLength() / KFHtPV.GetErrDecayLength(), // ldl
10411072
chi2geoNDF, chi2topoNDF,
10421073
KFHtPV.GetLifeTime(),
10431074
// V0
10441075
massV0, chi2massV0,
1045-
// daughter momenta
1076+
cpaFromKF(KFV0, kfpv),
1077+
// daughter momenta at vertex
10461078
kfpProton.GetPx(), kfpProton.GetPy(), kfpProton.GetPz(),
10471079
kfpPion.GetPx(), kfpPion.GetPy(), kfpPion.GetPz(),
10481080
kfpDeuteron.GetPx(), kfpDeuteron.GetPy(), kfpDeuteron.GetPz(),
1081+
// daughter momenta at inner wall of TPC
1082+
tpcInnerParamProton, tpcInnerParamPion, tpcInnerParamDeuteron,
10491083
// daughter DCAs KF
10501084
kfpProton.GetDistanceFromVertex(kfpv),
10511085
kfpPion.GetDistanceFromVertex(kfpv),
@@ -1060,10 +1094,6 @@ struct decay3bodyBuilder {
10601094
kfpProton.GetDistanceFromParticle(kfpDeuteron),
10611095
kfpPion.GetDistanceFromParticle(kfpDeuteron),
10621096
DCAvtxDaughters3D,
1063-
// daughter DCAs to PV in XY propagated with material
1064-
TrackPosDcaXY, TrackNegDcaXY, TrackBachDcaXY,
1065-
// daughter DCAs to PV in 3D propagated with material
1066-
TrackPosDca, TrackNegDca, TrackBachDca,
10671097
// daughter signs
10681098
kfpProton.GetQ(),
10691099
kfpPion.GetQ(),
@@ -1075,7 +1105,9 @@ struct decay3bodyBuilder {
10751105
dEdxProton,
10761106
dEdxPion,
10771107
dEdxDeuteron,
1078-
tofNSigmaDeuteron);
1108+
tofNSigmaDeuteron,
1109+
averageClusterSizeDeuteron,
1110+
trackBach.pidForTracking());
10791111
}
10801112
LOG(debug) << "Table filled.";
10811113

0 commit comments

Comments
 (0)