file: marx_hetg_plaw.txt

Use MARX to simulate an HETG observation of a powerlaw w/lines.

# - - - Setup and Run MARX

Roughly these steps are the same as for the ACIS example, except that
the HETG is inserted by using: GratingType="HETG".
  
It can be convenient to use pset to set the parameters in a
local par file, so that method is demonstrated and used here
(pset is provided with ciao or lheasoft.)

Set up the various marx parameters as desired.  Note that
the user can edit and re-do the following lines for new
or modified simulations.  As an example, the first simulation
here is for a POINT source; it is followed by a few additional
lines to do a second simulation with a DISK source.

In the working directory paste these sets of lines to the unix prompt:

- - - Simulation 1 - Point Source
# Make a par file for my marx simulation:
cp marx_orig.par mysim.par

# set the spectrum file to use:
pset mysim.par SpectrumType="FILE"
pset mysim.par SpectrumFile="plawflux_hetg.tbl"
pset mysim.par SourceFlux=-1

# Set other parameters of the simulation:
# Using 50 ks
pset mysim.par ExposureTime=50000
pset mysim.par OutputDir="hetg_plaw"
pset mysim.par DitherModel="INTERNAL"

# Use the HETG with ACIS-S:
pset mysim.par GratingType="HETG"
pset mysim.par DetectorType="ACIS-S"

# Some other parameters it can be useful to set:

# Date of observation (effects ACIS QE)
pset mysim.par TStart=2009.50
# Roll of the observation: 0 puts average dispersion along E -- W.
pset mysim.par Roll_Nom=0.0

# Source RA/DEC (degrees)
pset mysim.par SourceRA=250.000
pset mysim.par SourceDEC=-54.000
# Pointing RA/DEC (degrees)
pset mysim.par RA_Nom=250.000
pset mysim.par Dec_Nom=-54.000


# Finally, 
# Run the simulation:
marx @@mysim.par

   ... marx runs, ending with:
   Writing output to directory 'hetg_plaw' ...
   Total photons: 3495031, Total Photons detected: 255172, (efficiency: 0.073010)
     (efficiency this iteration  0.073986)  Total time: 50000.000492


# Create the fits event file and the aspect solution file:
marx2fits hetg_plaw hetg_plaw_evt2.fits
marxasp MarxDir="hetg_plaw" OutputFile="hetg_plaw_asol1.fits"
#
# all done with MARX simulation 1.
- - -

- - - Simulation 2 - Disk Source

# Do another simulation keeping most things the same as above
# by starting with the "mysim.par" as it is left from above
# but changing a few things to use a DISK source:

# New output dir for the Disk
pset mysim.par OutputDir="hetg_pldisk"

# Change the SourceType:
pset mysim.par SourceType="DISK"
# a thin disk with average radius ~ 2.0"
pset mysim.par S-DiskTheta0=1.7
pset mysim.par S-DiskTheta1=2.3

# Run the second simulation:
marx @@mysim.par

# Create the fits event file and the aspect solution file:
marx2fits hetg_pldisk hetg_pldisk_evt2.fits
marxasp MarxDir="hetg_pldisk" OutputFile="hetg_pldisk_asol1.fits"
#
# all done with MARX simulation 2.
- - -

We can look at the simulation output event files with ds9 to check that they
are as expected before continuuing with ciao-processing, e.g., 

[unix] ds9 hetg_pldisk_evt2.fits

Use a "log" scale to see the dispersed arms and streak as well as the
brighter zeroth-order image.  We can also use ds9 to record the center of the
disk (simulation 2) in X,Y coord.s:  4096.5, 4096.5, for further processing.


........
FYI: the routine marxcat allows simulations to be combined, e.g., we
could do the following to make a combination of the point and disk events:
[unix] marxcat hetg_plaw hetg_pldisk hetg_plboth
and then create fits and asol files as above.
This allows more complex spatial-spectral simulations to be done with marx.

-------