Skip to content

Commit 059f628

Browse files
Oussama3030Oussama3030
committed
Added additional MC information to study decay chains
1 parent 4442132 commit 059f628

4 files changed

Lines changed: 101 additions & 24 deletions

File tree

‎PWGLF/DataModel/LFSigmaTables.h‎

Lines changed: 6 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -790,6 +790,8 @@ DECLARE_SOA_COLUMN(PhotonMCPz, photonmcpz, float);
790790
DECLARE_SOA_COLUMN(IsPhotonPrimary, isPhotonPrimary, bool);
791791
DECLARE_SOA_COLUMN(PhotonPDGCode, photonPDGCode, int);
792792
DECLARE_SOA_COLUMN(PhotonPDGCodeMother, photonPDGCodeMother, int);
793+
DECLARE_SOA_COLUMN(PhotonPDGCodeGrandMother, photonPDGCodeGrandMother, int);
794+
DECLARE_SOA_COLUMN(PhotonGlobalIndexGrandMother, photonGlobalIndexGrandMother, int);
793795
DECLARE_SOA_COLUMN(PhotonIsCorrectlyAssoc, photonIsCorrectlyAssoc, bool);
794796

795797
DECLARE_SOA_COLUMN(KShortMCPx, kshortmcpx, float);
@@ -798,6 +800,8 @@ DECLARE_SOA_COLUMN(KShortMCPz, kshortmcpz, float);
798800
DECLARE_SOA_COLUMN(IsKShortPrimary, isKShortPrimary, bool);
799801
DECLARE_SOA_COLUMN(KShortPDGCode, kshortPDGCode, int);
800802
DECLARE_SOA_COLUMN(KShortPDGCodeMother, kshortPDGCodeMother, int);
803+
DECLARE_SOA_COLUMN(KShortPDGCodeGrandMother, kshortPDGCodeGrandMother, int);
804+
DECLARE_SOA_COLUMN(KShortGlobalIndexGrandMother, kshortGlobalIndexGrandMother, int);
801805
DECLARE_SOA_COLUMN(KShortIsCorrectlyAssoc, kshortIsCorrectlyAssoc, bool);
802806

