R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int lambda_decays () 
{
  PFireball* source_lambda = new PFireball (
   "Lambda", 2.0, 0.100, 0., 1., 0., 0., 0., 0., 0.);

  PDecayChannel* ch_FBall  = new PDecayChannel;
  ch_FBall ->AddChannel ( 1. , "Lambda" );

  PDecayManager* pdm = new PDecayManager;
  pdm->SetVerbose (1);

  /* To activate all the available decays of Lambda */

  pdm->SetDefault ("Lambda");

  /* Instead, to create your own decay list with your BR */
/*
  PDecayChannel* ch_lam_decays = new PDecayChannel;
  ch_lam_decays->AddChannel ( 0.639, "p", "pi-" );
  ch_lam_decays->AddChannel ( 0.358, "n", "pi0" );
  pdm->AddChannel ("Lambda", ch_lam_decays );
*/
  pdm->InitReaction ( source_lambda , ch_FBall ); 
  pdm->PrintReactionList() ;

  pdm->loop (10000, 0, "lambda_decays", 1, 0, 0, 1, 1);
  return 0;
}