# example: goscrut.sh "ccbar" 7.53 10000 "all"

echo "Using FAIRLIBDIR="$FAIRLIBDIR
echo "Press any key to continue, ctrl-c to abort now"
read -n 1 -s


prefix="test"
dec="DPM"
mom=7.53
nEvts=10000
detopt="12345"
dmode=0
rmode=0
pid="All"

if test "$1" != ""; then
  prefix=$1
fi

if test "$2" != ""; then
  mom=$2
fi

if test "$3" != ""; then
  nEvts=$3
fi

if test "$4" != ""; then
  detopt=$4
fi

if test "$5" != ""; then
  pid="$5"
fi

for dmode in 0 1 2 3
do
  for rmode in 0 1 2 3 4 5
  do
    if test $dmode -eq 0; then
      if test $rmode -eq 0; then
    	dec="DecFiles/etac_2pi_to_KKpi0_2pi.dec"
      fi
      if test $rmode -eq 1; then
    	dec="DecFiles/etac_2pi0_to_KKpi0_2pi0.dec"
      fi
      if test $rmode -eq 2; then
    	dec="DecFiles/etac_pi0_to_KKpi0_pi0.dec"
      fi
      if test $rmode -eq 3; then
    	dec="DecFiles/etac_gam_to_KKpi0_gam.dec"
      fi
      if test $rmode -eq 4; then
    	dec="DecFiles/etac_eta_to_KKpi0_eta.dec"
      fi
      if test $rmode -eq 5; then
    	dec="DecFiles/etac_2eta_to_KKpi0_2eta.dec"
      if test $rmode -eq 6; then
    	dec="DecFiles/etac_2K_to_KKpi0_2K.dec"
      fi
    fi

    if test $dmode -eq 1; then
      if test $rmode -eq 0; then
    	dec="DecFiles/etac_2pi_to_KsKpi_2pi.dec"
      fi
      if test $rmode -eq 1; then
    	dec="DecFiles/etac_2pi0_to_KsKpi_2pi0.dec"
      fi
      if test $rmode -eq 2; then
    	dec="DecFiles/etac_pi0_to_KsKpi_pi0.dec"
      fi
      if test $rmode -eq 3; then
    	dec="DecFiles/etac_gam_to_KsKpi_gam.dec"
      fi
      if test $rmode -eq 4; then
    	dec="DecFiles/etac_eta_to_KsKpi_eta.dec"
      fi
      if test $rmode -eq 5; then
    	dec="DecFiles/etac_2eta_to_KsKpi_2eta.dec"
      fi
      if test $rmode -eq 6; then
    	dec="DecFiles/etac_2K_to_KsKpi_2K.dec"
      fi
    fi

    if test $dmode -eq 2; then
      if test $rmode -eq 0; then
    	dec="DecFiles/jpsi_2pi_to_2e_2pi.dec"
      fi
      if test $rmode -eq 1; then
    	dec="DecFiles/jpsi_2pi0_to_2e_2pi0.dec"
      fi
      if test $rmode -eq 2; then
    	dec="DecFiles/jpsi_pi0_to_2e_pi0.dec"
      fi
      if test $rmode -eq 3; then
    	dec="DecFiles/jpsi_gam_to_2e_gam.dec"
      fi
      if test $rmode -eq 4; then
    	dec="DecFiles/jpsi_eta_to_2e_eta.dec"
      fi
      if test $rmode -eq 5; then
    	dec="DecFiles/jpsi_2eta_to_2e_2eta.dec"
      fi
      if test $rmode -eq 6; then
    	dec="DecFiles/jpsi_2K_to_2e_2K.dec"
      fi
    fi

    if test $dmode -eq 3; then
      if test $rmode -eq 0; then
    	dec="DecFiles/jpsi_2pi_to_2mu_2pi.dec"
      fi
      if test $rmode -eq 1; then
    	dec="DecFiles/jpsi_2pi0_to_2mu_2pi0.dec"
      fi
      if test $rmode -eq 2; then
    	dec="DecFiles/jpsi_pi0_to_2mu_pi0.dec"
      fi
      if test $rmode -eq 3; then
    	dec="DecFiles/jpsi_gam_to_2mu_gam.dec"
      fi
      if test $rmode -eq 4; then
    	dec="DecFiles/jpsi_eta_to_2mu_eta.dec"
      fi
      if test $rmode -eq 5; then
    	dec="DecFiles/jpsi_2eta_to_2mu_2eta.dec"
      fi
      if test $rmode -eq 6; then
    	dec="DecFiles/jpsi_2K_to_2mu_2K.dec"
      fi
    fi

#    echo "$prefix (prefix)"
#    echo "$dec (decay (BOX, DPM or decay file)"
#    echo "$mom (momentum)"
#    echo "$nEvts (nEvents)"
#    echo "$detopt (detector option (12345))"
#    echo "$dmode (Decay Mode: 0: etac->K+K-pi0, 1: etac->KsK+pi-(+cc.) 2: J/psi->e+e-, 3:J/psi->mu+mu-)"
#    echo "$rmode (Recoil mode: 0: pi+pi-, 1: pi0pi0, 2: pi0, 3: gam, 4: eta 5: eta eta)"
#    echo "$pid (PID selection)"

	echo "Doing: qsub job_scrut_all.sge $prefix $dec $mom $nEvts $detopt $dmode $rmode $pid"
	qsub job_scrut_all.sge $prefix $dec $mom $nEvts $detopt $dmode $rmode $pid

  done
done