hera=$VMCWORKDIR"/macro/ccbargsi"
_target=$hera"/data/"

prefix="dpm"
dec="DPM"
mom=4.5
nEvts=10000
res="pbarpSystem0"
detopt="12345"
dmode=0
rmode=0
pid="All"

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

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

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

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

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

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

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

# -----
if test "$dec" = "DecFiles/etac_2K_to_KKpi0_2K.dec"; then
  dmode=0
  rmode=6
fi
if test "$dec" = "DecFiles/etac_2K_to_KsKpi_2K.dec"; then
  dmode=1
  rmode=6
fi

if test "$dec" = "DecFiles/etac_2eta_to_KKpi0_2eta.dec"; then
  dmode=0
  rmode=5
fi
if test "$dec" = "DecFiles/etac_2eta_to_KsKpi_2eta.dec"; then
  dmode=1
  rmode=5
fi

if test "$dec" = "DecFiles/etac_2pi0_to_KKpi0_2pi0.dec"; then
  dmode=0
  rmode=1
fi
if test "$dec" = "DecFiles/etac_2pi0_to_KsKpi_2pi0.dec"; then
  dmode=1
  rmode=1
fi

if test "$dec" = "DecFiles/etac_2pi_to_KKpi0_2pi.dec"; then
  dmode=0
  rmode=0
fi
if test "$dec" = "DecFiles/etac_2pi_to_KsKpi_2pi.dec"; then
  dmode=1
  rmode=0
fi

if test "$dec" = "DecFiles/etac_eta_to_KKpi0_eta.dec"; then
  dmode=0
  rmode=4
fi
if test "$dec" = "DecFiles/etac_eta_to_KsKpi_eta.dec"; then
  dmode=1
  rmode=4
fi

if test "$dec" = "DecFiles/etac_gam_to_KKpi0_gam.dec"; then
  dmode=0
  rmode=3
fi
if test "$dec" = "DecFiles/etac_gam_to_KsKpi_gam.dec"; then
  dmode=1
  rmode=3
fi
if test "$dec" = "DecFiles/etac_pi0_to_KKpi0_pi0.dec"; then
  dmode=0
  rmode=2
fi
if test "$dec" = "DecFiles/etac_pi0_to_KsKpi_pi0.dec"; then
  dmode=1
  rmode=2
fi

# -----
if test "$dec" = "DecFiles/jpsi_2K_to_2e_2K.dec"; then
  dmode=2
  rmode=6
fi
if test "$dec" = "DecFiles/jpsi_2K_to_2mu_2K.dec"; then
  dmode=3
  rmode=6
fi

if test "$dec" = "DecFiles/jpsi_2eta_to_2e_2eta.dec"; then
  dmode=2
  rmode=5
fi
if test "$dec" = "DecFiles/jpsi_2eta_to_2mu_2eta.dec"; then
  dmode=3
  rmode=5
fi


if test "$dec" = "DecFiles/jpsi_2pi0_to_2e_2pi0.dec"; then
  dmode=2
  rmode=1
fi
if test "$dec" = "DecFiles/jpsi_2pi0_to_2mu_2pi0.dec"; then
  dmode=3
  rmode=1
fi

if test "$dec" = "DecFiles/jpsi_2pi_to_2e_2pi.dec"; then
  dmode=2
  rmode=0
fi
if test "$dec" = "DecFiles/jpsi_2pi_to_2mu_2pi.dec"; then
  dmode=3
  rmode=0
fi
if test "$dec" = "DecFiles/jpsi_eta_to_2e_eta.dec"; then
  dmode=2
  rmode=4
fi
if test "$dec" = "DecFiles/jpsi_eta_to_2mu_eta.dec"; then
  dmode=3
  rmode=4
fi

if test "$dec" = "DecFiles/jpsi_gam_to_2e_gam.dec"; then
  dmode=2
  rmode=3
fi
if test "$dec" = "DecFiles/jpsi_gam_to_2mu_gam.dec"; then
  dmode=3
  rmode=3
fi

if test "$dec" = "DecFiles/jpsi_pi0_to_2e_pi0.dec"; then
  dmode=2
  rmode=2
fi
if test "$dec" = "DecFiles/jpsi_pi0_to_2mu_pi0.dec"; then
  dmode=3
  rmode=2
fi

prefix=$prefix"_"$dmode"_"$rmode
outprefix=$_target$prefix

echo "$prefix (prefix)"
echo "$dec (decay (BOX, DPM or decay file)"
echo "$mom (momentum)"
echo "$nEvts (nEvents)"
echo "$res (resonance (for EvtGen))"
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 "Fast simulation..."
root -l -q -b $hera"/"simfast_opt.C\(\"$outprefix\",\"$dec\",$mom,$nEvts,\"$res\",11,\"$detopt\"\)  #&> $outprefix"_fast.log"
echo "Combinatorics and ntuples..."
root -l -q -b $hera"/"ana_etac_task.C\(\"$outprefix"_fast.root"\",$nEvts,$mom,$dmode,$rmode,\"$pid\"\)  #&> $outprefix"_ntp.log"
echo "Plots..."
root -l -q -b $hera"/"plotntuplesfromfile.C\(\"$outprefix"_fast_"$pid"_ntp.root"\"\) #&> $outprefix"_plots.log"