#!/bin/bash
#Script to run simplified study background vs. magnetic field on HIMster
#MOM=(15 8.9 11.91 4.06) #GeV
PANDAROOTmy=/home/karavdin/pandaRoot12/karavdina
PANDAROOT=/home/karavdin/pandaRoot12
cd ${PANDAROOTmy}
mkdir forward
mkdir backward
mkdir logs


MOM=(1.5)
MOMN=${#MOM[@]}
DMOM=(0 20)
DMOMN=${#DMOM[@]}
Qch=(-1 1)
QchN=${#Qch[@]}
for (( jq=0;jq<$QchN;jq++)); do
    charge=${Qch[${jq}]}
    for (( jp=0;jp<$MOMN;jp++)); do
	pbeam=${MOM[${jp}]}
	momStr=`echo $pbeam | sed -e '{ s/\./_/g; }'`
	for (( jdp=0;jdp<$DMOMN;jdp++)); do
	    dpbeam=${DMOM[${jdp}]}
	    cat >jobs/job_RK_${momStr}_${dpbeam}_${charge}.sh <<EOF
#!/bin/bash
#
#PBS -j oe 
#PBS -N RK${momStr}.${dpbeam}_${charge}
#PBS -o ${PANDAROOTmy}/logs/LOG_${momStr}_${dpbeam}_${charge}
#PBS -V
#PBS -l nodes=1:ppn=1,walltime=23:59:00
export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR
echo "Start Simulation for RK study with mom = ${pbeam}+/${dpbeam} and q=${charge}"
echo "Time of start: `date` "

	    cd ${PANDAROOTmy}
	    source ${PANDAROOT}/build/config.sh
	    time root -l -b -q RK_By.C+\(${pbeam},${dpbeam},${charge},true\) &> logs/${pbeam}_PLUS_${dpbeam}_q_${charge}_forward.log
	    time root -l -b -q RK_By.C+\(${pbeam},${dpbeam},${charge},false\) &> logs/${pbeam}_PLUS_${dpbeam}_q_${charge}_backward.log
echo "Time of end: `date` "
exit 0
EOF
### submit job to batch system
    	    iter=0
	    while [ $iter -lt 1 ];do 
		qsub jobs/job_RK_${momStr}_${dpbeam}_${charge}.sh
		done_test=$(echo $?)
		if [ "$done_test" -eq "0" ];then
		    echo "Job submitted!"
		    iter=1  
		else
		    echo "qsub did not work sleep for 60 seconds and then try again." 
		    sleep 60
		fi
	    done

	done 
    done
done
exit