Functions | Variables
plot_lowE_sig_spread_genie.C File Reference

Go to the source code of this file.

Functions

TFile fnoavg ("./ccqe_noavg_defkF.root","read")
 
TFile favgkFlo ("./ccqe_avg_kF200.root","read")
 
TFile favgkFhi ("./ccqe_avg_kF250.root","read")
 
TFile favgMalo ("./ccqe_avg_Mam15pct.root","read")
 
TFile favgMahi ("./ccqe_avg_Map15pct.root","read")
 
xsavg SetLineWidth (3)
 
xsavgkFlo SetLineColor (kGreen)
 
xsavgkFhi SetLineColor (kRed)
 
xsavg Draw ("L")
 
xsnoavg Draw ("LSAME")
 
 for (int i=0;i< nsteps;i++)
 
hframe2 GetXaxis () -> SetTitle("E_
 
hframe2 GetYaxis () -> SetTitle("#sigma spread %")
 
hframe2 Draw ()
 
grerrp_pF SetLineStyle (kDotted)
 
grerrp_kF SetLineStyle (kDashed)
 
legend SetFillColor (0)
 
legend SetBorderSize (0)
 
legend AddEntry (grerrp,"Combined error for #nu_{#mu} O^{16} CCQE","L")
 
legend AddEntry (grerrp_pF,"Effect of E_{#nu} smearing from the nucleon momentum distribution","L")
 
legend AddEntry (grerrp_kF,"k_{F} #pm 25 MeV (~10%) (effect of modifying Pauli-blocking term)","L")
 
legend AddEntry (grerrp_Ma,"M_{A} #pm 15%","L")
 

Variables

TDirectory * davg = (TDirectory *) favg -> Get("nu_mu_O16")
 
TDirectory * dnoavg = (TDirectory *) fnoavg -> Get("nu_mu_O16")
 
TDirectory * davgkFlo = (TDirectory *) favgkFlo -> Get("nu_mu_O16")
 
TDirectory * davgkFhi = (TDirectory *) favgkFhi -> Get("nu_mu_O16")
 
TDirectory * davgMalo = (TDirectory *) favgMalo -> Get("nu_mu_O16")
 
TDirectory * davgMahi = (TDirectory *) favgMahi -> Get("nu_mu_O16")
 
TGraph * xsavg = (TGraph*) davg -> Get("qel_cc_n")
 
TGraph * xsnoavg = (TGraph*) dnoavg -> Get("qel_cc_n")
 
TGraph * xsavgkFlo = (TGraph*) davgkFlo -> Get("qel_cc_n")
 
TGraph * xsavgkFhi = (TGraph*) davgkFhi -> Get("qel_cc_n")
 
TGraph * xsavgMalo = (TGraph*) davgMalo -> Get("qel_cc_n")
 
TGraph * xsavgMahi = (TGraph*) davgMahi -> Get("qel_cc_n")
 
TCanvas * c1 = new TCanvas()
 
TH1F * hframe1 = (TH1F*)c1->DrawFrame(0,0,1.8,10)
 
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_pF [nsteps]
 
double errn_pF [nsteps]
 
double errp_kF [nsteps]
 
double errn_kF [nsteps]
 
double errp_Ma [nsteps]
 
double errn_Ma [nsteps]
 
double enu [nsteps]
 
TGraph * grerrp_pF = new TGraph(nsteps, enu, errp_pF)
 
TGraph * grerrn_pF = new TGraph(nsteps, enu, errn_pF)
 
TGraph * grerrp_kF = new TGraph(nsteps, enu, errp_kF)
 
TGraph * grerrn_kF = new TGraph(nsteps, enu, errn_kF)
 
TGraph * grerrp_Ma = new TGraph(nsteps, enu, errp_Ma)
 
TGraph * grerrn_Ma = new TGraph(nsteps, enu, errn_Ma)
 
TGraph * grerrp = new TGraph(nsteps, enu, errp)
 
TGraph * grerrn = new TGraph(nsteps, enu, errn)
 
TCanvas * c2 = new TCanvas()
 
TH1F * hframe2 = (TH1F*)c2->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  ,
"Combined error for #nu_{#mu} O^{16} CCQE"  ,
"L"   
)
legend AddEntry ( grerrp_pF  ,
"Effect of E_{#nu} smearing from the nucleon momentum distribution"  ,
"L"   
)
legend AddEntry ( grerrp_kF  ,
"k_{F} #pm 25 MeV (~10%) (effect of modifying Pauli-blocking term)"  ,
"L"   
)
legend AddEntry ( grerrp_Ma  ,
"M_{A} #pm 15%"  ,
"L"   
)
xsavg Draw ( "L"  )
xsavgkFhi Draw ( "LSAME"  )
hframe2 Draw ( )
TFile favgkFhi ( "./ccqe_avg_kF250.root"  ,
"read"   
)
TFile favgkFlo ( "./ccqe_avg_kF200.root"  ,
"read"   
)
TFile favgMahi ( "./ccqe_avg_Map15pct.root"  ,
"read"   
)
TFile favgMalo ( "./ccqe_avg_Mam15pct.root"  ,
"read"   
)
TFile fnoavg ( "./ccqe_noavg_defkF.root"  ,
"read"   
)
for ( )

