# This script runs Lumi macros for the full Sim

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

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

#Declare array with beam energyes
#ENERGY=(1.5) #GeV
ENERGY=(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)"

	path="/home/karavdin/datastorage/TrackPerformanceStrip/"${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`
	cd /home/karavdin/pandaRoot12/macro/lmd
	source config.sh
	cat > ${path}/QA_jobs.log <<EOF
EOF
	cat >sum_strip_${SIMEVENTS}events_CA_${numTrks}trks_MinuitFit_${pbeam}.txt<<EOF
EOF
	cat>sum_strip_${SIMEVENTS}events_Follow_${numTrks}trks_MinuitFit_${pbeam}.txt<<EOF
EOF
	cat>sum_strip_${SIMEVENTS}events_CA_${numTrks}trks_KalmanFillter_${pbeam}.txt<<EOF
EOF
	cat>sum_strip_${SIMEVENTS}events_Follow_${numTrks}trks_KalmanFillter_${pbeam}.txt<<EOF
EOF
LEVEL=5
FIRSTCAMINUIT=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_0.root)
STANDARTCAMINUIT=$((${FIRSTCAMINUIT}*${LEVEL}/10))
FIRSTCAKAL=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_0.root)
STANDARTCAKAL=$((${FIRSTCAKAL}*${LEVEL}/10))
FIRSTFMINUIT=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_0.root)
STANDARTFMINUIT=$((${FIRSTFMINUIT}*${LEVEL}/10))
FIRSTFKAL=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_0.root)
STANDARTFKAL=$((${FIRSTFKAL}*${LEVEL}/10))
	for start in `seq 0 $step $last` ; do
	    if [ ! -r ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root ]; then
		echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: doesn't have output for CA and Minuit "  >> ${path}/QA_jobs.log 
	    else
		FILESIZE=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root)
		if [[ ${FILESIZE} -lt ${STANDARTCAMINUIT} ]]; then
		    echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: ouput for CA and Minuit  is too small"  >> ${path}/QA_jobs.log 
		else
		    echo -n "${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root "  >> sum_strip_${SIMEVENTS}events_CA_${numTrks}trks_MinuitFit_${pbeam}.txt
		fi
	    fi

	    if [ ! -r ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root ]; then
		echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: doesn't have output for Follow and Minuit "  >> ${path}/QA_jobs.log 
	    else
		FILESIZE=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root)
		if [[ ${FILESIZE} -lt ${STANDARTFMINUIT} ]]; then
		    echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: ouput for Follow and Minuit  is too small"  >> ${path}/QA_jobs.log 
		else
		    echo -n "${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}_GEANT4_${start}.root " >> sum_strip_${SIMEVENTS}events_Follow_${numTrks}trks_MinuitFit_${pbeam}.txt 
		fi
	    fi
	    if [ ! -r ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root ]; then
		echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: doesn't have output for CA and Kalman "  >> ${path}/QA_jobs.log 
	    else
		FILESIZE=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root)
		if [[ ${FILESIZE} -lt ${STANDARTCAKAL} ]]; then
		    echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: ouput for CA and Kalman  is too small"  >> ${path}/QA_jobs.log 
		else
		    echo -n "${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root " >> sum_strip_${SIMEVENTS}events_CA_${numTrks}trks_KalmanFillter_${pbeam}.txt
		fi
	    fi
	    if [ ! -r ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root ]; then
		echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: doesn't have output for Follow and Kalman "  >> ${path}/QA_jobs.log 
	    else
		FILESIZE=$(stat -c%s ${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root)
		if [[ ${FILESIZE} -lt ${STANDARTFKAL} ]]; then
		    echo "job_strip_${numTrks}trks_${SIMEVENTS}events_${start}.sh: ouput for Follow and Kalman  is too small"  >> ${path}/QA_jobs.log 
		else
		    echo -n "${path}/../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}_GEANT4_${start}.root " >> sum_strip_${SIMEVENTS}events_Follow_${numTrks}trks_KalmanFillter_${pbeam}.txt
		fi
	    fi
	done
	cat >/home/karavdin/myScripts/jobs/jobSumTrkRecPerfomanceStrip_${numTrks}trks_${momStr}mom.sh <<EOF
#!/bin/bash
#
#PBS -N lmdSUMstrip${momStr}
#PBS -j oe
#PBS -o LOGfilesumUPTrkRecPerformanceStrip${numTrks}trks${momStr}mom
#PBS -V
#PBS -l nodes=1:ppn=1,walltime=04:00:00
export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR
echo "Start analyse DPM Simulation with $nEvts events."
cd /home/karavdin/pandaRoot12/macro/lmd
source config.sh
### sum files
if [ -r ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}.root ]; then
rm ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}.root
fi
if [ -r ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}.root ]; then
rm ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}.root
fi
if [ -r ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}.root ]; then
rm ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}.root
fi
if [ -r ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}.root ]; then
rm ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}.root
fi
	${ROOTSYS}/bin/hadd -f ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_MinuitFit_${pbeam}.root ${filenames1} `cat sum_strip_${SIMEVENTS}events_CA_${numTrks}trks_MinuitFit_${pbeam}.txt`
	${ROOTSYS}/bin/hadd -f ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_CA_${numTrks}trks_KalmanFillter_${pbeam}.root ${filenames2} `cat sum_strip_${SIMEVENTS}events_CA_${numTrks}trks_KalmanFillter_${pbeam}.txt`
	${ROOTSYS}/bin/hadd -f ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_MinuitFit_${pbeam}.root ${filenames3} `cat sum_strip_${SIMEVENTS}events_Follow_${numTrks}trks_MinuitFit_${pbeam}.txt`
       ${ROOTSYS}/bin/hadd -f ${path}/../../Lumi_out_MC_and_REC_trks_matches_with_IDs0_Follow_${numTrks}trks_KalmanFillter_${pbeam}.root ${filenames4} `cat sum_strip_${SIMEVENTS}events_Follow_${numTrks}trks_KalmanFillter_${pbeam}.txt`
exit 0
EOF
qsub /home/karavdin/myScripts/jobs/jobSumTrkRecPerfomanceStrip_${numTrks}trks_${momStr}mom.sh
    done
done
exit