#!/bin/bash
################################################################
## script for running sim&ana for fast sim
## tests were done for h_c hadronic decays!
## 15/04/2014
## A.Karavdina
################################################################

pbeam=5.61
Nev=1000000
for opt in 'FwdSpec EmcBar Drc Dsc MvdGem' 
#'FwdSpec Drc Dsc MvdGem' 'EmcBar Drc Dsc MvdGem' 'FwdSpec EmcBar Drc MvdGem' 'FwdSpec EmcBar Drc Dsc'

#'FwdSpec BarrelSpec' 'EmcBar FwdTrk STT MvdGem' 'EmcBar FwdTrk'  'EmcBar STT MvdGem' 'EmcBar FwdTrk STT' 'EmcBar EmcFwCap FwdTrk STT MvdGem' 'EmcBar EmcFwCap EmcBwCap FwdTrk STT MvdGem' 'EmcBar EmcFwCap EmcBwCap EmcFwd FwdTrk STT MvdGem' 'EmcBar EmcFwd FwdTrk STT MvdGem' 'EmcFwCap EmcBwCap EmcFwd FwdTrk STT MvdGem' 'EmcFwd FwdTrk STT MvdGem'  'EmcFwCap EmcBwCap FwdTrk STT MvdGem'
do 

momStr=`echo $pbeam | sed -e '{ s/\./_/g; }'`
optStr=`echo $opt | sed -e '{ s/\ /_/g; }'`
mkdir -p results_pi+pi-3pi0_OPT/${optStr}/
echo ${opt}
method1=PHSP
method2=OmegaEta
# #root -l -b -q simfast_Veryopt.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}\"","\"hc_pipi3pi0_OmegaEta.dec\"",$pbeam,${Nev},"\"pbarpSystem0\"",88880,"\"${opt}\""\)
# #root -l -b -q simfast_opt.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}\"","\"hc_pipi3pi0_phsp.dec\"",$pbeam,${Nev},"\"pbarpSystem0\"",88880,"\"${opt}\""\)
# root -l -b -q simfast_opt.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}\"","\"hc_pipi3pi0_OmegaEta.dec\"",$pbeam,${Nev},"\"pbarpSystem0\"",88880,"\"${opt}\""\)
# root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}_fast.root\"",${Nev},$pbeam,"\"${method1}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
root -l -b -q  AnaAng.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
mv ${PWD}/plotRes.root results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method1}.root
mv ${PWD}/plotRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method1}.pdf
mv ${PWD}/plotAngRes.root results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_OmegaEta_${momStr}_${method1}.root
mv ${PWD}/plotAngRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_OmegaEta_${momStr}_${method1}.pdf
# root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}_fast.root\"",${Nev},$pbeam,"\"${method2}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
root -l -b -q  AnaAng.C\("\"${PWD}/hc_OmegaEta_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
mv ${PWD}/plotRes.root results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method2}.root
mv ${PWD}/plotRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method2}.pdf
mv ${PWD}/plotAngRes.root results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_OmegaEta_${momStr}_${method2}.root
mv ${PWD}/plotAngRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_OmegaEta_${momStr}_${method2}.pdf

# root -l -b -q simfast_opt.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}\"","\"hc_pipi3pi0_phsp.dec\"",$pbeam,${Nev},"\"pbarpSystem0\"",88880,"\"${opt}\""\)
# root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}_fast.root\"",${Nev},$pbeam,"\"${method1}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
root -l -b -q  AnaAng.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
mv ${PWD}/plotRes.root results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_PHSP_${momStr}_${method1}.root
mv ${PWD}/plotRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_PHSP_${momStr}_${method1}.pdf
mv ${PWD}/plotAngRes.root results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_PHSP_${momStr}_${method1}.root
mv ${PWD}/plotAngRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_PHSP_${momStr}_${method1}.pdf
# root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}_fast.root\"",${Nev},$pbeam,"\"${method2}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
root -l -b -q  AnaAng.C\("\"${PWD}/hc_PHSP_pi+pi-3pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
mv ${PWD}/plotRes.root results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_PHSP_${momStr}_${method2}.root
mv ${PWD}/plotRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotRes_sig_PHSP_${momStr}_${method2}.pdf
mv ${PWD}/plotAngRes.root results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_PHSP_${momStr}_${method2}.root
mv ${PWD}/plotAngRes.pdf results_pi+pi-3pi0_OPT/${optStr}/plotAngRes_sig_PHSP_${momStr}_${method2}.pdf

done
exit