# This script runs Lumi macros for the full Sim
################################################################
# Script for production and merge results of BkgStudies_MCmatch_withSignal_exe.C (Background\Signal study)
# !need all files from full reconstruction chain (results from GEANE, MC etc.)
# author: Anastasia Karavdina
# date: April, 2012
#
# to run it print and follow instuctions: 
# ./HIMsterDATA.sh
#!NB: Storagepath is a directory in ${path}
#change ${path} to right value in your case
################################################################

# reads number of Events, name of run and verbose Level
nEvts=10
nJobs=1
path="test"
verbose=0
mom=0

echo Enter number of Events:
read nEvts 
echo "Number of Events: $nEvts!"

echo Enter number of Jobs:
read nJobs
echo "Number of Jobs: $nJobs!"

step=`expr $nEvts / $nJobs`

echo Enter storage-folder:
read folder
echo "Storagepath: $folder!"

echo Enter verbosity Level:
read verbose
echo "Verbosity Level is: $verbose!"

echo Enter momentum from "1.5, 4.06, 8.9, 11.91, 15 [GeV/c]":
read mom
echo "Beam momentum is: $mom!"

momStr=`echo $mom | sed -e '{ s/\./_/g; }'`
# data directory. all data file go there and read from there, makes them easier to move around
path="/home/karavdin/datastorage/DPMoutput/all_"${momStr}"/"${folder}


# let's go

last=`expr $nEvts - 1`
##filenames=${path}/"names.log"
cat > ${path}/names.log <<EOF
EOF
cat > ${path}/namesMC.log <<EOF
EOF
cat > ${path}/QA_jobs.log <<EOF
EOF
FIRSTGEANE=$(stat -c%s ${path}/compMC_and_REC_withSignal0.root )
LEVEL=7
STANDARTGEANE=$((${FIRSTGEANE}*${LEVEL}/10))
FIRSTMC=$(stat -c%s ${path}/GenInfo_out_DPM_all_beamMom_${momStr}_0_inel_and_el.root )
STANDARTMC=$((${FIRSTMC}*${LEVEL}/10))
for start in `seq 0 $step $last`; do

if [ ! -r  ${path}/compMC_and_REC_withSignal${start}.root ]; then
echo "job_${start}.sh: ouput doesn't exist"  >> ${path}/QA_jobs.log 
fi
if [  -r  ${path}/compMC_and_REC_withSignal${start}.root ]; then
    echo -n "compMC_and_REC_withSignal${start}.root "  >> ${path}/names.log 
    FILESIZE=$(stat -c%s ${path}/compMC_and_REC_withSignal${start}.root )
    if [[ ${FILESIZE} -lt ${STANDARTGEANE} ]]; then
	echo "job_${start}.sh: ouput too small"  >> ${path}/QA_jobs.log 
    fi
fi
if [  -r  ${path}/GenInfo_out_DPM_all_beamMom_${momStr}_${start}_inel_and_el.root ]; then
    echo -n "GenInfo_out_DPM_all_beamMom_${momStr}_${start}_inel_and_el.root "  >> ${path}/namesMC.log 
    FILESIZE=$(stat -c%s ${path}/GenInfo_out_DPM_all_beamMom_${momStr}_${start}_inel_and_el.root )
    if [[ ${FILESIZE} -lt ${STANDARTMC} ]]; then
	echo "job_${start}.sh: ouput too small"  >> ${path}/QA_jobs.log 
    fi
fi

done


cat >/home/karavdin/myScripts/jobs/all_${momStr}/${folder}/jobHIMsum.sh <<EOF

#!/bin/bash
#
#PBS -N lumiJobSUM${momStr}
#PBS -j oe
#PBS -o ${path}/LOGfileSUM
#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
cd ${path}
ls -lah `cat ${path}/names.log` | awk '{ print $5}' &> ${path}/add_hists_QA.log
stat -c %s `cat ${path}/names.log` &> ${path}/add_hists_QA.log
stat -c %s `cat ${path}/namesMC.log` &> ${path}/add_histsMC_QA.log
${ROOTSYS}/bin/hadd -f sumAll.root ${filenames} `cat ${path}/names.log` &> ${path}/add_hists.log
if [ ! -r sum_MC_All.root ];then
rm sum_MC_All.root
fi
${ROOTSYS}/bin/hadd -f sum_MC_All.root ${filenamesMC} `cat ${path}/namesMC.log` &> ${path}/add_histsMC.log
cd /home/karavdin/pandaRoot12/build/bin
./sum_bkg_rec_mc_match -path ${path} -v 3 &> ${path}/summary_sumAll.log
./bkg_gen_info -path ${path} -mom ${mom} -fin 0
./bkg_gen_info -path ${path} -mom ${mom} -fin 1
exit 0
EOF
qsub /home/karavdin/myScripts/jobs/all_${momStr}/${folder}/jobHIMsum.sh
echo "Output files will be saved in $path"
#done