set -x
export LD_LIBRARY_PATH=.
base=two-beams
prob=$base/prob/matter
probext=0_0.8_1.0_0.04_5.0_2.5

do-p889 () {
for bl in 1522 1711 1840
do
# scale from 1km to baseline, scale from 0.5 MW to 1.5 MW
    sf=`echo "3.0 / ($bl * $bl)" | bc -l`

    out=$base/data/p889_${bl}_${probext}.ntuple
    rm -f $out
    ./collater \
               -s $sf \
               -f $base/tables/p889-flux.vec \
               -x "xsec/water.2.vec xsec/water.17.vec xsec/water.5.vec xsec/water.17.vec"  \
               -p "${prob}_numu_${bl}_${probext}.vec ${prob}_nue_${bl}_${probext}.vec ${prob}_antinumu_${bl}_${probext}.vec ${prob}_antinue_${bl}_${probext}.vec" \
               -e "0 10e9 10e6" \
	       -n $base/kumacs/nc_matrix.vec \
               -o $out
done
}

do-numi () {
for bl in 458 735 911
do
    case $bl in
    458) ang=5 ;;
    735) ang=4 ;;
    911) ang=3 ;;
    esac

# scale by from 1km to baseline, also scale by 10 due to change from
# .1GeV to .01GeV bins
    sf=`echo " 1.0 / ( 10.0 * $bl * $bl )" | bc -l`

    out=$base/data/numi_${bl}_${probext}.ntuple
    rm -f $out
    ./collater \
	       -s $sf \
               -f $base/tables/numi-flux-${ang}-1000.vec \
               -x "xsec/water.2.vec xsec/water.17.vec xsec/water.5.vec xsec/water.17.vec"  \
               -p "${prob}_numu_${bl}_${probext}.vec ${prob}_nue_${bl}_${probext}.vec ${prob}_antinumu_${bl}_${probext}.vec ${prob}_antinue_${bl}_${probext}.vec" \
               -e "0 10e9 10e6" \
	       -n $base/kumacs/nc_matrix.vec \
               -o $out
done
}

do-p889
do-numi