void UserDecayConfig() {
   cout << "Loading User Decay Config from macro"<< endl;  
   TDatabasePDG *db= TDatabasePDG::Instance();
   TParticlePDG *p=0;

   
     Int_t mode[6][3];                  
     Float_t bratio[6];
     Int_t AlphaPDG, He5PDG;
     // p= db->GetParticle("Alpha");
     // if(p) AlphaPDG=p->PdgCode();
     p=db->GetParticle("He5");
     if(p) He5PDG=p->PdgCode();
     cout<<" he5 "<<He5PDG<<" "<<AlphaPDG<<endl;
     
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio[kz] = 0.;
     mode[kz][0] = 0;
     mode[kz][1] = 0;
     mode[kz][2] = 0;
     //  cout << mode[kz][0] << " " << 	mode[kz][1] << " " << mode[kz][2] << endl;
     }
     bratio[0] = 100.;
     mode[0][0] =2212  ;
     mode[0][1] =1000020040  ;
     mode[0][2] =-211;
     
     
     
     gMC->SetDecayMode(He5PDG,bratio,mode);
     
     Int_t mode1[6][3];                  
     Float_t bratio1[6];
     Int_t He3PDG, H3LPDG;
     //p= db->GetParticle("He3");
     //if(p) He3PDG=p->PdgCode();
     p=db->GetParticle("H3L");
     if(p) H3LPDG=p->PdgCode();
     cout<<" he3 "<<He3PDG<<" "<<H3LPDG<<endl;
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio1[kz] = 0.;
     mode1[kz][0] = 0;
     mode1[kz][1] = 0; mode1[kz][2] = 0;
     }
     bratio1[0] = 100.;
     mode1[0][0] =1000020030;
     mode1[0][1] =-211  ;
     
     
     gMC->SetDecayMode(H3LPDG,bratio1,mode1);
     
     Int_t mode2[6][3];                  
     Float_t bratio2[6];
     Int_t H4LPDG;
     
     p=db->GetParticle("H4L");
     if(p) H4LPDG=p->PdgCode();
     
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio2[kz] = 0.;
     mode2[kz][0] = 0;
     mode2[kz][1] = 0; mode2[kz][2] = 0;
     }
     bratio2[0] = 100.;
     mode2[0][0] =1000020040;
     mode2[0][1] =-211;
     gMC->SetDecayMode(H4LPDG,bratio2,mode2);
     
     cout<<" h4L "<<H4LPDG<<" "<<endl;
     
     Int_t mode3[6][3];                  
     Float_t bratio3[6];
     Int_t He6LPDG, Li6PDG;
     //p= db->GetParticle("Ion_6_3");
     
     //if(p) Li6PDG=p->PdgCode();
     p=db->GetParticle("He6L");
     if(p) He6LPDG=p->PdgCode();
     
     cout<<" He6 "<<Li6PDG<<" "<<He6LPDG<<endl;
     
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio3[kz] = 0.;
     mode3[kz][0] = 0;
     mode3[kz][1] = 0; mode3[kz][2] = 0;
     }
     bratio3[0] = 100.;
     mode3[0][0] =1000030060;
     mode3[0][1] =-211  ;
     
     
     gMC->SetDecayMode(He6LPDG,bratio3,mode3);
     
     Int_t mode4[6][3];                  
     Float_t bratio4[6];
     Int_t He8LPDG, Li8PDG;
     //p= db->GetParticle("Ion_8_3");
     
     //if(p) Li8PDG=p->PdgCode();
     p=db->GetParticle("He8L");
     if(p) He8LPDG=p->PdgCode();
     cout<<" he8 "<<Li8PDG<<" "<<He8LPDG<<endl;  
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio4[kz] = 0.;
     mode4[kz][0] = 0;
     mode4[kz][1] = 0; mode4[kz][2] = 0;
     }
     bratio4[0] = 100.;
     mode4[0][0] =1000030080;//Li8PDG;
     mode4[0][1] =-211  ;
     
     
     gMC->SetDecayMode(He8LPDG,bratio4,mode4);
     
     Int_t mode5[6][3];                  
     Float_t bratio5[6];
     Int_t Be6PDG;
     p= db->GetParticle("Ion_6_4");
     if(p)Be6PDG=p->PdgCode();
     
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio5[kz] = 0.;
     mode5[kz][0] = 0;
     mode5[kz][1] = 0;mode5[kz][2] = 0;
     }
     bratio5[0] = 100.;
     mode5[0][0] =1000020040;
     mode5[0][1] =2212  ;
     mode5[0][2] =2212  ;
     cout<<" Be6 "<<Be6PDG<<" "<<He6LPDG<<endl;
     gMC->SetDecayMode(Be6PDG,bratio5,mode5);
     
     Int_t mode6[6][3];                  
     Float_t bratio6[6];
     Int_t Li6LPDG;
     
     if(p) Li6LPDG=p->PdgCode();
     p=db->GetParticle("Li6L");
     if(p) Li6LPDG=p->PdgCode();
     
     for (Int_t kz = 0; kz < 6; kz++) {
     bratio6[kz] = 0.;
     mode6[kz][0] = 0;
     mode6[kz][1] = 0; mode6[kz][2] = 0;
     }
     bratio6[0] = 100.;
     mode6[0][0] =Be6PDG;
     mode6[0][1] =-211  ;
     
     
     gMC->SetDecayMode(Li6LPDG,bratio6,mode6);
     cout<<" Li6 "<<Be6PDG<<" "<<Li6LPDG<<endl;
     Int_t mode7[6][3];                  
     Float_t bratio7[6];
     Int_t Li7LPDG, Be7PDG;
     // p=db->GetParticle("Ion_7_4");
     //if(p) Be7PDG=p->PdgCode();
     
     p=db->GetParticle("Li7L");
     if(p) Li7LPDG=p->PdgCode();
     
     for (Int_t kz = 0; kz < 6; kz++) {
      bratio7[kz] = 0.;
      mode7[kz][0] = 0;
      mode7[kz][1] = 0; mode7[kz][2] = 0;
      }
      bratio7[0] = 100.;
      mode7[0][0] =1000040070;//Be7PDG;
      mode7[0][1] =-211  ;
      
      
      gMC->SetDecayMode(Li7LPDG,bratio7,mode7);
      cout<<" Li7 "<<Be7PDG<<" "<<Li7LPDG<<endl;

      Int_t mode8[6][3];                  
      Float_t bratio8[6];
      Int_t Li8LPDG, Be8PDG;
      //p=db->GetParticle("Ion_8_4");
      //if(p) Be8PDG=p->PdgCode();
      
      p=db->GetParticle("Li8L");
      if(p) Li8LPDG=p->PdgCode();
      
      for (Int_t kz = 0; kz < 6; kz++) {
      bratio8[kz] = 0.;
      mode8[kz][0] = 0;
      mode8[kz][1] = 0; mode8[kz][2] = 0;
      }
      bratio8[0] = 100.;
      mode8[0][0] =1000040080;//Be8PDG;
      mode8[0][1] =-211  ;
      
      
      gMC->SetDecayMode(Li8LPDG,bratio8,mode8);
      cout<<" Be8 "<<Be8PDG<<" "<<Li8LPDG<<endl;
      
      Int_t mode9[6][3];                  
      Float_t bratio9[6];
      Int_t Li9LPDG, Be9PDG;
      //p=db->GetParticle("Ion_9_4");
      //if(p) Be9PDG=p->PdgCode();
      
      p=db->GetParticle("Li9L");
      if(p) Li9LPDG=p->PdgCode();
      
      for (Int_t kz = 0; kz < 6; kz++) {
      bratio9[kz] = 0.;
      mode9[kz][0] = 0;
      mode9[kz][1] = 0; mode9[kz][2] = 0;
      }
      bratio9[0] = 100.;
      mode9[0][0] =1000040090;//Be9PDG;
      mode9[0][1] =-211  ;
  
      
      gMC->SetDecayMode(Li9LPDG,bratio9,mode9);
      cout<<" Li9L "<< mode9[0][0]<<" "<<Li9LPDG<<endl;
      
      Int_t mode10[6][3];                  
      Float_t bratio10[6];
      Int_t Be7LPDG;
      //p=db->GetParticle("Ion_9_4");
      //if(p) Be9PDG=p->PdgCode();
      
      p=db->GetParticle("Be7L");
      if(p) Be7LPDG=p->PdgCode();
      
      for (Int_t kz = 0; kz < 6; kz++) {
      bratio10[kz] = 0.;
      mode10[kz][0] = 0;
      mode10[kz][1] = 0;
      mode10[kz][2] = 0;
      }
      bratio10[0] = 100.;
      mode10[0][0] =Be6PDG;
      mode10[0][1] =-211  ;
      mode10[0][1] =2212;
      
      
      gMC->SetDecayMode(Be7LPDG,bratio10,mode10);
      
      cout<<" Be7L "<<Be6PDG<<" "<<Be7LPDG<<endl;
      
      Int_t mode11[6][3];                  
      Float_t bratio11[6];
      Int_t Be8LPDG,B8PDG;
      //p=db->GetParticle("Ion_8_5");
      //if(p) B8PDG=p->PdgCode();
      
      p=db->GetParticle("Be8L");
      if(p) Be8LPDG=p->PdgCode();
      
      for (Int_t kz = 0; kz < 6; kz++) {
      bratio11[kz] = 0.;
      mode11[kz][0] = 0;
      mode11[kz][1] = 0;
      mode11[kz][2] = 0;
      }
      bratio11[0] = 100.;
      mode11[0][0] =1000050080;//B8PDG;
      mode11[0][1] =-211  ;
      //mode10[0][1] =2212;
      
      
      gMC->SetDecayMode(Be8LPDG,bratio11,mode11);
      
      cout<<" Be8L "<<B8PDG<<" "<<Be8LPDG<<endl;
     
  
   
   Int_t mode12[6][3];                  
   Float_t bratio12[6];
   Int_t B9PDG;
   p=db->GetParticle("Ion_9_5");
   if(p) B9PDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
     bratio12[kz] = 0.;
     mode12[kz][0] = 0;
     mode12[kz][1] = 0;
     mode12[kz][2] = 0;
   }
   bratio12[0] = 100.;
   mode12[0][0] =1000040080;//Be8PDG;
   mode12[0][1] =2212  ;
   //mode10[0][1] =2212;
   
   
   gMC->SetDecayMode(B9PDG,bratio12,mode12);
   
   cout<<" B9 "<<B9PDG<<" "<<endl;
   
   Int_t mode13[6][3];     Int_t Be9LPDG;
   Float_t bratio13[6];
   p=db->GetParticle("Be9L");
   if(p) Be9LPDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
     bratio13[kz] = 0.;
     mode13[kz][0] = 0;
     mode13[kz][1] = 0;
     mode13[kz][2] = 0;
   }
   bratio13[0] = 100.;
   mode13[0][0] =B9PDG;
   mode13[0][1] =-211  ;
   //mode10[0][1] =2212;
   
   
   gMC->SetDecayMode(Be9LPDG,bratio13,mode13);
   
   cout<<" Be9L "<<B9PDG<<" "<<Be9LPDG<<endl;


   Int_t mode14[6][3];                  
   Float_t bratio14[6];
   Int_t B10LPDG,C10PDG;
   //p=db->GetParticle("Ion_10_6");
   //if(p) C10PDG=p->PdgCode();
   p=db->GetParticle("B10L");
   if(p) B10LPDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio14[kz] = 0.;
      mode14[kz][0] = 0;
      mode14[kz][1] = 0;
       mode14[kz][2] = 0;
     }
   bratio14[0] = 100.;
   mode14[0][0] =1000060100;//C10PDG;
   mode14[0][1] =-211  ;
   //mode10[0][1] =2212;
  
 
   gMC->SetDecayMode(B10LPDG,bratio14,mode14);
 cout<<" B10L "<<C10PDG<<" "<<B10LPDG<<endl;

   Int_t mode15[6][3];                  
   Float_t bratio15[6];
   Int_t B11LPDG,C11PDG;
   //p=db->GetParticle("Ion_11_6");
   //if(p) C11PDG=p->PdgCode();
   p=db->GetParticle("B11L");
   if(p) B11LPDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio15[kz] = 0.;
      mode15[kz][0] = 0;
      mode15[kz][1] = 0;
       mode15[kz][2] = 0;
     }
   bratio15[0] = 100.;
   mode15[0][0] =1000060110;//C11PDG;
   mode15[0][1] =-211  ;
   //mode10[0][1] =2212;
  
 
   gMC->SetDecayMode(B11LPDG,bratio15,mode15);

 cout<<" B11L "<<C11PDG<<" "<<B11LPDG<<endl;

   Int_t mode16[6][3];                  
   Float_t bratio16[6];
   Int_t B12LPDG,C12PDG;
   // p=db->GetParticle("Ion_12_6");
   //if(p) C12PDG=p->PdgCode();
   p=db->GetParticle("B12L");
   if(p) B12LPDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio16[kz] = 0.;
      mode16[kz][0] = 0;
      mode16[kz][1] = 0;
      mode16[kz][2] = 0;
     }
   bratio16[0] = 100.;
   mode16[0][0] =1000060120;//C12PDG;
   mode16[0][1] =-211  ;
   //mode10[0][1] =2212;
  
 
   gMC->SetDecayMode(B12LPDG,bratio16,mode16);
