# file: marx_commands
#
# Use MARX to simulate Capella 1318 observation
#
# Do "source marx_commands" in a directory that already has:
#   Capella_HETG_08m10.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...
# do it first manually so comment it out here:
source marx_background


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

# where is marx:
set MARX_DIST_DIR = /nfs/wiwaxia/d4/ASC/src/marx-dist-devel

# 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 .
#
# Simulation Setup and Control
$MARX_DIST_DIR/bin/$ARCH/pset marx.par dNumRays=100000
$MARX_DIST_DIR/bin/$ARCH/pset marx.par ExposureTime=27000.0
$MARX_DIST_DIR/bin/$ARCH/pset marx.par OutputDir="marx_out"

# 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
# Select lower or upper line from these pairs:
##$MARX_DIST_DIR/bin/$ARCH/pset marx.par DetOffsetZ=0.00000
$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
$MARX_DIST_DIR/bin/$ARCH/pset marx.par SourceOffsetY=0.00000

# 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

# Grating Setup and Control
# set the exact HETG spacing...
# 3/3/00: simulation on web page was actually run with this
# rowland diameter - it is not the flight value, however.
# $MARX_DIST_DIR/bin/$ARCH/pset marx.par RowlandDiameter=8633.58
# Correct HETGS flight Rowland spacing is:
$MARX_DIST_DIR/bin/$ARCH/pset marx.par RowlandDiameter=8632.65
$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_DB0.35.fits
$MARX_DIST_DIR/bin/$ARCH/marx2fits marx_out Capella_offset.fits
# - - - - - - - - - - - - - - - - - -


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

rm -R total_out
##$MARX_DIST_DIR/bin/$ARCH/marxcat bkg_out marx_out total_out
$MARX_DIST_DIR/bin/$ARCH/marxcat bkg_out_offset marx_out total_out
##$MARX_DIST_DIR/bin/$ARCH/marx2fits total_out Capella_DB0.35_bkg.fits
$MARX_DIST_DIR/bin/$ARCH/marx2fits total_out Capella_offset_bkg.fits

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

unset MARX_DIST_DIR