From 62283a99ee5e856747714ea5c4f7ffcbdbf1e348 Mon Sep 17 00:00:00 2001
From: hristov <Peter.Hristov@cern.ch>
Date: Fri, 19 Dec 2014 17:25:29 +0100
Subject: [PATCH] ALIROOT-5728 Lambda(1520) is not decaying anymore when
 injected

---
 STEER/STEER/AliMC.cxx | 72 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 72 insertions(+)

diff --git a/STEER/STEER/AliMC.cxx b/STEER/STEER/AliMC.cxx
index 74b4a22d8d3..9afb60e2371 100644
--- a/STEER/STEER/AliMC.cxx
+++ b/STEER/STEER/AliMC.cxx
@@ -565,6 +565,78 @@ void  AliMC::AddParticles()
   
   TVirtualMC::GetMC()->SetDecayMode(225,bratio,mode);
 
+  // Lambda1520/Lambda1520bar
+
+  TVirtualMC::GetMC()->DefineParticle(3124, "Lambda1520", kPTNeutron, 1.5195 , 0.0, 4.22e-23,"Hadron", 0.0156, 3, -1, 0, 0, 0, 0, 3, 0, kTRUE);
+  TVirtualMC::GetMC()->DefineParticle(-3124, "Lambda1520bar", kPTNeutron, 1.5195 , 0.0, 4.22e-23,"Hadron", 0.0156, 3, -1, 0, 0, 0, 0, 3, 0, kTRUE);
+
+  // Lambda1520 decay modes
+
+  // L(1520) -> p K-
+  bratio[0] = 0.223547;
+  mode[0][0] = 2212;
+  mode[0][1] = -321;
+
+  // L(1520) -> n K0
+  bratio[1] = 0.223547;
+  mode[1][0] = 2112;
+  mode[1][1] = -311;
+
+  // L(1520) -> Sigma+ pi-
+  bratio[2] = 0.139096;
+  mode[2][0] = 3222;
+  mode[2][1] = -211;
+
+  // L(1520) -> Sigma0 pi0
+  bratio[3] = 0.139096;
+  mode[3][0] = 3212;
+  mode[3][1] = 111;
+
+  // L(1520) -> Sigma- pi+
+  bratio[4] = 0.139096;
+  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;
+
+  TVirtualMC::GetMC()->SetDecayMode(3124,bratio,mode);
+
+  // Lambda1520bar decay modes
+
+  // L(1520)bar -> p- K+
+  bratio[0] = 0.223547;
+  mode[0][0] = -2212;
+  mode[0][1] = 321;
+
+  // L(1520)bar -> nbar K0bar
+  bratio[1] = 0.223547;
+  mode[1][0] = -2112;
+  mode[1][1] = 311;
+
+  // L(1520)bar -> Sigmabar- pi+
+  bratio[2] = 0.139096;
+  mode[2][0] = -3222;
+  mode[2][1] = 211;
+
+  // L(1520)bar -> Sigma0bar pi0
+  bratio[3] = 0.139096;
+  mode[3][0] = -3212;
+  mode[3][1] = 111;
+
+  // L(1520)bar -> Sigmabar+ pi-
+  bratio[4] = 0.139096;
+  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;
+
+  TVirtualMC::GetMC()->SetDecayMode(-3124,bratio,mode);
 
   // --------------------------------------------------------------------
 }  
-- 
2.43.0