# file: narrow_commands
#
# Do a simulation with a bunch of narrow lines.
#
# Do "source marx_commands" in a directory that already has:
#   Narrow_lines.tbl
#
# Note that this file copies marx.par fresh from the marx distribution
# and then changes parameters in it as needed.
#-----------------------------------------------------------------


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

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

# one time copy the marx.par file here
##cp $MARX_DIST_DIR/marx.par dist_marx.par
# then change the Rowland spacing to 8637.00, save in mod_marx.par
# and in future do this ot start with a fresh .par file:
cp mod_marx.par 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=100000.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="Narrow_lines.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...
# Old standby value (gives zo FWHM 500um of 38.5 um):
##$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherBlur=0.35
# New value based on aspect rmd diameter of 0.2 instead of 0.34 arc sec:
# (gives FWHM E<2keV of 35.56 um)
## $MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherBlur=0.30
# New value based on interpolating to desired FWHM to get 37.9 um:
##$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherBlur=0.34
# gave 38.23 ...
# go a little smaller:
# And use RC = 8637.0 !
$MARX_DIST_DIR/bin/$ARCH/pset marx.par DitherBlur=0.33
# 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
$MARX_DIST_DIR/bin/$ARCH/marx2fits marx_out Narrow33_100ks.fits
# and create the asol file:
$MARX_DIST_DIR/bin/$ARCH/marxasp
cp acis_asol.fits Narrow33_100ks_asol.fits
cp marx.par Narrow33_100ks_marx.par
# - - - - - - - - - - - - - - - - - -

unset MARX_DIST_DIR