#!/bin/bash 
################################################################
## script for test cuts from rec as influence on determination of luminosity
## BOX data and DPM data should be generated before
## 26/03/2014 1st version
## A.Karavdina
################################################################
echo Enter momentum from "1.5, 4.06, 8.9, 11.91, 15 [GeV/c]":
read mom
echo "Beam mometum is: $mom!"
cd ${VMCWORKDIR}/macro/lmd/
pbeam=$mom
momStr=`echo $pbeam | sed -e '{ s/\./_/g; }'`

pathBOX=/panda/myResults/BOXoutpitMarch2014_pointBeam/mom_${momStr}/
pathDPM=/panda/myResults/DPMoutputMarch2014_pointBeam/mom_${momStr}/

#thcutVal=1.0

### preparation for the luminosity fit
cd ${VMCWORKDIR}/buildPanda/bin
#./createLumiFitData -t 0 -p ${pathDPM} -g 0.102459016393442626 ### 15 GeV/c, -g= 1./sigma_tot_el

#for thcutVal in 5000000 0.07 0.14 0.21 0.28 0.35 0.42 0.49 0.56 0.63 0.70
for thcutVal in 5000000
do 
  #  ./createLumiFitData -n 0  -m 1.5 -t a -p ${pathDPM} -c 61.919965  -a ${thcutVal} #data 1.5 GeV
    ./createLumiFitData -n 10000000  -m 15 -t a -p ${pathDPM} -c 9.76  -a ${thcutVal} # &> ${pathDPM}\log_dpm_data_create_${thcutValStr} #data 15 GeV 
    ./createLumiFitData -n 10000000 -m 15 -t er -p ${pathBOX} -a ${thcutVal} # &> ${pathDPM}\log_acc_res_data_create_${thcutValStr} #acc & res
    ./determineResolution -d ${pathBOX} #&> ${pathDPM}\log_determineResolution_${thcutValStr}
    ./plotResolutionParametrization -d  ${pathBOX} -l 0
    ### and now let's fit the luminosity
    ./runLumiFit -d ${pathDPM} -a ${pathBOX} &> ${pathDPM}\log_runLumiFit_${thcutValStr}
    ./plotLumiFitResults ${pathDPM}
    thcutValStr=`echo $thcutVal | sed -e '{ s/\./_/g; }'`
    mv fitresults_overview.pdf ${pathDPM}/fitresults_overview_Thcut_${thcutValStr}.pdf
    mv lumifit_reco.pdf ${pathDPM}/lumifit_reco_Thcut_${thcutValStr}.pdf
    mv resolution.pdf ${pathDPM}/resolution_Thcut_${thcutValStr}.pdf
    mv ${pathBOX}/resolution_params_1.root ${pathDPM}/resolution_params_1_${thcutValStr}.root
    mv ${pathBOX}/resolution_params_0.root ${pathDPM}/resolution_params_0_${thcutValStr}.root
    mv ${pathBOX}/lmd_acc_data.root ${pathDPM}/lmd_acc_data_${thcutValStr}.root
    mv ${pathDPM}/lmd_data.root ${pathDPM}/lmd_data_${thcutValStr}.root
done