Functions | Variables
plot_lowE_sig_spread_martinez10102329v1.C File Reference

Go to the source code of this file.

Functions

TSpline3 lo ("lo", xlo, ylo, 12,"0")
 
TSpline3 hi ("hi", xhi, yhi, 18,"0")
 
TSpline3 lo_no2 ("lo", xlo_no2p2h, ylo_no2p2h, 17,"0")
 
TSpline3 hi_no2 ("hi", xhi_no2p2h, yhi_no2p2h, 20,"0")
 
 for (int i=0;i< nsteps;i++)
 
grerrp SetLineStyle (kDashed)
 
hframe GetXaxis () -> SetTitle("E_
 
hframe GetYaxis () -> SetTitle("#sigma spread %")
 
hframe Draw ()
 
grerrp Draw ("L")
 
legend SetFillColor (0)
 
legend SetBorderSize (0)
 
legend AddEntry (grerrp,"Model spread (from Fig.1 in Martini arXiv:1010.2329v1","L")
 
legend AddEntry (grerrp_no2,"Model spread (as above but excluding RPA-2p2h)","L")
 

Variables

double ylo [12]
 
double xhi [18]
 
double yhi [18]
 
double xlo_no2p2h [17]
 
double ylo_no2p2h [17]
 
double xhi_no2p2h [20]
 
double yhi_no2p2h [20]
 
const int nsteps = 100
 
const double emin = 0.20
 
const double emax = 0.70
 
const double step = (emax-emin)/(nsteps-1)
 
double errp [nsteps]
 
double errn [nsteps]
 
double errp_no2 [nsteps]
 
double errn_no2 [nsteps]
 
double enu [nsteps]
 
TGraph * grerrp = new TGraph(nsteps, enu, errp)
 
TGraph * grerrn = new TGraph(nsteps, enu, errn)
 
TGraph * grerrp_no2 = new TGraph(nsteps, enu, errp_no2)
 
TGraph * grerrn_no2 = new TGraph(nsteps, enu, errn_no2)
 
TCanvas * c = new TCanvas()
 
TH1F * hframe = (TH1F*)c->DrawFrame(0.2,-50,0.7,50)
 
 GeV
 
TLegend * legend = new TLegend(0.6, 0.6, 0.9, 0.9)
 

Function Documentation

legend AddEntry ( grerrp  )
legend AddEntry ( grerrp_no2  ,
"Model spread (as above but excluding RPA-2p2h)"  ,
"L"   
)
hframe Draw ( )
grerrp Draw ( "L"  )
for ( )

Definition at line 175 of file plot_lowE_sig_spread_martinez10102329v1.C.

References MECModelEnuComparisons::i, step, and submit_syst::x.

175  {
176 
177  double x = emin + i*step;
178  double yh = hi->Eval(x);
179  double yl = lo->Eval(x);
180  double yc = 0.5*(yh+yl);
181  double yh_no2 = hi_no2->Eval(x);
182  double yl_no2 = lo_no2->Eval(x);
183  double yc_no2 = 0.5*(yh_no2+yl_no2);
184 
185  errp [i] = 100*(yh-yc)/yc;
186  errn [i] = 100*(yl-yc)/yc;
187  errp_no2[i] = 100*(yh_no2-yc_no2)/yc_no2;
188  errn_no2[i] = 100*(yl_no2-yc_no2)/yc_no2;
189  enu [i] = x;
190 }
TSpline3 lo("lo", xlo, ylo, 12,"0")
TSpline3 hi("hi", xhi, yhi, 18,"0")
TSpline3 lo_no2("lo", xlo_no2p2h, ylo_no2p2h, 17,"0")
TSpline3 hi_no2("hi", xhi_no2p2h, yhi_no2p2h, 20,"0")
hframe GetXaxis ( ) -> SetTitle("E_

Definition at line 202 of file plot_lowE_sig_spread_martinez10102329v1.C.

Referenced by fnex::BeamSyst::CalcBeamSystWeights(), fnex::CalibSyst::CalcCalibSystWeights(), fnex::GeniePCASyst::CalcGeniePCASystWeights(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), nueccinc_test::NueCCIncGlobalFitter::CheckSystsSize(), corrected_PE(), d4sigma_plot(), om::PlotClickHandler::DecodeClick(), ana::nueccinc::NueCCIncTemplateFitter::doFullFit(), evd::TQPad::Draw(), generic.tools.ImageMaker.ImageMaker::draw(), fnex::CorrectedSpectrum::DrawDataVsMC(), DrawFrame(), ana::CovMxSurface::DrawSurfacePoint(), energyResolution(), extract_hadnucleus_xsec(), FD_Data_PosComp(), FOMCalc(), ana::NumuCCIncPionTemplateFitter::getCorrelationMatrixTemplateBins(), ana::nueccinc::NueCCIncTemplateFitter::getCorrelationMatrixTemplateBins(), ana::NumuCCIncPionTemplateFitter::getCovarianceMatrixTemplateBins(), ana::nueccinc::NueCCIncTemplateFitter::getCovarianceMatrixTemplateBins(), GetHistVectors(), GetSpectToHist(), fnex::CorrectedSpectrum::GetStacksCanvasCopy(), main(), MakeCanvas(), MakeCanvasForData(), MakeCanvasForSim(), MakeCSV(), MakeNuEnergyPlots(), MakeOverlayCanvas(), MakePlots(), MakeSplitCans(), NewNearlinePlots(), nue_parameters_vs_fit_iteration(), fnex::ShifterAndWeighter::NueExtrapSig2017Weight(), nuint09_1pi3(), nuint09_coh3(), nuint09_qel4(), nuint09_qel5(), overlayFBplots(), SupernovaUtilities::ParameterisedDataToHistograms(), plot_datapredictions(), plot_predictions(), plot_quantile_boundaries_2020(), plot_rationoosc(), Plotting_DataAndPrediction(), PositionComparison(), RebinGenerators(), SetBin(), SetRange(), ana::NueExtrapSystSignalKin2017::Shift(), ana::NuWROSyst::Shift(), starPlot(), thresholdshadowing(), uncorrected_PE(), evd::THUnZoomable::UnZoom(), xsec_tot_uncert_optimization(), and xsec_uncertainty_per_bin().

hframe GetYaxis ( ) -> SetTitle("#sigma spread %")
TSpline3 hi ( "hi"  ,
xhi  ,
yhi  ,
18  ,
"0"   
)
TSpline3 hi_no2 ( "hi"  ,
xhi_no2p2h  ,
yhi_no2p2h  ,
20  ,
"0"   
)
TSpline3 lo ( "lo"  ,
xlo  ,
ylo  ,
12  ,
"0"   
)
TSpline3 lo_no2 ( "lo"  ,
xlo_no2p2h  ,
ylo_no2p2h  ,
17  ,
"0"   
)
legend SetBorderSize ( )
legend SetFillColor ( )
grerrp SetLineStyle ( kDashed  )

Variable Documentation

TCanvas* c = new TCanvas()

Definition at line 200 of file plot_lowE_sig_spread_martinez10102329v1.C.

const double emax = 0.70

Definition at line 166 of file plot_lowE_sig_spread_martinez10102329v1.C.

const double emin = 0.20

Definition at line 165 of file plot_lowE_sig_spread_martinez10102329v1.C.

double enu[nsteps]

Definition at line 173 of file plot_lowE_sig_spread_martinez10102329v1.C.

double errn[nsteps]

Definition at line 170 of file plot_lowE_sig_spread_martinez10102329v1.C.

double errn_no2[nsteps]

Definition at line 172 of file plot_lowE_sig_spread_martinez10102329v1.C.

double errp[nsteps]

Definition at line 169 of file plot_lowE_sig_spread_martinez10102329v1.C.

double errp_no2[nsteps]

Definition at line 171 of file plot_lowE_sig_spread_martinez10102329v1.C.

GeV

Definition at line 202 of file plot_lowE_sig_spread_martinez10102329v1.C.

TGraph* grerrn = new TGraph(nsteps, enu, errn)

Definition at line 193 of file plot_lowE_sig_spread_martinez10102329v1.C.

TGraph* grerrn_no2 = new TGraph(nsteps, enu, errn_no2)

Definition at line 195 of file plot_lowE_sig_spread_martinez10102329v1.C.

TGraph* grerrp = new TGraph(nsteps, enu, errp)

Definition at line 192 of file plot_lowE_sig_spread_martinez10102329v1.C.

TGraph* grerrp_no2 = new TGraph(nsteps, enu, errp_no2)

Definition at line 194 of file plot_lowE_sig_spread_martinez10102329v1.C.

TH1F* hframe = (TH1F*)c->DrawFrame(0.2,-50,0.7,50)

Definition at line 201 of file plot_lowE_sig_spread_martinez10102329v1.C.

Referenced by bwnorm(), main(), and plotTrainTest().

TLegend* legend = new TLegend(0.6, 0.6, 0.9, 0.9)

Definition at line 211 of file plot_lowE_sig_spread_martinez10102329v1.C.

const int nsteps = 100

Definition at line 164 of file plot_lowE_sig_spread_martinez10102329v1.C.

const double step = (emax-emin)/(nsteps-1)

Definition at line 167 of file plot_lowE_sig_spread_martinez10102329v1.C.

Referenced by for().

double xhi[18]
double xhi_no2p2h[20]
Initial value:
= {
0.14024,
0.17683,
0.19817,
0.21646,
0.24238,
0.26677,
0.29421,
0.32622,
0.35823,
0.38720,
0.42073,
0.45427,
0.49238,
0.52287,
0.55945,
0.60671,
0.64329,
0.68445,
0.71646,
0.74848,
}

Definition at line 114 of file plot_lowE_sig_spread_martinez10102329v1.C.

double xlo_no2p2h[17]
Initial value:
= {
0.15701,
0.18445,
0.20884,
0.24238,
0.28049,
0.31707,
0.35213,
0.39482,
0.42835,
0.46646,
0.51677,
0.56555,
0.60823,
0.65091,
0.69207,
0.72409,
0.74695,
}

Definition at line 76 of file plot_lowE_sig_spread_martinez10102329v1.C.

double yhi[18]
double yhi_no2p2h[20]
Initial value:
= {
0.08721,
0.46512,
0.75581,
1.19186,
1.65698,
2.09302,
2.67442,
3.28488,
3.72093,
4.21512,
4.65116,
5.11628,
5.52326,
5.78488,
6.13372,
6.45349,
6.68605,
6.88953,
7.03488,
7.12209,
}

Definition at line 136 of file plot_lowE_sig_spread_martinez10102329v1.C.

double ylo[12]
double ylo_no2p2h[17]
Initial value:
= {
0.02907,
0.26163,
0.55233,
1.04651,
1.62791,
2.20930,
2.70349,
3.28488,
3.77907,
4.15698,
4.70930,
5.17442,
5.49419,
5.72674,
5.95930,
6.10465,
6.25000,
}

Definition at line 95 of file plot_lowE_sig_spread_martinez10102329v1.C.