Update and modify Lambda1520 case

alibuild · Mar 19, 2024 · f6aa24b · f6aa24b
1 parent b3f52be
commit f6aa24b
Show file tree

Hide file tree

Showing 2 changed files with 62 additions and 23 deletions.
diff --git a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h
@@ -469,6 +469,15 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db)
     db->AddParticle("f2_1270", "f2_1270", 1.275, kFALSE,
                     0.185, 0, "Resonance", ionCode);
   }
+
+  // Lambda(1520)0
+  ionCode = 102134;
+  if (!db->GetParticle(ionCode)) {
+    db->AddParticle("Lambda_1520_0", "Lambda_1520_0", 1.5195, kFALSE, 0.0156, 0, "Resonance", ionCode);
+  }
+  if (!db->GetParticle(-ionCode)) {
+    db->AddParticle("AntiLambda_1520_0", "AntiLambda_1520_0", 1.5195, kFALSE, 0.0156, 0, "Resonance", -ionCode);
+  }
 
   // Xi-/+ (1820)
   ionCode = 123314;

diff --git a/Steer/src/O2MCApplication.cxx b/Steer/src/O2MCApplication.cxx
@@ -1123,7 +1123,7 @@ void addSpecialParticles()
   //  Int_t mode[6][3];
   //  Float_t bratio[6];
 
-  for (Int_t kz = 0; kz < 6; kz++) {
+  for (Int_t kz = 0; kz < 9; kz++) {
     bratio[kz] = 0.;
     mode[kz][0] = 0;
     mode[kz][1] = 0;
@@ -1137,76 +1137,106 @@ void addSpecialParticles()
 
   // Lambda1520/Lambda1520bar
 
-  TVirtualMC::GetMC()->DefineParticle(3124, "Lambda1520", kPTNeutron, 1.5195, 0.0, 4.22e-23, "Hadron", 0.0156, 3, -1, 0, 0, 0, 0, 0, 1, kTRUE);
-  TVirtualMC::GetMC()->DefineParticle(-3124, "Lambda1520bar", kPTNeutron, 1.5195, 0.0, 4.22e-23, "Hadron", 0.0156, 3, -1, 0, 0, 0, 0, 0, -1, kTRUE);
+  TVirtualMC::GetMC()->DefineParticle(102134, "Lambda1520", kPTNeutron, 1.5195, 0.0, 4.22e-23, "Hadron", 0.0156, 3, -1, 0, 0, 0, 0, 0, 1, kTRUE);
+  TVirtualMC::GetMC()->DefineParticle(-102134, "Lambda1520bar", kPTNeutron, 1.5195, 0.0, 4.22e-23, "Hadron", 0.0156, 3, -1, 0, 0, 0, 0, 0, -1, kTRUE);
 
   // Lambda1520 decay modes
 
   // L(1520) -> p K-
-  bratio[0] = 0.223547;
+  bratio[0] = 0.229944;
   mode[0][0] = 2212;
   mode[0][1] = -321;
 
   // L(1520) -> n K0
-  bratio[1] = 0.223547;
+  bratio[1] = 0.229944;
   mode[1][0] = 2112;
   mode[1][1] = -311;
 
   // L(1520) -> Sigma+ pi-
-  bratio[2] = 0.139096;
+  bratio[2] = 0.143076;
   mode[2][0] = 3222;
   mode[2][1] = -211;
 
   // L(1520) -> Sigma0 pi0
-  bratio[3] = 0.139096;
+  bratio[3] = 0.143076;
   mode[3][0] = 3212;
   mode[3][1] = 111;
 
   // L(1520) -> Sigma- pi+
-  bratio[4] = 0.139096;
+  bratio[4] = 0.143076;
   mode[4][0] = 3112;
   mode[4][1] = 211;
 
-  // The other decay modes are neglected
-  bratio[5] = 0.;
-  mode[5][0] = 0;
-  mode[5][1] = 0;
+  // L(1520) -> Sigma*- pi+
+  bratio[5] = 0.034066;
+  mode[5][0] = 3114;
+  mode[5][1] = 211;
+
+  // L(1520) -> Sigma*0 pi0
+  bratio[6] = 0.034066;
+  mode[6][0] = 3214;
+  mode[6][1] = 111;
+
+  // L(1520) -> Sigma*+ pi-
+  bratio[7] = 0.034066;
+  mode[7][0] = 3224;
+  mode[7][1] = -211;
+
+  // L(1520) -> Lambda gamma
+  bratio[8] = 0.008687;
+  mode[8][0] = 3122;
+  mode[8][1] = 22;
 
-  TVirtualMC::GetMC()->SetDecayMode(3124, bratio, mode);
+  TVirtualMC::GetMC()->SetDecayMode(102134, bratio, mode);
 
   // Lambda1520bar decay modes
 
   // L(1520)bar -> p- K+
-  bratio[0] = 0.223547;
+  bratio[0] = 0.229944;
   mode[0][0] = -2212;
   mode[0][1] = 321;
 
   // L(1520)bar -> nbar K0bar
-  bratio[1] = 0.223547;
+  bratio[1] = 0.229944;
   mode[1][0] = -2112;
   mode[1][1] = 311;
 
   // L(1520)bar -> Sigmabar- pi+
-  bratio[2] = 0.139096;
+  bratio[2] = 0.143076;
   mode[2][0] = -3222;
   mode[2][1] = 211;
 
   // L(1520)bar -> Sigma0bar pi0
-  bratio[3] = 0.139096;
+  bratio[3] = 0.143076;
   mode[3][0] = -3212;
   mode[3][1] = 111;
 
   // L(1520)bar -> Sigmabar+ pi-
-  bratio[4] = 0.139096;
+  bratio[4] = 0.143076;
   mode[4][0] = -3112;
   mode[4][1] = -211;
 
-  // The other decay modes are neglected
-  bratio[5] = 0.;
-  mode[5][0] = 0;
-  mode[5][1] = 0;
+  // L(1520)bar -> anti-Sigma*- pi-
+  bratio[5] = 0.034066;
+  mode[5][0] = -3114;
+  mode[5][1] = -211;
+
+  // L(1520)bar -> anti-Sigma*0 pi0
+  bratio[6] = 0.034066;
+  mode[6][0] = -3214;
+  mode[6][1] = 111;
+
+  // L(1520)bar -> anti-Sigma*+ pi+
+  bratio[7] = 0.034066;
+  mode[7][0] = -3224;
+  mode[7][1] = 211;
+
+  // L(1520)bar -> Anti-Lambda gamma
+  bratio[8] = 0.008687;
+  mode[8][0] = -3122;
+  mode[8][1] = 22;
 
-  TVirtualMC::GetMC()->SetDecayMode(-3124, bratio, mode);
+  TVirtualMC::GetMC()->SetDecayMode(-102134, bratio, mode);
 
   // --------------------------------------------------------------------