#!/bin/bash
#Script to run simplified study background vs. magnetic field on local PC
#MOM=(15 8.9 11.91 4.06) #GeV
PANDAROOTmy=/panda/karavdina/RungeKuttaSimplified
PANDAROOT=/panda/nov2013pandaroot/
cd ${PANDAROOTmy}
mkdir forward
mkdir backward
mkdir logs
nEv=2000
MOM=(15)
MOMN=${#MOM[@]}
#DMOM=(0)
DMOM=(-90 -80 -60 -40 -20 -10 10)
DMOMN=${#DMOM[@]}
DMOMN=${#DMOM[@]}
#Qch=(-3 -2 -1 1 2 3)
Qch=(-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}]}
	 
	    xycut=false
	    thcut=false

#	    th1=3 #mrad
#	    th2=8 #mrad
	    th1=0 #mrad
#	    th2=10 #mrad
#	    th2=785 #mrad =Pi/4.
	    th2=393 #mrad =Pi/8
#	    th2=1570 #mrad =Pi/2.
	    forwardAss=true
	    echo "Start Simulation with ${nEv} events for RK study with mom = ${pbeam}+/${dpbeam} and q=${charge}, th1=${th1}, th2=${th2} mrad"
	    echo "Time of start: `date` "
	    
	    cd ${PANDAROOTmy}
	    source ${PANDAROOT}/buildPanda/config.sh
	 
	    time root -l -b -q RK_By.C+\(${nEv},${pbeam},${dpbeam},${charge},${th1},${th2},true,${forwardAss},${xycut},${thcut},356,"\"${PANDAROOTmy}\""\) &> logs/${pbeam}_PLUS_${dpbeam}_q_${charge}_forward.log
		
	     	xycut2=true
	    	thcut2=false
		forwardAss2=true
	    	echo "back-propagation: isForward_for_backward=${forwardAss2}, cut_xy=${xycut2}, cut_thph=${thcut}"
	     	time root -l -b -q RK_By.C+\(${nEv},${pbeam},${dpbeam},${charge},${th1},${th2},false,${forwardAss2},${xycut2},${thcut2},356,"\"${PANDAROOTmy}\""\) &> logs/${pbeam}_PLUS_${dpbeam}_q_${charge}_backward.log
	
		echo "Time of end: `date` "
	done 
    done
done
exit