cout<<" B12L "<<C12PDG<<" "<<B12LPDG<<endl;

   Int_t mode36[6][3];                  
   Float_t bratio36[6];
   Int_t C12LPDG,N12PDG;
   //p=db->GetParticle("Ion_12_7");
   //if(p) N12PDG=p->PdgCode();
   p=db->GetParticle("C12L");
   if(p) C12LPDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio36[kz] = 0.;
      mode36[kz][0] = 0;
      mode36[kz][1] = 0;
      mode36[kz][2] = 0;
     }
   bratio36[0] = 100.;
   mode36[0][0] =1000070120;//N12PDG;
   mode36[0][1] =-211  ;
   //mode10[0][1] =2212;
  
 
   gMC->SetDecayMode(C12LPDG,bratio36,mode36);
cout<<" C12L "<<N12PDG<<" "<<C12LPDG<<endl;

   Int_t mode37[6][3];                  
   Float_t bratio37[6];
   Int_t C13LPDG,N13PDG;
   //p=db->GetParticle("Ion_13_7");
   //if(p) N13PDG=p->PdgCode();
   p=db->GetParticle("C13L");
   if(p) C13LPDG=p->PdgCode();
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio37[kz] = 0.;
      mode37[kz][0] = 0;
      mode37[kz][1] = 0;
      mode37[kz][2] = 0;
     }
   bratio37[0] = 100.;
   mode37[0][0] =1000070130;//N13PDG;
   mode37[0][1] =-211  ;
   //mode10[0][1] =2212;
  
 
   gMC->SetDecayMode(C13LPDG,bratio37,mode37);
