R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int eta_cocktail_simple ()
{
  PParticle p1 ("p", 2.2), p2 ("p"), comp = p1 + p2;

  PDecayManager* pdm = new PDecayManager;

  PDecayChannel* dc   = new PDecayChannel;
  dc->AddChannel (0.42, "p", "p", "eta" ); // direct eta production
  dc->AddChannel (0.58, "p", "NS11+" );    // or: first N*(1535)

  // The above decays will be connected with the composite particle by InitReaction.

  // Then, N*(1535) should decay

  PDecayChannel* nstar_decay = new PDecayChannel;
  nstar_decay->AddChannel (1.0 , "eta", "p");
  pdm->AddChannel ("NS11+", nstar_decay );

  pdm->InitReaction ( &comp , dc );
  pdm->Print();
  pdm->loop (100, 0, "eta_sample_simple", 1, 0, 0, 1, 1);

  return 0;
}