R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int eta_read_ascii () 
{
  PParticle eta ("eta"), pip ("pi+"), pim ("pi-"), pi0 ("pi0");
  PParticle* ar_eta_3pi[] = { &eta , &pip , &pim , &pi0 } ;
  PChannel* ch_eta_3pi  = new PChannel ( ar_eta_3pi , 3);

  PReaction* R = new PReaction (&ch_eta_3pi, "eta_decays", 1);

/* Or for all the decays of eta:
   Remove the commands above. Add instead:
     PReaction* R = new PReaction ("eta_decays");
     R->SetDecayAll (1.);
*/
  PProjector* input = new PProjector ();
  input->AddInputASCII ("eta_sample.txt"); 
  input->Do ("readline {@px @py @pz @m} ");
  input->Do ("myeta = P3M (px, py, pz, m) ");
  input->Do ("myeta->SetID (eta.pid)");  
  input->Do ("push (myeta)");

  R->AddPrologueBulk (input);
  R->Print();
  R->Loop(1000); 
  return 0;
}