| 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 | } |
|---|
