#!/bin/sh

#Needs to be set outside:
# number of events for each thread
#  LIT_NEVENTS=1000
# file id
#  XXXX=0001
# working directory
#  MAINDIRROOT=/hera/hades/user/ekrebs/cbm/dielectron/auau/centr
# Collision energy
#  LIT_ENERGY=25gev
# Embedded Pluto particles
#  PLUTOREAC=16
# MVD geometry
#  MVD=mvd4
# Field mape scale factor:
#  FIELDSCALE=1.0

cbmroot_version=JUL13
. /hera/cbm/users/$(whoami)/cbmroot/build_${cbmroot_version}/config.sh

export MYBUILDDIR=/hera/cbm/users/${whoami}/cbmroot

. ./pluto_table.sh

#Macro will be running via script
export LIT_SCRIPT=yes

#Reconstruction parameters
export LIT_GLOBAL_TRACKING_TYPE=branch
export LIT_STS_HITPRODUCER_TYPE=real
export LIT_TRD_HITPRODUCER_TYPE=smearing
export LIT_MUCH_HITPRODUCER_TYPE=advanced

#Normalization for efficiency
export LIT_NORM_STS_POINTS=4
export LIT_NORM_TRD_POINTS=8
export LIT_NORM_MUCH_POINTS=10
export LIT_NORM_TOF_POINTS=1
export LIT_NORM_TRD_HITS=8
export LIT_NORM_MUCH_HITS=10
export LIT_NORM_TOF_HITS=1


#Simulation parameters
#--------------------------------------------------

# If "yes" than UrQMD will be used as background
export LIT_URQMD=yes

# If "yes" PLUTO particles will be embedded
export LIT_PLUTO=yes

export LIT_UNIGEN=yes

PLUTOPARTICLE=${arr_dir[$PLUTOREAC]}
PLUTODECAY=${arr_decay[$PLUTOREAC]}

# Used for weight calculation in analysis macro
if [ "$arr_dalitz[$PLUTOREAC]" == 1 ] ; then
  export LIT_PLUTO_PARTICLE=${arr_pluto[$PLUTOREAC]}dalitz
else
  export LIT_PLUTO_PARTICLE=${arr_pluto[$PLUTOREAC]}
fi

# number of embedded electrons
export LIT_NOF_ELECTRONS=0

# number of embedded positrons
export LIT_NOF_POSITRONS=0

# If "yes" than primary electrons will be generated
export LIT_ELECTRONS=no

# number of embedded muons from FairBoxGenerator
export LIT_NOF_MUONS_PLUS=0
# number of embedded muons from FairBoxGenerator
export LIT_NOF_MUONS_MINUS=0
# number of embedded electrons from FairBoxGenerator
export LIT_NOF_ELECTRONS=0
# number of embedded positrons from FairBoxGenerator
export LIT_NOF_POSITRONS=0
# number of embedded pions from FairBoxGenerator
export LIT_NOF_PIONS_PLUS=0
# number of embedded pions from FairBoxGenerator
export LIT_NOF_PIONS_MINUS=0
# number of embedded J/Psi decayed to mu+ and mu-
export LIT_NOF_JPSI_TO_MUONS=0
# number of embedded J/Psi decayed to e+ and e-
export LIT_NOF_JPSI_TO_ELECTRONS=0
# number of generated Au ions
export LIT_NOF_AU_IONS=0

#---------------------------------------------------


#Geometries
#---------------------------------------------------
#export GEOMPATH=${MYBUILDDIR}/geometry
if [ ${MVD} = "nomvd" ]; then
   export LIT_MVD_GEOM=
elif [ ${MVD} = "mvd2" ]; then
   export LIT_MVD_GEOM=mvd/mvd_v07a.geo
elif [ ${MVD} = "mvd3" ]; then
   export LIT_MVD_GEOM=mvd/mvd_3_Nov_2012.geo
elif [ ${MVD} = "mvd3_051020" ]; then
   export LIT_MVD_GEOM=mvd/mvd3_5_10_20.geo
elif [ ${MVD} = "mvd4" ] ; then
   export LIT_MVD_GEOM=mvd/mvd_Nov_2012.geo
elif [ ${MVD} = "mvd4xl" ] ; then
   export LIT_MVD_GEOM=mvd/mvd_v10a.geo
elif [ ${MVD} = "mvd5" ]; then
   export LIT_MVD_GEOM=mvd/mvd_5_Nov_2012.geo
else
   echo "---- ERROR ----- ${MVD} is not a valid option for MVD geometry. -----"
   export LIT_MVD_GEOM=
fi

export LIT_CAVE_GEOM=cave.geo
export LIT_TARGET_GEOM=target/target_au_025mu.geo
export LIT_PIPE_GEOM=pipe/pipe_standard_short.geo
export LIT_SHIELD_GEOM=
export LIT_STS_GEOM=sts/sts_v12b.geo.root
export LIT_MUCH_GEOM=                      
export LIT_RICH_GEOM=rich/rich_v08a.geo
export LIT_TRD_GEOM=trd/trd_v13g.geo.root
export LIT_TOF_GEOM=tof/tof_v13b.root
export LIT_ECAL_GEOM=
export LIT_FIELD_MAP=field_v12b
export LIT_MAGNET_GEOM=magnet/magnet_v12a.geo
#---------------------------------------------------

