{
//look at higgs decays
Setup::decay_particle=-24;
Setup::mass_power=2;
Setup::mass_scale_on=true;

// Setup histograms
int n_bins=60;
double default_min_bin=0.0;
double default_max_bin=1.1;
Setup::SetHistogramDefaults(n_bins,default_min_bin,default_max_bin);

// Setup User Histograms
//Setup::UserTreeAnalysis = "RhoRhoUserTreeAnalysis";

// Description
Setup::gen1_desc_1="Pythia + Tauola Interface Test";
Setup::gen1_desc_2=" $W^- \\rightarrow \\pi^- \\nu_{\\tau} \\bar{\\nu_{\\tau}} $";
 Setup::gen1_desc_3="Using $Mass^2$. New";


if (Setup::stage==0)
    printf("Setup loaded from SETUP.C, ANALYSIS stage.\n");
else 
    printf("Setup loaded from SETUP.C, GENERATION stage %i.\n",Setup::stage);

Setup::SuppressDecay(111); // suppress pi0 decays

};