#!/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_BeamProf/mom_${momStr}/
pathBOX=/panda/myResults/BOXoutputMarch2014_BeamProf_newXYCut_newMCut_2/mom_${momStr}/
pathDPM=/panda/myResults/DPMoutputMarch2014_BeamProf_newXYCut_newMCut_2/mom_${momStr}/
pathBOX2=/panda/myResults/BOXoutputMarch2014_BeamProf_newXYCut/mom_${momStr}/
pathDPM2=/panda/myResults/DPMoutputMarch2014_BeamProf_newXYCut/mom_${momStr}/
pathBOX0=/panda/myResults/BOXoutpitMarch2014_BeamProf/mom_${momStr}/
pathDPM0=/panda/myResults/DPMoutputMarch2014_BeamProf/mom_${momStr}/

#pathDPM=/panda/myResults/DPMoutputMarch2014_BeamProf/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 0.7 1.1 1.4 1.75 2.1 2.45 2.8 3.15 3.5 4.2 7
#for thcutVal in 0.5 0.9 1.3 1.7 2.1 2.5 2.9 3.3 3.7 4.1 4.5 4.9 5.3 5.7 6.1 6.5 6.9 10000
#for thcutVal in 0.7 1.4 2.1 2.8 3.5 4.2 4.9 5.6 6.3 7.0 7.7 8.4 9.1 9.8 10.5 14 70 

#for thcutVal in 0.7
#for thcutVal in 0.7 1.4 2.1 2.8 3.5 
#for thcutVal in 5000000
#for thcutVal in -300.1
#for thcutVal in  0.7 1.4  
for thcutVal in -100
#1.7 1.9 2.3 2.5 4.9 5.6 6.3 7.0 14 70 
#for thcutVal in 0.7 1.4 2.1 2.8 3.5 70 5000000
#for thcutVal in  
do 
    ./createLumiFitData -n 0 -m 1.5 -t er -p ${pathBOX} -a ${thcutVal} #acc & res
    ./determineResolution -d ${pathBOX} &> ${pathDPM}\log_determineResolution
    ./plotResolutionParametrization -d  ${pathBOX} -l 0
    ./createLumiFitData -n 0  -m 1.5 -t a -p ${pathDPM} -c 61.919965  -a ${thcutVal} #data 1.5 GeV
    ### and now let's fit the luminosity
    ./runLumiFit -d ${pathDPM} -a ${pathBOX} &> ${pathDPM}\log_runLumiFit
    ./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

    ./createLumiFitData -n 0 -m 1.5 -t er -p ${pathBOX0} -a ${thcutVal} #acc & res
    ./determineResolution -d ${pathBOX0} #&> ${pathDPM0}\log_determineResolution
    ./plotResolutionParametrization -d  ${pathBOX0} -l 0
    ./createLumiFitData -n 0  -m 1.5 -t a -p ${pathDPM0} -c 61.919965  -a ${thcutVal} #data 1.5 GeV
    ### and now let's fit the luminosity
    ./runLumiFit -d ${pathDPM0} -a ${pathBOX0} #&> ${pathDPM0}\log_runLumiFit
    ./plotLumiFitResults ${pathDPM0}
    thcutValStr=`echo $thcutVal | sed -e '{ s/\./_/g; }'`
    mv fitresults_overview.pdf ${pathDPM0}/fitresults_overview_Thcut_${thcutValStr}.pdf
    mv lumifit_reco.pdf ${pathDPM0}/lumifit_reco_Thcut_${thcutValStr}.pdf
    mv resolution.pdf ${pathDPM0}/resolution_Thcut_${thcutValStr}.pdf
    mv ${pathBOX0}/resolution_params_1.root ${pathDPM0}/resolution_params_1_${thcutValStr}.root
    mv ${pathBOX0}/resolution_params_0.root ${pathDPM0}/resolution_params_0_${thcutValStr}.root
    mv ${pathBOX0}/lmd_acc_data.root ${pathDPM0}/lmd_acc_data_${thcutValStr}.root
    mv ${pathDPM0}/lmd_data.root ${pathDPM0}/lmd_data_${thcutValStr}.root

 ./createLumiFitData -n 0 -m 1.5 -t er -p ${pathBOX2} -a ${thcutVal} #acc & res
    ./determineResolution -d ${pathBOX2} #&> ${pathDPM2}\log_determineResolution
    ./plotResolutionParametrization -d  ${pathBOX2} -l 0
    ./createLumiFitData -n 0  -m 1.5 -t a -p ${pathDPM2} -c 61.919965  -a ${thcutVal} #data 1.5 GeV
    ### and now let's fit the luminosity
    ./runLumiFit -d ${pathDPM2} -a ${pathBOX2} #&> ${pathDPM2}\log_runLumiFit
    ./plotLumiFitResults ${pathDPM2}
    thcutValStr=`echo $thcutVal | sed -e '{ s/\./_/g; }'`
    mv fitresults_overview.pdf ${pathDPM2}/fitresults_overview_Thcut_${thcutValStr}.pdf
    mv lumifit_reco.pdf ${pathDPM2}/lumifit_reco_Thcut_${thcutValStr}.pdf
    mv resolution.pdf ${pathDPM2}/resolution_Thcut_${thcutValStr}.pdf
    mv ${pathBOX2}/resolution_params_1.root ${pathDPM2}/resolution_params_1_${thcutValStr}.root
    mv ${pathBOX2}/resolution_params_0.root ${pathDPM2}/resolution_params_0_${thcutValStr}.root
    mv ${pathBOX2}/lmd_acc_data.root ${pathDPM2}/lmd_acc_data_${thcutValStr}.root
    mv ${pathDPM2}/lmd_data.root ${pathDPM2}/lmd_data_${thcutValStr}.root

done