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  }
Xnue
TH1F * Xnue
Definition: Xsec_final.C:59
make_syst_table_plots.ibin
ibin
Definition: make_syst_table_plots.py:246
spline_nue
TSpline3 * spline_nue
Definition: Xsec_final.C:58
y_genie
double y_genie[100]
Definition: Xsec_final.C:61
x_genie
double x_genie[100]
Definition: Xsec_final.C:60
xnue
Definition: FillPIDs.h:16
h2f GetXaxis ( ) -> CenterTitle()

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

84 {#nu} (GeV)");
GeV
GeV
Definition: plot_lowE_sig_spread_genie.C:84
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.

Generated on Thu Dec 3 2020 05:15:05 for NOvA by   doxygen 1.8.11