#!/bin/tcsh
source ~/.login
cd /phenix/data07/hubert/pythia2/batch/
mkdir batch_$1
cd batch_$1
#
# 1) run pythia: -------------------------------------------------------
#
echo "First run Pythia: "
cp ../pp_c .
cp ../pp_c.in .
pp_c < pp_c.in > pp_c_$1.out
mv proot_xxx.root proot_$1.root
rm pp_c pp_c.in
#
# 2) convert to oscar file: --------------------------------------------
#
echo "Pythia done - start read_ntuple: "
cp ../read_ntuple.C .
ln -sf proot_$1.root pythia.root
root -b <<EOF
  .L /phenix/data07/hubert/pythia2/batch/read_ntuple.C
  read_ntuple t
  t.Loop()
EOF
mv ccbar_d_mu.oscar ccbar_d_mu_$1.oscar
rm read_ntuple.C
#
# now run pisa: -----------------------------------------------------------
#
echo "oscar file made - run PISA: "
cp ../pisa.kumac .
cp ../gffgo.dat .
cp ../phnx.par .
cp ../pisa.input .
ln -sf ccbar_d_mu_$1.oscar oscar.input
pisa <pisa.input > pisa_$1.out
rm pisa.kumac gffgo.dat phnx.par pisa.input fort.51 fort.66 
rm finishPisa svxPISA.par gintphnx.hbk last.kumac
pisaRootRead >& /dev/null
mv ancsvx.root bncsvx.root
rm anc*
mv bncsvx.root ancsvx_$1.root
rm PISAEvent.root
#
#------------------------------------------------------------------------
#
# echo "PISA done. "
# end