# This script runs Lumi macros for the full Sim

# reads number of Events, name of run and verbose Level
nEvts=10
nJobs=1
path="test"
verbose=0
mom=0
foilth=1
pitch=0
lowSmear=1
usefield=1

echo Enter number of Events:
read nEvts
echo "Number of Events: $nEvts!"

echo Enter number of Jobs:
read nJobs
echo "Number of Events: $nJobs!"

step=`expr $nEvts / $nJobs`

echo Enter storage-folder:
read folder
echo "Storagepath: $folder!"

echo Foil: 0=150mum, 1=300mum, 2=20mum
read foilth
echo "Use foil: $foilth!"

echo Pitch: 0=10, 1=50, 2=100, 3=200, 4=500, 5=1000 mum
read pitch
echo "Use Pitch: $pitch!"

echo Smearing: 0=low, 1=medium, 2=high
read lowSmear
echo "Use Smearing: $lowSmear!"

echo use Fields: 0=no, 1=yes
read usefield
echo "Use Field: $usefield!"

echo Enter verbosity Level:
read verbose
echo "Verbosity Level is: $verbose!"

for imom in `seq 0 1 3` ; do

if [ $imom -eq 0 ];
  then
   mom=2;
fi
if [ $imom -eq 1 ];
  then
   mom=16;
fi
if [ $imom -eq 2 ];
  then
   mom=10.1;
fi
if [ $imom -eq 3 ];
  then
   mom=6.2;
fi

momStr=`echo $mom | sed -e '{ s/\./_/g; }'`
input="DpmInput/el_100k_aida/el_"${momStr}"GeV.root"
path="dataDpm/"${folder}"/el_"${momStr}

# make data directory, if not there. all data file go there, makes them easier to move around

if [ ! -d $path ]; then 
  mkdir -p $path
fi

if [ ! -d jobs/${folder}/el_${momStr}/ ]; then 
  mkdir -p jobs/${folder}/el_${momStr}/
fi

# let's go

last=`expr $nEvts - 1`
for start in `seq 0 $step $last` ; do

cat >jobs/${folder}/el_${momStr}/job_${start}.sh <<EOF
#!/bin/bash
#
#PBS -N lumiJob.$start
#PBS -j oe
#PBS -o ${path}/0-job_${start}.log
#PBS -V
#PBS -l nodes=1:x86_64,walltime=04:00:00

export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR

echo "Start PandaRoot Simulation with $nEvts events."
root -l -q -b runLumi0SimDpm.C\($mom\,"\"${input}\"",$step,$start,"\"${path}\"",$verbose,$foilth,$usefield\) &> ${path}/0-sim_${start}.log
echo "Start Digitization."
root -l -q -b runLumi1Digi.C\($step,$start,"\"${path}\"",$verbose,$pitch,$lowSmear\) &> ${path}/1-digi_${start}.log
echo "Start Reco."
root -l -q -b runLumi2Reco.C\($step,$start,"\"${path}\"",$verbose\) &> ${path}/2-reco_${start}.log
echo "Start simple Trackfinding with 4 discs."
root -l -q -b runLumi3myTrackFit.C\($step,$start,"\"${path}\"",$verbose\) &> ${path}/3-track_${start}.log
#echo "Start simple Trackfinding with 3 discs."
#root -l -q -b runLumi3myTrackFit.C\($step,$start,"\"${path}\"",$verbose\,0.5,1\) &> ${path}/3-track_no4d_${start}.log

EOF


### submit job to batch system
qsub jobs/${folder}/el_${momStr}/job_${start}.sh


done
done