Functions | Variables
xf3_nu.C File Reference

Go to the source code of this file.

Functions

gStyle SetOptStat (0000)
 
ifstream base_neugen ("data/extract_sf.txt")
 
 for (int i=0;i< nfile;i++)
 
base_neugen close ()
 
ifstream lo_neugen ("data/lo_extract_sf.txt")
 
ifstream nutev ("data/final_nutev_xf3.txt", ios::in)
 
ifstream cdhsw ("data/final_cdhsw_xf3.txt", ios::in)
 
ifstream ccfr ("data/ccfr_f2xf3fe.txt", ios::in)
 
cout<<"MAKE THE PLOTS"<< endl;TCanvas c2("c2","", 0, 0, 1300, 900);c2-> Divide (2, 2)
 
c2 cd (1)
 
a1 SetAxisRange (0.01, 30,"y")
 
a1 Draw ()
 
a1 GetXaxis () -> SetTitle("Q2")
 
a1 GetYaxis () -> SetTitle("xF3")
 
gPad SetLogy ()
 
gPad SetLogx ()
 
gr_015 SetLineWidth (2)
 
gr_015 Draw ("C")
 
lo_gr_015 SetLineColor (2)
 
tex SetTextAlign (22)
 
tex SetTextSize (0.045)
 
tex DrawLatex (80, 8.0,"x = 0.015(20*xF3)")
 
tex DrawLatex (80, 2.22,"x = 0.045(5*xF3)")
 
tex DrawLatex (80, 0.520,"x = 0.080")
 
leg1 AddEntry (gr_015,"standard neugen xF3","l")
 
leg1 AddEntry (lo_gr_015,"+ no B-Y correction","l")
 
c2 cd (2)
 
tex DrawLatex (80, 9,"x = 0.125(5*xF3)")
 
c2 cd (3)
 
tex DrawLatex (80, 8,"x = 0.275(5*xF3)")
 
tex DrawLatex (80, 1.65161,"x = 0.350(2*xF3)")
 
c2 cd (4)
 
tex DrawLatex (55.9193, 1.914,"x = 0.550(5*xF3)")
 
tex DrawLatex (93.45, 0.312631,"x = 0.650(2*xF3)")
 
tex DrawLatex (95.565, 0.0456784,"x = 0.750")
 
c2 Print ("xf3_nu.eps")
 

Variables

const int nfile = 3672
 
Double_t row1_base [nfile]
 
Double_t row2_base [nfile]
 
Double_t row3_base [nfile]
 
Double_t row4_base [nfile]
 
Double_t row5_base [nfile]
 
Double_t lo_row1_base [nfile]
 
Double_t lo_row2_base [nfile]
 
Double_t lo_row3_base [nfile]
 
Double_t lo_row4_base [nfile]
 
Double_t lo_row5_base [nfile]
 
int nnutev = 75
 
Double_t nutev1 [nnutev]
 
Double_t nutev2 [nnutev]
 
Double_t nutev3 [nnutev]
 
Double_t nutev4 [nnutev]
 
int ncdhsw = 143
 
Double_t cdhsw1 [ncdhsw]
 
Double_t cdhsw2 [ncdhsw]
 
Double_t cdhsw3 [ncdhsw]
 
Double_t cdhsw4 [ncdhsw]
 
int nccfr = 116
 
Double_t ccfr1 [nccfr]
 
Double_t ccfr2 [nccfr]
 
Double_t ccfr3 [nccfr]
 
Double_t ccfr4 [nccfr]
 
Double_t ccfr5 [nccfr]
 
Double_t ccfr6 [nccfr]
 
Double_t ccfr7 [nccfr]
 
int num = 306
 
float x [12] = 0.015
 
Double_t q2 [12][num]
 
Double_t xf3 [12][num]
 
Int_t count [12]
 
Int_t num1 [12] ={20,5,1,5,2,1,5,2,1,5,2,1}
 
Double_t lo_q2 [12][num]
 
Double_t lo_xf3 [12][num]
 
Int_t lo_count [12]
 
Double_t nutev_q2 [12][num]
 
