f2_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++)
ifstream	lo_neugen ("data/lo_extract_sf.txt")
ifstream	newlo_neugen ("data/newlo_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)
TCanvas	c2 ("c2","", 0, 0, 1300, 950)
c2	Divide (2, 2)
c2	cd (1)
a1	SetAxisRange (0.05, 20,"y")
a1	Draw ()
a1	GetXaxis () -> SetTitle("Q2")
a1	GetYaxis () -> SetTitle("F2")
gPad	SetLogy ()
gPad	SetLogx ()
gr_015	SetLineWidth (2)
gr_015	Draw ("C")
lo_gr_015	SetLineColor (kRed)
lo_gr_015	SetLineStyle (2)
tex	SetTextAlign (22)
tex	SetTextSize (0.045)
tex	DrawLatex (28.4, 12.8,"x = 0.015(5*F2)")
tex	DrawLatex (171.55, 5.825,"x = 0.045(3*F2)")
tex	DrawLatex (74.85, 2.076,"x = 0.080(2*F2)")
leg1	AddEntry (gr_015,"standard neugen F2","l")
leg1	AddEntry (lo_gr_015,"turn off B-Y correction","l")
leg1	AddEntry (newlo_gr_015,"+ turn off q2 freezing","l")
c2	cd (2)
tex	DrawLatex (80, 9,"x = 0.125(5*F2)")
tex	DrawLatex (119.76, 4.23,"x = 0.175(3*F2)")
tex	DrawLatex (95.99, 1.32,"x = 0.225(2*F2)")
c2	cd (3)
tex	DrawLatex (79.106, 6.16,"x = 0.275(5*F2)")
tex	DrawLatex (202.5, 2.45,"x = 0.350(3*F2)")
tex	DrawLatex (80, 0.4545,"x = 0.450(2*F2)")
c2	cd (4)
tex	DrawLatex (107.223, 2.000,"x = 0.550(5*F2)")
tex	DrawLatex (252.645, 0.4062,"x = 0.650(3*F2)")
tex	DrawLatex (83.603, 0.0303,"x = 0.750(2*F2)")
c2	Print ("f2_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	row6_base [nfile]
Double_t	row7_base [nfile]
Double_t	row8_base [nfile]
Double_t	row9_base [nfile]
Double_t	row10_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]
Double_t	newlo_row1_base [nfile]
Double_t	newlo_row2_base [nfile]
Double_t	newlo_row3_base [nfile]
Double_t	newlo_row4_base [nfile]
Double_t	newlo_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	f2 [12][num]
Int_t	count [12]
Int_t	num1 [12] ={5,3,2,5,3,2,5,3,2,5,3,2}
Double_t	lo_q2 [12][num]
Double_t	lo_f2 [12][num]
Int_t	lo_count [12]
Double_t	newlo_q2 [12][num]
Double_t	newlo_f2 [12][num]
Int_t	newlo_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],f2[0])
TGraph *	lo_gr_015 = new TGraph(lo_count[0],lo_q2[0],lo_f2[0])
TGraph *	newlo_gr_015 = new TGraph(newlo_count[0],newlo_q2[0],newlo_f2[0])
TLatex *	tex = new TLatex()
TGraph *	gr_045 = new TGraph(count[1],q2[1],f2[1])
TGraph *	lo_gr_045 = new TGraph(lo_count[1],lo_q2[1],lo_f2[1])
TGraph *	newlo_gr_045 = new TGraph(newlo_count[1],newlo_q2[1],newlo_f2[1])
TGraph *	gr_080 = new TGraph(count[2],q2[2],f2[2])
TGraph *	lo_gr_080 = new TGraph(lo_count[2],lo_q2[2],lo_f2[2])
TGraph *	newlo_gr_080 = new TGraph(newlo_count[2],newlo_q2[2],newlo_f2[2])
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],f2[3])
TGraph *	lo_gr_125 = new TGraph(lo_count[3],lo_q2[3],lo_f2[3])
TGraph *	newlo_gr_125 = new TGraph(newlo_count[3],newlo_q2[3],newlo_f2[3])
TGraph *	gr_175 = new TGraph(count[4],q2[4],f2[4])
TGraph *	lo_gr_175 = new TGraph(lo_count[4],lo_q2[4],lo_f2[4])
TGraph *	newlo_gr_175 = new TGraph(newlo_count[4],newlo_q2[4],newlo_f2[4])
TGraph *	gr_225 = new TGraph(count[5],q2[5],f2[5])
TGraph *	lo_gr_225 = new TGraph(lo_count[5],lo_q2[5],lo_f2[5])
TGraph *	newlo_gr_225 = new TGraph(newlo_count[5],newlo_q2[5],newlo_f2[5])
TH1F *	a3 = new TH1F("a3","",10,0.005,2000)
TGraph *	gr_275 = new TGraph(count[6],q2[6],f2[6])
TGraph *	lo_gr_275 = new TGraph(lo_count[6],lo_q2[6],lo_f2[6])
TGraph *	newlo_gr_275 = new TGraph(newlo_count[6],newlo_q2[6],newlo_f2[6])
TGraph *	gr_350 = new TGraph(count[7],q2[7],f2[7])
TGraph *	lo_gr_350 = new TGraph(lo_count[7],lo_q2[7],lo_f2[7])
TGraph *	newlo_gr_350 = new TGraph(newlo_count[7],newlo_q2[7],newlo_f2[7])
TGraph *	gr_450 = new TGraph(count[8],q2[8],f2[8])
TGraph *	lo_gr_450 = new TGraph(lo_count[8],lo_q2[8],lo_f2[8])
TGraph *	newlo_gr_450 = new TGraph(newlo_count[8],newlo_q2[8],newlo_f2[8])
TH1F *	a4 = new TH1F("a4","",10,0.005,2000)
TGraph *	gr_550 = new TGraph(count[9],q2[9],f2[9])
TGraph *	lo_gr_550 = new TGraph(lo_count[9],lo_q2[9],lo_f2[9])
TGraph *	newlo_gr_550 = new TGraph(newlo_count[9],newlo_q2[9],newlo_f2[9])
TGraph *	gr_650 = new TGraph(count[10],q2[10],f2[10])
TGraph *	lo_gr_650 = new TGraph(lo_count[10],lo_q2[10],lo_f2[10])
TGraph *	newlo_gr_650 = new TGraph(newlo_count[10],newlo_q2[10],newlo_f2[10])
TGraph *	gr_750 = new TGraph(count[11],q2[11],f2[11])
TGraph *	lo_gr_750 = new TGraph(lo_count[11],lo_q2[11],lo_f2[11])
TGraph *	newlo_gr_750 = new TGraph(newlo_count[11],newlo_q2[11],newlo_f2[11])