803807
DECLARE_SOA_DYNAMIC_COLUMN(IsKStar, isKStar, //! IsSigma0
@@ -903,10 +907,10 @@ DECLARE_SOA_TABLE(KStarMCCores, "AOD", "KSTARMCCORES",
903907
kstarMCCore::MCradius, kstarMCCore::PDGCode, kstarMCCore::PDGCodeMother, kstarMCCore::MCprocess, kstarMCCore::IsProducedByGenerator,
904908

905909
kstarMCCore::PhotonMCPx, kstarMCCore::PhotonMCPy, kstarMCCore::PhotonMCPz,
906-
kstarMCCore::IsPhotonPrimary, kstarMCCore::PhotonPDGCode, kstarMCCore::PhotonPDGCodeMother, kstarMCCore::PhotonIsCorrectlyAssoc,
910+
kstarMCCore::IsPhotonPrimary, kstarMCCore::PhotonPDGCode, kstarMCCore::PhotonPDGCodeMother, kstarMCCore::PhotonPDGCodeGrandMother, kstarMCCore::PhotonGlobalIndexGrandMother, kstarMCCore::PhotonIsCorrectlyAssoc,
907911

908912
kstarMCCore::KShortMCPx, kstarMCCore::KShortMCPy, kstarMCCore::KShortMCPz,
909-
kstarMCCore::IsKShortPrimary, kstarMCCore::KShortPDGCode, kstarMCCore::KShortPDGCodeMother, kstarMCCore::KShortIsCorrectlyAssoc,
913+
kstarMCCore::IsKShortPrimary, kstarMCCore::KShortPDGCode, kstarMCCore::KShortPDGCodeMother, kstarMCCore::KShortPDGCodeGrandMother, kstarMCCore::KShortGlobalIndexGrandMother, kstarMCCore::KShortIsCorrectlyAssoc,
910914

911915
// Dynamic columns
912916
kstarMCCore::IsKStar<kstarMCCore::PDGCode>,

‎PWGLF/TableProducer/Strangeness/sigma0builder.cxx‎

Lines changed: 48 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -811,6 +811,10 @@ struct sigma0builder {
811811
int V02PDGCode = 0;
812812
int V01PDGCodeMother = 0;
813813
int V02PDGCodeMother = 0;
814+
int V01PDGCodeGrandMother = 0;
815+
int V02PDGCodeGrandMother = 0;
816+
int V01GlobalIndexGrandMother = 0;
817+
int V02GlobalIndexGrandMother = 0;
814818
int V0PairPDGCode = 0;
815819
int V0PairPDGCodeMother = 0;
816820
int V0PairMCProcess = -1;
@@ -956,6 +960,48 @@ struct sigma0builder {
956960
auto const& MCMother_v01 = MCMothersList_v01.front(); // First mother
957961
auto const& MCMother_v02 = MCMothersList_v02.front(); // First mother
958962

963+
// Add the grandmothers
964+
auto const& GrandMothersList_v01 = MCMother_v01.template mothers_as<aod::McParticles>();
965+
if (!GrandMothersList_v01.empty()) {
966+
MCinfo.V01PDGCodeGrandMother = GrandMothersList_v01.front().pdgCode();
967+
MCinfo.V01GlobalIndexGrandMother = GrandMothersList_v01.front().globalIndex();
968+
}
969+
970+
auto const& GrandMothersList_v02 = MCMother_v02.template mothers_as<aod::McParticles>();
971+
if (!GrandMothersList_v02.empty()) {
972+
MCinfo.V02PDGCodeGrandMother = GrandMothersList_v02.front().pdgCode();
973+
MCinfo.V02GlobalIndexGrandMother = GrandMothersList_v02.front().globalIndex();
974+
}
975+
976+
// check grandmothers and fill histograms
977+
int kShortMotherCode = 0;
978+
int photonMotherCode = 0;
979+
if ((std::abs(MCParticle_v01.pdgCode()) == PDG_t::kGamma) && (std::abs(MCParticle_v02.pdgCode()) == PDG_t::kK0Short) && (!fIsKStar)) {
980+
981+
kShortMotherCode = MCMother_v02.pdgCode();
982+
983+
// If the KShort mother is a (anti)Kaon, use the grandmother instead
984+
if (std::abs(kShortMotherCode) == PDG_t::kK0 ||
985+
std::abs(kShortMotherCode) == PDG_t::kK0Long ||
986+
std::abs(kShortMotherCode) == PDG_t::kKPlus) {
987+
auto const& kShortGrandMothers = MCMother_v02.template mothers_as<aod::McParticles>();
988+
if (!kShortGrandMothers.empty()) {
989+
kShortMotherCode = kShortGrandMothers.front().pdgCode();
990+
}
991+
}
992+
993+
photonMotherCode = MCMother_v01.pdgCode();
994+
// If the photon mother is a pi0, climb to grandmother
995+
if (std::abs(photonMotherCode) == PDG_t::kPi0) {
996+
auto const& photonGrandMothers = MCMother_v01.template mothers_as<aod::McParticles>();
997+
if (!photonGrandMothers.empty()) {
998+
photonMotherCode = photonGrandMothers.front().pdgCode();
999+
}
1000+
}
1001+
1002+
histos.fill(HIST("MCQA/h2dTrueDaughtersMatrix"), kShortMotherCode, photonMotherCode);
1003+
}
1004+
9591005
if (MCMother_v01.globalIndex() == MCMother_v02.globalIndex()) { // Is it the same mother?
9601006

9611007
MCinfo.fV0PairProducedByGenerator = MCMother_v01.producedByGenerator();
@@ -2561,9 +2607,9 @@ struct sigma0builder {
25612607

25622608
kstarmccores(kstarMCInfo.V0PairMCRadius, kstarMCInfo.V0PairPDGCode, kstarMCInfo.V0PairPDGCodeMother, kstarMCInfo.V0PairMCProcess, kstarMCInfo.fV0PairProducedByGenerator,
25632609
kstarMCInfo.V01MCpx, kstarMCInfo.V01MCpy, kstarMCInfo.V01MCpz,
2564-
kstarMCInfo.fIsV01Primary, kstarMCInfo.V01PDGCode, kstarMCInfo.V01PDGCodeMother, kstarMCInfo.fIsV01CorrectlyAssign,
2610+
kstarMCInfo.fIsV01Primary, kstarMCInfo.V01PDGCode, kstarMCInfo.V01PDGCodeMother, kstarMCInfo.V01PDGCodeGrandMother, kstarMCInfo.V01GlobalIndexGrandMother, kstarMCInfo.fIsV01CorrectlyAssign,
25652611
kstarMCInfo.V02MCpx, kstarMCInfo.V02MCpy, kstarMCInfo.V02MCpz,
2566-
kstarMCInfo.fIsV02Primary, kstarMCInfo.V02PDGCode, kstarMCInfo.V02PDGCodeMother, kstarMCInfo.fIsV02CorrectlyAssign);
2612+
kstarMCInfo.fIsV02Primary, kstarMCInfo.V02PDGCode, kstarMCInfo.V02PDGCodeMother, kstarMCInfo.V02PDGCodeGrandMother, kstarMCInfo.V02GlobalIndexGrandMother, kstarMCInfo.fIsV02CorrectlyAssign);
25672613
}
25682614

25692615
// KStar -> stracollisions link

‎PWGLF/Tasks/Resonances/k892hadronphoton.cxx‎

Lines changed: 44 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -383,11 +383,9 @@ struct k892hadronphoton {
383383
histos.add(histodir + "/MC/KStar/hDCAPairDauVsPt", "hDCAPairDauVsPt", kTH2D, {axisDCAdau, axisPt});
384384

385385
// 1/pT Resolution:
386-
if (fillResoQAhistos && histodir == "BeforeSel") {
387-
386+
if (fillResoQAhistos) {
388387
histos.add(histodir + "/MC/Reso/h2dKShortPtResolution", "h2dKShortPtResolution", kTH2D, {axisInvPt, axisDeltaPt});
389-
histos.add(histodir + "/MC/Reso/h3dKShortPtResoVsTPCCR", "h3dKShortPtResoVsTPCCR", kTH3D, {axisInvPt, axisDeltaPt, axisTPCrows});
390-
histos.add(histodir + "/MC/Reso/h3dKShortPtResoVsTPCCR", "h3dKShortPtResoVsTPCCR", kTH3D, {axisInvPt, axisDeltaPt, axisTPCrows});
388+
histos.add(histodir + "/MC/Reso/h2dGammaPtResolution", "h2dGammaPtResolution", kTH2D, {axisInvPt, axisDeltaPt});
391389
histos.add(histodir + "/MC/Reso/h2dKStarPtResolution", "h2dKStarPtResolution", kTH2D, {axisInvPt, axisDeltaPt});
392390
histos.add(histodir + "/MC/Reso/h2dKStarRadiusResolution", "h2dKStarRadiusResolution", kTH2D, {axisPt, axisDeltaPt});
393391
}
@@ -396,10 +394,16 @@ struct k892hadronphoton {
396394
if (fillBkgQAhistos) {
397395
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_All", "h2dPtVsMassKStar_All", kTH2D, {axisPt, axisKStarMass});
398396
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_TrueDaughters", "h2dPtVsMassKStar_TrueDaughters", kTH2D, {axisPt, axisKStarMass});
399-
histos.add(histodir + "/MC/BkgStudy/h2dTrueDaughtersMatrix", "h2dTrueDaughtersMatrix", kTHnSparseD, {{10001, -5000.5f, +5000.5f}, {10001, -5000.5f, +5000.5f}});
397+
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_PhotonOmega", "h2dPtVsMassKStar_PhotonOmega", kTH2D, {axisPt, axisKStarMass});
398+
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_PhotonRho", "h2dPtVsMassKStar_PhotonRho", kTH2D, {axisPt, axisKStarMass});
399+
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_PhotonEta", "h2dPtVsMassKStar_PhotonEta", kTH2D, {axisPt, axisKStarMass});
400+
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_KShortKCharged", "h2dPtVsMassKStar_KShortKCharged", kTH2D, {axisPt, axisKStarMass});
401+
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_KStarPionKaon", "h2dPtVsMassKStar_KStarPionKaon", kTH2D, {axisPt, axisKStarMass});
400402
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_TrueGammaFakeKShort", "h2dPtVsMassKStar_TrueGammaFakeKShort", kTH2D, {axisPt, axisKStarMass});
401403
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_FakeGammaTrueKShort", "h2dPtVsMassKStar_FakeGammaTrueKShort", kTH2D, {axisPt, axisKStarMass});
402404
histos.add(histodir + "/MC/BkgStudy/h2dPtVsMassKStar_FakeDaughters", "h2dPtVsMassKStar_FakeDaughters", kTH2D, {axisPt, axisKStarMass});
405+
histos.add(histodir + "/MC/BkgStudy/h2dTrueDaughtersMatrix", "h2dTrueDaughtersMatrix", kTHnSparseD, {{10001, -5000.5f, +5000.5f}, {10001, -5000.5f, +5000.5f}});
406+
histos.add(histodir + "/MC/BkgStudy/h2dTrueDaughtersMatrixGrandMother", "h2dTrueDaughtersMatrixGrandMother", kTHnSparseD, {{10001, -5000.5f, +5000.5f}, {10001, -5000.5f, +5000.5f}});
403407
}
404408
}
405409
}
@@ -761,36 +765,36 @@ struct k892hadronphoton {
761765
return trkCode;
762766
}
763767

764-
template <typename TKStarObject>
768+
template <int mode, typename TKStarObject>
765769
void getResolution(TKStarObject const& kstar)
766770
{
771+
// Check whether it is before or after selections
772+
static constexpr std::string_view MainDir[] = {"BeforeSel", "AfterSel"};
767773

768774
//_______________________________________
769775
// Gamma MC association
770776
if (std::abs(kstar.photonPDGCode()) == PDG_t::kGamma) {
771777
if (kstar.photonmcpt() > 0) {
772-
histos.fill(HIST("BeforeSel/MC/Reso/h3dGammaPtResoVsTPCCR"), 1.f / kstar.kshortmcpt(), 1.f / kstar.kshortPt() - 1.f / kstar.kshortmcpt(), -1 * kstar.photonNegTPCCrossedRows()); // 1/pT resolution
773-
histos.fill(HIST("BeforeSel/MC/Reso/h3dGammaPtResoVsTPCCR"), 1.f / kstar.kshortmcpt(), 1.f / kstar.kshortPt() - 1.f / kstar.kshortmcpt(), kstar.photonPosTPCCrossedRows()); // 1/pT resolution
774-
histos.fill(HIST("BeforeSel/MC/Reso/h2dGammaPtResolution"), 1.f / kstar.photonmcpt(), 1.f / kstar.photonPt() - 1.f / kstar.photonmcpt()); // pT resolution
778+
// histos.fill(HIST(MainDir[mode]) + HIST("/MC/Reso/h2dGammaPtResolution"), 1.f / kstar.photonmcpt(), kstar.photonPt() - kstar.photonmcpt());
779+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/Reso/h2dGammaPtResolution"), 1.f / kstar.photonmcpt(), (kstar.photonPt() - kstar.photonmcpt()) / kstar.photonmcpt());
775780
}
776781
}
777782

778783
//_______________________________________
779784
// KShort MC association
780785
if (std::abs(kstar.kshortPDGCode()) == PDG_t::kK0Short) {
781786
if (kstar.kshortmcpt() > 0) {
782-
histos.fill(HIST("BeforeSel/MC/Reso/h2dKShortPtResolution"), 1.f / kstar.kshortmcpt(), 1.f / kstar.kshortPt() - 1.f / kstar.kshortmcpt()); // 1/pT resolution
783-
histos.fill(HIST("BeforeSel/MC/Reso/h3dKShortPtResoVsTPCCR"), 1.f / kstar.kshortmcpt(), 1.f / kstar.kshortPt() - 1.f / kstar.kshortmcpt(), -1 * kstar.kshortNegTPCCrossedRows()); // 1/pT resolution
784-
histos.fill(HIST("BeforeSel/MC/Reso/h3dKShortPtResoVsTPCCR"), 1.f / kstar.kshortmcpt(), 1.f / kstar.kshortPt() - 1.f / kstar.kshortmcpt(), kstar.kshortPosTPCCrossedRows()); // 1/pT resolution
787+
// histos.fill(HIST(MainDir[mode]) + HIST("/MC/Reso/h2dKShortPtResolution"), 1.f / kstar.kshortmcpt(), kstar.kshortPt() - kstar.kshortmcpt());
788+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/Reso/h2dKShortPtResolution"), 1.f / kstar.kshortmcpt(), (kstar.kshortPt() - kstar.kshortmcpt()) / kstar.kshortmcpt());
785789
}
786790
}
787791

788792
//_______________________________________
789793
// KStar MC association
790794
if (kstar.isKStar()) {
791-
histos.fill(HIST("BeforeSel/MC/Reso/h2dKStarRadiusResolution"), kstar.mcpt(), kstar.radius() - kstar.mcradius()); // pT resolution
795+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/Reso/h2dKStarRadiusResolution"), kstar.mcpt(), kstar.radius() - kstar.mcradius()); // pT resolution
792796
if (kstar.mcpt() > 0)
793-
histos.fill(HIST("BeforeSel/MC/Reso/h2dKStarPtResolution"), 1.f / kstar.mcpt(), 1.f / kstar.pt() - 1.f / kstar.mcpt()); // pT resolution
797+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/Reso/h2dKStarPtResolution"), 1.f / kstar.mcpt(), (1.f / kstar.pt() - 1.f / kstar.mcpt()) / (1.f / kstar.mcpt())); // pT resolution
794798
}
795799
}
796800

@@ -804,18 +808,40 @@ struct k892hadronphoton {
804808
bool fIsKStar = kstar.isKStar();
805809
int photonPDGCode = kstar.photonPDGCode();
806810
int photonPDGCodeMother = kstar.photonPDGCodeMother();
811+
int photonPDGCodeGrandMother = kstar.photonPDGCodeGrandMother();
812+
int photonGlobalIndexGrandMother = kstar.photonGlobalIndexGrandMother();
807813
int kshortPDGCode = kstar.kshortPDGCode();
808814
int kshortPDGCodeMother = kstar.kshortPDGCodeMother();
815+
int kshortPDGCodeGrandMother = kstar.kshortPDGCodeGrandMother();
816+
int kshortGlobalIndexGrandMother = kstar.kshortGlobalIndexGrandMother();
809817
float kstarpT = kstar.pt();
810818
float kstarMass = kstar.kstarMass();
811819

812820
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_All"), kstarpT, kstarMass);
813821

814822
//_______________________________________
815-
// Real Gamma x Real KShort - but not from the same kstar!
816-
if ((!fIsKStar)) { //(std::abs(photonPDGCode) == PDG_t::kGamma) && (std::abs(KShortPDGCode) == PDG_t::kK0Short) &&
823+
// Real Gamma x Real KShort - but not direct decays from the same kstar!
824+
if ((!fIsKStar) && (std::abs(photonPDGCode) == PDG_t::kGamma) && (std::abs(kshortPDGCode) == PDG_t::kK0Short)) {
817825
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_TrueDaughters"), kstarpT, kstarMass);
818826
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dTrueDaughtersMatrix"), kshortPDGCodeMother, photonPDGCodeMother);
827+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dTrueDaughtersMatrixGrandMother"), kshortPDGCodeGrandMother, photonPDGCodeGrandMother);
828+
829+
// coming from the same kstar (coupled), but via a different decay channel (K* -> K0 pi0 -> K0s gamma)
830+
if (std::abs(kshortPDGCodeGrandMother) == o2::constants::physics::Pdg::kK0Star892 && std::abs(photonPDGCodeGrandMother) == o2::constants::physics::Pdg::kK0Star892 && (photonGlobalIndexGrandMother == kshortGlobalIndexGrandMother)) // K*(892)0
831+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_KStarPionKaon"), kstarpT, kstarMass);
832+
833+
// Break down of different photon / kshort sources
834+
if (std::abs(photonPDGCodeGrandMother) == o2::constants::physics::Pdg::kOmega) // omega(782)
835+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_PhotonOmega"), kstarpT, kstarMass);
836+
837+
if (std::abs(photonPDGCodeGrandMother) == PDG_t::kRho770Plus) // rho+-(770)
838+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_PhotonRho"), kstarpT, kstarMass);
839+
840+
if (std::abs(photonPDGCodeMother) == o2::constants::physics::Pdg::kEta) // eta
841+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_PhotonEta"), kstarpT, kstarMass);
842+
843+
if (std::abs(kshortPDGCodeGrandMother) == o2::constants::physics::Pdg::kKPlusStar892) // K*(892)+-
844+
histos.fill(HIST(MainDir[mode]) + HIST("/MC/BkgStudy/h2dPtVsMassKStar_KShortKCharged"), kstarpT, kstarMass);
819845
}
820846

821847
//_______________________________________
@@ -978,8 +1004,8 @@ struct k892hadronphoton {
9781004

9791005
//_______________________________________
9801006
// pT resolution histos
981-
if ((mode == 0) && fillResoQAhistos)
982-
getResolution(kstar);
1007+
if (fillResoQAhistos)
1008+
getResolution<mode>(kstar);
9831009
}
9841010
}
9851011
}

‎PWGLF/Tasks/Resonances/k892hadronphotonBkg.cxx‎

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -90,6 +90,7 @@ struct k892hadronphotonBkg {
9090
Configurable<float> kstarMaxRap{"kstarMaxRap", 0.5f, "Max |y(K*)|"};
9191
Configurable<int> nBkgRot{"nBkgRot", 3, "Rotations per pair (rotational bkg)"};
9292
Configurable<int> rotationalCut{"rotationalCut", 10, "theta band: [pi - pi/cut, pi + pi/cut]"};
93+
Configurable<float> rotationalFactor{"rotationalFactor", 1.f, "Factor to scale the angle of rotation (rotationalFactor * PI)"};
9394
} kstarBkgConfig;
9495

9596
ConfigurableAxis axisVertexMixBkg{"axisVertexMixBkg", {VARIABLE_WIDTH, -10.f, -8.f, -6.f, -4.f, -2.f, 0.f, 2.f, 4.f, 6.f, 8.f, 10.f}, "z-vertex bins for mixing"};
@@ -595,8 +596,8 @@ struct k892hadronphotonBkg {
595596
0.0);
596597

597598
for (int irot = 0; irot < kstarBkgConfig.nBkgRot; ++irot) {
598-
float theta = rotRng.Uniform(o2::constants::math::PI - o2::constants::math::PI / kstarBkgConfig.rotationalCut,
599-
o2::constants::math::PI + o2::constants::math::PI / kstarBkgConfig.rotationalCut);
599+
float theta = rotRng.Uniform(kstarBkgConfig.rotationalFactor * o2::constants::math::PI - o2::constants::math::PI / kstarBkgConfig.rotationalCut,
600+
kstarBkgConfig.rotationalFactor * o2::constants::math::PI + o2::constants::math::PI / kstarBkgConfig.rotationalCut);
600601

601602
ROOT::Math::PtEtaPhiMVector kRot(kshort.pt(), kshort.eta(), kshort.phi() + theta, o2::constants::physics::MassK0Short);
602603

0 commit comments

Comments
 (0)