R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

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

                                  /* further flags - optional */

                1, 		  /* 1 = save all particles, 0 = only tracked  */
		0,                /* Unused                                    */
                0,                /* 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;
}