Skip to content

Updating hypernuclei information in O2 Databases #13750

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 4 commits into from
Nov 29, 2024
Merged

Updating hypernuclei information in O2 Databases #13750

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
f40931a
Updating hypernuclei
jditzelnew Nov 28, 2024
d015e9f
Updating hypernuclei
jditzelnew Nov 28, 2024
c14d1a5
Please consider the following formatting changes
alibuild Nov 28, 2024
2a33a08
Merge pull request #4 from alibuild/alibot-cleanup-13750
jditzelnew Nov 29, 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
54 changes: 22 additions & 32 deletions DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -235,62 +235,38 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db)
//Hyper nuclei and exotica
ionCode = 1010010030;
if (!db->GetParticle(ionCode)) {
db->AddParticle("HyperTriton", "HyperTriton", 2.99131, kFALSE,
db->AddParticle("HyperTriton", "HyperTriton", 2.991134, kFALSE,
2.5e-15, 3, "Ion", ionCode);
}

ionCode = -1010010030;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperTriton", "AntiHyperTriton", 2.99131, kFALSE,
db->AddParticle("AntiHyperTriton", "AntiHyperTriton", 2.991134, kFALSE,
2.5e-15, 3, "Ion", ionCode);
}

//hyper hydrogen 4 ground state
ionCode = 1010010040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("Hyperhydrog4", "Hyperhydrog4", 3.9226, kFALSE,
db->AddParticle("Hyperhydrog4", "Hyperhydrog4", 3.922434, kFALSE,
2.5e-15, 3, "Ion", ionCode);
}
//anti hyper hydrogen 4 ground state
ionCode = -1010010040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperhydrog4", "AntiHyperhydrog4", 3.9226, kFALSE,
2.5e-15, 3, "Ion", ionCode);
}
//hyper hydrogen 4 excited state
ionCode = 1010010041;
if (!db->GetParticle(ionCode)) {
db->AddParticle("Hyperhydrog4*", "Hyperhydrog4*", 3.9237, kFALSE,
2.5e-15, 3, "Ion", ionCode);
}
//anti hyper hydrogen 4 excited state
ionCode = -1010010041;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperhydrog4*", "AntiHyperhydrog4*", 3.9237, kFALSE,
db->AddParticle("AntiHyperhydrog4", "AntiHyperhydrog4", 3.922434, kFALSE,
2.5e-15, 3, "Ion", ionCode);
}
//hyper helium 4 ground state
ionCode = 1010020040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("Hyperhelium4", "Hyperhelium4", 3.9217, kFALSE,
db->AddParticle("Hyperhelium4", "Hyperhelium4", 3.921728, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}
//anti hyper helium 4 ground state
ionCode = -1010020040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperhelium4", "AntiHyperhelium4", 3.9217, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}
//hyper helium 4 excited state
ionCode = 1010020041;
if (!db->GetParticle(ionCode)) {
db->AddParticle("Hyperhelium4*", "Hyperhelium4*", 3.9231, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}
//anti hyper helium 4 excited state
ionCode = -1010020041;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperhelium4*", "AntiHyperhelium4*", 3.9231, kFALSE,
db->AddParticle("AntiHyperhelium4", "AntiHyperhelium4", 3.921728, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}

Expand All @@ -309,13 +285,13 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db)

ionCode = 1010020050;
if (!db->GetParticle(ionCode)) {
db->AddParticle("Hyperhelium5", "Hyperhelium5", 4.841, kFALSE,
db->AddParticle("Hyperhelium5", "Hyperhelium5", 4.839961, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}

ionCode = -1010020050;
if (!db->GetParticle(ionCode)) {
db->AddParticle("AntiHyperhelium5", "AntiHyperhelium5", 4.841, kFALSE,
db->AddParticle("AntiHyperhelium5", "AntiHyperhelium5", 4.839961, kFALSE,
2.5e-15, 6, "Ion", ionCode);
}

Expand All @@ -331,6 +307,20 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db)
2.5e-15, 6, "Ion", ionCode);
}

// 4-Xi-He
ionCode = 1120020040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("4XiHe", "4XiHe", 4.128, kFALSE, 4.04e-15, 3, "Ion", ionCode);
db->AddAntiParticle("Anti4XiHe", -ionCode);
}

// 4-Xi-H
ionCode = 1120010040;
if (!db->GetParticle(ionCode)) {
db->AddParticle("4XiH", "4XiH", 4.128, kFALSE, 4.04e-15, 3, "Ion", ionCode);
db->AddAntiParticle("Anti4XiH", -ionCode);
}

