# 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

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

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

step=`expr $nEvts / $nJobs`

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

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

echo Enter momentum from "1.5, 4.06, 8.9, 11.91, 15 [GeV/c]":
read mom
echo "Beam mometum is: $mom!"

momStr=`echo $mom | sed -e '{ s/\./_/g; }'`

path="/home/karavdin/datastorage/DPMoutput/all_"${momStr}"/"${folder}
# 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
# echo "Output files will be saved in $path"
if [ ! -d jobs/all_${momStr}/${folder}/ ]; then 
  mkdir -p jobs/all_${momStr}/${folder}/
fi

# let's go

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

cat >jobs/all_${momStr}/${folder}/job_${start}.sh <<EOF
#!/bin/bash
#
#PBS -N lmdDPMJob.$start
#PBS -j oe
##PBS -o ${path}/0-job_${start}.log
#PBS -o ${path}/testLOGDPMfile${start}
#PBS -V
#PBS -l nodes=1:ppn=1,walltime=23:00:00
export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR
echo "Start DPM Simulation with $nEvts events."

cd /home/karavdin/pandaRoot12/macro/lmd
source config.sh


# # # # [OLD FORMAT]
# #cd /home/karavdin/pandaRoot12/pgenerators/DpmEvtGen
# #./DPMGen -s $RANDOM -m 1 -p ${mom} \-n ${step} \-f DPM_all_${momStr}_${start}.root &> ${path}/0-sim_${start}.log

## # ##with DPMDirect
## # time root -l -b -q runLumi0DPMDirectbkg.C\($step,$start,$mom,"\"${path}\"",$verbose,$RANDOM\) &> ${path}/0a-sim_${start}.log

# # # #[NEW FORMAT]
# cd /home/karavdin/pandaRoot12/pgenerators/DpmEvtGenOct2012
# ## ./DPMGen -s $RANDOM -m ${mom} -e 1 -n ${step} -t 0.1 -f DPM_all_${momStr}_${start}.root &> ${path}/0-sim_${start}.log 
# ./DPMGen -s 1947 -m ${mom} -e 1 -n ${step} -t 0.1 -f DPM_all_${momStr}_${start}.root &> ${path}/0-sim_${start}.log 
#  mv DPM_all_${momStr}_${start}.root  ${path}/.
#  cd /home/karavdin/pandaRoot12/macro/lmd
#  time root -l -b -q runLumi0DPMbkg.C\($step,$start,$mom,"\"${path}\"",0,"\"${path}/DPM_all_${momStr}_${start}.root\""\) &> ${path}/0a-sim_${start}.log
# cd /home/karavdin/pandaRoot12/build/bin
# ./dpm_out_check -path ${path}  -s ${start} -n ${step} -p ${momStr} &> ${path}/0b-sim_${start}.log

cd /home/karavdin/pandaRoot12/macro/lmd
## time root -l -b -q runLumi1Digi.C\($step,$start,"\"${path}\"",$verbose\) &> ${path}/1-digi_${start}.log
## time root -l -b -q runLumi2Reco.C\($step,$start,"\"${path}\"",$verbose\) &> ${path}/2-reco_${start}.log
## time root -l -b -q runLumi3Finder.C\($step,$start,"\"${path}\"",$verbose,"\"CA\"",false\) &> ${path}/3-track_search${start}.log
## time root -l -b -q runLumi4Fitter.C\($step,$start,"\"${path}\"",$verbose\) &> ${path}/4-track_fit${start}.log
## time root -l -b -q runLumi5Geane.C\($step,$mom,$start,"\"${path}\"",$verbose\) &> ${path}/5-track_back_propagation${start}.log

 root -l -b -q runLumiPixel1Digi.C\(${step},${start},"\"${path}\"",$verbose\)  &> ${path}/1-digi_${start}.log
 ### addMS=false
 root -l -b -q runLumiPixel2Reco.C\(${step},${start},"\"${path}\"",$verbose,false\)  &> ${path}/2-reco_${start}.log
 ## CA & Kalman
 root -l -b -q runLumiPixel2bHitMerge.C\(${step},${start},"\"${path}\""\)  &>  ${path}/2b-reco_${start}.log
 root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",$verbose,"\"CA\"",true,true\)  &> ${path}/3-track_search${start}.log
 root -l -b -q runLumi4KalmanFitter.C\(${step},${start},"\"${path}\"",$verbose,1,true\)  &> ${path}/4-track_fit${start}.log
 root -l -b -q runLumi5Geane.C\(${step},${mom},${start},"\"${path}\"",$verbose,true\) &> ${path}/5-track_back_propagation${start}.log
 cd /home/karavdin/pandaRoot12/build/bin
./bkg_rec_mc_match -s ${start} -n ${step} -mom  ${mom} -path ${path} -v  $verbose >> ${path}/bkgStudy_ALL_${start}.log
cd /home/karavdin/pandaRoot12/macro/lmd
exit 0
EOF
### submit job to batch system
qsub jobs/all_${momStr}/${folder}/job_${start}.sh
echo "Output files will be saved in $path"
done