cout<<" C13L "<<N13PDG<<" "<<C12LPDG<<endl;

 

  Int_t mode25[6][3];                  
   Float_t bratio25[6];
   Int_t He7LLPDG;
   p=db->GetParticle("He7LL");
   if(p) He7LLPDG=p->PdgCode();
  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio25[kz] = 0.;
      mode25[kz][0] = 0;
      mode25[kz][1] = 0;mode25[kz][2] = 0;
      
     }
   bratio25[0] = 100.;
   mode25[0][0] =Li7LPDG;
   mode25[0][1] =-211  ;
  
  
 gMC->SetDecayMode(He7LLPDG,bratio25,mode25);
 cout<<" He7LL "<<Li7LPDG<<" "<<He7LLPDG<<endl;

 Int_t mode26[6][3];                  
   Float_t bratio26[6];
   Int_t He6LLPDG;
   p=db->GetParticle("He6LL");
   if(p) He6LLPDG=p->PdgCode();
  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio26[kz] = 0.;
      mode26[kz][0] = 0;
      mode26[kz][1] = 0;mode26[kz][2] = 0;
      
     }
   bratio26[0] = 100.;
   mode26[0][0] =Li6LPDG;
   mode26[0][1] =-211;
  
  
 gMC->SetDecayMode(He6LLPDG,bratio26,mode26);
 cout<<" He6LL "<<Li6LPDG<<" "<<He6LLPDG<<endl;


 Int_t mode24[6][3];                  
   Float_t bratio24[6];
   Int_t Li7LLPDG;
   p=db->GetParticle("Li7LL");
   if(p) Li7LLPDG=p->PdgCode();
  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio24[kz] = 0.;
      mode24[kz][0] = 0;
      mode24[kz][1] = 0;mode24[kz][2] = 0;
      
     }
   bratio24[0] = 100.;
   mode24[0][0] =Be7LPDG;
   mode24[0][1] =-211  ;
  
  
 gMC->SetDecayMode(Li7LLPDG,bratio24,mode24);
 cout<<" Li7LL "<<Be7LPDG<<" "<<Li7LLPDG<<endl;

 Int_t mode23[6][3];                  
   Float_t bratio23[6];
   Int_t Li8LLPDG;
   p=db->GetParticle("Li8LL");
   if(p) Li8LLPDG=p->PdgCode();  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio23[kz] = 0.;
      mode23[kz][0] = 0;
      mode23[kz][1] = 0; mode23[kz][2] = 0;
      
     }
   bratio23[0] = 100.;
   mode23[0][0] =Be8LPDG;
   mode23[0][1] =-211  ;
  
  
 gMC->SetDecayMode(Li8LLPDG,bratio23,mode23);
 cout<<" Li8LL "<<Be8LPDG<<" "<<Li8LLPDG<<endl;

 
