Functions | Variables
Xsec_final.C File Reference

Go to the source code of this file.

Functions

gargamelle SetTitle ("Gargamelle #nu_{e} CC data")
 
gargamelle SetMarkerStyle (21)
 
gargamelle SetMarkerColor (kViolet-1)
 
gargamelle SetLineColor (kViolet-1)
 
gargamelle SetLineWidth (2)
 
t2k SetMarkerStyle (24)
 
t2k SetLineColor (4)
 
t2k SetMarkerColor (4)
 
nova SetLineColor (2)
 
nova SetLineWidth (3)
 
nova_stat SetMarkerStyle (8)
 
nova_stat SetMarkerSize (1.2)
 
 for (int ibin=1;ibin< Xnue->GetNbinsX()+1;++ibin)
 
band SetFillColor (3)
 
band SetLineColor (3)
 
genie SetLineStyle (2)
 
genie SetLineColor (1)
 
gStyle SetOptStat (0)
 
h2f SetXTitle ("Electron Neutrino Energy (GeV)")
 
h2f SetYTitle ("#nu_{e} CC #it{#sigma} (#times10^{-39} cm^{2}/nucleon)")
 
h2f GetXaxis () -> CenterTitle()
 
h2f GetYaxis () -> CenterTitle()
 
h2f Draw ()
 
prelim SetTextColor (kBlue+1)
 
prelim SetNDC ()
 
prelim SetTextSize (2/30.)
 
prelim SetTextAlign (32)
 
band Draw ("E3 same")
 
gargamelle Draw ("P same")
 
nova Draw ("P z same")
 
gStyle SetEndErrorSize (4)
 
genie Draw ("c same")
 
leg SetLineColor (0)
 
leg SetFillColor (0)
 
leg AddEntry (genie,"GENIE 2.8.0 #pm 1#sigma on C12","fl")
 
leg AddEntry (gargamelle,"Gargamelle #nu_{e} data","pl")
 
leg AddEntry (t2k,"T2K #nu_{e} data","pl")
 
leg AddEntry (nova_stat,"NOvA #nu_{e} data","pl")
 
leg AddEntry (nova_stat,"stat. uncert.","E")
 
leg AddEntry (nova,"stat. + syst. uncert.","E")
 
leg Draw ("same")
 

Variables

double y [6] = {13.53,17.8427,24.5361,33.1036,37.9249,36.5359}
 
double exl [6] = {0.917728,0.98571,0.91773,0.95133,1.46119,2.10699}
 
double exh [6] = {1.08768,0.98571,0.98571,0.91773,1.52993,3.94362}
 
double eyl [6] = {2.2105,3.9125,4.5932,5.9587,8.5096,6.9779}
 
double eyh [6] = {3.576,3.4038,4.5971,5.4423,8.3394,6.4654}
 
TGraphAsymmErrors * gargamelle = new TGraphAsymmErrors(6,x,y,exl,exh,eyl,eyh)
 
double x_t2k [1] ={1.343}
 
double y_t2k [1] ={11.1}
 
double exl_t2k [1] ={0.918}
 
double exh_t2k [1] ={1.782}
 
double eyl_t2k [1] ={2.1}
 
double eyh_t2k [1] ={2.1}
 
TGraphAsymmErrors * t2k =new TGraphAsymmErrors(1,x_t2k,y_t2k,exl_t2k,exh_t2k,eyl_t2k,eyh_t2k)
 
double x_nova [4] ={1.25,1.75,2.25,2.75}
 
double y_nova [4] ={22.0921,28.8235,28.16,37.6431}
 
double ex_nova [4] ={0.25,0.25,0.25,0.25}
 
double ey_nova [4] ={6.58498,5.99768,5.53289,6.66383}
 
TGraphErrors * nova =new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova)
 
double ey_nova_stat [4] ={4.88965,3.12006,2.64657,2.54683}
 
TGraphErrors * nova_stat =new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova_stat)
 
TFile * fgenie =new TFile("xsec_genie.root")
 
TGraph * gr_genie =(TGraph*)fgenie->Get("nu_e_C12/tot_cc")
 
TSpline3 * spline_nue =new TSpline3("mysplinenue",gr_genie)
 
TH1F * Xnue =new TH1F("Xnue","",100,0,10)
 
double x_genie [100] ={0.}
 
double y_genie [100] ={0.}
 
TFile * fband =new TFile("genie_band.root")
 
TH1F * hplus =(TH1F*)fband->Get("ccnue_1p_total")
 
TH1F * hminus =(TH1F*)fband->Get("ccnue_1m_total")
 
double x_band [49] ={0.}
 
double y_band [49] ={0.}
 
double exl_band [49] ={0.}
 
