Go to the source code of this file.
|
| my_reaction | Do ("pp_miss = [p + p] - ([p,1] + [p,2]); total_miss = [p + p] - ([p,1] + [p,2] + [pi+] + [pi-])") |
| my_reaction | Do (histo1,"_x= pp_miss->M2()") |
| my_reaction | Do (histo2,"_x= pp_miss->P()") |
| my_reaction | Do (histo3,"_x= ([p,1]->Theta() * 180.)/TMath::Pi()") |
| my_reaction | Do (histo1_cut,"if pp_miss->P() > 0.3; _x= pp_miss->M2()") |
| my_reaction | Do (histo2_cut,"if pp_miss->P() > 0.3; _x= pp_miss->P()") |
| my_reaction | Do (histo2_cut2,"#mom = 0.; if pp_miss->P() > 0.3 && pp_miss->P() < 0.6 ; #mom = 1.") |
| my_reaction | Print () |
| cout<< my_reaction.Loop(50000)<< " events recorded"<< endl;TCanvas *c1=new TCanvas("c1","Canvas");c1-> | SetLogy (1) |
| c1 | Divide (2) |
| c1 | cd (1) |
| histo1 | Draw ("") |
| histo1_cut | Draw ("same") |
| histo1_cut | SetLineColor (2) |
| c1 | cd (2) |
| histo2_cut2 | SetLineColor (3) |
|
| TH1F * | histo1_cut = new TH1F ("histo1_cut","p p missing mass",100,0.1,1.0) |
| TH1F * | histo2 = new TH1F ("histo2","p p missing momentum",100,0.1,1.0) |
| TH1F * | histo2_cut = new TH1F ("histo2_cut","p p missing momentum",100,0.1,1.0) |
| TH1F * | histo2_cut2 = new TH1F ("histo2_cut2","p p missing momentum",100,0.1,1.0) |
| TH1F * | histo3 = new TH1F ("histo3","Theta of the first proton",100,0,180) |
| PReaction | my_reaction ("3.5","p","p","p p w [pi+ pi- pi0]") |
◆ cd() [1/2]
◆ cd() [2/2]
◆ Divide()
◆ Do() [1/7]
◆ Do() [2/7]
◆ Do() [3/7]
◆ Do() [4/7]
◆ Do() [5/7]
◆ Do() [6/7]
◆ Do() [7/7]
◆ Draw() [1/2]
◆ Draw() [2/2]
◆ Print()
◆ SetLineColor() [1/2]
◆ SetLineColor() [2/2]
◆ SetLogy()
◆ histo1_cut
| TH1F* histo1_cut = new TH1F ("histo1_cut","p p missing mass",100,0.1,1.0) |
◆ histo2
| TH1F* histo2 = new TH1F ("histo2","p p missing momentum",100,0.1,1.0) |
◆ histo2_cut
| TH1F* histo2_cut = new TH1F ("histo2_cut","p p missing momentum",100,0.1,1.0) |
◆ histo2_cut2
| TH1F* histo2_cut2 = new TH1F ("histo2_cut2","p p missing momentum",100,0.1,1.0) |
◆ histo3
| TH1F* histo3 = new TH1F ("histo3","Theta of the first proton",100,0,180) |
◆ my_reaction
| PReaction my_reaction("3.5","p","p","p p w [pi+ pi- pi0]") |
( |
"3.5" | , |
|
|
"p" | , |
|
|
"p" | , |
|
|
"p p w " | [pi+pi- pi0] ) |