Int_t m[6][3];                  
Float_t b[6];
Int_t Li9LLPDG;
p=db->GetParticle("Li9LL");
if(p) Li9LLPDG=p->PdgCode();


for (Int_t kz = 0; kz < 6; kz++) {
  b[kz] = 0.;
  m[kz][0] = 0;
  m[kz][1] = 0;m[kz][2] = 0;
  
}
b[0] = 100.;
m[0][0] =Be9LPDG;
m[0][1] =-211;

 
 gMC->SetDecayMode(Li9LLPDG,b,m);
 cout<<" Li9LL "<<Be9LPDG<<" "<<Li9LLPDG<<endl;




Int_t mode21[6][3];                  
 Float_t bratio21[6];
 Int_t Be9LLPDG;
 p=db->GetParticle("Be9LL");
 if(p) Be9LLPDG=p->PdgCode();
 
 
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio21[kz] = 0.;
   mode21[kz][0] = 0;
   mode21[kz][1] = 0;
   mode21[kz][2] = 0;
 }
 bratio21[0] = 100.;
 mode21[0][0] =Be9LPDG;
 mode21[0][1] =-211  ;
 mode21[0][1] =2212;
 
 gMC->SetDecayMode(Be9LLPDG,bratio21,mode21);
 cout<<" Be9LL "<<Be9LPDG<<" "<<Be9LLPDG<<endl;

 
 
 Int_t mode20[6][3];                  
 Float_t bratio20[6];
 Int_t Be10LLPDG;
 p=db->GetParticle("Be10LL");
 if(p) Be10LLPDG=p->PdgCode();
 
 
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio20[kz] = 0.;
   mode20[kz][0] = 0;
   mode20[kz][1] = 0;
   mode20[kz][2] = 0;
 }
 bratio20[0] = 100.;
 mode20[0][0] =B10LPDG;
 mode20[0][1] =-211  ;
 //mode10[0][1] =2212;
 
 
 gMC->SetDecayMode(Be10LLPDG,bratio20,mode20);
 cout<<" Be10LL "<<B10LPDG<<" "<<Be10LLPDG<<endl;
 
 Int_t mode19[6][3];                  
 Float_t bratio19[6];
 Int_t Be11LLPDG;
 p=db->GetParticle("Be11LL");
 if(p) Be11LLPDG=p->PdgCode();
 
 
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio19[kz] = 0.;
   mode19[kz][0] = 0;
   mode19[kz][1] = 0;
   mode19[kz][2] = 0;
 }
 bratio19[0] = 100.;
 mode19[0][0] =B11LPDG;
 mode19[0][1] =-211  ;
 //mode10[0][1] =2212;
 
 
 gMC->SetDecayMode(Be11LLPDG,bratio19,mode19);
 cout<<" Be11LL "<<B11LPDG<<" "<<Be11LLPDG<<endl;
 
 Int_t mode18[6][3];                  
 Float_t bratio18[6];
 Int_t Be12LLPDG;
 p=db->GetParticle("Be12LL");
 if(p) Be12LLPDG=p->PdgCode();
 
 
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio18[kz] = 0.;
   mode18[kz][0] = 0;
   mode18[kz][1] = 0;
   mode18[kz][2] = 0;
 }
 bratio18[0] = 100.;
 mode18[0][0] =B12LPDG;
 mode18[0][1] =-211  ;
 //mode10[0][1] =2212;
 
 
 gMC->SetDecayMode(Be12LLPDG,bratio18,mode18);
 cout<<" Be12LL "<<B12LPDG<<" "<<Be12LLPDG<<endl;
 
 
 Int_t mode17[6][3];                  
 Float_t bratio17[6];
 Int_t B12LLPDG;
 
 p=db->GetParticle("B12LL");
 if(p) B12LLPDG=p->PdgCode();
 
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio17[kz] = 0.;
   mode17[kz][0] = 0;
   mode17[kz][1] = 0;
   mode17[kz][2] = 0;
 }
 bratio17[0] = 100.;
 mode17[0][0] =C12LPDG;
 mode17[0][1] =-211  ;
 //mode10[0][1] =2212;
 
 
 gMC->SetDecayMode(B12LLPDG,bratio17,mode17);
 cout<<" B12LL "<<C12LPDG<<" "<<B12LLPDG<<endl;
 
 
 Int_t mode27[6][3];                  
 Float_t bratio27[6];
 Int_t B11LLPDG;
 p=db->GetParticle("B11LL");
 if(p) B11LLPDG=p->PdgCode();
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio27[kz] = 0.;
   mode27[kz][0] = 0;
   mode27[kz][1] = 0;
   mode27[kz][2] = 0;
   
     }
 bratio27[0] = 100.;
 mode27[0][0] =B10LPDG;
 mode27[0][1] =-211;
 mode27[0][2] =2212;
 gMC->SetDecayMode(B11LLPDG,bratio27,mode27);
 cout<<" B11LL "<<B10LPDG<<" "<<B11LLPDG<<endl;
 
 Int_t mode28[6][3];                  
 Float_t bratio28[6];
 Int_t B13LLPDG;
 p=db->GetParticle("B13LL");
 if(p) B13LLPDG=p->PdgCode();
 
 
 for (Int_t kz = 0; kz < 6; kz++) {
   bratio28[kz] = 0.;
   mode28[kz][0] = 0;
   mode28[kz][1] = 0;mode28[kz][2] = 0;
   
 }
 bratio28[0] = 100.;
 mode28[0][0] =C13LPDG;
 mode28[0][1] =-211;
 
 
 gMC->SetDecayMode(B13LLPDG,bratio28,mode28);
 cout<<" B13LL "<<C13LPDG<<" "<<B13LLPDG<<endl;

 
 /*
   Int_t mode27[6][3];                  
   Float_t bratio27[6];
   Int_t H5LLPDG;
   p=db->GetParticle("H5LL");
   if(p) H5LLPDG=p->PdgCode();
  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio27[kz] = 0.;
      mode27[kz][0] = 0;
      mode27[kz][1] = 0;
      
     }
   bratio27[0] = 100.;
   mode27[0][0] =He5PDG;
   mode27[0][1] =-211  ;
  
  
 gMC->SetDecayMode(H5LLPDG,bratio27,mode27);

 Int_t mode28[6][3];                  
   Float_t bratio28[6];
   Int_t He5LLPDG;
   p=db->GetParticle("He5LL");
   if(p) He5LLPDG=p->PdgCode();
  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio28[kz] = 0.;
      mode28[kz][0] = 0;
      mode28[kz][1] = 0;mode28[kz][2] = 0;
      
     }
   bratio28[0] = 100.;
   mode28[0][0] =He5PDG;
   mode28[0][1] =-211  ;
   mode28[0][1] =2212;
  
 gMC->SetDecayMode(He5LLPDG,bratio28,mode28);

 Int_t mode29[6][3];                  
   Float_t bratio29[6];
   Int_t Li10LPDG;
   p=db->GetParticle("Li10L");
   if(p) Li10LPDG=p->PdgCode();
  
   
   for (Int_t kz = 0; kz < 6; kz++) {
      bratio29[kz] = 0.;
      mode29[kz][0] = 0;
      mode29[kz][1] = 0;
      
     }
   bratio29[0] = 100.;
   mode29[0][0] =He6LPDG;
   mode29[0][1] =1000020040  ;
   
  
   gMC->SetDecayMode(Li10LPDG,bratio29,mode29);
*/



    
}