#!/bin/bash
################################################################
## script for running sim&ana for fast sim
## tests were done for h_c hadronic decays!
## 15/04/2014
## A.Karavdina
################################################################
#for pbeam in 5.3119 5.4119 5.5119 5.7119 5.8119 5.9119 5.6119
#for pbeam in 4.6 4.9 5.2 5.5 5.612 5.7 6.0 6.3 6.6 5.3 5.4 5.65 5.8
#for pbeam in 5.1 5.3 5.5 5.612 5.7 5.9 6.1
#for pbeam in 5.55 5.61 5.65 5.7 5.75 5.8 
#pbeam=5.61
#for pbeam in 5.58219 5.6035 5.60617 5.6075 5.60883 5.61017 5.6115 5.61417 5.63553 5.87802
for pbeam in 5.61
do 
momStr=`echo $pbeam | sed -e '{ s/\./_/g; }'`
#A
opt='FwdSpec EmcBar Drc Dsc MvdGem'
#B
#opt='FwdSpec Drc Dsc MvdGem'
#C
#opt='EmcBar Drc Dsc MvdGem'
#D
#opt='FwdSpec EmcBar Drc MvdGem'
#E
#opt='FwdSpec EmcBar Drc Dsc'
optStr=`echo $opt | sed -e '{ s/\ /_/g; }'`
echo ${opt}
method1=PHSP
method2=OmegaEta
root -l -b -q simfast_opt.C\("\"${PWD}/hc_phsp_2pi+2pi-pi0_${momStr}_${optStr}\"","\"hc_5pi_phsp.dec\"",$pbeam,100000,"\"pbarpSystem0\"",88880,"\"${opt}\""\)
# root -l -b -q  ana_hc_task.C\("\"${PWD}/pbarp${ppbar}_hc_phsp_2pi+2pi-pi0_${momStr}_${optStr}_fast.root\"",100000,$pbeam,"\"${method1}\""\)
# root -l -b -q  AnaAna.C\("\"${PWD}/pbarp${ppbar}_hc_phsp_2pi+2pi-pi0_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
# mv plotRes.root results/plotRes_pbarp${ppbar}_sig_${optStr}_${momStr}_${method1}.root
# mv plotRes.pdf results/plotRes_pbarp${ppbar}_sig_${optStr}_${momStr}_${method1}.pdf
root -l -b -q  ana_hc_task.C\("\"${PWD}/hc_phsp_2pi+2pi-pi0_${momStr}_${optStr}_fast.root\"",100000,$pbeam,"\"${method2}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/hc_phsp_2pi+2pi-pi0_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
mv plotRes.root results/plotRes_phsp_noCuts_sig_${optStr}_${momStr}_${method2}.root
mv plotRes.pdf results/plotRes_phsp_noCuts_sig_${optStr}_${momStr}_${method2}.pdf


root -l -b -q simfast_opt.C\("\"${PWD}/dpm_${momStr}_${optStr}\"","\"DPM\"",$pbeam,100000,"\"dpm\"",11,"\"${opt}\""\)
# root -l -b -q  ana_hc_task.C\("\"${PWD}/dpm_${momStr}_${optStr}_fast.root\"",100000,$pbeam,"\"${method1}\""\)
# root -l -b -q  AnaAna.C\("\"${PWD}/dpm_${momStr}_${optStr}_fast.root_${method1}_taskana.root\"","\"${PWD}\""\)
# mv plotRes.root results/plotRes_bkg_${optStr}_${momStr}_${method1}.root
# mv plotRes.pdf results/plotRes_bkg_${optStr}_${momStr}_${method1}.pdf
root -l -b -q  ana_hc_task.C\("\"${PWD}/dpm_${momStr}_${optStr}_fast.root\"",100000,$pbeam,"\"${method2}\""\)
root -l -b -q  AnaAna.C\("\"${PWD}/dpm_${momStr}_${optStr}_fast.root_${method2}_taskana.root\"","\"${PWD}\""\)
mv plotRes.root results/plotRes_bkg_noCuts_${optStr}_${momStr}_${method2}.root
mv plotRes.pdf results/plotRes_bkg_noCuts_${optStr}_${momStr}_${method2}.pdf

done
exit