# This script runs Lumi macros for the full Sim

#-- fix params of simulation --------------------------------------------------
#numer of events used for resolution study
SIMEVENTS=1000000
nJobs=100
#SIMEVENTS=1000
#nJobs=1

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

#Declare array with beam energyes
#ENERGY=(1.5 4.06 8.9 11.91) #GeV
#ENERGY=(1.5 15) #GeV
ENERGY=(1.5 4.06 8.9 11.91) #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)"
##	echo ${T} ": Now we'll do test with "${SIMEVENTS} "events for track performance with "${numTrks} "trks per event and" ${pbeam} " GeV/c">> logTrkPerformanceRUN
	path="/home/karavdin/datastorage/TrackPerformancePixel/"${SIMEVENTS}"/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 ">> logTrkPerformancePIXEL

cat >jobs/mom_${momStr}/job_pixel_${TRKS[${jt}]}trks_${SIMEVENTS}events_${start}.sh <<EOF
#!/bin/bash
#
#PBS -N lmdPixel.$start
#PBS -j oe
##PBS -o ${path}/0-job_${start}.log
#PBS -o ${path}/../../testLOGfileTrkPerfomancePixel${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

    ##mergedHits=true
    ##addMS=true
    root -l -b -q runLumiPixel0SimBox.C\(${step},${start},${pbeam},"\"${path}\"",0,-2212,${numTrks},$RANDOM\) &>  ${path}/log${start}Pixel_SimBox_${pbeam}GeV
    root -l -b -q runLumiPixel1Digi.C\(${step},${start},"\"${path}\""\)  &>  ${path}/log${start}Pixel_Digi_${pbeam}GeV
    root -l -b -q runLumiPixel2Reco.C\(${step},${start},"\"${path}\"",0,true\)  &> ${path}/log${start}Pixel_HitRec_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumiPixel2bHitMerge.C\(${step},${start},"\"${path}\""\)  &> ${path}/log${start}Pixel_HitMerge_${numTrks}trks_mergedHits_${pbeam}GeV

    ##CA & Minuit
    root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",0,"\"CA\"",true,true\)  &> ${path}/log${start}Pixel_CA_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumiPixel4Fitter.C\(${step},${start},"\"${path}\"",0,true\)  &> ${path}/log${start}Pixel_MinuitFit_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\"",0,true\)  &> ${path}/log${start}Pixel_GEANE_${numTrks}trks_mergedHits_${pbeam}GeV
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 1 -mh 1 -v 0 -path ${path} &> ${path}/log${start}PIXEL_CA_Minuit_mc_rec_match_${pbeam}GeV
   
      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_mergedHits_true_addMS_true_MinuitFit_${pbeam}_GEANT4_${start}.root

    ###Trk-Follow & Minuit
    root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",0,"\"Follow\"",true,true\)  &> ${path}/log${start}Pixel_Follow_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumiPixel4Fitter.C\(${step},${start},"\"${path}\"",0,true\)
    root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\"",0,true\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 1  -mh 1 -v 0 -path ${path}
    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_mergedHits_true_addMS_true_MinuitFit_${pbeam}_GEANT4_${start}.root

    ###addMS=false
    root -l -b -q runLumiPixel2Reco.C\(${step},${start},"\"${path}\"",0,false\)
    ## CA & Kalman
    root -l -b -q runLumiPixel2bHitMerge.C\(${step},${start},"\"${path}\""\)
    root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",0,"\"CA\"",true,true\)
    time root -l -b -q runLumi4KalmanFitter.C\(${step},${start},"\"${path}\"",0,1,true\)  &> ${path}/log${start}Pixel_KalmanFit_${numTrks}trks_mergedHits_${pbeam}GeV
    time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\"",0,true\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 1  -mh 1 -v 0 -path ${path}  &> ${path}/log${start}PIXEL_Follow_Minuit_mc_rec_match_${pbeam}GeV
  
    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_mergedHits_true_addMS_false_KalmanFillter_${pbeam}_GEANT4_${start}.root

    ##Trk-Follow & Kalman
    root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",0,"\"Follow\"",true,true\)
    time root -l -b -q runLumi4KalmanFitter.C\(${step},${start},"\"${path}\"",0,1,true\)
    time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\"",0,true\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 1  -mh 1 -v 0 -path ${path}
    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_mergedHits_true_addMS_false_KalmanFillter_${pbeam}_GEANT4_${start}.root

  # #   # ##don't merge hits
# #    ## mergedHits=false
 
# # ##    addMS=true
# #     root -l -b -q runLumiPixel2Reco.C\(${step},${start},"\"${path}\"",0,true\)
# #     ##CA & Minuit
# #     root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",0,"\"CA\"",true,false\)  &> ${path}/log${start}Pixel_CA_${numTrks}trks_singleHits_${pbeam}GeV
# #     root -l -b -q runLumiPixel4Fitter.C\(${step},${start},"\"${path}\"",0,false\)
# #     root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\"",0,false\)
# #     cd ${VMCWORKDIR}/build/bin
# #     ./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 1 -mh 0 -v 0 -path ${path}

# #     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_mergedHits_false_addMS_true_MinuitFit_${pbeam}_GEANT4_${start}.root

# #     ##addMS=false
# #     root -l -b -q runLumiPixel2Reco.C\(${step},${start},"\"${path}\"",0,false\)
# #     ## CA & Kalman
# #     root -l -b -q runLumiPixel3Finder.C\(${step},${start},"\"${path}\"",0,"\"CA\"",true,false\)
# #     time root -l -b -q runLumi4KalmanFitter.C\(${step},${start},"\"${path}\"",0,1,false\)
# #     time root -l -b -q runLumi5Geane.C\(${step},${pbeam},${start},"\"${path}\"",0,false\)
# #     cd ${VMCWORKDIR}/build/bin
# #     ./rec_mc_match -s ${start} -n ${step} -t ${numTrks} -mom ${pbeam} -npx 1 -mh 0 -v 0 -path ${path}
# #     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_mergedHits_false_addMS_false_KalmanFillter_${pbeam}_GEANT4_${start}.root
exit 0
EOF
### submit job to batch system
qsub jobs/mom_${momStr}/job_pixel_${TRKS[${jt}]}trks_${SIMEVENTS}events_${start}.sh 
echo "Output files will be saved in $path"
	 done
     done
done
exit