proton_Li7_macro.C

Go to the documentation of this file.

//TITLE Example for a p + 7Li reaction
 
{
 
    makeDistributionManager()->Exec("nucleus_fermi:proton");
 
    PReaction my_reaction("3.5","p","7Li","p (n) phi [K+ K-] (6Li)","delme");
    //both participant and spectator must be in brackets "( )"
    //the last particle is the spectator
 
    my_reaction.Print();
 
    //check momentum distribution of K+
    TH1F * histo1 = new TH1F ("histo1","K+ momentum",100,0,3);
    my_reaction.Do(histo1,"_x = [K+]->P()");
 
    //mass of phi
    TH1F * histo2 = new TH1F ("histo2","phi mass",100,0.9,1.2);
    my_reaction.Do(histo2,"_x = [phi]->M()");
 
    my_reaction.Loop(1000);
 
    histo2->Draw();
 
 
    //same at much lower energy
    //slightly below threshold, this is really time consuming
    PReaction my_reaction2("2.5","p","7Li","p (n) phi [K+ K-] (6Li)","delme");
    TH1F * histo2a = new TH1F ("histo2a","phi mass",100,0.9,1.2);
    my_reaction2.Do(histo2a,"_x = [phi]->M()");
    my_reaction2.Loop(1000);
    histo2a->Draw("same");
 
}