Double_t nutev_xf3 [12][num]
 
Double_t nutev_q2err [12][num]
 
Double_t nutev_err [12][num]
 
Int_t nutev_count [12]
 
Double_t cdhsw_q2 [12][num]
 
Double_t cdhsw_xf3 [12][num]
 
Double_t cdhsw_q2err [12][num]
 
Double_t cdhsw_err [12][num]
 
Int_t cdhsw_count [12]
 
Double_t ccfr_q2 [12][num]
 
Double_t ccfr_xf3 [12][num]
 
Double_t ccfr_q2err [12][num]
 
Double_t ccfr_err [12][num]
 
Int_t ccfr_count [12] = count_sp1
 
Int_t count_sp1 =0
 
Int_t count_sp2 =0
 
Double_t temp_q21 [num]
 
Double_t temp_sp1 [num]
 
Double_t temp_errsp1 [num]
 
Double_t temp_q22 [num]
 
Double_t temp_sp2 [num]
 
Double_t temp_errsp2 [num]
 
TH1F * a1 = new TH1F("a1","",10,0.005,2000)
 
TGraph * gr_015 = new TGraph(count[0],q2[0],xf3[0])
 
TGraph * lo_gr_015 = new TGraph(count[0],lo_q2[0],lo_xf3[0])
 
TGraph * gr_045 = new TGraph(count[1],q2[1],xf3[1])
 
TGraph * lo_gr_045 = new TGraph(count[0],lo_q2[1],lo_xf3[1])
 
TGraph * gr_080 = new TGraph(count[2],q2[2],xf3[2])
 
TGraph * lo_gr_080 = new TGraph(count[2],lo_q2[2],lo_xf3[2])
 
TLatex * tex = new TLatex()
 
TLegend * leg1 = new TLegend(0.356,0.262,0.81668,0.541984)
 
TH1F * a2 = new TH1F("a2","",10,0.005,2000)
 
TGraph * gr_125 = new TGraph(count[3],q2[3],xf3[3])
 
TGraph * lo_gr_125 = new TGraph(count[3],lo_q2[3],lo_xf3[3])
 
TGraph * gr_175 = new TGraph(count[4],q2[4],xf3[4])
 
TGraph * lo_gr_175 = new TGraph(count[4],lo_q2[4],lo_xf3[4])
 
TGraph * gr_225 = new TGraph(count[5],q2[5],xf3[5])
 
TGraph * lo_gr_225 = new TGraph(count[5],lo_q2[5],lo_xf3[5])
 
TH1F * a3 = new TH1F("a3","",10,0.005,2000)
 
TGraph * gr_275 = new TGraph(count[6],q2[6],xf3[6])
 
TGraph * lo_gr_275 = new TGraph(count[6],lo_q2[6],lo_xf3[6])
 
TGraph * gr_350 = new TGraph(count[7],q2[7],xf3[7])
 
TGraph * lo_gr_350 = new TGraph(count[7],lo_q2[7],lo_xf3[7])
 
TGraph * gr_450 = new TGraph(count[8],q2[8],xf3[8])
 
TGraph * lo_gr_450 = new TGraph(count[8],lo_q2[8],lo_xf3[8])
 
TH1F * a4 = new TH1F("a4","",10,0.005,2000)
 
TGraph * gr_550 = new TGraph(count[9],q2[9],xf3[9])
 
TGraph * lo_gr_550 = new TGraph(count[9],lo_q2[9],lo_xf3[9])
 
TGraph * gr_650 = new TGraph(count[10],q2[10],xf3[10])
 
TGraph * lo_gr_650 = new TGraph(count[10],lo_q2[10],lo_xf3[10])
 
TGraph * gr_750 = new TGraph(count[11],q2[11],xf3[11])
 
TGraph * lo_gr_750 = new TGraph(count[11],lo_q2[11],lo_xf3[11])
 

Function Documentation

leg1 AddEntry ( gr_015  ,
"standard neugen xF3"  ,
"l"   
)
leg1 AddEntry ( lo_gr_015  ,
"+ no B-Y correction"  ,
"l"   
)
ifstream base_neugen ( "data/extract_sf.txt"  )