// hyper helium 4 sigma
ionCode = 1110020040;
if (!db->GetParticle(ionCode)) {
Expand Down
173 changes: 143 additions & 30 deletions Steer/src/O2MCApplication.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -264,45 +264,44 @@ void addSpecialParticles()
LOG(info) << "Adding custom particles to VMC";

//Hypertriton
TVirtualMC::GetMC()->DefineParticle(1010010030, "HyperTriton", kPTHadron, 2.99131, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 3, kFALSE);
TVirtualMC::GetMC()->DefineParticle(1010010030, "HyperTriton", kPTHadron, 2.991134, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 3, kFALSE);
//Anti-Hypertriton
TVirtualMC::GetMC()->DefineParticle(-1010010030, "AntiHyperTriton", kPTHadron, 2.99131, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 3, kFALSE);
TVirtualMC::GetMC()->DefineParticle(-1010010030, "AntiHyperTriton", kPTHadron, 2.991134, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 3, kFALSE);

//Hyper hydrogen 4 ground state
TVirtualMC::GetMC()->DefineParticle(1010010040, "Hyperhydrog4", kPTHadron, 3.9226, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
TVirtualMC::GetMC()->DefineParticle(1010010040, "Hyperhydrog4", kPTHadron, 3.922434, 1.0, 2.08e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
//Anti-Hyper hydrogen 4 ground state
TVirtualMC::GetMC()->DefineParticle(-1010010040, "AntiHyperhydrog4", kPTHadron, 3.9226, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

//Hyper hydrogen 4 excited state
TVirtualMC::GetMC()->DefineParticle(1010010041, "Hyperhydrog4*", kPTHadron, 3.9237, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
//Anti-Hyper hydrogen 4 excited state
TVirtualMC::GetMC()->DefineParticle(-1010010041, "AntiHyperhydrog4*", kPTHadron, 3.9237, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
TVirtualMC::GetMC()->DefineParticle(-1010010040, "AntiHyperhydrog4", kPTHadron, 3.922434, 1.0, 2.08e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

//Hyper helium 4 ground state
TVirtualMC::GetMC()->DefineParticle(1010020040, "Hyperhelium4", kPTHadron, 3.9217, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
TVirtualMC::GetMC()->DefineParticle(1010020040, "Hyperhelium4", kPTHadron, 3.921728, 2.0, 2.50e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
//Anti-Hyper helium 4 ground state
TVirtualMC::GetMC()->DefineParticle(-1010020040, "AntiHyperhelium4", kPTHadron, 3.9217, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

//Hyper helium 4 excited state
TVirtualMC::GetMC()->DefineParticle(1010020041, "Hyperhelium4*", kPTHadron, 3.9231, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
//Anti-Hyper helium 4 excited state
TVirtualMC::GetMC()->DefineParticle(-1010020041, "AntiHyperhelium4*", kPTHadron, 3.9231, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
TVirtualMC::GetMC()->DefineParticle(-1010020040, "AntiHyperhelium4", kPTHadron, 3.921728, 2.0, 2.50e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

// Lithium 4 ground state
TVirtualMC::GetMC()->DefineParticle(1000030040, "Lithium4", kPTHadron, 3.7513, 3.0, 9.1e-23, "Ion", 0.003, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
// Anti-Lithium 4 ground state
TVirtualMC::GetMC()->DefineParticle(-1000030040, "AntiLithium4", kPTHadron, 3.7513, 3.0, 9.1e-23, "Ion", 0.003, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);

//Hyper helium 5
TVirtualMC::GetMC()->DefineParticle(1010020050, "Hyperhelium5", kPTHadron, 4.841, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);
TVirtualMC::GetMC()->DefineParticle(1010020050, "Hyperhelium5", kPTHadron, 4.839961, 2.0, 2.74e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);
//Anti-Hyper helium 5
TVirtualMC::GetMC()->DefineParticle(-1010020050, "AntiHyperhelium5", kPTHadron, 4.841, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);
TVirtualMC::GetMC()->DefineParticle(-1010020050, "AntiHyperhelium5", kPTHadron, 4.839961, 2.0, 2.74e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);

//Double Hyper hydrogen 4
TVirtualMC::GetMC()->DefineParticle(1020010040, "DoubleHyperhydrogen4", kPTHadron, 4.106, 1.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
//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);

// 4Xi(-)H
TVirtualMC::GetMC()->DefineParticle(1120010040, "4XiH", kPTHadron, 4.128, 1.0, 1.639e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
// Anti-4Xi(-)H
TVirtualMC::GetMC()->DefineParticle(-1120010040, "Anti4XiH", kPTHadron, 4.128, 1.0, 1.639e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
// 4Xi(-)He
TVirtualMC::GetMC()->DefineParticle(1120020040, "4XiHe", kPTHadron, 4.128, 1.0, 1.639e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
// Anti-4Xi(-)He
TVirtualMC::GetMC()->DefineParticle(-1120020040, "Anti4XiHe", kPTHadron, 4.128, 1.0, 1.639e-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
Expand Down Expand Up @@ -586,8 +585,6 @@ void addSpecialParticles()
mode3[1][2] = -211; // negative pion

TVirtualMC::GetMC()->SetDecayMode(1010010040, bratio3, mode3);
//Decay for the excited state (after em transition)
TVirtualMC::GetMC()->SetDecayMode(1010010041, bratio3, mode3);

// Define the 2- and 3-body phase space decay for the Hyper Hydrogen 4
Int_t amode3[6][3];
Expand All @@ -608,8 +605,6 @@ void addSpecialParticles()
amode3[1][2] = 211; // positive pion

TVirtualMC::GetMC()->SetDecayMode(-1010010040, abratio3, amode3);
//Decay for the excited state (after em transition)
TVirtualMC::GetMC()->SetDecayMode(-1010010041, abratio3, amode3);

// Define the 3-body phase space decay for the Hyper Helium 4
Int_t mode4[6][3];
Expand All @@ -621,14 +616,16 @@ void addSpecialParticles()
mode4[kz][1] = 0;
mode4[kz][2] = 0;
}
bratio4[0] = 100.;
bratio4[0] = 50.;
mode4[0][0] = 1000020030; // Helium3
mode4[0][1] = -211; // negative pion
mode4[0][2] = 2212; // proton

bratio4[1] = 50.;
mode4[1][0] = 1000030040; // lithium-4
mode4[1][1] = -211; // negative pion

TVirtualMC::GetMC()->SetDecayMode(1010020040, bratio4, mode4);
//Decay for the excited state (after em transition)
TVirtualMC::GetMC()->SetDecayMode(1010020041, bratio4, mode4);

// Define the 2-body phase space decay for the Anti-Hyper Helium 4
Int_t amode4[6][3];
Expand All @@ -640,14 +637,16 @@ void addSpecialParticles()
amode4[kz][1] = 0;
amode4[kz][2] = 0;
}
abratio4[0] = 100.;
abratio4[0] = 50.;
amode4[0][0] = -1000020030; // anti-Helium 3
amode4[0][1] = 211; // positive pion
amode4[0][2] = -2212; // anti proton

abratio4[1] = 50.;
amode4[1][0] = -1000030040; // antilithium-4
amode4[1][1] = 211; // positive pion

TVirtualMC::GetMC()->SetDecayMode(-1010020040, abratio4, amode4);
//Decay for the excited state (after em transition)
TVirtualMC::GetMC()->SetDecayMode(-1010020041, abratio4, amode4);

// Define the 2-body phase space decay for the Lithium 4
Int_t model4[6][3];
Expand Down Expand Up @@ -733,10 +732,15 @@ void addSpecialParticles()
mode42[kz][1] = 0;
mode42[kz][2] = 0;
}
bratio42[0] = 100.;
bratio42[0] = 50.;
mode42[0][0] = 1010020040; // Hyper-Helium4
mode42[0][1] = -211; // negative pion

bratio42[1] = 50.;
mode42[1][0] = 1010010030; // Hypertriton
mode42[1][1] = 2212; // proton
mode42[1][2] = -211; // negative pion

TVirtualMC::GetMC()->SetDecayMode(1020010040, bratio42, mode42);

// Define the 2-body phase space decay for the Anti Double Hyper Hydrogen 4
Expand All @@ -749,12 +753,121 @@ void addSpecialParticles()
amode42[kz][1] = 0;
amode42[kz][2] = 0;
}
abratio42[0] = 100.;
abratio42[0] = 50.;
amode42[0][0] = -1010020040; // anti-Hyper-Helium 4
amode42[0][1] = 211; // positive pion

abratio42[1] = 50.;
amode42[1][0] = -1010010030; // anti-Hypertriton
amode42[1][1] = -2212; // antiproton
amode42[1][2] = 211; // positive pion

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

// Define the decay for the 4Xi(-)He
Int_t mode4XiHe[6][3];
Float_t bratio4XiHe[6];

for (Int_t kz = 0; kz < 6; kz++) {
bratio4XiHe[kz] = 0.;
mode4XiHe[kz][0] = 0;
mode4XiHe[kz][1] = 0;
mode4XiHe[kz][2] = 0;
}
bratio4XiHe[0] = 33.;
mode4XiHe[0][0] = 1010020040; // HyperHelium4
mode4XiHe[0][1] = -211; // negative pion

bratio4XiHe[1] = 33.;
mode4XiHe[1][0] = 3122; // lambda
mode4XiHe[1][1] = 1000020030; // helium-3
mode4XiHe[1][2] = -211; // negative pion

bratio4XiHe[2] = 33.;
mode4XiHe[2][0] = 1000030040; // lithium-4
mode4XiHe[2][1] = -211; // negative pion
mode4XiHe[2][2] = -211; // negative pion

TVirtualMC::GetMC()->SetDecayMode(1120020040, bratio4XiHe, mode4XiHe);

// Define the decay for the Anti-4Xi(-)He
Int_t amode4XiHe[6][3];
Float_t abratio4XiHe[6];

for (Int_t kz = 0; kz < 6; kz++) {
abratio4XiHe[kz] = 0.;
amode4XiHe[kz][0] = 0;
amode4XiHe[kz][1] = 0;
amode4XiHe[kz][2] = 0;
}
abratio4XiHe[0] = 33.;
amode4XiHe[0][0] = -1010020040; // antiHyperHelium-4
amode4XiHe[0][1] = 211; // positive pion

abratio4XiHe[1] = 33.;
amode4XiHe[1][0] = -3122; // antilambda
amode4XiHe[1][1] = -1000020030; // antihelium-3
amode4XiHe[1][2] = 211; // positive pion

abratio4XiHe[2] = 33.;
amode4XiHe[2][0] = -1000030040; // antilithium-4
amode4XiHe[2][1] = 211; // positive pion
amode4XiHe[2][2] = 211; // positive pion

TVirtualMC::GetMC()->SetDecayMode(-1120020040, abratio4XiHe, amode4XiHe);

// Define the decay for the 4Xi(-)H
Int_t mode4XiH[6][3];
Float_t bratio4XiH[6];

for (Int_t kz = 0; kz < 6; kz++) {
bratio4XiH[kz] = 0.;
mode4XiH[kz][0] = 0;
mode4XiH[kz][1] = 0;
mode4XiH[kz][2] = 0;
}
bratio4XiH[0] = 33.;
mode4XiH[0][0] = 1010010040; // HyperHydrogen4
mode4XiH[0][1] = -211; // negative pion

bratio4XiH[1] = 33.;
mode4XiH[1][0] = 3122; // lambda
mode4XiH[1][1] = 1000010030; // triton
mode4XiH[1][2] = -211; // negative pion

bratio4XiH[2] = 33.;
mode4XiH[2][0] = 1000020040; // alpha
mode4XiH[2][1] = -211; // negative pion
mode4XiH[2][2] = -211; // negative pion

TVirtualMC::GetMC()->SetDecayMode(1120010040, bratio4XiH, mode4XiH);

// Define the decay for the Anti-4Xi(-)H
Int_t amode4XiH[6][3];
Float_t abratio4XiH[6];

for (Int_t kz = 0; kz < 6; kz++) {
abratio4XiH[kz] = 0.;
amode4XiH[kz][0] = 0;
amode4XiH[kz][1] = 0;
amode4XiH[kz][2] = 0;
}
abratio4XiH[0] = 33.;
amode4XiH[0][0] = -1010010040; // antiHyperHydrogen-4
amode4XiH[0][1] = 211; // positive pion

abratio4XiH[1] = 33.;
amode4XiH[1][0] = -3122; // antilambda
amode4XiH[1][1] = -1000010030; // antitriton
amode4XiH[1][2] = 211; // positive pion

abratio4XiH[2] = 33.;
amode4XiH[2][0] = -1000020040; // antialpha
amode4XiH[2][1] = 211; // positive pion
amode4XiH[2][2] = 211; // positive pion

TVirtualMC::GetMC()->SetDecayMode(-1120010040, abratio4XiH, amode4XiH);

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