Skip to content

Add new particle: Hyper Helium4Sigma #13642

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
sawenzel merged 3 commits into from
Oct 31, 2024
Merged

Add new particle: Hyper Helium4Sigma #13642

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
b256f49
Add new particle: Hyper He4Sigma
wang-yuanzhe Oct 30, 2024
e7d2a5f
Please consider the following formatting changes
alibuild Oct 30, 2024
b37a220
Merge pull request #4 from alibuild/alibot-cleanup-13642
wang-yuanzhe Oct 30, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 4 additions & 1 deletion Common/Constants/include/CommonConstants/PhysicsConstants.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,7 +75,8 @@ enum Pdg {
kHyperTriton = 1010010030,
kHyperHydrogen4 = 1010010040,
kHyperHelium4 = 1010020040,
kHyperHelium5 = 1010020050
kHyperHelium5 = 1010020050,
kHyperHelium4Sigma = 1110020040
};

/// \brief Declarations of masses for additional particles
Expand Down Expand Up @@ -124,6 +125,7 @@ constexpr double MassHyperTriton = 2.99131;
constexpr double MassHyperHydrogen4 = 3.9226;
constexpr double MassHyperHelium4 = 3.9217;
constexpr double MassHyperHelium5 = 4.841;
constexpr double MassHyperHelium4Sigma = 3.995;

/// \brief Declarations of masses for particles in ROOT PDG_t
constexpr double MassDown = 0.00467;
Expand Down Expand Up @@ -194,6 +196,7 @@ constexpr double MassKaonNeutral = MassK0;
constexpr double MassLambda = MassLambda0;
constexpr double MassHyperhydrog4 = MassHyperHydrogen4;
constexpr double MassHyperhelium4 = MassHyperHelium4;
constexpr double MassHyperhelium4sigma = MassHyperHelium4Sigma;

// Light speed
constexpr float LightSpeedCm2S = 299792458.e2; // C in cm/s
Expand Down
1 change: 1 addition & 0 deletions Common/Constants/include/CommonConstants/make_pdg_header.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -131,6 +131,7 @@ class Pdg(Enum):
kHyperHydrogen4 = 1010010040
kHyperHelium4 = 1010020040
kHyperHelium5 = 1010020050
kHyperHelium4Sigma = 1110020040


dbPdg = ROOT.o2.O2DatabasePDG
Expand Down
13 changes: 13 additions & 0 deletions DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -331,6 +331,19 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db)
2.5e-15, 6, "Ion", ionCode);
}

// hyper helium 4 sigma
ionCode = 1110020040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("Hyperhelium4sigma", "Hyperhelium4sigma", 3.995, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}
// anti-hyper helium 4 sigma
ionCode = -1110020040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperhelium4sigma", "AntiHyperhelium4sigma", 3.995, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}

ionCode = 1010000020;
if (!db->GetParticle(ionCode)) {
db->AddParticle("LambdaNeutron", "LambdaNeutron", 2.054, kFALSE,
Expand Down
53 changes: 53 additions & 0 deletions Steer/src/O2MCApplication.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -303,6 +303,11 @@ void addSpecialParticles()
//Double Anti-Hyper hydrogen 4
TVirtualMC::GetMC()->DefineParticle(-1020010040, "DoubleAntiHyperhydrogen4", kPTHadron, 4.106, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

// Hyper helium 4 sigma
TVirtualMC::GetMC()->DefineParticle(1110020040, "Hyperhelium4sigma", kPTHadron, 3.995, 2.0, 8.018e-11, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
// Anti-Hyper helium 4 sigma
TVirtualMC::GetMC()->DefineParticle(-1110020040, "AntiHyperhelium4sigma", kPTHadron, 3.995, 2.0, 8.018e-11, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

//Lambda-Neutron
TVirtualMC::GetMC()->DefineParticle(1010000020, "LambdaNeutron", kPTNeutron, 2.054, 0.0, 2.632e-10, "Hadron", 0.0, 0, 1, 0, 0, 0, 0, 0, 2, kFALSE);

Expand Down Expand Up @@ -748,6 +753,54 @@ void addSpecialParticles()

TVirtualMC::GetMC()->SetDecayMode(-1020010040, abratio42, amode42);

// Define the 2- and 3-body phase space decay for the Hyper Helium 4 sigma
Int_t mode4s[6][3];
Float_t bratio4s[6];

for (Int_t kz = 0; kz < 6; kz++) {
bratio4s[kz] = 0.;
mode4s[kz][0] = 0;
mode4s[kz][1] = 0;
mode4s[kz][2] = 0;
}
bratio4s[0] = 20.;
mode4s[0][0] = 1000020040; // Helium4
mode4s[0][1] = 111; // pion0
bratio4s[1] = 40.;
mode4s[1][0] = 1000010030; // tritium
mode4s[1][2] = 2212; // proton
mode4s[1][1] = 111; // pion0
bratio4s[2] = 40.;
mode4s[2][0] = 1000010030; // tritium
mode4s[2][2] = 2212; // pion+
mode4s[2][1] = 2112; // neutron

TVirtualMC::GetMC()->SetDecayMode(1110020040, bratio4s, mode4s);

// Define the 2- and 3-body phase space decay for the Anti Hyper Helium 4 sigma
Int_t amode4s[6][3];
Float_t abratio4s[6];

for (Int_t kz = 0; kz < 6; kz++) {
abratio4s[kz] = 0.;
amode4s[kz][0] = 0;
amode4s[kz][1] = 0;
amode4s[kz][2] = 0;
}
abratio4s[0] = 50.;
amode4s[0][0] = -1000020040; // anti-Helium4
amode4s[0][1] = 111; // pion0
abratio4s[1] = 50.;
amode4s[1][0] = -1000010030; // anti-tritium
amode4s[1][2] = -2212; // anti-proton
amode4s[1][1] = 111; // pion0
abratio4s[2] = 50.;
amode4s[2][0] = -1000010030; // anti-tritium
amode4s[2][2] = -211; // pion-
amode4s[2][1] = -2112; // anti-neutron

TVirtualMC::GetMC()->SetDecayMode(-1110020040, abratio4s, amode4s);

// Define the 2-body phase space decay for the Lambda-neutron boundstate
Int_t mode1[6][3];
Float_t bratio1[6];
Expand Down