f2_nu.C
Go to the documentation of this file.
1 //
2 // Data vs GENIE neutrino F2 comparisons
3 //
4 // Contributed to GENIE/ValidationTools by Debdatta Bhattacharya (Pittsburgh Univ.)
5 // Jan 15h, 2009
6 //
7 // Revisions:
8 //
9 
10 
11 {
12 gROOT->Reset();
13 gStyle->SetOptStat(0000);
14 
15 const int nfile = 3672;
16 
17 // original f2,xf3 from neugen
18 ifstream base_neugen("data/extract_sf.txt");
19 
20 Double_t row1_base[nfile];
21 Double_t row2_base[nfile];
22 Double_t row3_base[nfile];
23 Double_t row4_base[nfile];
24 Double_t row5_base[nfile];
25 Double_t row6_base[nfile];
26 Double_t row7_base[nfile];
27 Double_t row8_base[nfile];
28 Double_t row9_base[nfile];
29 Double_t row10_base[nfile];
30 
31 for(int i=0;i<nfile;i++){
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];
33  //cout<<" "<<row1_base[i]<<" "<<row2_base[i]<<" "<<row3_base[i]<<" "<<row4_base[i]<<" "<<row5_base[i]<<endl;
34 }
35 
36 ifstream lo_neugen("data/lo_extract_sf.txt");
37 
38 Double_t lo_row1_base[nfile];
39 Double_t lo_row2_base[nfile];
40 Double_t lo_row3_base[nfile];
41 Double_t lo_row4_base[nfile];
42 Double_t lo_row5_base[nfile];
43 
44 for(int i=0;i<nfile;i++){
45  lo_neugen>>lo_row1_base[i]>>lo_row2_base[i]>>lo_row3_base[i]>>lo_row4_base[i]>>lo_row5_base[i];
46  //cout<<" "<<row1_base[i]<<" "<<row2_base[i]<<" "<<row3_base[i]<<" "<<row4_base[i]<<" "<<row5_base[i]<<endl;
47 }
48 
49 
50 ifstream newlo_neugen("data/newlo_extract_sf.txt");
51 
57 
58 for(int i=0;i<nfile;i++){
59  newlo_neugen>>newlo_row1_base[i]>>newlo_row2_base[i]>>newlo_row3_base[i]>>newlo_row4_base[i]>>newlo_row5_base[i];
60  //cout<<" "<<row1_base[i]<<" "<<row2_base[i]<<" "<<row3_base[i]<<" "<<row4_base[i]<<" "<<row5_base[i]<<endl;
61 }
62 
63 
64 //
65 // nutev
66 //
67 
68 ifstream nutev("data/final_nutev_xf3.txt",ios::in);
69 
70 int nnutev = 75;
71 Double_t nutev1[nnutev];
72 Double_t nutev2[nnutev];
73 Double_t nutev3[nnutev];
74 Double_t nutev4[nnutev];
75 
76 for(int i=0;i<nnutev;i++){
77  nutev>>nutev1[i]>>nutev2[i]>>nutev3[i]>>nutev4[i];
78  //cout<<" "<<nutev1[i]<<" "<<nutev2[i]<<" "<<nutev3[i]<<" "<<nutev4[i]<<endl;
79 }
80 
81 //
82 // cdhsw
83 //
84 
85 ifstream cdhsw("data/final_cdhsw_xf3.txt",ios::in);
86 
87 int ncdhsw = 143;
88 Double_t cdhsw1[ncdhsw];
89 Double_t cdhsw2[ncdhsw];
90 Double_t cdhsw3[ncdhsw];
91 Double_t cdhsw4[ncdhsw];
92 
93 for(int i=0;i<ncdhsw;i++){
94  cdhsw>>cdhsw1[i]>>cdhsw2[i]>>cdhsw3[i]>>cdhsw4[i];
95  //cout<<" "<<cdhsw1[i]<<" "<<cdhsw2[i]<<" "<<cdhsw3[i]<<" "<<cdhsw4[i]<<endl;
96 }
97 
98 //
99 // ccfr
100 //
101 
102 ifstream ccfr("data/ccfr_f2xf3fe.txt",ios::in);
103 
104 int nccfr = 116;
105 Double_t ccfr1[nccfr];
106 Double_t ccfr2[nccfr];
107 Double_t ccfr3[nccfr];
108 Double_t ccfr4[nccfr];
109 Double_t ccfr5[nccfr];
110 Double_t ccfr6[nccfr];
111 Double_t ccfr7[nccfr];
112 
113 for(int i=0;i<nccfr;i++){
114  ccfr>>ccfr1[i]>>ccfr2[i]>>ccfr3[i]>>ccfr4[i]>>ccfr5[i]>>ccfr6[i]>>ccfr7[i];
115  //cout<<" "<<ccfr1[i]<<" "<<ccfr2[i]<<" "<<ccfr3[i]<<" "<<ccfr4[i]<<endl;
116 }
117 
118 
119 int num = 306;
120 float x[12];
121 
122 x[0] = 0.015;
123 x[1] = 0.045;
124 x[2] = 0.080;
125 x[3] = 0.125;
126 x[4] = 0.175;
127 x[5] = 0.225;
128 x[6] = 0.275;
129 x[7] = 0.350;
130 x[8] = 0.450;
131 x[9] = 0.550;
132 x[10] = 0.650;
133 x[11] = 0.750;
134 
135 // base neugen sf
136 
137 Double_t q2[12][num];
138 Double_t f2[12][num];
139 
140 Int_t count[12];
141 
142 Int_t num1[12]={5,3,2,5,3,2,5,3,2,5,3,2};
143 
144 for(int k=0;k<12;k++){
145  count[k]=0;
146 }
147 
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){
151 
152  q2[l][count[l]] = row2_base[k];
153  f2[l][count[l]] = num1[l]*row4_base[k];
154 
155  count[l]++;
156  }
157  }
158 }
159 
160 
161 // lo neugen sf
162 
163 Double_t lo_q2[12][num];
164 Double_t lo_f2[12][num];
165 
166 Int_t lo_count[12];
167 
168 for(int k=0;k<12;k++){
169  lo_count[k]=0;
170 }
171 
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){
175 
176  lo_q2[l][lo_count[l]] = lo_row2_base[k];
177  lo_f2[l][lo_count[l]] = num1[l]*lo_row4_base[k];
178 
179  lo_count[l]++;
180  }
181  }
182 }
183 
184 
185 // newlo neugen sf
186 
187 Double_t newlo_q2[12][num];
188 Double_t newlo_f2[12][num];
189 
190 Int_t newlo_count[12];
191 
192 for(int k=0;k<12;k++){
193  newlo_count[k]=0;
194 }
195 
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){
199 
200  newlo_q2[l][newlo_count[l]] = newlo_row2_base[k];
201  newlo_f2[l][newlo_count[l]] = num1[l]*newlo_row4_base[k];
202 
203  newlo_count[l]++;
204  }
205  }
206 }
207 
208 
209 
210 Double_t nutev_q2[12][num];
211 Double_t nutev_xf3[12][num];
212 Double_t nutev_q2err[12][num];
213 Double_t nutev_err[12][num];
214 Int_t nutev_count[12];
215 
216 for(int k=0;k<12;k++){
217  nutev_count[k]=0;
218 }
219 
220 for(int k=0;k<nnutev;k++){
221  for(int l=0;l<12;l++){
222  if(fabs(nutev1[k]-x[l])<.001){
223 
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];
228 
229  nutev_count[l]++;
230  }
231  }
232 }
233 
234 //cdhsw
235 Double_t cdhsw_q2[12][num];
236 Double_t cdhsw_xf3[12][num];
237 Double_t cdhsw_q2err[12][num];
238 Double_t cdhsw_err[12][num];
239 Int_t cdhsw_count[12];
240 
241 for(int k=0;k<12;k++){
242  cdhsw_count[k]=0;
243 }
244 
245 for(int k=0;k<ncdhsw;k++){
246  for(int l=0;l<12;l++){
247  if(fabs(cdhsw1[k]-x[l])<.001){
248 
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];
253  cdhsw_count[l]++;
254  }
255  }
256 }
257 
258 
259 //ccfr
260 Double_t ccfr_q2[12][num];
261 Double_t ccfr_xf3[12][num];
262 Double_t ccfr_q2err[12][num];
263 Double_t ccfr_err[12][num];
264 Int_t ccfr_count[12];
265 
266 for(int k=0;k<12;k++){
267  ccfr_count[k]=0;
268 }
269 
270 for(int k=0;k<nccfr;k++){
271  for(int l=0;l<12;l++){
272  if(fabs(ccfr1[k]-x[l])<.001){
273 
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];
278  ccfr_count[l]++;
279  }
280  }
281 }
282 
283 Int_t count_sp1=0;
284 Int_t count_sp2=0;
285 
286 Double_t temp_q21[num];
287 Double_t temp_sp1[num];
288 Double_t temp_errsp1[num];
289 
290 Double_t temp_q22[num];
291 Double_t temp_sp2[num];
292 Double_t temp_errsp2[num];
293 
294 //special ccfr case - x=0.015
295 for(int k=0;k<nccfr;k++){
296 
297  if(fabs(ccfr1[k]-0.0125)<.001){
298 
299  temp_q21[count_sp1]= ccfr2[k];
300  temp_sp1[count_sp1]= ccfr5[k];
301  temp_errsp1[count_sp1]= ccfr6[k];
302 
303  //cout<<temp_sp1[count_sp1]<<endl;
304  count_sp1++;
305  }
306 
307  if(fabs(ccfr1[k]-0.0175)<.001){
308 
309  temp_q22[count_sp2]= ccfr2[k];
310  temp_sp2[count_sp2]= ccfr5[k];
311  temp_errsp2[count_sp2]= ccfr6[k];
312 
313  //cout<<temp_sp2[count_sp2]<<endl;
314  count_sp2++;
315  }
316 }
317 
318 ccfr_count[0] = count_sp1;
319 
320 for(int k=0;k<ccfr_count[0];k++){
321 
322  ccfr_q2[0][k] = temp_q21[k];
323  ccfr_q2err[0][k] = 0;
324 
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;
327 }
328 
329 
330 count_sp1=0;
331 count_sp2=0;
332 
333 //special ccfr case - x=0.0425
334 for(int k=0;k<nccfr;k++){
335 
336  if(fabs(ccfr1[k]-0.035)<.001){
337 
338  temp_q21[count_sp1]= ccfr2[k];
339  temp_sp1[count_sp1]= ccfr5[k];
340  temp_errsp1[count_sp1]= ccfr6[k];
341 
342  //cout<<temp_sp1[count_sp1]<<endl;
343  count_sp1++;
344  }
345 
346  if(fabs(ccfr1[k]-0.050)<.001){
347 
348  temp_q22[count_sp2]= ccfr2[k];
349  temp_sp2[count_sp2]= ccfr5[k];
350  temp_errsp2[count_sp2]= ccfr6[k];
351 
352  //cout<<temp_sp2[count_sp2]<<endl;
353  count_sp2++;
354  }
355 }
356 
357 ccfr_count[1] = count_sp1;
358 
359 for(int k=0;k<ccfr_count[1];k++){
360 
361  ccfr_q2[1][k] = temp_q21[k];
362  ccfr_q2err[1][k] = 0;
363 
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;
366 
367  //cout<<ccfr_q2[1][k]<<" "<<ccfr_xf3[1][k]<<endl;
368 }
369 
370 count_sp1=0;
371 count_sp2=0;
372 
373 //special ccfr case - x=0.0425
374 for(int k=0;k<nccfr;k++){
375 
376  if(fabs(ccfr1[k]-0.070)<.001){
377 
378  temp_q21[count_sp1]= ccfr2[k];
379  temp_sp1[count_sp1]= ccfr5[k];
380  temp_errsp1[count_sp1]= ccfr6[k];
381 
382  //cout<<temp_sp1[count_sp1]<<endl;
383  count_sp1++;
384  }
385 
386  if(fabs(ccfr1[k]-0.090)<.001){
387 
388  temp_q22[count_sp2]= ccfr2[k];
389  temp_sp2[count_sp2]= ccfr5[k];
390  temp_errsp2[count_sp2]= ccfr6[k];
391 
392  //cout<<temp_sp2[count_sp2]<<endl;
393  count_sp2++;
394  }
395 }
396 
397 ccfr_count[2] = count_sp1;
398 
399 for(int k=0;k<ccfr_count[2];k++){
400 
401  ccfr_q2[2][k] = temp_q21[k];
402  ccfr_q2err[2][k] = 0;
403 
404  ccfr_xf3[2][k] = ((temp_sp1[k]+temp_sp2[k]))/2;
405  ccfr_err[2][k] = ((temp_errsp1[k]+temp_errsp2[k]))/2;
406 
407  //cout<<ccfr_q2[2][k]<<" "<<ccfr_xf3[2][k]<<endl;
408 }
409 
410 
411 count_sp1=0;
412 count_sp2=0;
413 
414 //special ccfr case - x=0.125
415 for(int k=0;k<nccfr;k++){
416 
417  if(fabs(ccfr1[k]-0.110)<.001){
418 
419  temp_q21[count_sp1]= ccfr2[k];
420  temp_sp1[count_sp1]= ccfr5[k];
421  temp_errsp1[count_sp1]= ccfr6[k];
422 
423  //cout<<temp_sp1[count_sp1]<<endl;
424  count_sp1++;
425  }
426 
427  if(fabs(ccfr1[k]-0.140)<.001){
428 
429  temp_q22[count_sp2]= ccfr2[k];
430  temp_sp2[count_sp2]= ccfr5[k];
431  temp_errsp2[count_sp2]= ccfr6[k];
432 
433  //cout<<temp_sp2[count_sp2]<<endl;
434  count_sp2++;
435  }
436 }
437 
438 ccfr_count[3] = count_sp1;
439 
440 for(int k=0;k<ccfr_count[3];k++){
441 
442  ccfr_q2[3][k] = temp_q21[k];
443  ccfr_q2err[3][k] = 0;
444 
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;
447 
448  //cout<<ccfr_q2[2][k]<<" "<<ccfr_xf3[2][k]<<endl;
449 }
450 
451 
452 count_sp1=0;
453 count_sp2=0;
454 
455 //special ccfr case - x=0.180
456 for(int k=0;k<nccfr;k++){
457 
458  if(fabs(ccfr1[k]-0.180)<.001){
459 
460  ccfr_q2[4][ccfr_count[4]] = ccfr2[k];
461  ccfr_q2err[4][ccfr_count[4]] = 0;
462 
463  ccfr_xf3[4][ccfr_count[4]] = 2*ccfr5[k];
464  ccfr_err[4][ccfr_count[4]] = 2*ccfr6[k];
465 
466  ccfr_count[4]++;
467  }
468 }
469 
470 
471 
472 TCanvas c2("c2","",0,0,1300,950);
473 c2->Divide(2,2);
474 
475 c2->cd(1);
476 TH1F *a1 = new TH1F("a1","",10,0.005,2000);
477 a1->SetAxisRange(0.05,20,"y");
478 a1->Draw();
479 a1->GetXaxis()->SetTitle("Q2");
480 a1->GetYaxis()->SetTitle("F2");
481 gPad->SetLogy();
482 gPad->SetLogx();
483 
484 TGraph *gr_015 = new TGraph(count[0],q2[0],f2[0]);
485 gr_015->SetLineWidth(2);
486 gr_015->Draw("C");
487 
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");
493 
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");
498 
499 TLatex *tex = new TLatex();
500 tex->SetTextAlign(22); //centers the text on the chosen coordinates
501 tex->SetTextSize(0.045);//a reasonably large size
502 tex->DrawLatex(28.4, 12.8, "x = 0.015(5*F2)");
503 
504 TGraph *gr_045 = new TGraph(count[1],q2[1],f2[1]);
505 gr_045->SetLineWidth(2);
506 gr_045->Draw("C");
507 
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");
513 
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");
518 
519 tex->DrawLatex(171.55, 5.825, "x = 0.045(3*F2)");
520 
521 TGraph *gr_080 = new TGraph(count[2],q2[2],f2[2]);
522 gr_080->SetLineWidth(2);
523 gr_080->Draw("C");
524 
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");
530 
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");
535 
536 tex->DrawLatex(74.85, 2.076, "x = 0.080(2*F2)");
537 
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");
542 leg1->Draw();
543 
544 c2->cd(2);
545 TH1F *a2 = new TH1F("a2","",10,0.005,2000);
546 a2->SetAxisRange(0.05,20,"y");
547 a2->Draw();
548 a2->GetXaxis()->SetTitle("Q2");
549 a2->GetYaxis()->SetTitle("F2");
550 gPad->SetLogy();
551 gPad->SetLogx();
552 
553 TGraph *gr_125 = new TGraph(count[3],q2[3],f2[3]);
554 gr_125->SetLineWidth(2);
555 gr_125->Draw("C");
556 
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");
562 
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");
567 
568 tex->DrawLatex(80, 9, "x = 0.125(5*F2)");
569 
570 TGraph *gr_175 = new TGraph(count[4],q2[4],f2[4]);
571 gr_175->SetLineWidth(2);
572 gr_175->Draw("C");
573 
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");
579 
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");
584 
585 tex->DrawLatex(119.76, 4.23, "x = 0.175(3*F2)");
586 
587 TGraph *gr_225 = new TGraph(count[5],q2[5],f2[5]);
588 gr_225->SetLineWidth(2);
589 gr_225->Draw("C");
590 
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");
596 
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");
601 
602 tex->DrawLatex(95.99, 1.32, "x = 0.225(2*F2)");
603 
604 
605 c2->cd(3);
606 TH1F *a3 = new TH1F("a3","",10,0.005,2000);
607 a3->SetAxisRange(0.05,20,"y");
608 a3->Draw();
609 a3->GetXaxis()->SetTitle("Q2");
610 a3->GetYaxis()->SetTitle("F2");
611 gPad->SetLogy();
612 gPad->SetLogx();
613 
614 TGraph *gr_275 = new TGraph(count[6],q2[6],f2[6]);
615 gr_275->SetLineWidth(2);
616 gr_275->Draw("C");
617 
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");
623 
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");
628 
629 tex->DrawLatex(79.106,6.16, "x = 0.275(5*F2)");
630 
631 TGraph *gr_350 = new TGraph(count[7],q2[7],f2[7]);
632 gr_350->SetLineWidth(2);
633 gr_350->Draw("C");
634 
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");
640 
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");
645 
646 tex->DrawLatex(202.5, 2.45, "x = 0.350(3*F2)");
647 
648 TGraph *gr_450 = new TGraph(count[8],q2[8],f2[8]);
649 gr_450->SetLineWidth(2);
650 gr_450->Draw("C");
651 
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");
657 
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");
662 
663 tex->DrawLatex(80, 0.4545, "x = 0.450(2*F2)");
664 
665 c2->cd(4);
666 TH1F *a4 = new TH1F("a4","",10,0.005,2000);
667 a4->SetAxisRange(0.01,20,"y");
668 a4->Draw();
669 a4->GetXaxis()->SetTitle("Q2");
670 a4->GetYaxis()->SetTitle("F2");
671 gPad->SetLogy();
672 gPad->SetLogx();
673 
674 TGraph *gr_550 = new TGraph(count[9],q2[9],f2[9]);
675 gr_550->SetLineWidth(2);
676 gr_550->Draw("C");
677 
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");
683 
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");
688 
689 tex->DrawLatex(107.223, 2.000, "x = 0.550(5*F2)");
690 
691 TGraph *gr_650 = new TGraph(count[10],q2[10],f2[10]);
692 gr_650->SetLineWidth(2);
693 gr_650->Draw("C");
694 
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");
700 
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");
705 
706 tex->DrawLatex(252.645, 0.4062, "x = 0.650(3*F2)");
707 
708 TGraph *gr_750 = new TGraph(count[11],q2[11],f2[11]);
709 gr_750->SetLineWidth(2);
710 gr_750->Draw("C");
711 
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");
717 
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");
722 
723 tex->DrawLatex(83.603, 0.0303, "x = 0.750(2*F2)");
724 
725 
726 c2->Print("f2_nu.eps");
727 
728 
729 }
730 
731 
732 
733 /*TGraphErrors *grnutev_015 = new TGraphErrors(nutev_count[0],nutev_q2[0],nutev_xf3[0],nutev_q2err[0],nutev_err[0]);
734 grnutev_015->Draw("P*");
735 grnutev_015->SetMarkerColor(kMagenta);
736 grnutev_015->SetMarkerStyle(20);
737 grnutev_015->SetMarkerSize(0.9);
738 grnutev_015->SetLineColor(kMagenta);
739 TGraphErrors *grcdhsw_015 = new TGraphErrors(cdhsw_count[0],cdhsw_q2[0],cdhsw_xf3[0],cdhsw_q2err[0],cdhsw_err[0]);
740 grcdhsw_015->Draw("P*");
741 grcdhsw_015->SetMarkerColor(kGreen);
742 grcdhsw_015->SetLineColor(kGreen);
743 TGraphErrors *grccfr_015 = new TGraphErrors(ccfr_count[0],ccfr_q2[0],ccfr_xf3[0],ccfr_q2err[0],ccfr_err[0]);
744 grccfr_015->Draw("P*");
745 grccfr_015->SetMarkerColor(kMagenta);
746 grccfr_015->SetLineColor(kMagenta);*/
747 
748 
749 /*TGraphErrors *grnutev_045 = new TGraphErrors(nutev_count[1],nutev_q2[1],nutev_xf3[1],nutev_q2err[1],nutev_q2err[1]);
750 grnutev_045->Draw("P*");
751 grnutev_045->SetMarkerColor(kMagenta);
752 grnutev_045->SetMarkerStyle(20);
753 grnutev_045->SetMarkerSize(0.9);
754 grnutev_045->SetLineColor(kMagenta);
755 TGraphErrors *grcdhsw_045 = new TGraphErrors(cdhsw_count[1],cdhsw_q2[1],cdhsw_xf3[1],cdhsw_q2err[1],cdhsw_q2err[1]);
756 grcdhsw_045->Draw("P*");
757 grcdhsw_045->SetMarkerColor(kGreen);
758 grcdhsw_045->SetLineColor(kGreen);
759 TGraphErrors *grccfr_045 = new TGraphErrors(ccfr_count[1],ccfr_q2[1],ccfr_xf3[1],ccfr_q2err[1],ccfr_err[1]);
760 grccfr_045->Draw("P*");
761 grccfr_045->SetMarkerColor(kMagenta);
762 grccfr_045->SetLineColor(kMagenta);*/
763 
764 /*TGraphErrors *grnutev_080 = new TGraphErrors(nutev_count[2],nutev_q2[2],nutev_xf3[2],nutev_q2err[2],nutev_q2err[2]);
765 grnutev_080->Draw("P*");
766 grnutev_080->SetMarkerColor(kMagenta);
767 grnutev_080->SetMarkerStyle(20);
768 grnutev_080->SetMarkerSize(0.9);
769 grnutev_080->SetLineColor(kMagenta);
770 TGraphErrors *grcdhsw_080 = new TGraphErrors(cdhsw_count[2],cdhsw_q2[2],cdhsw_xf3[2],cdhsw_q2err[2],cdhsw_q2err[2]);
771 grcdhsw_080->Draw("P*");
772 grcdhsw_080->SetMarkerColor(kGreen);
773 grcdhsw_080->SetLineColor(kGreen);
774 TGraphErrors *grccfr_080 = new TGraphErrors(ccfr_count[2],ccfr_q2[2],ccfr_xf3[2],ccfr_q2err[2],ccfr_err[2]);
775 grccfr_080->Draw("P*");
776 grccfr_080->SetMarkerColor(kMagenta);
777 grccfr_080->SetLineColor(kMagenta);*/
778 
779 /*TGraphErrors *grnutev_125 = new TGraphErrors(nutev_count[3],nutev_q2[3],nutev_xf3[3],nutev_q2err[3],nutev_q2err[3]);
780 grnutev_125->Draw("P*");
781 grnutev_125->SetMarkerColor(kMagenta);
782 grnutev_125->SetMarkerStyle(20);
783 grnutev_125->SetMarkerSize(0.9);
784 grnutev_125->SetLineColor(kMagenta);
785 TGraphErrors *grcdhsw_125 = new TGraphErrors(cdhsw_count[3],cdhsw_q2[3],cdhsw_xf3[3],cdhsw_q2err[3],cdhsw_q2err[3]);
786 grcdhsw_125->Draw("P*");
787 grcdhsw_125->SetMarkerColor(kGreen);
788 grcdhsw_125->SetLineColor(kGreen);
789 TGraphErrors *grccfr_125 = new TGraphErrors(ccfr_count[3],ccfr_q2[3],ccfr_xf3[3],ccfr_q2err[3],ccfr_q2err[3]);
790 grccfr_125->Draw("P*");
791 grccfr_125->SetMarkerColor(kMagenta);
792 grccfr_125->SetLineColor(kMagenta);*/
793 
794 /*TGraphErrors *grnutev_175 = new TGraphErrors(nutev_count[4],nutev_q2[4],nutev_xf3[4],nutev_q2err[4],nutev_q2err[4]);
795 grnutev_175->Draw("P*");
796 grnutev_175->SetMarkerColor(kMagenta);
797 grnutev_175->SetMarkerStyle(20);
798 grnutev_175->SetMarkerSize(0.9);
799 grnutev_175->SetLineColor(kMagenta);
800 TGraphErrors *grcdhsw_175 = new TGraphErrors(cdhsw_count[4],cdhsw_q2[4],cdhsw_xf3[4],cdhsw_q2err[4],cdhsw_q2err[4]);
801 grcdhsw_175->Draw("P*");
802 grcdhsw_175->SetMarkerColor(kGreen);
803 grcdhsw_175->SetLineColor(kGreen);
804 TGraphErrors *grccfr_175 = new TGraphErrors(ccfr_count[4],ccfr_q2[4],ccfr_xf3[4],ccfr_q2err[4],ccfr_q2err[4]);
805 grccfr_175->Draw("P*");
806 grccfr_175->SetMarkerColor(kMagenta);
807 grccfr_175->SetLineColor(kMagenta);*/
808 
809 
810 /*TGraphErrors *grnutev_225 = new TGraphErrors(nutev_count[5],nutev_q2[5],nutev_xf3[5],nutev_q2err[5],nutev_q2err[5]);
811 grnutev_225->Draw("P*");
812 grnutev_225->SetMarkerColor(kMagenta);
813 grnutev_225->SetMarkerStyle(20);
814 grnutev_225->SetMarkerSize(0.9);
815 grnutev_225->SetLineColor(kMagenta);
816 TGraphErrors *grcdhsw_225 = new TGraphErrors(cdhsw_count[5],cdhsw_q2[5],cdhsw_xf3[5],cdhsw_q2err[5],cdhsw_q2err[5]);
817 grcdhsw_225->Draw("P*");
818 grcdhsw_225->SetMarkerColor(kGreen);
819 grcdhsw_225->SetLineColor(kGreen);
820 TGraphErrors *grccfr_225 = new TGraphErrors(ccfr_count[5],ccfr_q2[5],ccfr_xf3[5],ccfr_q2err[5],ccfr_q2err[5]);
821 grccfr_225->Draw("P*");
822 grccfr_225->SetMarkerColor(kMagenta);
823 grccfr_225->SetLineColor(kMagenta);*/
824 
825 
826 /*TGraphErrors *grnutev_275 = new TGraphErrors(nutev_count[6],nutev_q2[6],nutev_xf3[6],nutev_q2err[6],nutev_q2err[6]);
827 grnutev_275->Draw("P*");
828 grnutev_275->SetMarkerColor(kMagenta);
829 grnutev_275->SetMarkerStyle(20);
830 grnutev_275->SetMarkerSize(0.9);
831 grnutev_275->SetLineColor(kMagenta);
832 TGraphErrors *grcdhsw_275 = new TGraphErrors(cdhsw_count[6],cdhsw_q2[6],cdhsw_xf3[6],cdhsw_q2err[6],cdhsw_q2err[6]);
833 grcdhsw_275->Draw("P*");
834 grcdhsw_275->SetMarkerColor(kGreen);
835 grcdhsw_275->SetLineColor(kGreen);
836 TGraphErrors *grccfr_275 = new TGraphErrors(ccfr_count[6],ccfr_q2[6],ccfr_xf3[6],ccfr_q2err[6],ccfr_q2err[6]);
837 grccfr_275->Draw("P*");
838 grccfr_275->SetMarkerColor(kMagenta);
839 grccfr_275->SetLineColor(kMagenta);*/
840 
841 /*TGraphErrors *grnutev_350 = new TGraphErrors(nutev_count[7],nutev_q2[7],nutev_xf3[7],nutev_q2err[7],nutev_q2err[7]);
842 grnutev_350->Draw("P*");
843 grnutev_350->SetMarkerColor(kMagenta);
844 grnutev_350->SetMarkerStyle(20);
845 grnutev_350->SetMarkerSize(0.9);
846 grnutev_350->SetLineColor(kMagenta);
847 TGraphErrors *grcdhsw_350 = new TGraphErrors(cdhsw_count[7],cdhsw_q2[7],cdhsw_xf3[7],cdhsw_q2err[7],cdhsw_q2err[7]);
848 grcdhsw_350->Draw("P*");
849 grcdhsw_350->SetMarkerColor(kGreen);
850 grcdhsw_350->SetLineColor(kGreen);
851 TGraphErrors *grccfr_350 = new TGraphErrors(ccfr_count[7],ccfr_q2[7],ccfr_xf3[7],ccfr_q2err[7],ccfr_q2err[7]);
852 grccfr_350->Draw("P*");
853 grccfr_350->SetMarkerColor(kMagenta);
854 grccfr_350->SetLineColor(kMagenta);*/
855 
856 /*TGraphErrors *grnutev_450 = new TGraphErrors(nutev_count[8],nutev_q2[8],nutev_xf3[8],nutev_q2err[8],nutev_q2err[8]);
857 grnutev_450->Draw("P*");
858 grnutev_450->SetMarkerColor(kMagenta);
859 grnutev_450->SetMarkerStyle(20);
860 grnutev_450->SetMarkerSize(0.9);
861 grnutev_450->SetLineColor(kMagenta);
862 TGraphErrors *grcdhsw_450 = new TGraphErrors(cdhsw_count[8],cdhsw_q2[8],cdhsw_xf3[8],cdhsw_q2err[8],cdhsw_q2err[8]);
863 grcdhsw_450->Draw("P*");
864 grcdhsw_450->SetMarkerColor(kGreen);
865 grcdhsw_450->SetLineColor(kGreen);
866 TGraphErrors *grccfr_450 = new TGraphErrors(ccfr_count[8],ccfr_q2[8],ccfr_xf3[8],ccfr_q2err[8],ccfr_q2err[8]);
867 grccfr_450->Draw("P*");
868 grccfr_450->SetMarkerColor(kMagenta);
869 grccfr_450->SetLineColor(kMagenta);*/
870 
871 /*TGraphErrors *grnutev_550 = new TGraphErrors(nutev_count[9],nutev_q2[9],nutev_xf3[9],nutev_q2err[9],nutev_q2err[9]);
872 grnutev_550->Draw("P*");
873 grnutev_550->SetMarkerColor(kMagenta);
874 grnutev_550->SetMarkerStyle(20);
875 grnutev_550->SetMarkerSize(0.9);
876 grnutev_550->SetLineColor(kMagenta);
877 TGraphErrors *grcdhsw_550 = new TGraphErrors(cdhsw_count[9],cdhsw_q2[9],cdhsw_xf3[9],cdhsw_q2err[9],cdhsw_q2err[9]);
878 grcdhsw_550->Draw("P*");
879 grcdhsw_550->SetMarkerColor(kGreen);
880 grcdhsw_550->SetLineColor(kGreen);
881 TGraphErrors *grccfr_550 = new TGraphErrors(ccfr_count[9],ccfr_q2[9],ccfr_xf3[9],ccfr_q2err[9],ccfr_q2err[9]);
882 grccfr_550->Draw("P*");
883 grccfr_550->SetMarkerColor(kMagenta);
884 grccfr_550->SetLineColor(kMagenta);*/
885 
886 
887 /*TGraphErrors *grnutev_650 = new TGraphErrors(nutev_count[10],nutev_q2[10],nutev_xf3[10],nutev_q2err[10],nutev_q2err[10]);
888 grnutev_650->Draw("P*");
889 grnutev_650->SetMarkerColor(kMagenta);
890 grnutev_650->SetMarkerStyle(20);
891 grnutev_650->SetMarkerSize(0.9);
892 grnutev_650->SetLineColor(kMagenta);
893 TGraphErrors *grcdhsw_650 = new TGraphErrors(cdhsw_count[10],cdhsw_q2[10],cdhsw_xf3[10],cdhsw_q2err[10],cdhsw_q2err[10]);
894 grcdhsw_650->Draw("P*");
895 grcdhsw_650->SetMarkerColor(kGreen);
896 grcdhsw_650->SetLineColor(kGreen);
897 TGraphErrors *grccfr_650 = new TGraphErrors(ccfr_count[10],ccfr_q2[10],ccfr_xf3[10],ccfr_q2err[10],ccfr_q2err[10]);
898 grccfr_650->Draw("P*");
899 grccfr_650->SetMarkerColor(kMagenta);
900 grccfr_650->SetLineColor(kMagenta);*/
901 
902 
903 /*TGraphErrors *grnutev_750 = new TGraphErrors(nutev_count[11],nutev_q2[11],nutev_xf3[11],nutev_q2err[11],nutev_q2err[11]);
904 grnutev_750->Draw("P*");
905 grnutev_750->SetMarkerColor(kMagenta);
906 grnutev_750->SetMarkerStyle(20);
907 grnutev_750->SetMarkerSize(0.9);
908 grnutev_750->SetLineColor(kMagenta);
909 TGraphErrors *grccfr_750 = new TGraphErrors(ccfr_count[11],ccfr_q2[11],ccfr_xf3[11],ccfr_q2err[11],ccfr_q2err[11]);
910 grccfr_750->Draw("P*");
911 grccfr_750->SetMarkerColor(kMagenta);
912 grccfr_750->SetLineColor(kMagenta);*/
Double_t nutev4[nnutev]
Definition: f2_nu.C:74
Double_t cdhsw_xf3[12][num]
Definition: f2_nu.C:236
TH1F * a3
Definition: f2_nu.C:606
Double_t row8_base[nfile]
Definition: f2_nu.C:27
TGraph * newlo_gr_015
Definition: f2_nu.C:494
Double_t row7_base[nfile]
Definition: f2_nu.C:26
TLegend * leg1
Definition: f2_nu.C:538
TH1F * a2
Definition: f2_nu.C:545
Double_t newlo_row3_base[nfile]
Definition: f2_nu.C:54
fvar< T > fabs(const fvar< T > &x)
Definition: fabs.hpp:15
Double_t cdhsw_q2[12][num]
Definition: f2_nu.C:235
TGraph * lo_gr_750
Definition: f2_nu.C:712
TGraph * gr_045
Definition: f2_nu.C:504
TGraph * lo_gr_015
Definition: f2_nu.C:488
Double_t newlo_row1_base[nfile]
Definition: f2_nu.C:52
Double_t nutev1[nnutev]
Definition: f2_nu.C:71
Double_t ccfr7[nccfr]
Definition: f2_nu.C:111
TGraph * gr_275
Definition: f2_nu.C:614
Int_t count_sp1
Definition: f2_nu.C:283
Double_t newlo_row5_base[nfile]
Definition: f2_nu.C:56
TGraph * gr_550
Definition: f2_nu.C:674
TGraph * newlo_gr_045
Definition: f2_nu.C:514
Double_t temp_sp1[num]
Definition: f2_nu.C:287
Double_t row3_base[nfile]
Definition: f2_nu.C:22
Int_t cdhsw_count[12]
Definition: f2_nu.C:239
TGraph * newlo_gr_550
Definition: f2_nu.C:684
TGraph * lo_gr_350
Definition: f2_nu.C:635
Int_t count[12]
Definition: f2_nu.C:140
Double_t lo_f2[12][num]
Definition: f2_nu.C:164
TGraph * lo_gr_045
Definition: f2_nu.C:508
Double_t newlo_f2[12][num]
Definition: f2_nu.C:188
Int_t nutev_count[12]
Definition: f2_nu.C:214
Double_t lo_row3_base[nfile]
Definition: f2_nu.C:40
TGraph * gr_650
Definition: f2_nu.C:691
TGraph * lo_gr_175
Definition: f2_nu.C:574
TGraph * newlo_gr_650
Definition: f2_nu.C:701
TGraph * lo_gr_650
Definition: f2_nu.C:695
Double_t ccfr2[nccfr]
Definition: f2_nu.C:106
Double_t row1_base[nfile]
Definition: f2_nu.C:20
Int_t num1[12]
Definition: f2_nu.C:142
Double_t cdhsw_err[12][num]
Definition: f2_nu.C:238
TGraph * newlo_gr_450
Definition: f2_nu.C:658
const int nfile
Definition: f2_nu.C:15
TGraph * newlo_gr_225
Definition: f2_nu.C:597
TGraph * gr_015
Definition: f2_nu.C:484
Double_t q2[12][num]
Definition: f2_nu.C:137
Double_t ccfr5[nccfr]
Definition: f2_nu.C:109
ifstream newlo_neugen("data/newlo_extract_sf.txt")
int ncdhsw
Definition: f2_nu.C:87
int nccfr
Definition: f2_nu.C:104
TGraph * gr_175
Definition: f2_nu.C:570
Double_t ccfr6[nccfr]
Definition: f2_nu.C:110
TGraph * gr_350
Definition: f2_nu.C:631
Double_t lo_row4_base[nfile]
Definition: f2_nu.C:41
TGraph * gr_225
Definition: f2_nu.C:587
Double_t ccfr_q2err[12][num]
Definition: f2_nu.C:262
Double_t nutev_err[12][num]
Definition: f2_nu.C:213
Double_t row9_base[nfile]
Definition: f2_nu.C:28
TGraph * lo_gr_125
Definition: f2_nu.C:557
TGraph * newlo_gr_750
Definition: f2_nu.C:718
Double_t f2[12][num]
Definition: f2_nu.C:138
TGraph * newlo_gr_275
Definition: f2_nu.C:624
Int_t lo_count[12]
Definition: f2_nu.C:166
TH1F * a1
Definition: f2_nu.C:476
Double_t cdhsw_q2err[12][num]
Definition: f2_nu.C:237
Double_t nutev_q2err[12][num]
Definition: f2_nu.C:212
float x[12]
Definition: f2_nu.C:120
TGraph * newlo_gr_350
Definition: f2_nu.C:641
Double_t row2_base[nfile]
Definition: f2_nu.C:21
Double_t ccfr3[nccfr]
Definition: f2_nu.C:107
Double_t ccfr_xf3[12][num]
Definition: f2_nu.C:261
Double_t newlo_q2[12][num]
Definition: f2_nu.C:187
TGraph * newlo_gr_080
Definition: f2_nu.C:531
Double_t nutev_xf3[12][num]
Definition: f2_nu.C:211
Double_t cdhsw3[ncdhsw]
Definition: f2_nu.C:90
Double_t temp_q21[num]
Definition: f2_nu.C:286
Double_t row5_base[nfile]
Definition: f2_nu.C:24
int nnutev
Definition: f2_nu.C:70
TGraph * gr_080
Definition: f2_nu.C:521
Double_t lo_row2_base[nfile]
Definition: f2_nu.C:39
Double_t temp_errsp2[num]
Definition: f2_nu.C:292
Double_t temp_q22[num]
Definition: f2_nu.C:290
TGraph * lo_gr_550
Definition: f2_nu.C:678
Double_t lo_q2[12][num]
Definition: f2_nu.C:163
ifstream cdhsw("data/final_cdhsw_xf3.txt", ios::in)
ifstream nutev("data/final_nutev_xf3.txt", ios::in)
TGraph * newlo_gr_125
Definition: f2_nu.C:563
ifstream in
Definition: comparison.C:7
Double_t nutev2[nnutev]
Definition: f2_nu.C:72
Double_t lo_row1_base[nfile]
Definition: f2_nu.C:38
Double_t ccfr_err[12][num]
Definition: f2_nu.C:263
Double_t ccfr4[nccfr]
Definition: f2_nu.C:108
Double_t nutev_q2[12][num]
Definition: f2_nu.C:210
int num
Definition: f2_nu.C:119
Double_t ccfr1[nccfr]
Definition: f2_nu.C:105
ifstream ccfr("data/ccfr_f2xf3fe.txt", ios::in)
TGraph * lo_gr_450
Definition: f2_nu.C:652
TGraph * newlo_gr_175
Definition: f2_nu.C:580
TLatex * tex
Definition: f2_nu.C:499
Int_t newlo_count[12]
Definition: f2_nu.C:190
TGraph * gr_125
Definition: f2_nu.C:553
TH1F * a4
Definition: f2_nu.C:666
Double_t temp_sp2[num]
Definition: f2_nu.C:291
Double_t row4_base[nfile]
Definition: f2_nu.C:23
Int_t ccfr_count[12]
Definition: f2_nu.C:264
TCanvas c2("c2","", 0, 0, 1300, 950)
Double_t row6_base[nfile]
Definition: f2_nu.C:25
Double_t ccfr_q2[12][num]
Definition: f2_nu.C:260
Double_t temp_errsp1[num]
Definition: f2_nu.C:288
TGraph * gr_450
Definition: f2_nu.C:648
TGraph * lo_gr_225
Definition: f2_nu.C:591
Double_t lo_row5_base[nfile]
Definition: f2_nu.C:42
Double_t cdhsw1[ncdhsw]
Definition: f2_nu.C:88
Int_t count_sp2
Definition: f2_nu.C:284
Double_t row10_base[nfile]
Definition: f2_nu.C:29
Double_t cdhsw4[ncdhsw]
Definition: f2_nu.C:91
Double_t newlo_row2_base[nfile]
Definition: f2_nu.C:53
Double_t cdhsw2[ncdhsw]
Definition: f2_nu.C:89
TGraph * gr_750
Definition: f2_nu.C:708
Double_t newlo_row4_base[nfile]
Definition: f2_nu.C:55
Double_t nutev3[nnutev]
Definition: f2_nu.C:73
TGraph * lo_gr_275
Definition: f2_nu.C:618
ifstream lo_neugen("data/lo_extract_sf.txt")
ifstream base_neugen("data/extract_sf.txt")
TGraph * lo_gr_080
Definition: f2_nu.C:525