#Digitizers
export LIT_STS_DIGI=${VMCWORKDIR}/parameters/sts/sts_v12b_std.digi.par
export LIT_TOF_DIGI=${VMCWORKDIR}/parameters/tof/tof_v13b.digi.par
export LIT_TRD_DIGI=${VMCWORKDIR}/parameters/trd/trd_v13g.digi.par
export LIT_MUCH_DIGI=${VMCWORKDIR}/parameters/much/much_v12b.digi.root
export LIT_STS_MAT_BUDGET_FILE=${VMCWORKDIR}/parameters/sts/sts_matbudget_v12b.root

#Field map
export LIT_FIELD_MAP_SCALE=${FIELDSCALE}
export LIT_FIELD_Z=50.0
export LIT_FIELD_SYM_TYPE=3

#If "yes" DELTA electrons will be embedded
#export LIT_DELTA=no
#export DELTA_FILE=/lustre/cbm/user/ebelolap/aug11/sep12/${LIT_ENERGY}/${FIELDDIR}/deltasource/mc.delta.root
if [ ${LIT_DELTA} = "yes" ]; then
    RECO_DIR_ADD=${RECO_DIR_ADD}_delta_${LIT_INTERACTION_RATE}khz
fi

export LIT_PION_MISIDENTIFICATION_LEVEL=-1.0
export LIT_USE_MC_MOMENTUM=no


# Set directories
MAINDIR=${MAINDIRROOT}   # mainroot set outside
	
FIELDDIR=${LIT_FIELD_MAP_SCALE}field_${LIT_FIELD_Z}z

MVDDIR=${MVD}

MAINDIR=${MAINDIR}/${FIELDDIR}_${MVDDIR}

if [[ "${LIT_PLUTO}" = "yes" && "${LIT_URQMD}" = "no" ]] ; then
  EVTGENDIR=pluto_only/${PLUTOPARTICLE}/${PLUTODECAY}
fi
if [[ "${LIT_PLUTO}" = "no" && "${LIT_URQMD}" = "yes" ]] ; then
  EVTGENDIR=urqmd_only
fi
if [[ "${LIT_PLUTO}" = "yes" && "${LIT_URQMD}" = "yes" ]] ; then
  EVTGENDIR=${PLUTOPARTICLE}/${PLUTODECAY}
fi

MAINDIR=${MAINDIR}/${EVTGENDIR}

export LIT_WORKDIR=${MAINDIR}

export LIT_MC_FILE=${MAINDIR}/sim/mc.${XXXX}.root
export LIT_PAR_FILE=${MAINDIR}/par/param.${XXXX}.root
export LIT_GLOBAL_RECO_FILE=${MAINDIR}/reco${RECO_DIR_ADD}/reco.${XXXX}.root
export LIT_MATCH_FILE=${MAINDIR}/reco${RECO_DIR_ADD}/match.${XXXX}.root
export LIT_GLOBAL_HITS_FILE=${MAINDIR}/reco/hits.${XXXX}.root
export LIT_GLOBAL_TRACKS_FILE=${MAINDIR}/reco/tracks.${XXXX}.root
export LIT_RESULT_DIR=${MAINDIR}/qa
export LIT_ANALYSIS_FILE=${MAINDIR}/analysis${RECO_DIR_ADD}/analysis.${XXXX}.root
export LIT_DILEP_ANALYSIS_FILE=${MAINDIR}/analysis.mc_momentum.${XXXX}.root 
export LIT_MC_DIEL_ANALYSIS_FILE=${MAINDIR}/analysis${RECO_DIR_ADD}/analysis.mc.diel.${XXXX}.root 
export LIT_RECO_DIEL_ANALYSIS_FILE=${MAINDIR}/analysis${RECO_DIR_ADD}/analysis.reco.diel.${XXXX}.root 
export LIT_GAMMA_DIEL_ANALYSIS_FILE=${MAINDIR}/analysis${RECO_DIR_ADD}/analysis.gamma.diel.${XXXX}.root 
export LIT_URQMD_FILE=/hera/cbm/prod/gen/urqmd/auau/${LIT_ENERGY}/centr/urqmd.auau.${LIT_ENERGY}.centr.0${XXXX}.root
export LIT_PLUTO_FILE=/hera/cbm/users/ekrebs/pluto/cktA/${LIT_ENERGY}/${PLUTOPARTICLE}/${PLUTODECAY}/pluto.auau.${LIT_ENERGY}.${PLUTOPARTICLE}.${PLUTODECAY}.${XXXX}.root
export LIT_DELTA_FILE=/hera/cbm/users/ekrebs/prod/delta/delta.mc.root
#export LIT_DELTA_FILE=/hera/cbm/users/ekrebs/prod/delta/auau/${FIELDDIR}/mvd.mc.delta.${XXXX}.root
#---------------------------------------------------

# Make missing directories
if [ ! -d ${MAINDIR}/reco${RECO_DIR_ADD} ]
then
   mkdir -p  ${MAINDIR}/reco${RECO_DIR_ADD}
fi

if [ ! -d ${MAINDIR}/sim ]
then
   mkdir -p  ${MAINDIR}/sim
fi

if [ ! -d ${MAINDIR}/par ]
then
   mkdir -p  ${MAINDIR}/par
fi

if [ ! -d ${MAINDIR}/analysis${RECO_DIR_ADD} ]
then
   mkdir -p  ${MAINDIR}/analysis${RECO_DIR_ADD}
fi

#if [ ! -d ${LIT_RESULT_DIR} ]
#then
#   mkdir -p  ${LIT_RESULT_DIR}
#fi
#---------------------------------------------------