1 | { |
---|
2 | |
---|
3 | TCanvas* c1 = new TCanvas("c1","SPL Flux",100,100,600,600); |
---|
4 | c1->SetLogy(); |
---|
5 | |
---|
6 | TLegend* leg = new TLegend(0.20,0.20,0.45,0.55); |
---|
7 | leg->SetTextFont(gStyle->GetTextFont()); |
---|
8 | leg->SetTextSize(0.04); |
---|
9 | |
---|
10 | |
---|
11 | |
---|
12 | TH1D* hdum = new TH1D("hdum","Rates 5yrs (+)",5,0.08,1.08); |
---|
13 | hdum->SetMaximum(1e3); |
---|
14 | hdum->SetMinimum(1e-3); |
---|
15 | hdum->SetXTitle("E_{#nu} (GeV)"); |
---|
16 | hdum->SetYTitle("Number of Events"); |
---|
17 | hdum->Draw(); |
---|
18 | |
---|
19 | TTree T1("t1","SPL Rate for 5yrs +"); |
---|
20 | T1.ReadFile("./rate5+.dat","ene/F:numuCC/F:nueB/F:nubare/F:pi0/F:nue1d/F:nue3d/F"); |
---|
21 | |
---|
22 | T1.Draw("ene>>hnueB(5,0.08,1.08)","nueB","SAME"); |
---|
23 | |
---|
24 | T1.SetLineColor(6); |
---|
25 | T1.Draw("ene>>hnumuCC(5,0.08,1.08)","numuCC","SAME"); |
---|
26 | |
---|
27 | T1.SetLineColor(7); |
---|
28 | T1.Draw("ene>>hnubare(5,0.08,1.08)","nubare","SAME"); |
---|
29 | |
---|
30 | T1.SetLineColor(3); |
---|
31 | T1.Draw("ene>>hpi0(5,0.08,1.08)","pi0","SAME"); |
---|
32 | |
---|
33 | T1.SetLineColor(4); |
---|
34 | T1.Draw("ene>>hnue3d(5,0.08,1.08)","nue3d","SAME"); |
---|
35 | |
---|
36 | T1.SetLineColor(4); |
---|
37 | T1.SetLineStyle(2); |
---|
38 | T1.Draw("ene>>hnue1d(5,0.08,1.08)","nue1d","SAME"); |
---|
39 | |
---|
40 | |
---|
41 | leg->AddEntry(hnue3d,"#nu_{e} (#theta_{13} = 3^{o})","l"); |
---|
42 | leg->AddEntry(hnue1d,"#nu_{e} (#theta_{13} = 1^{o})","l"); |
---|
43 | leg->AddEntry(hnueB,"#nu_{e} Beam","l"); |
---|
44 | leg->AddEntry(hpi0,"#pi^{0} (NC)","l"); |
---|
45 | leg->AddEntry(hnumuCC,"#mu missId ","l"); |
---|
46 | leg->AddEntry(hnubare,"#bar{#nu_{e}}","l"); |
---|
47 | leg->Draw(); |
---|
48 | |
---|
49 | } |
---|