file: casa_sim.txt

working directory:  /nfs/spectra/d1/dd/Science/CasA/marx_070301


Make simulations of a second HETG Cas A observation at
a different roll, etc.


- - - FITS image for input to marx
Use the first observation, acisf01046N002_evt2.fits,
to create a FITS image for marx input.
Filter this image on just the Si (H- and He- like) range:
ENERGY  1750 to 2070  (based on histogram of data)
and select just the zeroth-order image:
X 3700 to 4500 ; Y  3800 to 4420.
Use ds9 to create a FITS image file by saving the 
frame (1024x1024, bin=1):  casa_zo_si.fits


- - - Input spectrum for the image:
Use (for now) the spectrum model for r1 (from '/Spect_all/):
casa_r1_spec.dat


- - - Setup marx
source /nfs/cxc/a1/setup/marx-setup-devel
cp /nfs/cxc/a1/src/marx/marx-dist-devel/marx.par marx_orig.par



- - - Do a 1046-like simulation:
Use Ideal HRMA to avoid adding more blur to the image.
Reduce DitherBlur for same reason.
RA and DEC Nom are from the observation and 
RA DEC Source are the same as Nom to match fits image to real data.


cp marx_orig.par marx.par


pset marx.par TStart=2001.4
pset marx.par GratingType="HETG"

pset marx.par dNumRays=100000
pset marx.par ExposureTime=0.0
pset marx.par HRMA_Ideal="yes"
pset marx.par DitherBlur=0.10

pset marx.par SourceRA=3.5087057410749E+02
pset marx.par SourceDEC=5.8815806203148E+01
pset marx.par RA_Nom=3.5087057410749E+02
pset marx.par Dec_Nom=5.8815806203148E+01

pset marx.par Roll_Nom=8.5955941473876E+01

pset marx.par SpectrumType="FILE"
pset marx.par SourceType="IMAGE"

pset marx.par S-ImageFile="casa_zo_si.fits"
pset marx.par SpectrumFile="casa_r1_spec.dat"
pset marx.par OutputDir="casa_sim1046"

pset marx.par NumRays=5000000

marx

marx2fits casa_sim1046 casa_sim1046.fits


- - -
Compare real and simulated zeroth-order images

Can use evt2img on the simulated and real fits files
to compare them, etc.

evt2img /nfs/spectra/d1/dd/Science/CasA/Data/1046/acisf01046N002_E500toE2300evt2.fits
  Using: Red=0.500016-1.47002, Green=1.47002-1.82002, Blue=1.82002-2.29002

Do similar evts--> image "by hand" using hydra/S-Lang routine
and compare the counts in the 1.75 to 2.00 keV range.

hydra hy3d_isisrc
hydra> ()=evalfile("casa_images");

hydra> casa_images(0,[3]);       % Real data
   Number in Si-peak  : 629.611 x10^3

hydra> casa_images(1,[3]);       % 5Mray simulation
   Number in Si-peak  : 137.613 x10^3

-->  So 5 Mrays is ~  (137.6/629.6) * 70.0 ks = 15.3 ks equiv Si peak.