Function Documentation

leg1 AddEntry	(	gr_015	,
		"standard neugen F2"	,
		"l"
	)

leg1 AddEntry	(	lo_gr_015	,
		"turn off B-Y correction"	,
		"l"
	)

leg1 AddEntry	(	newlo_gr_015	,
		"+ turn off q2 freezing"	,
		"l"
	)

ifstream base_neugen

(

"data/extract_sf.txt"

)

TCanvas c2	(	"c2"	,
		""	,
		0	,
		0	,
		1300	,
		950
	)

ifstream ccfr	(	"data/ccfr_f2xf3fe.txt"	,
		ios::in
	)

c2 cd

(

1

)

c2 cd

(

2

)

c2 cd

(

3

)

c2 cd

(

4

)

ifstream cdhsw	(	"data/final_cdhsw_xf3.txt"	,
		ios::in
	)

c2 Divide	(	2	,
		2
	)

fBaryonPT2pdf Draw

(

)

lo_gr_750 Draw

(

"C"

)

tex DrawLatex	(	28.	4,
		12.	8
	)

tex DrawLatex	(	171.	55,
		5.	825
	)

tex DrawLatex	(	74.	85,
		2.	076
	)

tex DrawLatex	(	80	,
		9
	)

tex DrawLatex	(	119.	76,
		4.	23
	)

tex DrawLatex	(	95.	99,
		1.	32
	)

tex DrawLatex	(	79.	106,
		6.	16
	)

tex DrawLatex	(	202.	5,
		2.	45
	)

tex DrawLatex	(	80	,
		0.	4545
	)

tex DrawLatex	(	107.	223,
		2.	000
	)

