# Edit the following lists
baselines="295"
#dmsols="8.2e-5"			# kamland hep-ex/0406035 + solar
#dmsols="7.9e-5"
dmsols="0.0"
dmatms="2.5e-3"
#sstt12s="0.816"			# kamland hep-ex/0406035 + solar
sstt12s="0.824"
sstt23s="1.0"
sstt13s="0.02"
#cps="0 45 90 135 180 225 270 -45 -90 -135 -180 -225 -270"
cps="45"
densities="2.3"
neutrinos="+2 -2"

# Edit the energy range
#energyrange="0.01 10 0.01"
energyrange="0.001 1.5 0.001"

# Maybe edit output locations
table_dir="./tables"
info_dir="./info"
plot_dir="./plots"

# For now, any arg is taken to mean force the prob vector file to be
# regenerated.
force="$1"


# leave the rest alone

for outputdirs in $table_dir $info_dir $plot_dir ; do
    if [ ! -d $outputdirs ] ; then
	mkdir -p $outputdirs
    fi
done

# Loop over parameters
for neutrino  in $neutrinos; do
for baseline  in $baselines; do
for dmsol     in $dmsols; do
for dmatm     in $dmatms; do
for sstt12    in $sstt12s; do
for sstt23    in $sstt23s; do
for sstt13    in $sstt13s; do
for cp        in $cps; do
for density   in $densities; do

    nuname="numu"
    case $neutrino in
	+1) nuname="nue";;
	+2) nuname="numu";;
	+3) nuname="nutau";;
	-1) nuname="anue";;
	-2) nuname="anumu";;
	-3) nuname="anutau";;
    esac

    base=${nuname}_${baseline}_${dmsol}_${dmatm}_${sstt12}_${sstt23}_${sstt13}_${cp}_${density}

    table_file="$table_dir/${base}.vec"

    # Don't regenerate.
    if [ -n "$force" -o ! -f $table_file ] ; then
	info_file="$info_dir/${base}.info"
	args="-n $neutrino -b $baseline -s $dmsol -a $dmatm -S $sstt12 $sstt23 $sstt13 -d $cp -D $density -E $energyrange"

	# do it
	echo "genprob $args > $table_file 2> $info_file"
	genprob $args > $table_file 2> $info_file
    fi

    # plot it
    plot_file=${plot_dir}/${base}.ps
    pargs="$nuname $baseline $dmsol $dmatm $sstt12 $sstt23 $sstt13 $cp $density $energyrange"
    pcmd="\n\ngen $table_file $plot_file $pargs\nquit\n"
    echo $pcmd
    echo -e $pcmd | paw

done;done;done;done;done;done;done;done;done;