# This script runs Lumi macros for the full Sim

#-- fix params of simulation --------------------------------------------------
#numer of events used for resolution study
SIMEVENTS=1000000

#Declare array with number of tracks per event
TRKS=(1 2 3 4 5 10 15)
#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 15) #GeV
#ENERGY=(1.5) #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}
	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
cat >jobs/mom_${momStr}/job_${TRKS[${jt}]}trks_${SIMEVENTS}events.sh <<EOF
#!/bin/bash
#
 #PBS -N lumiJob.$start
#PBS -j oe
##PBS -o ${path}/0-job_${start}.log
#PBS -o testLOGfile
#PBS -V
#PBS -l nodes=1:ppn=1,walltime=06:00:00
export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR
cd /home/karavdin/pandaRoot12/macro/lmd
source config.sh
rm ${path}/FairRunInfo*
rm ${path}/Lumi_MC_0.root
rm ${path}/Lumi_Params_0.root
rm ${path}/Lumi_digi_0.root
rm ${path}/Lumi_reco_0.root
rm ${path}/Lumi_TCand_0.root
rm ${path}/Lumi_Track_0.root
rm ${path}/Lumi_Geane_0.root

    ##mergedHits=true
    ##addMS=true
    root -l -b -q runLumiPixel0SimBox.C\(${SIMEVENTS},${pbeam},"\"${path}\"",0,-2212,${numTrks}\) &>  ${path}/logPixel_SimBox_${pbeam}GeV
    root -l -b -q runLumiPixel1Digi.C\(${SIMEVENTS},0,"\"${path}\""\)  &>  ${path}/logPixel_Digi_${pbeam}GeV
    root -l -b -q runLumiPixel2Reco.C\(${SIMEVENTS},0,"\"${path}\"",0,true\)  &> ${path}/logPixel_HitRec_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumiPixel2bHitMerge.C\(${SIMEVENTS},0,"\"${path}\""\)  &> ${path}/logPixel_HitMerge_${numTrks}trks_mergedHits_${pbeam}GeV

    ##CA & Minuit
    root -l -b -q runLumiPixel3Finder.C\(${SIMEVENTS},0,"\"${path}\"",0,"\"CA\"",true,true\)  &> ${path}/logPixel_CA_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumiPixel4Fitter.C\(${SIMEVENTS},0,"\"${path}\"",0,true\)  &> ${path}/logPixel_MinuitFit_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumi5Geane.C\(${SIMEVENTS},${pbeam},0,"\"${path}\"",0,true\)  &> ${path}/logPixel_GEANE_${numTrks}trks_mergedHits_${pbeam}GeV
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s 0 -n ${SIMEVENTS} -t ${numTrks} -mom ${pbeam} -npx 1 -mh 1 -v 0 -path ${path} &> ${path}/logPIXEL_CA_Minuit_mc_rec_match_${pbeam}GeV
   

   cd ${VMCWORKDIR}/macro/lmd
   mv ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0.root ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_mergedHits_true_addMS_true_MinuitFit_${pbeam}.root

    # ##Trk-Follow & Minuit
    root -l -b -q runLumiPixel3Finder.C\(${SIMEVENTS},0,"\"${path}\"",0,"\"Follow\"",true,true\)  &> ${path}/logPixel_Follow_${numTrks}trks_mergedHits_${pbeam}GeV
    root -l -b -q runLumiPixel4Fitter.C\(${SIMEVENTS},0,"\"${path}\"",0,true\)
    root -l -b -q runLumi5Geane.C\(${SIMEVENTS},${pbeam},0,"\"${path}\"",0,true\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s 0 -n ${SIMEVENTS} -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_IDs0.root ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_mergedHits_true_addMS_true_MinuitFit_${pbeam}.root

    ###addMS=false
    root -l -b -q runLumiPixel2Reco.C\(${SIMEVENTS},0,"\"${path}\"",0,false\)
    ## CA & Kalman
    root -l -b -q runLumiPixel2bHitMerge.C\(${SIMEVENTS},0,"\"${path}\""\)
    root -l -b -q runLumiPixel3Finder.C\(${SIMEVENTS},0,"\"${path}\"",0,"\"CA\"",true,true\)
    time root -l -b -q runLumi4KalmanFitter.C\(${SIMEVENTS},0,"\"${path}\"",0,1,true\)  &> ${path}/logPixel_KalmanFit_${numTrks}trks_mergedHits_${pbeam}GeV
    time root -l -b -q runLumi5Geane.C\(${SIMEVENTS},${pbeam},0,"\"${path}\"",0,true\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s 0 -n ${SIMEVENTS} -t ${numTrks} -mom ${pbeam} -npx 1  -mh 1 -v 0 -path ${path}  &> ${path}/logPIXEL_Follow_Minuit_mc_rec_match_${pbeam}GeV
  
    cd ${VMCWORKDIR}/macro/lmd
    mv ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0.root ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_mergedHits_true_addMS_false_KalmanFillter_${pbeam}.root

    ##Trk-Follow & Kalman
    root -l -b -q runLumiPixel3Finder.C\(${SIMEVENTS},0,"\"${path}\"",0,"\"Follow\"",true,true\)
    time root -l -b -q runLumi4KalmanFitter.C\(${SIMEVENTS},0,"\"${path}\"",0,1,true\)
    time root -l -b -q runLumi5Geane.C\(${SIMEVENTS},${pbeam},0,"\"${path}\"",0,true\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s 0 -n ${SIMEVENTS} -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_IDs0.root ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_mergedHits_true_addMS_false_KalmanFillter_${pbeam}.root

    # ##don't merge hits
   ## mergedHits=false
 
##    addMS=true
    root -l -b -q runLumiPixel2Reco.C\(${SIMEVENTS},0,"\"${path}\"",0,true\)
    ##CA & Minuit
    root -l -b -q runLumiPixel3Finder.C\(${SIMEVENTS},0,"\"${path}\"",0,"\"CA\"",true,false\)  &> ${path}/logPixel_CA_${numTrks}trks_singleHits_${pbeam}GeV
    root -l -b -q runLumiPixel4Fitter.C\(${SIMEVENTS},0,"\"${path}\"",0,false\)
    root -l -b -q runLumi5Geane.C\(${SIMEVENTS},${pbeam},0,"\"${path}\"",0,false\)
    cd ${VMCWORKDIR}/build/bin
    ./rec_mc_match -s 0 -n ${SIMEVENTS} -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_IDs0.root ${path}/Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_mergedHits_false_addMS_true_MinuitFit_${pbeam}.root

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