{ makeDistributionManager()->Exec("elementary"); //Activate the Krivoruchenko Dalitz 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 ND13+","delta_bold",0,0,0,0); //Attach a control histogram TH1F * pp_sum = new TH1F ("pp_sum","pp DiLepton mass (coherent sum)",100,1.,2.3); pp_sum->Sumw2(); TH1F * delta = new TH1F ("delta","D mass ",100,1.0,2.3); delta->Sumw2(); // my_reaction->Do(pp_sum,"_x=[dilepton]->M()"); my_reaction->Do(delta,"_x=[ND13+]->M()"); my_reaction->Print(); my_reaction->Preheating(100); my_reaction->Loop(100000); // */ //=================================== PReaction *my_reaction2 = new PReaction("3.5","p","p","p ND13+ [dilepton [e+ e-] p]","delta_dalitz",0,0,0,0); // PReaction *my_reaction2 = new PReaction("3.5","p","p","p D+ [dilepton [e+ e-] p]","delta_dalitz",1,0,0,1); //PReaction *my_reaction2 = new PReaction("3.5","p","p","p ND13+","delta_mass",1,0,0,1); //Attach a control histogram TH1F *pp_sum2 = new TH1F ("pp_sum2","pp DiLepton mass (coherent sum)",100,0.0,1.2); pp_sum2->Sumw2(); TH1F *delta2 = new TH1F ("delta2","D mass ",100,1.0,2.3); TH1F *ratio = new TH1F ("ratio","ratio mass ",100,1.0,2.3); delta2->Sumw2(); //my_reaction2->Do("_w = _w * {ND13+_dalitz}->GetBR([ND13+]->M());"); my_reaction2->Do(pp_sum2,"_x=[dilepton]->M()"); my_reaction2->Do(delta2,"_x=[ND13+]->M()"); my_reaction2->Print(); my_reaction2->Preheating(100); my_reaction2->Loop(100000); ratio = (TH1F*)delta2->Clone(); delta2->Divide(delta); TCanvas *n1 = new TCanvas(); delta->Draw(""); TCanvas *n2 = new TCanvas(); delta2->Draw(""); TCanvas *n3 = new TCanvas(); ratio->Draw(""); //pp_sum2->Draw(""); //pp_sum2->Draw(""); }