R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int pp_pppi0_gg () 
{

    PReaction *my_reaction = new PReaction ("0.2797", /* Ekin of beam particle */

						      /* Caution: 0.2797  = momentum
								 "0.2797" = Ekin ! */

				"p", "p", 	      /* Substrates */
				"p p pi0 [g g]",      /* Scheme of production & decay */
				"pp_pppi0_gg", 	      /* Output TTree filename */
				1, 0, 0, 0 );
    my_reaction->Print ();

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

    return 0;
}