double exh_band [49] ={0.}
 
double eyl_band [49] ={0.}
 
double eyh_band [49] ={0.}
 
TH1F * Xnue2 =new TH1F("Xnue2","",50,0,10)
 
TGraphAsymmErrors * band =new TGraphAsymmErrors(49,x_band,y_band,exl_band,exh_band,eyl_band,eyh_band)
 
TGraph * genie =new TGraph(100,x_genie,y_genie)
 
TCanvas * c1 =new TCanvas("Xsec_Enue_withNOvA","",700,600)
 
TH2F * h2f =new TH2F("h2f","",10,0,10,10,0,59)
 
TLatex * prelim = new TLatex(.9, .95, "NO#nuA Preliminary")
 
TLegend * leg =new TLegend(0.6,0.6,0.89,0.89)
 

Function Documentation

leg AddEntry ( genie  ,
"GENIE 2.8.0 #pm 1#sigma on C12"  ,
"fl"   
)
leg AddEntry ( gargamelle  ,
"Gargamelle #nu_{e} data ,
"pl"   
)
leg AddEntry ( t2k  ,
"T2K #nu_{e} data ,
"pl"   
)
leg AddEntry ( nova_stat  ,
"NOvA #nu_{e} data ,
"pl"   
)
leg AddEntry ( nova_stat  ,
"stat. uncert."  ,
"E"   
)
leg AddEntry ( nova  ,
"stat. + syst. uncert."  ,
"E"   
)
h2f Draw ( )
band Draw ( "E3 same"  )
gargamelle Draw ( "P same"  )
nova Draw ( "P z same"  )
genie Draw ( "c same"  )
leg Draw ( "same"  )
for ( int  ibin = 1; ibin<Xnue->GetNbinsX()+1; ++ibin)

Definition at line 63 of file Xsec_final.C.

References make_syst_table_plots::ibin.

63  {
64  double binctr=Xnue->GetBinCenter(ibin);
65  double xnue=(spline_nue->Eval(binctr)/12.)*10.0;
66  Xnue->SetBinContent(ibin,xnue);
67 
68  x_genie[ibin-1]=binctr;
69  y_genie[ibin-1]=xnue;
70 
71  }
TH1F * Xnue
Definition: Xsec_final.C:59
TSpline3 * spline_nue
Definition: Xsec_final.C:58
double y_genie[100]
Definition: Xsec_final.C:61
double x_genie[100]
Definition: Xsec_final.C:60
Definition: FillPIDs.h:16
h2f GetXaxis ( ) -> CenterTitle()

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