Referenced by for().

ifstream ccfr ( "data/ccfr_f2xf3fe.txt"  ,
ios::in   
)
c2 cd ( )
c2 cd ( )
c2 cd ( )
c2 cd ( )
ifstream cdhsw ( "data/final_cdhsw_xf3.txt"  ,
ios::in   
)
base_neugen close ( )
cout<<"MAKE THE PLOTS"<<endl;TCanvas c2("c2","",0,0,1300,900);c2-> Divide ( ,
 
)
a1 Draw ( )
gr_015 Draw ( "C"  )
tex DrawLatex ( 80  ,
8.  0 
)
tex DrawLatex ( 80  ,
2.  22 
)
tex DrawLatex ( 80  ,
0.  520 
)
tex DrawLatex ( 80  ,
 
)
tex DrawLatex ( 80  ,
 
)
tex DrawLatex ( 80  ,
1.  65161 
)
tex DrawLatex ( 55.  9193,
1.  914 
)
tex DrawLatex ( 93.  45,
0.  312631 
)
tex DrawLatex ( 95.  565,
0.  0456784 
)
for ( )

Definition at line 26 of file xf3_nu.C.

References base_neugen(), MECModelEnuComparisons::i, and lo_neugen().

26  {
28  //cout<<" "<<row1_base[i]<<" "<<row2_base[i]<<" "<<row3_base[i]<<" "<<row4_base[i]<<" "<<row5_base[i]<<endl;
29 }
ifstream base_neugen("data/extract_sf.txt")
Double_t row3_base[nfile]
Definition: xf3_nu.C:22
Double_t row1_base[nfile]
Definition: xf3_nu.C:20
Double_t row5_base[nfile]
Definition: xf3_nu.C:24
Double_t row2_base[nfile]
Definition: xf3_nu.C:21
Double_t row4_base[nfile]
Definition: xf3_nu.C:23
a1 GetXaxis ( ) -> SetTitle("Q2")

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

