diff --git a/Common/Constants/include/CommonConstants/PhysicsConstants.h b/Common/Constants/include/CommonConstants/PhysicsConstants.h index 5f169f799eb5b..71f978aa3d1f6 100644 --- a/Common/Constants/include/CommonConstants/PhysicsConstants.h +++ b/Common/Constants/include/CommonConstants/PhysicsConstants.h @@ -51,7 +51,10 @@ enum Pdg { kDSBar = -431, kDSStar = 433, kDS1 = 10433, + kDS1Star2700 = 30433, + kDS1Star2860 = 40433, kDS2Star = 435, + kDS3Star2860 = 437, kDStar = 413, kDStar0 = 423, kChiC1 = 20443, @@ -72,6 +75,10 @@ enum Pdg { kXiCCPlusPlus = 4422, kXiCPlus = 4232, kXiC0 = 4132, + kXiC3055Plus = 4325, + kXiC3080Plus = 4326, + kXiC3055_0 = 4315, + kXiC3080_0 = 4316, kDeuteron = 1000010020, kTriton = 1000010030, kHelium3 = 1000020030, @@ -104,7 +111,10 @@ constexpr double MassDS = 1.96835; constexpr double MassDSBar = 1.96835; constexpr double MassDSStar = 2.1122; constexpr double MassDS1 = 2.53511; +constexpr double MassDS1Star2700 = 2.714; +constexpr double MassDS1Star2860 = 2.859; constexpr double MassDS2Star = 2.5691; +constexpr double MassDS3Star2860 = 2.86; constexpr double MassDStar = 2.01026; constexpr double MassDStar0 = 2.00685; constexpr double MassChiC1 = 3.51067; @@ -125,6 +135,10 @@ constexpr double MassXiB0 = 5.7919; constexpr double MassXiCCPlusPlus = 3.62155; constexpr double MassXiCPlus = 2.46771; constexpr double MassXiC0 = 2.47044; +constexpr double MassXiC3055Plus = 3.0559; +constexpr double MassXiC3080Plus = 3.0772; +constexpr double MassXiC3055_0 = 3.059; +constexpr double MassXiC3080_0 = 3.0799; constexpr double MassDeuteron = 1.87561294257; constexpr double MassTriton = 2.80892113298; constexpr double MassHelium3 = 2.80839160743; diff --git a/Common/Constants/include/CommonConstants/make_pdg_header.py b/Common/Constants/include/CommonConstants/make_pdg_header.py index 4fdfd052ff613..9f8e7ca6ec6da 100755 --- a/Common/Constants/include/CommonConstants/make_pdg_header.py +++ b/Common/Constants/include/CommonConstants/make_pdg_header.py @@ -106,7 +106,10 @@ class Pdg(Enum): kDSBar = -431 kDSStar = 433 kDS1 = 10433 + kDS1Star2700 = 30433 + kDS1Star2860 = 40433 kDS2Star = 435 + kDS3Star2860 = 437 kDStar = 413 kDStar0 = 423 kChiC1 = 20443 @@ -127,6 +130,10 @@ class Pdg(Enum): kXiCCPlusPlus = 4422 kXiCPlus = 4232 kXiC0 = 4132 + kXiC3055Plus = 4325 + kXiC3080Plus = 4326 + kXiC3055_0 = 4315 + kXiC3080_0 = 4316 kDeuteron = 1000010020 kTriton = 1000010030 kHelium3 = 1000020030 diff --git a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h index 0ee5c06a74411..ec8e6db889167 100644 --- a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h +++ b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h @@ -640,10 +640,47 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db) db->AddParticle("Anti-Theta_c_3100", "Anti-Theta_c_3100", 3.099, kFALSE, 83.e-6, 0, "Resonance", ionCode); // same width as D*+ (83 keV) } - if (!db->GetParticle(-ionCode)) { - db->AddParticle("Theta_c_3100", "Theta_c_3100", 3.099, kFALSE, - 83.e-6, 0, "Resonance", -ionCode); // same width as D*+ (83 keV) + db->AddAntiParticle("Theta_c_3100", -ionCode); + + // Charm resonances not present in PYTHIA (consistent with https://github.com/AliceO2Group/O2DPG/blob/master/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_DResoTrigger.cfg) + // Mesons + ionCode = 30433; + if (!db->GetParticle(ionCode)) { + db->AddParticle("D*_s1_Plus_2700", "D*_s1_Plus_2700", 2.714, false, 0.122, 3, "Resonance", ionCode); + } + db->AddAntiParticle("D*_s1_Minus_2700", -ionCode); + ionCode = 40433; + if (!db->GetParticle(ionCode)) { + db->AddParticle("D*_s1_Plus_2860", "D*_s1_Plus_2860", 2.859, false, 0.160, 3, "Resonance", ionCode); + } + db->AddAntiParticle("D*_s1_Minus_2860", -ionCode); + ionCode = 437; + if (!db->GetParticle(ionCode)) { + db->AddParticle("D*_s3_Plus_2860", "D*_s3_Plus_2860", 2.860, false, 0.053, 3, "Resonance", ionCode); + } + db->AddAntiParticle("D*_s3_Minus_2860", -ionCode); + + // Baryons + ionCode = 4325; + if (!db->GetParticle(ionCode)) { + db->AddParticle("Xi_c_Plus_3055", "Xi_c_Plus_3055", 3.0559, false, 0.0078, 3, "Resonance", ionCode); + } + db->AddAntiParticle("Xi_c_Minus_3055", -ionCode); + ionCode = 4326; + if (!db->GetParticle(ionCode)) { + db->AddParticle("Xi_c_Plus_3080", "Xi_c_Plus_3080", 3.0772, false, 0.0036, 3, "Resonance", ionCode); + } + db->AddAntiParticle("Xi_c_Minus_3080", -ionCode); + ionCode = 4315; + if (!db->GetParticle(ionCode)) { + db->AddParticle("Xi_c_0_3055", "Xi_c_0_3055", 3.0590, false, 0.0064, 0, "Resonance", ionCode); + } + db->AddAntiParticle("Anti-Xi_c_0_3055", -ionCode); + ionCode = 4316; + if (!db->GetParticle(ionCode)) { + db->AddParticle("Xi_c_0_3080", "Xi_c_0_3080", 3.0799, false, 0.0056, 0, "Resonance", ionCode); } + db->AddAntiParticle("Anti-Xi_c_0_3080", -ionCode); // d*(2380) - dibaryon resonance