Definition at line 50 of file plot_lowE_sig_spread_genie.C.

References MECModelEnuComparisons::i, ana::Sqrt(), step, and submit_syst::x.

50  {
51 
52  double x = emin + i*step;
53  double yavg = xsavg->Eval(x);
54  double ynoavg = xsnoavg->Eval(x);
55  double yavgnoavgc = 0.5*(yavg+ynoavg);
56  double yavgkFlo = xsavgkFlo->Eval(x);
57  double yavgkFhi = xsavgkFhi->Eval(x);
58  double yavgMalo = xsavgMalo->Eval(x);
59  double yavgMahi = xsavgMahi->Eval(x);
60 
61  enu [i] = x;
62 
63  errp_pF[i] = 100*( TMath::Max(yavg,ynoavg) - yavgnoavgc)/yavgnoavgc;
64  errn_pF[i] = 100*( TMath::Min(yavg,ynoavg) - yavgnoavgc)/yavgnoavgc;
65  errp_kF[i] = 100*( TMath::Max(yavgkFhi,yavgkFlo) - yavg)/yavg;
66  errn_kF[i] = 100*( TMath::Min(yavgkFhi,yavgkFlo) - yavg)/yavg;
67  errp_Ma[i] = 100*( TMath::Max(yavgMahi,yavgMalo) - yavg)/yavg;
68  errn_Ma[i] = 100*( TMath::Min(yavgMahi,yavgMalo) - yavg)/yavg;
71 }
TGraph * xsavg
double enu[nsteps]
double errp_Ma[nsteps]
TGraph * xsavgkFhi
const double emin
double errp_pF[nsteps]
double errn[nsteps]
const double step
TGraph * xsnoavg
Var Sqrt(const Var &v)
Use to take sqrt of a var.
Definition: Var.cxx:326
double errp_kF[nsteps]
TGraph * xsavgMahi
TGraph * xsavgMalo
TGraph * xsavgkFlo
double errp[nsteps]
double errn_kF[nsteps]
double errn_Ma[nsteps]
double errn_pF[nsteps]
hframe2 GetXaxis ( ) -> SetTitle("E_

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

84 {#nu} (GeV)");
hframe2 GetYaxis ( ) -> SetTitle("#sigma spread %")
legend SetBorderSize ( )
legend SetFillColor ( )
xsavgkFlo SetLineColor ( kGreen  )
xsavgkFhi SetLineColor ( kRed  )
grerrn_pF SetLineStyle ( kDotted  )
grerrp_kF SetLineStyle ( kDashed  )
xsavg SetLineWidth ( )

Variable Documentation

TCanvas* c1 = new TCanvas()

Definition at line 23 of file plot_lowE_sig_spread_genie.C.

TCanvas* c2 = new TCanvas()

Definition at line 82 of file plot_lowE_sig_spread_genie.C.

TDirectory* davg = (TDirectory *) favg -> Get("nu_mu_O16")

Definition at line 9 of file plot_lowE_sig_spread_genie.C.

TDirectory* davgkFhi = (TDirectory *) favgkFhi -> Get("nu_mu_O16")

Definition at line 12 of file plot_lowE_sig_spread_genie.C.

TDirectory* davgkFlo = (TDirectory *) favgkFlo -> Get("nu_mu_O16")

Definition at line 11 of file plot_lowE_sig_spread_genie.C.

TDirectory* davgMahi = (TDirectory *) favgMahi -> Get("nu_mu_O16")

Definition at line 14 of file plot_lowE_sig_spread_genie.C.

TDirectory* davgMalo = (TDirectory *) favgMalo -> Get("nu_mu_O16")

Definition at line 13 of file plot_lowE_sig_spread_genie.C.

TDirectory* dnoavg = (TDirectory *) fnoavg -> Get("nu_mu_O16")

Definition at line 10 of file plot_lowE_sig_spread_genie.C.

const double emax = 0.70
const double emin = 0.20
double enu[nsteps]

Definition at line 48 of file plot_lowE_sig_spread_genie.C.

double errn[nsteps]

Definition at line 41 of file plot_lowE_sig_spread_genie.C.

double errn_kF[nsteps]

Definition at line 45 of file plot_lowE_sig_spread_genie.C.

double errn_Ma[nsteps]

Definition at line 47 of file plot_lowE_sig_spread_genie.C.

double errn_pF[nsteps]

Definition at line 43 of file plot_lowE_sig_spread_genie.C.

double errp[nsteps]

Definition at line 40 of file plot_lowE_sig_spread_genie.C.

double errp_kF[nsteps]

Definition at line 44 of file plot_lowE_sig_spread_genie.C.

double errp_Ma[nsteps]

Definition at line 46 of file plot_lowE_sig_spread_genie.C.

double errp_pF[nsteps]

Definition at line 42 of file plot_lowE_sig_spread_genie.C.

GeV

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

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

Definition at line 80 of file plot_lowE_sig_spread_genie.C.

TGraph* grerrn_kF = new TGraph(nsteps, enu, errn_kF)

Definition at line 76 of file plot_lowE_sig_spread_genie.C.

TGraph* grerrn_Ma = new TGraph(nsteps, enu, errn_Ma)

Definition at line 78 of file plot_lowE_sig_spread_genie.C.

TGraph* grerrn_pF = new TGraph(nsteps, enu, errn_pF)

Definition at line 74 of file plot_lowE_sig_spread_genie.C.

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

Definition at line 79 of file plot_lowE_sig_spread_genie.C.

TGraph* grerrp_kF = new TGraph(nsteps, enu, errp_kF)

Definition at line 75 of file plot_lowE_sig_spread_genie.C.

TGraph* grerrp_Ma = new TGraph(nsteps, enu, errp_Ma)

Definition at line 77 of file plot_lowE_sig_spread_genie.C.

TGraph* grerrp_pF = new TGraph(nsteps, enu, errp_pF)

Definition at line 73 of file plot_lowE_sig_spread_genie.C.

TH1F* hframe1 = (TH1F*)c1->DrawFrame(0,0,1.8,10)

Definition at line 24 of file plot_lowE_sig_spread_genie.C.

TH1F* hframe2 = (TH1F*)c2->DrawFrame(0.2,-50,0.7,50)

Definition at line 83 of file plot_lowE_sig_spread_genie.C.

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

Definition at line 104 of file plot_lowE_sig_spread_genie.C.

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

Definition at line 38 of file plot_lowE_sig_spread_genie.C.

Referenced by for().

TGraph* xsavg = (TGraph*) davg -> Get("qel_cc_n")

Definition at line 16 of file plot_lowE_sig_spread_genie.C.

TGraph* xsavgkFhi = (TGraph*) davgkFhi -> Get("qel_cc_n")

Definition at line 19 of file plot_lowE_sig_spread_genie.C.

TGraph* xsavgkFlo = (TGraph*) davgkFlo -> Get("qel_cc_n")

Definition at line 18 of file plot_lowE_sig_spread_genie.C.

TGraph* xsavgMahi = (TGraph*) davgMahi -> Get("qel_cc_n")

Definition at line 21 of file plot_lowE_sig_spread_genie.C.

TGraph* xsavgMalo = (TGraph*) davgMalo -> Get("qel_cc_n")

Definition at line 20 of file plot_lowE_sig_spread_genie.C.

TGraph* xsnoavg = (TGraph*) dnoavg -> Get("qel_cc_n")

Definition at line 17 of file plot_lowE_sig_spread_genie.C.