4 double x[6] = {1.46513,2.4752,3.44589,4.44633,5.42126,7.93998};
5 double y[6] = {13.53,17.8427,24.5361,33.1036,37.9249,36.5359};
6 double exl[6] = {0.917728,0.98571,0.91773,0.95133,1.46119,2.10699};
7 double exh[6] = {1.08768,0.98571,0.98571,0.91773,1.52993,3.94362};
8 double eyl[6] = {2.2105,3.9125,4.5932,5.9587,8.5096,6.9779};
9 double eyh[6] = {3.576,3.4038,4.5971,5.4423,8.3394,6.4654};
11 TGraphAsymmErrors*
gargamelle =
new TGraphAsymmErrors(6,x,y,exl,exh,eyl,eyh);
12 gargamelle->SetTitle(
"Gargamelle #nu_{e} CC data");
13 gargamelle->SetMarkerStyle(21);
14 gargamelle->SetMarkerColor(
kViolet-1);
15 gargamelle->SetLineColor(
kViolet-1);
16 gargamelle->SetLineWidth(2);
30 TGraphAsymmErrors *
t2k=
new TGraphAsymmErrors(1,x_t2k,y_t2k,exl_t2k,exh_t2k,eyl_t2k,eyh_t2k);
31 t2k->SetMarkerStyle(24);
33 t2k->SetMarkerColor(4);
37 double x_nova[4]={1.25,1.75,2.25,2.75};
38 double y_nova[4]={22.0921,28.8235,28.16,37.6431};
40 double ey_nova[4]={6.58498,5.99768,5.53289,6.66383};
42 TGraphErrors *
nova=
new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova);
45 nova->SetLineColor(2);
47 nova->SetLineWidth(3);
50 TGraphErrors *
nova_stat=
new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova_stat);
51 nova_stat->SetMarkerStyle(8);
52 nova_stat->SetMarkerSize(1.2);
53 nova_stat->SetLineWidth(3);
56 TFile *
fgenie=
new TFile(
"xsec_genie.root");
57 TGraph *
gr_genie=(TGraph*)fgenie->Get(
"nu_e_C12/tot_cc");
58 TSpline3*
spline_nue=
new TSpline3(
"mysplinenue",gr_genie);
59 TH1F *
Xnue=
new TH1F(
"Xnue",
"",100,0,10);
64 double binctr=Xnue->GetBinCenter(
ibin);
65 double xnue=(spline_nue->Eval(binctr)/12.)*10.0;
66 Xnue->SetBinContent(
ibin,xnue);
68 x_genie[
ibin-1]=binctr;
74 TFile *
fband=
new TFile(
"genie_band.root");
75 TH1F *
hplus=(TH1F*)fband->Get(
"ccnue_1p_total");
76 TH1F *
hminus=(TH1F*)fband->Get(
"ccnue_1m_total");
83 TH1F *
Xnue2=
new TH1F(
"Xnue2",
"",50,0,10);
85 double binctr=Xnue2->GetBinCenter(
ibin);
87 double xnue2=(spline_nue->Eval(binctr)/12.)*10.0;
88 Xnue2->SetBinContent(
ibin,xnue2);
91 x_band[
ibin-2]=binctr;
95 eyl_band[
ibin-2]=xnue2*(hminus->GetBinContent(
ibin));
96 eyh_band[
ibin-2]=xnue2*(hplus->GetBinContent(
ibin));
100 TGraphAsymmErrors *
band=
new TGraphAsymmErrors(49,x_band,y_band,exl_band,exh_band,eyl_band,eyh_band);
101 band->SetFillColor(3);
102 band->SetLineColor(3);
116 TGraph *
genie=
new TGraph(100,x_genie,y_genie);
117 genie->SetLineStyle(2);
118 genie->SetLineWidth(3);
119 genie->SetLineColor(1);
120 genie->SetFillColor(3);
122 gStyle->SetOptStat(0);
124 TCanvas *
c1=
new TCanvas(
"Xsec_Enue_withNOvA",
"",700,600);
125 TH2F *
h2f=
new TH2F(
"h2f",
"",10,0,10,10,0,59);
127 h2f->SetXTitle(
"Electron Neutrino Energy (GeV)");
128 h2f->SetYTitle(
"#nu_{e} CC #it{#sigma} (#times10^{-39} cm^{2}/nucleon)");
129 h2f->GetXaxis()->CenterTitle();
130 h2f->GetYaxis()->CenterTitle();
133 TLatex*
prelim =
new TLatex(.9, .95,
"NO#nuA Preliminary");
134 prelim->SetTextColor(
kBlue+1);
136 prelim->SetTextSize(2/30.);
137 prelim->SetTextAlign(32);
148 band->Draw(
"E3 same");
150 gargamelle->Draw(
"P same");
152 nova->Draw(
"P z same");
154 gStyle->SetEndErrorSize(4);
156 nova_stat->Draw(
"P same");
158 genie->Draw(
"c same");
163 TLegend *
leg=
new TLegend(0.6,0.6,0.89,0.89);
164 leg->SetLineColor(0);
165 leg->SetFillColor(0);
166 leg->AddEntry(genie,
"GENIE 2.8.0 #pm 1#sigma on C12",
"fl");
167 leg->AddEntry(gargamelle,
"Gargamelle #nu_{e} data",
"pl");
168 leg->AddEntry(t2k,
"T2K #nu_{e} data",
"pl");
169 leg->AddEntry(nova_stat,
"NOvA #nu_{e} data",
"pl");
170 leg->AddEntry(nova_stat,
"stat. uncert.",
"E");
171 leg->AddEntry(nova,
"stat. + syst. uncert.",
"E");
THE MAIN GENIE PROJECT NAMESPACE
TGraphAsymmErrors * gargamelle
Hold drift constants from current run.