From f45d41bc175159762411fb4b642f4e8d4a6757bd Mon Sep 17 00:00:00 2001 From: fmazzasc Date: Tue, 25 Jul 2023 17:41:01 +0200 Subject: [PATCH 1/2] Adding (anti-)lithium4 particle --- .../SimulationDataFormat/O2DatabasePDG.h | 13 ++++++ Steer/src/O2MCApplication.cxx | 41 +++++++++++++++++++ 2 files changed, 54 insertions(+) diff --git a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h index 1882e1e31903f..4c9c78b93ed11 100644 --- a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h +++ b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h @@ -296,6 +296,19 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db) 2.5e-15, 6, "Ion", ionCode); } + //Lithium 4 ground state + ionCode = 1000030040; + if (!db->GetParticle(ionCode)) { + db->AddParticle("Lithium4", "Lithium4", 3.74976, kFALSE, + 0.005, 9, "Ion", ionCode); + } + //anti Lithium 4 ground state + ionCode = -1000030040; + if (!db->GetParticle(ionCode)) { + db->AddParticle("AntiLithium4", "AntiLithium4", 3.74976, kFALSE, + 0.005, 9, "Ion", ionCode); + } + ionCode = 1010020050; if (!db->GetParticle(ionCode)) { db->AddParticle("Hyperhelium5", "Hyperhelium5", 4.841, kFALSE, diff --git a/Steer/src/O2MCApplication.cxx b/Steer/src/O2MCApplication.cxx index 64f79d4c46ac9..5795ccb4927bc 100644 --- a/Steer/src/O2MCApplication.cxx +++ b/Steer/src/O2MCApplication.cxx @@ -264,6 +264,11 @@ void addSpecialParticles() //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); + //Lithium 4 ground state + TVirtualMC::GetMC()->DefineParticle(1000030040, "Lithium4", kPTHadron, 3.74976, 3.0, 9.1e-23, "Ion", 0.005, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE); + //Anti-Lithium 4 ground state + TVirtualMC::GetMC()->DefineParticle(-1000030040, "AntiLithium4", kPTHadron, 3.74976, 3.0, 9.1e-23, "Ion", 0.005, 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); //Anti-Hyper helium 5 @@ -597,6 +602,42 @@ void addSpecialParticles() //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]; + Float_t bratiol4[6]; + + for (Int_t kz = 0; kz < 6; kz++) { + bratiol4[kz] = 0.; + model4[kz][0] = 0; + model4[kz][1] = 0; + model4[kz][2] = 0; + } + bratiol4[0] = 100.; + model4[0][0] = 1000020030; // Helium3 + model4[0][1] = 2212; // proton + + TVirtualMC::GetMC()->SetDecayMode(1000030040, bratiol4, model4); + + + // Define the 2-body phase space decay for the Anti-Lithium 4 + Int_t amodel4[6][3]; + Float_t abratiol4[6]; + + for (Int_t kz = 0; kz < 6; kz++) { + abratiol4[kz] = 0.; + amodel4[kz][0] = 0; + amodel4[kz][1] = 0; + amodel4[kz][2] = 0; + } + abratiol4[0] = 100.; + amodel4[0][0] = -1000020030; // Anti-Helium3 + amodel4[0][1] = -2212; // Anti-proton + + TVirtualMC::GetMC()->SetDecayMode(-1000030040, abratiol4, amodel4); + + // Define the 3-body phase space decay for the Hyper Helium 5 Int_t mode41[6][3]; Float_t bratio41[6]; From dd03ad085dcbde9620135ddb7e17145f042386dd Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Tue, 25 Jul 2023 15:43:02 +0000 Subject: [PATCH 2/2] Please consider the following formatting changes --- .../include/SimulationDataFormat/O2DatabasePDG.h | 4 ++-- Steer/src/O2MCApplication.cxx | 8 ++------ 2 files changed, 4 insertions(+), 8 deletions(-) diff --git a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h index 4c9c78b93ed11..3bd92b4e11664 100644 --- a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h +++ b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h @@ -296,13 +296,13 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db) 2.5e-15, 6, "Ion", ionCode); } - //Lithium 4 ground state + // Lithium 4 ground state ionCode = 1000030040; if (!db->GetParticle(ionCode)) { db->AddParticle("Lithium4", "Lithium4", 3.74976, kFALSE, 0.005, 9, "Ion", ionCode); } - //anti Lithium 4 ground state + // anti Lithium 4 ground state ionCode = -1000030040; if (!db->GetParticle(ionCode)) { db->AddParticle("AntiLithium4", "AntiLithium4", 3.74976, kFALSE, diff --git a/Steer/src/O2MCApplication.cxx b/Steer/src/O2MCApplication.cxx index 5795ccb4927bc..118374751ee69 100644 --- a/Steer/src/O2MCApplication.cxx +++ b/Steer/src/O2MCApplication.cxx @@ -264,9 +264,9 @@ void addSpecialParticles() //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); - //Lithium 4 ground state + // Lithium 4 ground state TVirtualMC::GetMC()->DefineParticle(1000030040, "Lithium4", kPTHadron, 3.74976, 3.0, 9.1e-23, "Ion", 0.005, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE); - //Anti-Lithium 4 ground state + // Anti-Lithium 4 ground state TVirtualMC::GetMC()->DefineParticle(-1000030040, "AntiLithium4", kPTHadron, 3.74976, 3.0, 9.1e-23, "Ion", 0.005, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE); //Hyper helium 5 @@ -602,8 +602,6 @@ void addSpecialParticles() //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]; Float_t bratiol4[6]; @@ -620,7 +618,6 @@ void addSpecialParticles() TVirtualMC::GetMC()->SetDecayMode(1000030040, bratiol4, model4); - // Define the 2-body phase space decay for the Anti-Lithium 4 Int_t amodel4[6][3]; Float_t abratiol4[6]; @@ -637,7 +634,6 @@ void addSpecialParticles() TVirtualMC::GetMC()->SetDecayMode(-1000030040, abratiol4, amodel4); - // Define the 3-body phase space decay for the Hyper Helium 5 Int_t mode41[6][3]; Float_t bratio41[6];