UserDecay.C

Go to the documentation of this file.

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();
   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] =2112  ;
   mode[0][1] =AlphaPDG  ;
   
 /*  bratio[1] = 50.;
   mode[1][0] =2212  ;
   mode[1][1] =AlphaPDG  ;
    
  */
   gMC->SetDecayMode(He5PDG,bratio,mode);
    

    
}