R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int pip_mupnu_TNtuple () 
{ 
  TFile *file_TNtuple = new TFile ("pip_mupnu_tntuple.root", "RECREATE");

  TNtuple *my_ntuple = new TNtuple (
			"ntu",  "Monitor ntuple", 
			"thelab:philab:plab:pt:ylab"  );

  PParticle pip ("pi+", 0.1 );                    /* Source particle, Ekin [GeV] */

  PReaction *my_reaction = new PReaction ( &pip, "mu+ nu", "pip_mupnu", 1, 0, 0, 0);

  my_reaction->Print ();

  my_reaction->Do ("thelab = [mu+]->Theta() * 180./TMath::Pi()" );
  my_reaction->Do ("philab = [mu+]->Phi() * 180./TMath::Pi()" );
  my_reaction->Do ("plab = [mu+]->P()" );
  my_reaction->Do ("pt = [mu+]->Pt()" );
  my_reaction->Do ("ylab = [mu+]->Rapidity()" );

  my_reaction->Do ("#th = 1; if thelab<2 ; #th = 0");
  my_reaction->Do ("#ph = 1; if philab > 30. || phi < -30. ; #ph = 0");

  my_reaction->Output (my_ntuple);

  my_reaction->Loop (10000, 1);		    /* No of events, status @ end */

  file_TNtuple->Write ();
  file_TNtuple->Close ();
  return 0;
}