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

do-p889 () {
for bl in 350 1720 2540
do
# scale from 1km to baseline
    sf=`echo "1.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/tables/nc_matrix.vec \
               -o $out
done
}


do-p889