R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int pip_mupnu_epnunu () 
{
  PParticle pip ("pi+", 1.0 );                    /* Source particle, Ekin [GeV] */
    
  PReaction *my_reaction = new PReaction (
                &pip ,                 /* Source particle */
	        "nu mu+ [e+ nu nu]",   /* Decay scheme of source particle */
        	"pip_mupnu_epnunu" ,   /* Output TTree filename */

                1, 		  /* 1 = save all particles, 0 = only tracked  */
		0,                /* Unused                                    */
                1,                /* 1 = generate production vertex, 0 = don't */
		0                 /* 1 = generate text output file , 0 = don't */                
               );
    my_reaction->Print ();

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

    return 0;
}