84 {#nu} (GeV)");
h2f GetYaxis ( ) -> CenterTitle()
gStyle SetEndErrorSize ( )
band SetFillColor ( )
leg SetFillColor ( )
gargamelle SetLineColor ( kViolet 1)
t2k SetLineColor ( )
nova SetLineColor ( )
band SetLineColor ( )
genie SetLineColor ( )
leg SetLineColor ( )
genie SetLineStyle ( )
gargamelle SetLineWidth ( )
nova SetLineWidth ( )
gargamelle SetMarkerColor ( kViolet 1)
t2k SetMarkerColor ( )
nova_stat SetMarkerSize ( 1.  2)
gargamelle SetMarkerStyle ( 21  )
t2k SetMarkerStyle ( 24  )
nova_stat SetMarkerStyle ( )
prelim SetNDC ( )
gStyle SetOptStat ( )
prelim SetTextAlign ( 32  )
prelim SetTextColor ( kBlue 1)

Referenced by Save1DContour().

prelim SetTextSize ( 2/  30.)
gargamelle SetTitle ( "Gargamelle #nu_{e} CC data )
h2f SetXTitle ( "Electron Neutrino Energy (GeV)"  )
h2f SetYTitle ( "#nu_{e} CC #it{#sigma} (#times10^{-39} cm^{2}/nucleon)"  )

Variable Documentation

TGraphAsymmErrors* band =new TGraphAsymmErrors(49,x_band,y_band,exl_band,exh_band,eyl_band,eyh_band)

Definition at line 100 of file Xsec_final.C.

Referenced by sensitivity_plot().

TCanvas* c1 =new TCanvas("Xsec_Enue_withNOvA","",700,600)

Definition at line 124 of file Xsec_final.C.

double ex_nova[4] ={0.25,0.25,0.25,0.25}

Definition at line 39 of file Xsec_final.C.

double exh[6] = {1.08768,0.98571,0.98571,0.91773,1.52993,3.94362}

Definition at line 7 of file Xsec_final.C.

double exh_band[49] ={0.}

Definition at line 80 of file Xsec_final.C.

double exh_t2k[1] ={1.782}

Definition at line 26 of file Xsec_final.C.

double exl[6] = {0.917728,0.98571,0.91773,0.95133,1.46119,2.10699}

Definition at line 6 of file Xsec_final.C.

double exl_band[49] ={0.}

Definition at line 79 of file Xsec_final.C.

double exl_t2k[1] ={0.918}

Definition at line 25 of file Xsec_final.C.

double ey_nova[4] ={6.58498,5.99768,5.53289,6.66383}

Definition at line 40 of file Xsec_final.C.

double ey_nova_stat[4] ={4.88965,3.12006,2.64657,2.54683}

Definition at line 49 of file Xsec_final.C.

double eyh[6] = {3.576,3.4038,4.5971,5.4423,8.3394,6.4654}

Definition at line 9 of file Xsec_final.C.

double eyh_band[49] ={0.}

Definition at line 82 of file Xsec_final.C.

double eyh_t2k[1] ={2.1}

Definition at line 28 of file Xsec_final.C.

double eyl[6] = {2.2105,3.9125,4.5932,5.9587,8.5096,6.9779}

Definition at line 8 of file Xsec_final.C.

double eyl_band[49] ={0.}

Definition at line 81 of file Xsec_final.C.

double eyl_t2k[1] ={2.1}

Definition at line 27 of file Xsec_final.C.

TFile* fband =new TFile("genie_band.root")

Definition at line 74 of file Xsec_final.C.

TFile* fgenie =new TFile("xsec_genie.root")

Definition at line 56 of file Xsec_final.C.

TGraphAsymmErrors* gargamelle = new TGraphAsymmErrors(6,x,y,exl,exh,eyl,eyh)

Definition at line 11 of file Xsec_final.C.

TGraph* genie =new TGraph(100,x_genie,y_genie)
TGraph* gr_genie =(TGraph*)fgenie->Get("nu_e_C12/tot_cc")

Definition at line 57 of file Xsec_final.C.

TH2F* h2f =new TH2F("h2f","",10,0,10,10,0,59)

Definition at line 125 of file Xsec_final.C.

TH1F* hminus =(TH1F*)fband->Get("ccnue_1m_total")

Definition at line 76 of file Xsec_final.C.

Referenced by ppfx_smooth_weights_make(), SavePCAForCAFAna(), and test_beam_errorband().

TH1F* hplus =(TH1F*)fband->Get("ccnue_1p_total")

Definition at line 75 of file Xsec_final.C.

Referenced by SavePCAForCAFAna(), and test_beam_errorband().

TLegend* leg =new TLegend(0.6,0.6,0.89,0.89)

Definition at line 163 of file Xsec_final.C.

TGraphErrors* nova =new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova)

Definition at line 42 of file Xsec_final.C.

TGraphErrors* nova_stat =new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova_stat)

Definition at line 50 of file Xsec_final.C.

TLatex* prelim = new TLatex(.9, .95, "NO#nuA Preliminary")
TSpline3* spline_nue =new TSpline3("mysplinenue",gr_genie)

Definition at line 58 of file Xsec_final.C.

TGraphAsymmErrors* t2k =new TGraphAsymmErrors(1,x_t2k,y_t2k,exl_t2k,exh_t2k,eyl_t2k,eyh_t2k)

Definition at line 30 of file Xsec_final.C.

Referenced by nova_t2k_comparison().

double x_band[49] ={0.}

Definition at line 77 of file Xsec_final.C.

double x_genie[100] ={0.}

Definition at line 60 of file Xsec_final.C.

double x_nova[4] ={1.25,1.75,2.25,2.75}

Definition at line 37 of file Xsec_final.C.

Referenced by testbeam::BeamlineEventDisplay::DrawWCTrack().

double x_t2k[1] ={1.343}

Definition at line 23 of file Xsec_final.C.

TH1F* Xnue =new TH1F("Xnue","",100,0,10)

Definition at line 59 of file Xsec_final.C.

TH1F* Xnue2 =new TH1F("Xnue2","",50,0,10)

Definition at line 83 of file Xsec_final.C.

double y[6] = {13.53,17.8427,24.5361,33.1036,37.9249,36.5359}

Definition at line 5 of file Xsec_final.C.

double y_band[49] ={0.}

Definition at line 78 of file Xsec_final.C.

double y_genie[100] ={0.}

Definition at line 61 of file Xsec_final.C.

double y_nova[4] ={22.0921,28.8235,28.16,37.6431}

Definition at line 38 of file Xsec_final.C.

Referenced by testbeam::BeamlineEventDisplay::DrawWCTrack().

double y_t2k[1] ={11.1}

Definition at line 24 of file Xsec_final.C.