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