the flux files for Miller_82 and Barish_77 are the same.
for barish_77 only the 1horn is put into the flux
for the Miller_82 the flux is 0.75*1horn + 0.25*2horn

q2_nu_Baker_81.dat - one that Experiment reads
q2_nu_Baker_81_events.dat - one read by events 
q2_nu_Baker_81_original_errors.dat - the original one typed in, can be read by 
                          "errors"
                       note there is a 0.26 mm off set from zero which is 
                       not in
q2_nu_Baker_81_raw.dat

flux_spline_area_Kitagaki_83_fig_6.dat is the spline fit from figure 6 pg 440

to create the spline file using curves we copy table 1 of 
Kit_82 and run the awk script
awk '{print $1" " 10.0*log($2/0.05)/log(100.0/0.05)}' flux_Kitagaki_83_spline_table_1_orig.dat
or in terms of the code
awk '{print $1 " " yruler*log(yval/flux_low)/log(flux_high/flux_low)
x scale is 1  yruler=10.0, flux_high=100.0 flux_low=0.05


The output of the program for data flux
pg 440 of Kitagaki is
5 6.91824 10 131.433
10 14.4654 20 156.403
20 24.5283 30 69.5739
30 34.5912 40 30.643
40 45.283 50 12.906
50 59.7484 70 7.00364
70 84.9057 109.937 1.97226
109.937 133.333 200 0.598382

flux_Kitagaki_83_spline_table_1.dat
is used as the flux to create a 
spline flux
 14.8 66.8998 <- not actually spline flux
 24.4 30.7    <- not actually spline flux
 30 17.6
 40 9.43001
 60 2.76001
 80 1.19
 100 0.671002
 120 0.415999
 140 0.259
 160 0.157
 180 0.091
 200 0.051

to create the splinearea flux we
    cout<<" 30.   40. "<<ExpB.IntegFlux(30.,40.)/10.<<endl;
    cout<<" 40.   60. "<<ExpB.IntegFlux(40.,60.)/20.0<<endl;
    cout<<" 60.   80. "<<ExpB.IntegFlux(60.,80.)/20.0<<endl;
    cout<<" 80.   100. "<<ExpB.IntegFlux(80.,100.)/20.0<<endl;
    cout<<" 100.   120. "<<ExpB.IntegFlux(100.,120.)/20.0<<endl;
    cout<<" 120.   140. "<<ExpB.IntegFlux(120.,140.)/20.0<<endl;
    cout<<" 140.   160. "<<ExpB.IntegFlux(140.,160.)/20.0<<endl;
    cout<<" 160.   180. "<<ExpB.IntegFlux(160.,180.)/20.0<<endl;
    cout<<" 180.   200. "<<ExpB.IntegFlux(180.,200.)/20.0<<endl;
below is the output , but need to add the first 3 entries

  5.    10.0  56.21 <- scaling from Fig 6 in KIt_83
 10.    20.   66.9 <- from Table 1
 20.   30.   30.7 <- from table 1
 30.   40.   12.5016 
 40.   60.   5.7492
 60.   80.   1.67766
 80.   100.  0.927946
 100.   120. 0.522319
 120.   140. 0.334291
 140.   160. 0.203771
 160.   180. 0.121376
 180.   200. 0.0702249

to create 
flux_spline_area_Kitagaki_83_tab_1.dat
   0     5.000000    29.520214     3.558638     0.000000
   1    10.000000   -18.479144    13.158510    -0.479994
   2    20.000000   272.068875   -15.896292     0.246376
   3    30.000000    81.412197    -3.185847     0.034536
   4    40.000000    35.078721    -0.869173     0.005577
   5    60.000000    33.092200    -0.802956     0.005025
   6    80.000000    -3.970151     0.123603    -0.000766
   7   100.000000     9.857231    -0.152945     0.000617
   8   120.000000    -0.058908     0.012324    -0.000071
   9   140.000000     3.293915    -0.035573     0.000100
  10   160.000000     1.604283    -0.014453     0.000034
  11   180.000000     0.513517    -0.002333     0.000000
  12   200.000000     0.000000     0.000000     0.000000

Hardly any difference between
flux_spline_area_Kitagaki_83_tab_1.dat 
flux_spline_area_Kitagaki_83_fig_6.dat 
in the fits
also not much difference between 
flux_spline_area_Kitagaki_83_tab_1.dat and
flux_Kitagaki_83.dat  ( just uses the points of Fig 6, 1983) 

q2_nu_Barish_77.dat will be made into two files
q2_nu_Barish_77_with_overflow.dat -  puts in overflow correctly
                                  - but may not be good for plots
q2_nu_Barish_77_normal.dat - like the plots in the papers
