|
| 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") |
| |
|
| 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]) |
| |
1.8.11