84 {#nu} (GeV)");
a1 GetYaxis ( ) -> SetTitle("xF3")
ifstream lo_neugen ( "data/lo_extract_sf.txt"  )

Referenced by for().

ifstream nutev ( "data/final_nutev_xf3.txt"  ,
ios::in   
)
c2 Print ( "xf3_nu.eps"  )
a4 SetAxisRange ( 0.  01,
30  ,
"y"   
)
lo_gr_015 SetLineColor ( )
gr_015 SetLineWidth ( )
gPad SetLogx ( )
gPad SetLogy ( )
gStyle SetOptStat ( 0000  )
tex SetTextAlign ( 22  )
tex SetTextSize ( 0.  045)

Variable Documentation

TH1F* a1 = new TH1F("a1","",10,0.005,2000)

Definition at line 426 of file xf3_nu.C.

TH1F* a2 = new TH1F("a2","",10,0.005,2000)

Definition at line 477 of file xf3_nu.C.

TH1F* a3 = new TH1F("a3","",10,0.005,2000)

Definition at line 517 of file xf3_nu.C.

TH1F* a4 = new TH1F("a4","",10,0.005,2000)

Definition at line 557 of file xf3_nu.C.

Double_t ccfr1[nccfr]

Definition at line 81 of file xf3_nu.C.

Double_t ccfr2[nccfr]

Definition at line 82 of file xf3_nu.C.

Double_t ccfr3[nccfr]

Definition at line 83 of file xf3_nu.C.

Double_t ccfr4[nccfr]

Definition at line 84 of file xf3_nu.C.

Double_t ccfr5[nccfr]

Definition at line 85 of file xf3_nu.C.

Double_t ccfr6[nccfr]

Definition at line 86 of file xf3_nu.C.

Double_t ccfr7[nccfr]

Definition at line 87 of file xf3_nu.C.

ccfr_count[3] = count_sp1

Definition at line 214 of file xf3_nu.C.

Double_t ccfr_err[12][num]

Definition at line 213 of file xf3_nu.C.

Double_t ccfr_q2[12][num]

Definition at line 210 of file xf3_nu.C.

Double_t ccfr_q2err[12][num]

Definition at line 212 of file xf3_nu.C.

Double_t ccfr_xf3[12][num]

Definition at line 211 of file xf3_nu.C.

Double_t cdhsw1[ncdhsw]

Definition at line 66 of file xf3_nu.C.

Double_t cdhsw2[ncdhsw]

Definition at line 67 of file xf3_nu.C.

Double_t cdhsw3[ncdhsw]

Definition at line 68 of file xf3_nu.C.

Double_t cdhsw4[ncdhsw]

Definition at line 69 of file xf3_nu.C.

Int_t cdhsw_count[12]

Definition at line 189 of file xf3_nu.C.

Double_t cdhsw_err[12][num]

Definition at line 188 of file xf3_nu.C.

Double_t cdhsw_q2[12][num]

Definition at line 185 of file xf3_nu.C.

Double_t cdhsw_q2err[12][num]

Definition at line 187 of file xf3_nu.C.

Double_t cdhsw_xf3[12][num]

Definition at line 186 of file xf3_nu.C.

Int_t count[12]

Definition at line 115 of file xf3_nu.C.

count_sp1 =0

Definition at line 233 of file xf3_nu.C.

count_sp2 =0

Definition at line 234 of file xf3_nu.C.

TGraph* gr_015 = new TGraph(count[0],q2[0],xf3[0])

Definition at line 434 of file xf3_nu.C.

TGraph* gr_045 = new TGraph(count[1],q2[1],xf3[1])

Definition at line 443 of file xf3_nu.C.

TGraph* gr_080 = new TGraph(count[2],q2[2],xf3[2])

Definition at line 452 of file xf3_nu.C.

TGraph* gr_125 = new TGraph(count[3],q2[3],xf3[3])

Definition at line 485 of file xf3_nu.C.

TGraph* gr_175 = new TGraph(count[4],q2[4],xf3[4])

Definition at line 494 of file xf3_nu.C.

TGraph* gr_225 = new TGraph(count[5],q2[5],xf3[5])

Definition at line 503 of file xf3_nu.C.

TGraph* gr_275 = new TGraph(count[6],q2[6],xf3[6])

Definition at line 525 of file xf3_nu.C.

TGraph* gr_350 = new TGraph(count[7],q2[7],xf3[7])

Definition at line 534 of file xf3_nu.C.

TGraph* gr_450 = new TGraph(count[8],q2[8],xf3[8])

Definition at line 543 of file xf3_nu.C.

TGraph* gr_550 = new TGraph(count[9],q2[9],xf3[9])

Definition at line 565 of file xf3_nu.C.

TGraph* gr_650 = new TGraph(count[10],q2[10],xf3[10])

Definition at line 575 of file xf3_nu.C.

TGraph* gr_750 = new TGraph(count[11],q2[11],xf3[11])

Definition at line 585 of file xf3_nu.C.

TLegend* leg1 = new TLegend(0.356,0.262,0.81668,0.541984)

Definition at line 469 of file xf3_nu.C.

Int_t lo_count[12]

Definition at line 141 of file xf3_nu.C.

TGraph* lo_gr_015 = new TGraph(count[0],lo_q2[0],lo_xf3[0])

Definition at line 438 of file xf3_nu.C.

TGraph* lo_gr_045 = new TGraph(count[0],lo_q2[1],lo_xf3[1])

Definition at line 447 of file xf3_nu.C.

TGraph* lo_gr_080 = new TGraph(count[2],lo_q2[2],lo_xf3[2])

Definition at line 456 of file xf3_nu.C.

TGraph* lo_gr_125 = new TGraph(count[3],lo_q2[3],lo_xf3[3])

Definition at line 489 of file xf3_nu.C.

TGraph* lo_gr_175 = new TGraph(count[4],lo_q2[4],lo_xf3[4])

Definition at line 498 of file xf3_nu.C.

TGraph* lo_gr_225 = new TGraph(count[5],lo_q2[5],lo_xf3[5])

Definition at line 507 of file xf3_nu.C.

TGraph* lo_gr_275 = new TGraph(count[6],lo_q2[6],lo_xf3[6])

Definition at line 529 of file xf3_nu.C.

TGraph* lo_gr_350 = new TGraph(count[7],lo_q2[7],lo_xf3[7])

Definition at line 538 of file xf3_nu.C.

TGraph* lo_gr_450 = new TGraph(count[8],lo_q2[8],lo_xf3[8])

Definition at line 547 of file xf3_nu.C.

TGraph* lo_gr_550 = new TGraph(count[9],lo_q2[9],lo_xf3[9])

Definition at line 569 of file xf3_nu.C.

TGraph* lo_gr_650 = new TGraph(count[10],lo_q2[10],lo_xf3[10])

Definition at line 579 of file xf3_nu.C.

TGraph* lo_gr_750 = new TGraph(count[11],lo_q2[11],lo_xf3[11])

Definition at line 589 of file xf3_nu.C.

Double_t lo_q2[12][num]

Definition at line 138 of file xf3_nu.C.

Double_t lo_row1_base[nfile]

Definition at line 35 of file xf3_nu.C.

Double_t lo_row2_base[nfile]

Definition at line 36 of file xf3_nu.C.

Double_t lo_row3_base[nfile]

Definition at line 37 of file xf3_nu.C.

Double_t lo_row4_base[nfile]

Definition at line 38 of file xf3_nu.C.

Double_t lo_row5_base[nfile]

Definition at line 39 of file xf3_nu.C.

Double_t lo_xf3[12][num]

Definition at line 139 of file xf3_nu.C.

int nccfr = 116

Definition at line 80 of file xf3_nu.C.

int ncdhsw = 143

Definition at line 65 of file xf3_nu.C.

const int nfile = 3672

Definition at line 16 of file xf3_nu.C.

int nnutev = 75

Definition at line 50 of file xf3_nu.C.

int num = 306

Definition at line 96 of file xf3_nu.C.

Int_t num1[12] ={20,5,1,5,2,1,5,2,1,5,2,1}

Definition at line 116 of file xf3_nu.C.

Double_t nutev1[nnutev]

Definition at line 51 of file xf3_nu.C.

Double_t nutev2[nnutev]

Definition at line 52 of file xf3_nu.C.

Double_t nutev3[nnutev]

Definition at line 53 of file xf3_nu.C.

Double_t nutev4[nnutev]

Definition at line 54 of file xf3_nu.C.

Int_t nutev_count[12]

Definition at line 164 of file xf3_nu.C.

Double_t nutev_err[12][num]

Definition at line 163 of file xf3_nu.C.

Double_t nutev_q2[12][num]

Definition at line 160 of file xf3_nu.C.

Double_t nutev_q2err[12][num]

Definition at line 162 of file xf3_nu.C.

Double_t nutev_xf3[12][num]

Definition at line 161 of file xf3_nu.C.

Double_t q2[12][num]

Definition at line 112 of file xf3_nu.C.

Double_t row1_base[nfile]

Definition at line 20 of file xf3_nu.C.

Double_t row2_base[nfile]

Definition at line 21 of file xf3_nu.C.

Double_t row3_base[nfile]

Definition at line 22 of file xf3_nu.C.

Double_t row4_base[nfile]

Definition at line 23 of file xf3_nu.C.

Double_t row5_base[nfile]

Definition at line 24 of file xf3_nu.C.

Double_t temp_errsp1[num]

Definition at line 238 of file xf3_nu.C.

Double_t temp_errsp2[num]

Definition at line 242 of file xf3_nu.C.

Double_t temp_q21[num]

Definition at line 236 of file xf3_nu.C.

Double_t temp_q22[num]

Definition at line 240 of file xf3_nu.C.

Double_t temp_sp1[num]

Definition at line 237 of file xf3_nu.C.

Double_t temp_sp2[num]

Definition at line 241 of file xf3_nu.C.

TLatex* tex = new TLatex()

Definition at line 461 of file xf3_nu.C.

x[11] = 0.015

Definition at line 97 of file xf3_nu.C.

Double_t xf3[12][num]

Definition at line 113 of file xf3_nu.C.