#!/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=100000
for opt in '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_2pi+2pi-pi0_OPT/${optStr}/
echo ${opt}
method1=PHSP
method2=OmegaEta
root -l -b -q simfast_Veryopt.C\("\"${PWD}/hc_OmegaEta_2pi+2pi-pi0_${momStr}_${optStr}\"","\"hc_5pi_OmegaEta_Ang.dec\"",$pbeam,${Nev},"\"pbarpSystem0\"",88880,"\"${opt}\""\)
root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_OmegaEta_2pi+2pi-pi0_${momStr}_${optStr}_fast.root\"",${Nev},$pbeam,"\"${method1}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_OmegaEta_2pi+2pi-pi0_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
mv ${PWD}/plotRes.root results_2pi+2pi-pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method1}.root
mv ${PWD}/plotRes.pdf results_2pi+2pi-pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method1}.pdf
root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_OmegaEta_2pi+2pi-pi0_${momStr}_${optStr}_fast.root\"",${Nev},$pbeam,"\"${method2}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_OmegaEta_2pi+2pi-pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
root -l -b -q  AnaAng.C\("\"${PWD}/hc_OmegaEta_2pi+2pi-pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
mv ${PWD}/plotRes.root results_2pi+2pi-pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method2}.root
mv ${PWD}/plotRes.pdf results_2pi+2pi-pi0_OPT/${optStr}/plotRes_sig_OmegaEta_${momStr}_${method2}.pdf
mv ${PWD}/plotAngRes.root results_2pi+2pi-pi0_OPT/${optStr}/plotAngRes_sig_OmegaEta_${momStr}_${method2}.root
mv ${PWD}/plotAngRes.pdf results_2pi+2pi-pi0_OPT/${optStr}/plotAngRes_sig_OmegaEta_${momStr}_${method2}.pdf
done
exit