R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int pi0_gg_long () 
{
    //   pi0 (Tb = 0.100)
    PParticle pi0 ("pi0", 0.100);

    //    Decay   pi0 ---> gamma+gamma

    PParticle g1 ("g") , g2 ("g");

    PParticle* parttable_gg[] = { &pi0, &g1 , &g2 };

    PChannel* pi0_ggdecay = new PChannel ( parttable_gg , 2 , 1 );

    //    Constructing the full reaction

    PChannel* table_of_channels[] = { pi0_ggdecay };

    PReaction* my_reaction = new PReaction 
               ( table_of_channels , "myoutfile_pi0_gg", 1 , 1 , 0 , 0 , 1 );
/*
   o Table of channels
   o Output file name
   o Number of channels in table
   o Options:
      f0: 0=Tracked, 1=All particles in ROOT output file_name.root
      f1: unused
      f2: 0=Do not calculate, 1=Calculate production vertices for product particles
      f3: 0=Do not produce, 1=Produce GEANT ascii output file_name.evt (tracked particles always)
*/

    my_reaction->Print ();

    //    Generate events!

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

    return 0;
}