tex DrawLatex	(	252.	645,
		0.	4062
	)

tex DrawLatex	(	83.	603,
		0.	0303
	)

for

(

)

Definition at line 31 of file f2_nu.C.

                        {
  base_neugen>>row1_base[i]>>row2_base[i]>>row3_base[i]>>row4_base[i]>>row5_base[i]>>row6_base[i]>>row7_base[i]>>row8_base[i]>>row9_base[i]>>row10_base[i];
  //cout<<"  "<<row1_base[i]<<"  "<<row2_base[i]<<"  "<<row3_base[i]<<"  "<<row4_base[i]<<"  "<<row5_base[i]<<endl;
}

a4 GetXaxis

(

)

-> SetTitle("Q2")

Definition at line 84 of file plot_lowE_sig_spread_genie.C.

{#nu} (GeV)");

a4 GetYaxis

(

)

-> SetTitle("F2")

ifstream lo_neugen

(

"data/lo_extract_sf.txt"

)

ifstream newlo_neugen

(

"data/newlo_extract_sf.txt"

)

ifstream nutev	(	"data/final_nutev_xf3.txt"	,
		ios::in
	)

c2 Print

(

"f2_nu.eps"

)

a4 SetAxisRange	(	0.	05,
		20	,
		"y"
	)

newlo_gr_750 SetLineColor

(

kRed

)

lo_gr_750 SetLineStyle

(

2

)

lo_gr_750 SetLineWidth

(

2

)

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 476 of file f2_nu.C.

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

Definition at line 545 of file f2_nu.C.

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

Definition at line 606 of file f2_nu.C.

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

Definition at line 666 of file f2_nu.C.

Double_t ccfr1[nccfr]

Definition at line 105 of file f2_nu.C.

Double_t ccfr2[nccfr]

Definition at line 106 of file f2_nu.C.

Double_t ccfr3[nccfr]

Definition at line 107 of file f2_nu.C.

Double_t ccfr4[nccfr]

Definition at line 108 of file f2_nu.C.

Double_t ccfr5[nccfr]

Definition at line 109 of file f2_nu.C.

Double_t ccfr6[nccfr]

Definition at line 110 of file f2_nu.C.

Double_t ccfr7[nccfr]

Definition at line 111 of file f2_nu.C.

ccfr_count[3] = count_sp1

Definition at line 264 of file f2_nu.C.

Double_t ccfr_err[12][num]

Definition at line 263 of file f2_nu.C.

Double_t ccfr_q2[12][num]

Definition at line 260 of file f2_nu.C.

Double_t ccfr_q2err[12][num]

Definition at line 262 of file f2_nu.C.

Double_t ccfr_xf3[12][num]

Definition at line 261 of file f2_nu.C.

Double_t cdhsw1[ncdhsw]

Definition at line 88 of file f2_nu.C.

Double_t cdhsw2[ncdhsw]

Definition at line 89 of file f2_nu.C.

Double_t cdhsw3[ncdhsw]

Definition at line 90 of file f2_nu.C.

Double_t cdhsw4[ncdhsw]

Definition at line 91 of file f2_nu.C.

Int_t cdhsw_count[12]

Definition at line 239 of file f2_nu.C.

Double_t cdhsw_err[12][num]

Definition at line 238 of file f2_nu.C.

Double_t cdhsw_q2[12][num]

Definition at line 235 of file f2_nu.C.

Double_t cdhsw_q2err[12][num]

Definition at line 237 of file f2_nu.C.

Double_t cdhsw_xf3[12][num]

Definition at line 236 of file f2_nu.C.

Int_t count[12]

Definition at line 140 of file f2_nu.C.

count_sp1 =0

Definition at line 283 of file f2_nu.C.

count_sp2 =0

Definition at line 284 of file f2_nu.C.

Double_t f2[12][num]

Definition at line 138 of file f2_nu.C.

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

Definition at line 484 of file f2_nu.C.

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

Definition at line 504 of file f2_nu.C.

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

Definition at line 521 of file f2_nu.C.

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

Definition at line 553 of file f2_nu.C.

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

Definition at line 570 of file f2_nu.C.

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

Definition at line 587 of file f2_nu.C.

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

Definition at line 614 of file f2_nu.C.

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

Definition at line 631 of file f2_nu.C.

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

Definition at line 648 of file f2_nu.C.

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

Definition at line 674 of file f2_nu.C.

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

Definition at line 691 of file f2_nu.C.

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

Definition at line 708 of file f2_nu.C.

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

Definition at line 538 of file f2_nu.C.

Int_t lo_count[12]

Definition at line 166 of file f2_nu.C.

Double_t lo_f2[12][num]

Definition at line 164 of file f2_nu.C.

TGraph* lo_gr_015 = new TGraph(lo_count[0],lo_q2[0],lo_f2[0])

Definition at line 488 of file f2_nu.C.

TGraph* lo_gr_045 = new TGraph(lo_count[1],lo_q2[1],lo_f2[1])

Definition at line 508 of file f2_nu.C.

TGraph* lo_gr_080 = new TGraph(lo_count[2],lo_q2[2],lo_f2[2])

Definition at line 525 of file f2_nu.C.

TGraph* lo_gr_125 = new TGraph(lo_count[3],lo_q2[3],lo_f2[3])

Definition at line 557 of file f2_nu.C.

TGraph* lo_gr_175 = new TGraph(lo_count[4],lo_q2[4],lo_f2[4])

Definition at line 574 of file f2_nu.C.

TGraph* lo_gr_225 = new TGraph(lo_count[5],lo_q2[5],lo_f2[5])

Definition at line 591 of file f2_nu.C.

TGraph* lo_gr_275 = new TGraph(lo_count[6],lo_q2[6],lo_f2[6])

Definition at line 618 of file f2_nu.C.

TGraph* lo_gr_350 = new TGraph(lo_count[7],lo_q2[7],lo_f2[7])

Definition at line 635 of file f2_nu.C.

TGraph* lo_gr_450 = new TGraph(lo_count[8],lo_q2[8],lo_f2[8])

Definition at line 652 of file f2_nu.C.

TGraph* lo_gr_550 = new TGraph(lo_count[9],lo_q2[9],lo_f2[9])

Definition at line 678 of file f2_nu.C.

TGraph* lo_gr_650 = new TGraph(lo_count[10],lo_q2[10],lo_f2[10])

Definition at line 695 of file f2_nu.C.

TGraph* lo_gr_750 = new TGraph(lo_count[11],lo_q2[11],lo_f2[11])

Definition at line 712 of file f2_nu.C.

Double_t lo_q2[12][num]

Definition at line 163 of file f2_nu.C.

Double_t lo_row1_base[nfile]

Definition at line 38 of file f2_nu.C.

Double_t lo_row2_base[nfile]

Definition at line 39 of file f2_nu.C.

Double_t lo_row3_base[nfile]

Definition at line 40 of file f2_nu.C.

Double_t lo_row4_base[nfile]

Definition at line 41 of file f2_nu.C.

Double_t lo_row5_base[nfile]

Definition at line 42 of file f2_nu.C.

int nccfr = 116

Definition at line 104 of file f2_nu.C.

int ncdhsw = 143

Definition at line 87 of file f2_nu.C.

Int_t newlo_count[12]

Definition at line 190 of file f2_nu.C.

Double_t newlo_f2[12][num]

Definition at line 188 of file f2_nu.C.

TGraph* newlo_gr_015 = new TGraph(newlo_count[0],newlo_q2[0],newlo_f2[0])

Definition at line 494 of file f2_nu.C.

TGraph* newlo_gr_045 = new TGraph(newlo_count[1],newlo_q2[1],newlo_f2[1])

Definition at line 514 of file f2_nu.C.

TGraph* newlo_gr_080 = new TGraph(newlo_count[2],newlo_q2[2],newlo_f2[2])

Definition at line 531 of file f2_nu.C.

TGraph* newlo_gr_125 = new TGraph(newlo_count[3],newlo_q2[3],newlo_f2[3])

Definition at line 563 of file f2_nu.C.

TGraph* newlo_gr_175 = new TGraph(newlo_count[4],newlo_q2[4],newlo_f2[4])

Definition at line 580 of file f2_nu.C.

TGraph* newlo_gr_225 = new TGraph(newlo_count[5],newlo_q2[5],newlo_f2[5])

Definition at line 597 of file f2_nu.C.

TGraph* newlo_gr_275 = new TGraph(newlo_count[6],newlo_q2[6],newlo_f2[6])

Definition at line 624 of file f2_nu.C.

TGraph* newlo_gr_350 = new TGraph(newlo_count[7],newlo_q2[7],newlo_f2[7])

Definition at line 641 of file f2_nu.C.

TGraph* newlo_gr_450 = new TGraph(newlo_count[8],newlo_q2[8],newlo_f2[8])

Definition at line 658 of file f2_nu.C.

TGraph* newlo_gr_550 = new TGraph(newlo_count[9],newlo_q2[9],newlo_f2[9])

Definition at line 684 of file f2_nu.C.

TGraph* newlo_gr_650 = new TGraph(newlo_count[10],newlo_q2[10],newlo_f2[10])

Definition at line 701 of file f2_nu.C.

TGraph* newlo_gr_750 = new TGraph(newlo_count[11],newlo_q2[11],newlo_f2[11])

Definition at line 718 of file f2_nu.C.

Double_t newlo_q2[12][num]

Definition at line 187 of file f2_nu.C.

Double_t newlo_row1_base[nfile]

Definition at line 52 of file f2_nu.C.

Double_t newlo_row2_base[nfile]

Definition at line 53 of file f2_nu.C.

Double_t newlo_row3_base[nfile]

Definition at line 54 of file f2_nu.C.

Double_t newlo_row4_base[nfile]

Definition at line 55 of file f2_nu.C.

Double_t newlo_row5_base[nfile]

Definition at line 56 of file f2_nu.C.

const int nfile = 3672

Definition at line 15 of file f2_nu.C.

int nnutev = 75

Definition at line 70 of file f2_nu.C.

int num = 306

Definition at line 119 of file f2_nu.C.

Int_t num1[12] ={5,3,2,5,3,2,5,3,2,5,3,2}

Definition at line 142 of file f2_nu.C.

Double_t nutev1[nnutev]

Definition at line 71 of file f2_nu.C.

Double_t nutev2[nnutev]

Definition at line 72 of file f2_nu.C.

Double_t nutev3[nnutev]

Definition at line 73 of file f2_nu.C.

Double_t nutev4[nnutev]

Definition at line 74 of file f2_nu.C.

Int_t nutev_count[12]

Definition at line 214 of file f2_nu.C.

Double_t nutev_err[12][num]

Definition at line 213 of file f2_nu.C.

Double_t nutev_q2[12][num]

Definition at line 210 of file f2_nu.C.

Double_t nutev_q2err[12][num]

Definition at line 212 of file f2_nu.C.

Double_t nutev_xf3[12][num]

Definition at line 211 of file f2_nu.C.

Double_t q2[12][num]

Definition at line 137 of file f2_nu.C.

Double_t row10_base[nfile]

Definition at line 29 of file f2_nu.C.

Double_t row1_base[nfile]

Definition at line 20 of file f2_nu.C.

Double_t row2_base[nfile]

Definition at line 21 of file f2_nu.C.

Double_t row3_base[nfile]

Definition at line 22 of file f2_nu.C.

Double_t row4_base[nfile]

Definition at line 23 of file f2_nu.C.

Double_t row5_base[nfile]

Definition at line 24 of file f2_nu.C.

Double_t row6_base[nfile]

Definition at line 25 of file f2_nu.C.

Double_t row7_base[nfile]

Definition at line 26 of file f2_nu.C.

Double_t row8_base[nfile]

Definition at line 27 of file f2_nu.C.

Double_t row9_base[nfile]

Definition at line 28 of file f2_nu.C.

Double_t temp_errsp1[num]

Definition at line 288 of file f2_nu.C.

Double_t temp_errsp2[num]

Definition at line 292 of file f2_nu.C.

Double_t temp_q21[num]

Definition at line 286 of file f2_nu.C.

Double_t temp_q22[num]

Definition at line 290 of file f2_nu.C.

Double_t temp_sp1[num]

Definition at line 287 of file f2_nu.C.

Double_t temp_sp2[num]

Definition at line 291 of file f2_nu.C.

TLatex* tex = new TLatex()

Definition at line 499 of file f2_nu.C.

x[11] = 0.015

Definition at line 120 of file f2_nu.C.