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.

                                                     {
    double binctr=Xnue->GetBinCenter(ibin);
    double xnue=(spline_nue->Eval(binctr)/12.)*10.0;
    Xnue->SetBinContent(ibin,xnue);

    x_genie[ibin-1]=binctr;
    y_genie[ibin-1]=xnue;
    
   }

h2f GetXaxis

(

)

-> CenterTitle()

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

{#nu} (GeV)");

h2f GetYaxis

(

)

-> CenterTitle()

gStyle SetEndErrorSize

(

4

)

band SetFillColor

(

3

)

leg SetFillColor

(

0

)

gargamelle SetLineColor

(

kViolet-

1

)

t2k SetLineColor

(

4

)

nova SetLineColor

(

2

)

band SetLineColor

(

3

)

genie SetLineColor

(

1

)

leg SetLineColor

(

0

)

genie SetLineStyle

(

2

)

gargamelle SetLineWidth

(

2

)

nova SetLineWidth

(

3

)

gargamelle SetMarkerColor

(

kViolet-

1

)

t2k SetMarkerColor

(

4

)

nova_stat SetMarkerSize

(

1.

2

)

gargamelle SetMarkerStyle

(

21

)

t2k SetMarkerStyle

(

24

)

nova_stat SetMarkerStyle

(

8

)

prelim SetNDC

(

)

gStyle SetOptStat

(

0

)

prelim SetTextAlign

(

32

)

prelim SetTextColor

(

kBlue+

1

)

Referenced by Save1DContour().

prelim SetTextSize

(

2/

30.

)

Referenced by canMan::drawOverlay(), MakeSplitCans(), plot_uncertainty(), and PlotDeltaSigmaWithSigma().

gargamelle SetTitle

(

"Gargamelle #nu_{e} CC data"

)

Referenced by plot_validation_datamc.Spectrum::Draw(), generic.tools.ImageMaker.ImageMaker::draw(), plot_validation_datamc.Ratio::Draw(), plot_validation_datamc_2018.Spectrum::Draw(), plot_validation_datamc_2018.Ratio::Draw(), generic.controllers.MakeCanvases.HistogramExtractor::ExtractHistsFromDirectory(), GetHistVectors(), Canvas.Canvas::label_axes(), generic.tools.ImageMaker.ImageMaker::label_axes(), main(), make_muonid_opt(), MakeCanvasForData(), MakeCanvasForSim(), MakeNuEnergyPlots(), MakeOverlayCanvas(), MakePlots(), make_reco_miniprod5_plots::MakePlots(), make_reco_plots::MakePlots(), validation::mcvalidate(), plot_quantile_boundaries_2020(), PlotCovariance(), PlotGenRatios(), PlotPreSelectionPlots(), PlotSelectionPlots(), make_mec_shifts_plots::PrepHist(), make_xsec_tuning_stage_plots::PrepHist(), validation::recovalidate(), and make_true_q0q3_plots::SetTitleQ0().

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)

Definition at line 116 of file Xsec_final.C.

Referenced by genie_syst_pca(), ana::TrivialCrossSectionAnalysis::operator=(), and ana::TrivialCrossSectionAnalysis::TrivialCrossSectionAnalysis().

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")

Definition at line 133 of file Xsec_final.C.

Referenced by Beam(), BiCountCanvas(), BiProbCanvas(), FakeData(), myPreliminary(), NDSimulation(), NDSimulation1(), pidPlotProng(), Prelim(), Preliminary(), preliminary(), Preliminary2(), PreliminarySide(), PrintCaption(), PrintChiSq(), PrintChiSq2(), PrintPreliminary(), PrintSimulation(), Simulation(), simulation(), numu::Simulation(), and SimulationSide().

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.