#!/bin/bash
if [ $# != 2 ]; then
    echo "usage: genr8.bash <Ebeam_GeV> <tslope_GeV>"
    exit 1
fi
if [ $2 = 0 ]; then
    Ebeam=$1
    tslope=0.001
else
    Ebeam=$1
    tslope=$2
fi
sed "s/<Ebeam>/$Ebeam/" p_pbar_proton.template |\
sed "s/<tslope>/$tslope/" > p_pbar_proton.input
genr8 -Agenr8_output.txt -M10000000 -r10000 < p_pbar_proton.input
genr8_2_hddm -V"0 0 50 80" genr8_output.txt
mv genr8_output.hddm genr8_${Ebeam}slope${tslope}.hddm
rm -f genr8_output.txt