# file: marx_commands
#
# Use MARX to simulate a Capella observation (e.g. like Obsid 1318)
#
# Do "source marx_commands" in a directory that already has:
#   Capella_HETG_08m10.tbl
# and (only needed if background is simulated too):
#   ACIS_background.tbl
#
# Note that this file copies marx.par fresh from the marx distribution
# and then changes parameters in it as needed.
#-----------------------------------------------------------------


# - - - - - - - - - - - - - - - - - -
# - - - - - - - - - - - - - - - - - -
# BACKGROUND
#
# run a new background for this simulation...
# (Could do it first manually and comment it out here)
source marx_background


# - - - - - - - - - - - - - - - - - -
# - - - - - - - - - - - - - - - - - -
# CAPELLA
#

# Get access to MARX:
# where is marx:
set MARX_DIST_DIR = /nfs/wiwaxia/d4/ASC/src/marx-dist
# where is marx data:
setenv MARX_DATA_DIR $MARX_DIST_DIR/marx/data

# copy the marx.par file here
cp $MARX_DIST_DIR/marx.par .
# and the marxasp file:
cp $MARX_DIST_DIR/marxasp.par .

#
# Simulation Setup and Control
$MARX_DIST_DIR/bin/$ARCH/pset marx.par dNumRays=100000
$MARX_DIST_DIR/bin/$ARCH/pset marx.par ExposureTime=30000.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par OutputDir="marx_out"
# Same output file for marxasp:
$MARX_DIST_DIR/bin/$ARCH/pset marxasp.par MarxDir="marx_out"
$MARX_DIST_DIR/bin/$ARCH/pset marxasp.par OutputFile="acis_asol.fits"

# Science Instrument Setup and Control
$MARX_DIST_DIR/bin/$ARCH/pset marx.par GratingType="HETG"
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DetectorType="ACIS-S"
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DetIdeal="no"
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DetOffsetX=0.00000
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DetOffsetY=0.00000
# Offset to 3.0 mm is nominal recommendation to reduce CTI
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DetOffsetZ=3.00000

# Source Spectral Parameters
# for Capella, flux expected to be 0.1105
# adjust it a bit to better agree w/data (total dispersed +/-1 MEG events)
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SourceFlux=0.077
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SpectrumType="FILE"
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SpectrumFile="Capella_HETG_08m10.tbl"

# Source Spatial Parameters
# on-axis point source at infinity
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SourceType="POINT"
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SourceDistance=0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SourceOffsetZ=0.00000
# Offset point to move away from quadrant boundary...
# into node 1 of S3.
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SourceOffsetY=-0.33000

# Dither ON !
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherModel="INTERNAL"
# Add additional blur to get accurate Sky X,Y blur...
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherBlur=0.35
# Nominal HETGS dither parameters:
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherAmp_RA=8.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherAmp_Dec=8.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherAmp_Roll=0.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherPeriod_RA=1000.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherPeriod_Dec=707.0
# Source location on sky
$MARX_DIST_DIR/bin/$ARCH/pset marx.par Pointing_RA=79.2
$MARX_DIST_DIR/bin/$ARCH/pset marx.par Pointing_Dec=46.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherRoll=98.6

# Grating Setup and Control
$MARX_DIST_DIR/bin/$ARCH/pset marx.par UseGratingEffFiles="yes"

# - - - - - - - - - - - - - - - - - -
# finally, execute marx:
$MARX_DIST_DIR/bin/$ARCH/marx @@marx.par
# and do marx2fits on just this part of simulation
$MARX_DIST_DIR/bin/$ARCH/marx2fits marx_out Capella_nobkg.fits
# and create the asol file:
$MARX_DIST_DIR/bin/$ARCH/marxasp
# - - - - - - - - - - - - - - - - - -


# - - - - - - - - - - - - - - - - - -
# - - - - - - - - - - - - - - - - - -
# ASSEMBLE the COMPLETE OUTPUT FITS FILE
# cat the source and background files together,
# and do marx2fits

# remove the final output directory...
rm -R total_out
# combine the simulation with the background simulation
$MARX_DIST_DIR/bin/$ARCH/marxcat marx_out bkg_out total_out
# convert to FITS file:
$MARX_DIST_DIR/bin/$ARCH/marx2fits total_out Capella_wbkg.fits

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

unset MARX_DIST_DIR