makeDistributionManager()->Exec("elementary");
    makeDistributionManager()->Exec("dalitz_mod: krivoruchenko");
    makeDistributionManager()->Exec("dalitz_mod: static_br_thresh=0.100 ; flat_generator");
    
    PReaction *my_reaction = 
	new PReaction("3.5","p","p","p D+",NULL,1,0,0,0);

    TH1F * pp_sum = 
	new TH1F ("pp_sum","pp DiLepton mass (coherent sum)",100,0.,0.6);
    pp_sum->Sumw2();
    TH1F * delta = 
	new TH1F ("delta","D mass ",100,0.9,1.6);
    delta->Sumw2();


    my_reaction->Do(pp_sum,"_x=[dilepton]->M()");
    my_reaction->Do(delta,"_x=[D+]->M()");
    my_reaction->Print();

    my_reaction->Preheating(100);
    my_reaction->Loop(100000);

    PReaction *my_reaction2 = 
	new PReaction("3.5","p","p","p D+ [dilepton p]",NULL,1,0,0,0);

    TH1F * pp_sum2 = 
	new TH1F ("pp_sum2","pp DiLepton mass (coherent sum)",100,0.,0.6);
    pp_sum2->Sumw2();
    TH1F * delta2 = 
	new TH1F ("delta2","D mass ",100,0.9,1.6);
    delta2->Sumw2();

    my_reaction2->Do(pp_sum2,"_x=[dilepton]->M()");
    my_reaction2->Do(delta2,"_x=[D+]->M()");
    my_reaction2->Print();

    my_reaction2->Preheating(100);
    my_reaction2->Loop(100000);

    //    delta2->Divide(delta);
    delta->Draw("");
    delta2->Draw("same");[/QUOTE]