# This script runs Lumi macros for the full Sim

#-- fix params of simulation --------------------------------------------------
#numer of events used for resolution study
SIMEVENTS=1000000
nJobs=100
step=`expr $SIMEVENTS / $nJobs`

#Declare array with number of tracks per event
TRKS=(1)
#TRKS=(1 2 3 4)
# get number of elements in the array
TRKSN=${#TRKS[@]}

#Declare array with beam energyes
#ENERGY=(1.5 4.06 8.9 11.91 15) #GeV
ENERGY=(1.5 4.06 8.9 11.91 15) #GeV
#ENERGY=(15) #GeV
ENERGYN=${#ENERGY[@]}
#---------------------------------------------------------------------------------------

# reads number of Events, name of run and verbose Level
path="test"

# let's go
for (( jt=0;jt<$TRKSN;jt++)); do
    for (( jp=0;jp<$ENERGYN;jp++)); do
	momStr=`echo  ${ENERGY[${jp}]} | sed -e '{ s/\./_/g; }'`
	pbeam=${ENERGY[${jp}]} 
	numTrks=${TRKS[${jt}]}
##	T="$(date)"

       path="/home/karavdin/datastorage/TrackPerformanceStrip/"${SIMEVENTS}"/mom_"${momStr}"/"${numTrks}"trks"
	##path="/home/karavdin/datastorage/TrackPerformanceStrip/"${SIMEVENTS}"/GEANT3/mom_"${momStr}"/"${numTrks}"trks"
	if [ ! -d $path ]; then 
	    mkdir -p $path
	fi
 #echo "Output files will be saved in $path"
	if [ ! -d jobs/mom_${momStr} ]; then 
	    mkdir -p jobs/mom_${momStr}
	fi
	last=`expr $SIMEVENTS - 1`
	for start in `seq 0 $step $last` ; do
##	    echo ${T} ": Now we'll do test with "${step} "events (start at "${start}") for track performance with "${numTrks} "trks per event and" ${pbeam} " GeV/c">> logTrkPerformanceSTRIP
	    cat >jobs/mom_${momStr}/job_strip_${TRKS[${jt}]}trks_${SIMEVENTS}events_${start}.sh <<EOF
#!/bin/bash
#
 #PBS -N lmdStrip.$start
#PBS -j oe
##PBS -o ${path}/0-job_${start}.log
#PBS -o ${path}/../../testLOGfileTrkPerfomanceStrip${TRKS[${jt}]}trks${momStr}mom${start}start
#PBS -V
##PBS -l nodes=1:ppn=1,walltime=13:55:00,pmem=2gb,vmem=2gb
#PBS -l nodes=1:ppn=1,walltime=04:00:00
export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR
cd /home/karavdin/pandaRoot12/macro/lmd
source config.sh
time root -l -b -q runLumi0SimBox.C\(${step},${start},${pbeam},"\"${path}\"",0,-2212,${numTrks},$RANDOM\)  &> ${path}/log${start}Strip_Sim_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi1Digi.C\(${step},${start},"\"${path}\"",0\)  &> ${path}/log${start}Strip_Digi_${numTrks}trks_${pbeam}GeV

####addMS=true
time root -l -b -q runLumi2Reco.C\(${step},${start},"\"${path}\"",0,true\)  &> ${path}/log${start}Strip_HitRec_${numTrks}trks_${pbeam}GeV

time root -l -b -q runLumi3Finder.C\(${step},${start},"\"${path}\"",0,"\"CA\"",true\)  &> ${path}/log${start}Strip_CA_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi4Fitter.C\(${step},${start},"\"${path}\""\)  &> ${path}/log${start}Strip_MinuitFit_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\""\)  &> ${path}/log${start}Strip_GEANE_${numTrks}trks_${pbeam}GeV
cd ${VMCWORKDIR}/build/bin
./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 0 -v 0 -path ${path}  &> ${path}/log${start}STRIP_CA_Minuit_mc_rec_match

cd ${VMCWORKDIR}/macro/lmd/
mv ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs${start}.root ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root


time root -l -b -q runLumi3Finder.C\(${step},${start},"\"${path}\"",0,"\"Follow\"",true\)  &> ${path}/log${start}Strip_Follow_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi4Fitter.C\(${step},${start},"\"${path}\""\)
time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\""\)
cd ${VMCWORKDIR}/build/bin
./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 0 -v 0 -path ${path}  &> ${path}/log${start}STRIP_Follow_Minuit_mc_rec_match

cd ${VMCWORKDIR}/macro/lmd/
mv ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs${start}.root ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root

##addMS=false
time root -l -b -q runLumi2Reco.C\(${step},${start},"\"${path}\"",0,false\)

time root -l -b -q runLumi3Finder.C\(${step},${start},"\"${path}\"",0,"\"CA\"",true\) &> ${path}/log${start}Strip_CA2_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi4KalmanFitter.C\(${step},${start},"\"${path}\"",0,0\)  &> ${path}/log${start}Strip_CA2_KalmanFit_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\""\)
 
cd ${VMCWORKDIR}/build/bin
./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 0 -v 3 -path ${path} &> ${path}/log${start}STRIP_CA_Kalman_mc_rec_match
cd ${VMCWORKDIR}/macro/lmd
mv ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs${start}.root ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root

time root -l -b -q runLumi3Finder.C\(${step},${start},"\"${path}\"",0,"\"Follow\"",true\) &> ${path}/log${start}Strip_Follow2_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi4KalmanFitter.C\(${step},${start},"\"${path}\"",0,0\) &> ${path}/log${start}Strip_Follow2_KalmanFit_${numTrks}trks_${pbeam}GeV
time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\""\)
 
cd ${VMCWORKDIR}/build/bin
./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 0 -v 0 -path ${path}  &> ${path}/log${start}STRIP_Follow_Kalman_mc_rec_match
cd ${VMCWORKDIR}/macro/lmd
mv ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs${start}.root ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root 
exit 0
EOF
### submit job to batch system
	    qsub jobs/mom_${momStr}/job_strip_${TRKS[${jt}]}trks_${SIMEVENTS}events_${start}.sh 
	    echo "Output files will be saved in $path"
	done
    done
done
exit