13 gStyle->SetOptStat(0000);
32 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];
36 ifstream
lo_neugen(
"data/lo_extract_sf.txt");
45 lo_neugen>>lo_row1_base[
i]>>lo_row2_base[
i]>>lo_row3_base[
i]>>lo_row4_base[
i]>>lo_row5_base[
i];
59 newlo_neugen>>newlo_row1_base[
i]>>newlo_row2_base[
i]>>newlo_row3_base[
i]>>newlo_row4_base[
i]>>newlo_row5_base[
i];
77 nutev>>nutev1[
i]>>nutev2[
i]>>nutev3[
i]>>nutev4[
i];
94 cdhsw>>cdhsw1[
i]>>cdhsw2[
i]>>cdhsw3[
i]>>cdhsw4[
i];
114 ccfr>>ccfr1[
i]>>ccfr2[
i]>>ccfr3[
i]>>ccfr4[
i]>>ccfr5[
i]>>ccfr6[
i]>>ccfr7[
i];
142 Int_t
num1[12]={5,3,2,5,3,2,5,3,2,5,3,2};
144 for(
int k=0;k<12;k++){
148 for(
int k=0;k<
nfile;k++){
149 for(
int l=0;
l<12;
l++){
150 if(
fabs(row1_base[k]-x[
l])<.001){
152 q2[
l][count[
l]] = row2_base[k];
153 f2[
l][count[
l]] = num1[
l]*row4_base[k];
168 for(
int k=0;k<12;k++){
172 for(
int k=0;k<
nfile;k++){
173 for(
int l=0;
l<12;
l++){
174 if(
fabs(lo_row1_base[k]-x[
l])<.001){
176 lo_q2[
l][lo_count[
l]] = lo_row2_base[k];
177 lo_f2[
l][lo_count[
l]] = num1[
l]*lo_row4_base[k];
192 for(
int k=0;k<12;k++){
196 for(
int k=0;k<
nfile;k++){
197 for(
int l=0;
l<12;
l++){
198 if(
fabs(newlo_row1_base[k]-x[
l])<.001){
200 newlo_q2[
l][newlo_count[
l]] = newlo_row2_base[k];
201 newlo_f2[
l][newlo_count[
l]] = num1[
l]*newlo_row4_base[k];
216 for(
int k=0;k<12;k++){
220 for(
int k=0;k<
nnutev;k++){
221 for(
int l=0;
l<12;
l++){
222 if(
fabs(nutev1[k]-x[
l])<.001){
224 nutev_q2[
l][nutev_count[
l]] = nutev2[k];
225 nutev_q2err[
l][nutev_count[
l]] = 0;
226 nutev_xf3[
l][nutev_count[
l]] = num1[
l]*nutev3[k];
227 nutev_err[
l][nutev_count[
l]] = num1[
l]*nutev4[k];
241 for(
int k=0;k<12;k++){
245 for(
int k=0;k<
ncdhsw;k++){
246 for(
int l=0;
l<12;
l++){
247 if(
fabs(cdhsw1[k]-x[
l])<.001){
249 cdhsw_q2[
l][cdhsw_count[
l]] = cdhsw2[k];
250 cdhsw_q2err[
l][cdhsw_count[
l]] = 0;
251 cdhsw_xf3[
l][cdhsw_count[
l]] = num1[
l]*cdhsw3[k];
252 cdhsw_err[
l][cdhsw_count[
l]] = num1[
l]*cdhsw4[k];
266 for(
int k=0;k<12;k++){
270 for(
int k=0;k<
nccfr;k++){
271 for(
int l=0;
l<12;
l++){
272 if(
fabs(ccfr1[k]-x[
l])<.001){
274 ccfr_q2[
l][ccfr_count[
l]] = ccfr2[k];
275 ccfr_q2err[
l][ccfr_count[
l]] = 0;
276 ccfr_xf3[
l][ccfr_count[
l]] = num1[
l]*ccfr5[k];
277 ccfr_err[
l][ccfr_count[
l]] = num1[
l]*ccfr6[k];
295 for(
int k=0;k<
nccfr;k++){
297 if(
fabs(ccfr1[k]-0.0125)<.001){
307 if(
fabs(ccfr1[k]-0.0175)<.001){
320 for(
int k=0;k<ccfr_count[0];k++){
322 ccfr_q2[0][k] = temp_q21[k];
323 ccfr_q2err[0][k] = 0;
325 ccfr_xf3[0][k] = (20*(temp_sp1[k]+temp_sp2[k]))/2;
326 ccfr_err[0][k] = (20*(temp_errsp1[k]+temp_errsp2[k]))/2;
334 for(
int k=0;k<
nccfr;k++){
336 if(
fabs(ccfr1[k]-0.035)<.001){
346 if(
fabs(ccfr1[k]-0.050)<.001){
359 for(
int k=0;k<ccfr_count[1];k++){
361 ccfr_q2[1][k] = temp_q21[k];
362 ccfr_q2err[1][k] = 0;
364 ccfr_xf3[1][k] = (5*(temp_sp1[k]+temp_sp2[k]))/2;
365 ccfr_err[1][k] = (5*(temp_errsp1[k]+temp_errsp2[k]))/2;
374 for(
int k=0;k<
nccfr;k++){
376 if(
fabs(ccfr1[k]-0.070)<.001){
386 if(
fabs(ccfr1[k]-0.090)<.001){
399 for(
int k=0;k<ccfr_count[2];k++){
401 ccfr_q2[2][k] = temp_q21[k];
402 ccfr_q2err[2][k] = 0;
404 ccfr_xf3[2][k] = ((temp_sp1[k]+temp_sp2[k]))/2;
405 ccfr_err[2][k] = ((temp_errsp1[k]+temp_errsp2[k]))/2;
415 for(
int k=0;k<
nccfr;k++){
417 if(
fabs(ccfr1[k]-0.110)<.001){
427 if(
fabs(ccfr1[k]-0.140)<.001){
440 for(
int k=0;k<ccfr_count[3];k++){
442 ccfr_q2[3][k] = temp_q21[k];
443 ccfr_q2err[3][k] = 0;
445 ccfr_xf3[3][k] = (5*(temp_sp1[k]+temp_sp2[k]))/2;
446 ccfr_err[3][k] = (5*(temp_errsp1[k]+temp_errsp2[k]))/2;
456 for(
int k=0;k<
nccfr;k++){
458 if(
fabs(ccfr1[k]-0.180)<.001){
460 ccfr_q2[4][ccfr_count[4]] = ccfr2[k];
461 ccfr_q2err[4][ccfr_count[4]] = 0;
463 ccfr_xf3[4][ccfr_count[4]] = 2*ccfr5[k];
464 ccfr_err[4][ccfr_count[4]] = 2*ccfr6[k];
472 TCanvas
c2(
"c2",
"",0,0,1300,950);
476 TH1F *
a1 =
new TH1F(
"a1",
"",10,0.005,2000);
477 a1->SetAxisRange(0.05,20,
"y");
479 a1->GetXaxis()->SetTitle(
"Q2");
480 a1->GetYaxis()->SetTitle(
"F2");
484 TGraph *
gr_015 =
new TGraph(count[0],q2[0],f2[0]);
485 gr_015->SetLineWidth(2);
488 TGraph *
lo_gr_015 =
new TGraph(lo_count[0],lo_q2[0],lo_f2[0]);
489 lo_gr_015->SetLineWidth(2);
490 lo_gr_015->SetLineColor(
kRed);
491 lo_gr_015->SetLineStyle(2);
492 lo_gr_015->Draw(
"C");
494 TGraph *
newlo_gr_015 =
new TGraph(newlo_count[0],newlo_q2[0],newlo_f2[0]);
495 newlo_gr_015->SetLineWidth(2);
496 newlo_gr_015->SetLineColor(
kRed);
497 newlo_gr_015->Draw(
"C");
499 TLatex *
tex =
new TLatex();
500 tex->SetTextAlign(22);
501 tex->SetTextSize(0.045);
502 tex->DrawLatex(28.4, 12.8,
"x = 0.015(5*F2)");
504 TGraph *
gr_045 =
new TGraph(count[1],q2[1],f2[1]);
505 gr_045->SetLineWidth(2);
508 TGraph *
lo_gr_045 =
new TGraph(lo_count[1],lo_q2[1],lo_f2[1]);
509 lo_gr_045->SetLineWidth(2);
510 lo_gr_045->SetLineColor(
kRed);
511 lo_gr_045->SetLineStyle(2);
512 lo_gr_045->Draw(
"C");
514 TGraph *
newlo_gr_045 =
new TGraph(newlo_count[1],newlo_q2[1],newlo_f2[1]);
515 newlo_gr_045->SetLineWidth(2);
516 newlo_gr_045->SetLineColor(
kRed);
517 newlo_gr_045->Draw(
"C");
519 tex->DrawLatex(171.55, 5.825,
"x = 0.045(3*F2)");
521 TGraph *
gr_080 =
new TGraph(count[2],q2[2],f2[2]);
522 gr_080->SetLineWidth(2);
525 TGraph *
lo_gr_080 =
new TGraph(lo_count[2],lo_q2[2],lo_f2[2]);
526 lo_gr_080->SetLineWidth(2);
527 lo_gr_080->SetLineColor(
kRed);
528 lo_gr_080->SetLineStyle(2);
529 lo_gr_080->Draw(
"C");
531 TGraph *
newlo_gr_080 =
new TGraph(newlo_count[2],newlo_q2[2],newlo_f2[2]);
532 newlo_gr_080->SetLineWidth(2);
533 newlo_gr_080->SetLineColor(
kRed);
534 newlo_gr_080->Draw(
"C");
536 tex->DrawLatex(74.85, 2.076,
"x = 0.080(2*F2)");
538 TLegend *
leg1 =
new TLegend(0.356,0.262,0.81668,0.541984);
539 leg1->AddEntry(gr_015,
"standard neugen F2",
"l");
540 leg1->AddEntry(lo_gr_015,
"turn off B-Y correction",
"l");
541 leg1->AddEntry(newlo_gr_015,
"+ turn off q2 freezing",
"l");
545 TH1F *
a2 =
new TH1F(
"a2",
"",10,0.005,2000);
546 a2->SetAxisRange(0.05,20,
"y");
548 a2->GetXaxis()->SetTitle(
"Q2");
549 a2->GetYaxis()->SetTitle(
"F2");
553 TGraph *
gr_125 =
new TGraph(count[3],q2[3],f2[3]);
554 gr_125->SetLineWidth(2);
557 TGraph *
lo_gr_125 =
new TGraph(lo_count[3],lo_q2[3],lo_f2[3]);
558 lo_gr_125->SetLineWidth(2);
559 lo_gr_125->SetLineColor(
kRed);
560 lo_gr_125->SetLineStyle(2);
561 lo_gr_125->Draw(
"C");
563 TGraph *
newlo_gr_125 =
new TGraph(newlo_count[3],newlo_q2[3],newlo_f2[3]);
564 newlo_gr_125->SetLineWidth(2);
565 newlo_gr_125->SetLineColor(
kRed);
566 newlo_gr_125->Draw(
"C");
568 tex->DrawLatex(80, 9,
"x = 0.125(5*F2)");
570 TGraph *
gr_175 =
new TGraph(count[4],q2[4],f2[4]);
571 gr_175->SetLineWidth(2);
574 TGraph *
lo_gr_175 =
new TGraph(lo_count[4],lo_q2[4],lo_f2[4]);
575 lo_gr_175->SetLineWidth(2);
576 lo_gr_175->SetLineColor(
kRed);
577 lo_gr_175->SetLineStyle(2);
578 lo_gr_175->Draw(
"C");
580 TGraph *
newlo_gr_175 =
new TGraph(newlo_count[4],newlo_q2[4],newlo_f2[4]);
581 newlo_gr_175->SetLineWidth(2);
582 newlo_gr_175->SetLineColor(
kRed);
583 newlo_gr_175->Draw(
"C");
585 tex->DrawLatex(119.76, 4.23,
"x = 0.175(3*F2)");
587 TGraph *
gr_225 =
new TGraph(count[5],q2[5],f2[5]);
588 gr_225->SetLineWidth(2);
591 TGraph *
lo_gr_225 =
new TGraph(lo_count[5],lo_q2[5],lo_f2[5]);
592 lo_gr_225->SetLineWidth(2);
593 lo_gr_225->SetLineColor(
kRed);
594 lo_gr_225->SetLineStyle(2);
595 lo_gr_225->Draw(
"C");
597 TGraph *
newlo_gr_225 =
new TGraph(newlo_count[5],newlo_q2[5],newlo_f2[5]);
598 newlo_gr_225->SetLineWidth(2);
599 newlo_gr_225->SetLineColor(
kRed);
600 newlo_gr_225->Draw(
"C");
602 tex->DrawLatex(95.99, 1.32,
"x = 0.225(2*F2)");
606 TH1F *
a3 =
new TH1F(
"a3",
"",10,0.005,2000);
607 a3->SetAxisRange(0.05,20,
"y");
609 a3->GetXaxis()->SetTitle(
"Q2");
610 a3->GetYaxis()->SetTitle(
"F2");
614 TGraph *
gr_275 =
new TGraph(count[6],q2[6],f2[6]);
615 gr_275->SetLineWidth(2);
618 TGraph *
lo_gr_275 =
new TGraph(lo_count[6],lo_q2[6],lo_f2[6]);
619 lo_gr_275->SetLineWidth(2);
620 lo_gr_275->SetLineColor(
kRed);
621 lo_gr_275->SetLineStyle(2);
622 lo_gr_275->Draw(
"C");
624 TGraph *
newlo_gr_275 =
new TGraph(newlo_count[6],newlo_q2[6],newlo_f2[6]);
625 newlo_gr_275->SetLineWidth(2);
626 newlo_gr_275->SetLineColor(
kRed);
627 newlo_gr_275->Draw(
"C");
629 tex->DrawLatex(79.106,6.16,
"x = 0.275(5*F2)");
631 TGraph *
gr_350 =
new TGraph(count[7],q2[7],f2[7]);
632 gr_350->SetLineWidth(2);
635 TGraph *
lo_gr_350 =
new TGraph(lo_count[7],lo_q2[7],lo_f2[7]);
636 lo_gr_350->SetLineWidth(2);
637 lo_gr_350->SetLineColor(
kRed);
638 lo_gr_350->SetLineStyle(2);
639 lo_gr_350->Draw(
"C");
641 TGraph *
newlo_gr_350 =
new TGraph(newlo_count[7],newlo_q2[7],newlo_f2[7]);
642 newlo_gr_350->SetLineWidth(2);
643 newlo_gr_350->SetLineColor(
kRed);
644 newlo_gr_350->Draw(
"C");
646 tex->DrawLatex(202.5, 2.45,
"x = 0.350(3*F2)");
648 TGraph *
gr_450 =
new TGraph(count[8],q2[8],f2[8]);
649 gr_450->SetLineWidth(2);
652 TGraph *
lo_gr_450 =
new TGraph(lo_count[8],lo_q2[8],lo_f2[8]);
653 lo_gr_450->SetLineWidth(2);
654 lo_gr_450->SetLineColor(
kRed);
655 lo_gr_450->SetLineStyle(2);
656 lo_gr_450->Draw(
"C");
658 TGraph *
newlo_gr_450 =
new TGraph(newlo_count[8],newlo_q2[8],newlo_f2[8]);
659 newlo_gr_450->SetLineWidth(2);
660 newlo_gr_450->SetLineColor(
kRed);
661 newlo_gr_450->Draw(
"C");
663 tex->DrawLatex(80, 0.4545,
"x = 0.450(2*F2)");
666 TH1F *
a4 =
new TH1F(
"a4",
"",10,0.005,2000);
667 a4->SetAxisRange(0.01,20,
"y");
669 a4->GetXaxis()->SetTitle(
"Q2");
670 a4->GetYaxis()->SetTitle(
"F2");
674 TGraph *
gr_550 =
new TGraph(count[9],q2[9],f2[9]);
675 gr_550->SetLineWidth(2);
678 TGraph *
lo_gr_550 =
new TGraph(lo_count[9],lo_q2[9],lo_f2[9]);
679 lo_gr_550->SetLineWidth(2);
680 lo_gr_550->SetLineColor(
kRed);
681 lo_gr_550->SetLineStyle(2);
682 lo_gr_550->Draw(
"C");
684 TGraph *
newlo_gr_550 =
new TGraph(newlo_count[9],newlo_q2[9],newlo_f2[9]);
685 newlo_gr_550->SetLineWidth(2);
686 newlo_gr_550->SetLineColor(
kRed);
687 newlo_gr_550->Draw(
"C");
689 tex->DrawLatex(107.223, 2.000,
"x = 0.550(5*F2)");
691 TGraph *
gr_650 =
new TGraph(count[10],q2[10],f2[10]);
692 gr_650->SetLineWidth(2);
695 TGraph *
lo_gr_650 =
new TGraph(lo_count[10],lo_q2[10],lo_f2[10]);
696 lo_gr_650->SetLineWidth(2);
697 lo_gr_650->SetLineColor(
kRed);
698 lo_gr_650->SetLineStyle(2);
699 lo_gr_650->Draw(
"C");
701 TGraph *
newlo_gr_650 =
new TGraph(newlo_count[10],newlo_q2[10],newlo_f2[10]);
702 newlo_gr_650->SetLineWidth(2);
703 newlo_gr_650->SetLineColor(
kRed);
704 newlo_gr_650->Draw(
"C");
706 tex->DrawLatex(252.645, 0.4062,
"x = 0.650(3*F2)");
708 TGraph *
gr_750 =
new TGraph(count[11],q2[11],f2[11]);
709 gr_750->SetLineWidth(2);
712 TGraph *
lo_gr_750 =
new TGraph(lo_count[11],lo_q2[11],lo_f2[11]);
713 lo_gr_750->SetLineWidth(2);
714 lo_gr_750->SetLineColor(
kRed);
715 lo_gr_750->SetLineStyle(2);
716 lo_gr_750->Draw(
"C");
718 TGraph *
newlo_gr_750 =
new TGraph(newlo_count[11],newlo_q2[11],newlo_f2[11]);
719 newlo_gr_750->SetLineWidth(2);
720 newlo_gr_750->SetLineColor(
kRed);
721 newlo_gr_750->Draw(
"C");
723 tex->DrawLatex(83.603, 0.0303,
"x = 0.750(2*F2)");
726 c2->Print(
"f2_nu.eps");
Double_t cdhsw_xf3[12][num]
Double_t row8_base[nfile]
Double_t row7_base[nfile]
Double_t newlo_row3_base[nfile]
fvar< T > fabs(const fvar< T > &x)
Double_t cdhsw_q2[12][num]
Double_t newlo_row1_base[nfile]
Double_t newlo_row5_base[nfile]
Double_t row3_base[nfile]
Double_t newlo_f2[12][num]
Double_t lo_row3_base[nfile]
Double_t row1_base[nfile]
Double_t cdhsw_err[12][num]
ifstream newlo_neugen("data/newlo_extract_sf.txt")
Double_t lo_row4_base[nfile]
Double_t ccfr_q2err[12][num]
Double_t nutev_err[12][num]
Double_t row9_base[nfile]
Double_t cdhsw_q2err[12][num]
Double_t nutev_q2err[12][num]
Double_t row2_base[nfile]
Double_t ccfr_xf3[12][num]
Double_t newlo_q2[12][num]
Double_t nutev_xf3[12][num]
Double_t row5_base[nfile]
Double_t lo_row2_base[nfile]
Double_t temp_errsp2[num]
ifstream cdhsw("data/final_cdhsw_xf3.txt", ios::in)
ifstream nutev("data/final_nutev_xf3.txt", ios::in)
Double_t lo_row1_base[nfile]
Double_t ccfr_err[12][num]
Double_t nutev_q2[12][num]
ifstream ccfr("data/ccfr_f2xf3fe.txt", ios::in)
Double_t row4_base[nfile]
TCanvas c2("c2","", 0, 0, 1300, 950)
Double_t row6_base[nfile]
Double_t ccfr_q2[12][num]
Double_t temp_errsp1[num]
Double_t lo_row5_base[nfile]
Double_t row10_base[nfile]
Double_t newlo_row2_base[nfile]
Double_t newlo_row4_base[nfile]
ifstream lo_neugen("data/lo_extract_sf.txt")
ifstream base_neugen("data/extract_sf.txt")