#!/bin/bash
# HIMster script ----------------------------------------------------------------------
#part0: produce MC data (ideal)
# ver 0.0 [20/12/2013] 
# author: A.Karavdina
# ---------------------------------------------------------------------------------------
nDataSamp=100
numEv=50000 #number of events in one sample
#-----------------------------------------------------------------------------------------

verbose=0
PANDAROOT=/home/karavdin/pandaroot13
# #----------------------------------------------------------------------
#pbeam=1.5
echo Enter momentum from "1.5, 4.06, 8.9, 11.91, 15 [GeV/c]":
read pbeam
echo "Beam mometum is: $pbeam!"
momStr=`echo $pbeam | sed -e '{ s/\./_/g; }'`
pathG=/home/karavdin/datastorage/AlignmentLMDpixel/MCdataJuly2014/BOX/mom_${momStr}
mkdir -p ${pathG}
numTrks=1

# ### pick up half solenoid maps for 1.5 GeV/c [tmp, should be removed one day!] ###
# source /home/karavdin/pandaroot13/buildPanda/config.sh
# solsw=3 #GeV, switch between half and full solenoid: below half of solenoid
# if [[ $(echo "if (${pbeam} > ${solsw}) 1 else 0" | bc) -eq 1 ]]; then
#     echo "FULL solenoid field maps"
#     ln -sf ${VMCWORKDIR}/input/solenoid_fc/* ${VMCWORKDIR}/input/.
# else
#     echo "HALF solenoid field maps"
#     ln -sf ${VMCWORKDIR}/input/solenoid_hc/* ${VMCWORKDIR}/input/.
# fi
# ######################################################
if [ ! -d jobs/align/mom_${momStr}/ ]; then 
    mkdir -p jobs/align/mom_${momStr}/
fi
for(( jn=0;jn<${nDataSamp};jn++));do #diff data samples for statistic study
    path=${pathG}/Sample${jn}/
    if [ ! -r ${path} ]; then
	mkdir -p ${path}
    fi
    path_loc=/local/scratch/karavdina/align_MCdata/mom_${momStr}/Sample_${jn}/
	#	path_loc=${path}
    cat >jobs/align/mom_${momStr}/job_MCbox_sample${jn}.sh <<EOF
#!/bin/bash
#
#PBS -j oe 
#PBS -N Align_MC${jn}
#PBS -o ${path}/LOGhimsterFile_part0
#PBS -V
#PBS -l nodes=1:ppn=1,walltime=07:59:00
export PATH=\$PBS_O_PATH
cd \$PBS_O_WORKDIR
echo "Start MC data (BOX) simulation for Alignment study"
echo "Time of start: `date` "
       if [ ! -r ${path_loc} ]; then 
        mkdir -p ${path_loc}
       fi

     cd ${PANDAROOT}/macro/lmd
     source ${PANDAROOT}/buildPanda/config.sh
     root -l -b -q runLumiPixel0SimBox.C\(${numEv},0,"\"${path_loc}\"",0,-2212,${pbeam},${numTrks},$RANDOM\)  &> ${path_loc}/0-sim.log

   mv ${path_loc}/Lumi_* ${path}/.
   ##mv ${path_loc}/0-sim.log ${path}/.

rm -r ${path_loc}

echo "Time of end: `date` "
exit 0
EOF

### submit job to batch system
    		iter=0
		while [ $iter -lt 1 ];do 
		    qsub jobs/align/mom_${momStr}/job_MCbox_sample${jn}.sh
		    done_test=$(echo $?)
		    if [ "$done_test" -eq "0" ];then
			echo "Output files will be saved in $path"
			iter=1
		    else
			echo "qsub did not work sleep for 60 seconds and then try again." 
			sleep 60
		    fi
		done
	    done