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