Research & Science Home ESA Public Web Site Sci-Tech Portal      XMM-Newton Public Web Site XMM-Newton Sci-Tech Portal
Astrophysics Missions Planetary Exploration Missions Solar Terrestrial Science Missions Fundamental Physics Missions Science Faculty

Creation of EPIC ESAS spectral analysis files for a cluster radial profile

Notice
Due to programmed maintenance, XMM-Newton web service, Science Archive and some SOC computing services may be unavailable from Friday, 24 May, 14:30 until Sunday, 26 May, 20:00 (GMT). During this period no ToO requests can be processed.

Introduction

This thread describes how to create EPIC spectral and response files for a complicated multi-region analysis combining data from the three instruments.

This thread uses data from the observation of Abell 1795, ObsID 0097820101, which is also used as the Spectral Analysis example where a script and output files are provided.

Expected Outcome

The final outcome of this thread is the spectral analysis of the radial profile of the cluster of galaxies Abell 1795 with ten annuli.

SAS Tasks to be Used

Prerequisites

Useful Links

  • This thread makes use of the image display software ds9 and the spectral fitting tool Xspec
  • A full guide to the use of the ESAS software in SAS can be found here or here
  • Cluster analysis paper - Snowden et al. 2008, A&A, 478, 615

Caveats



Procedure

This thread contains a step-by-step recipe to create EPIC background and source spectra and response files from all three instruments.
  1. Set up your SAS environment (following the SAS Startup Thread)

  2. Create cleaned (filtered for soft proton excesses) MOS and pn (including OOT processing) event files for your observation:

     epchain
     epchain withoutoftime=true
     emchain
     pn-filter
     mos-filter

    Note that mos-filter indicates which CCDs are operating in an anomalous mode (the one marked by **** in the example below). These CCDs should be excluded in downstream processing.

     1S003
         Limit>1.5 CCD=2 Hardness=5.100 Uncertainty=0.882
         Limit>1.5 CCD=3 Hardness=5.607 Uncertainty=1.150
         Limit>2.5 CCD=4 Hardness=3.110 Uncertainty=0.418
         Limit>2.5 CCD=5 Hardness=1.526 Uncertainty=0.149 ****
         Limit>1.5 CCD=6 Hardness=4.000 Uncertainty=0.716
         Limit>1.5 CCD=7 Hardness=5.629 Uncertainty=1.032
     2S004
         Limit>2.5 CCD=2 Hardness=2.639 Uncertainty=0.365
         Limit>1.5 CCD=3 Hardness=4.419 Uncertainty=0.746
         Limit>1.5 CCD=4 Hardness=6.875 Uncertainty=1.301
         Limit>2.5 CCD=5 Hardness=3.982 Uncertainty=0.590
         Limit>1.5 CCD=6 Hardness=4.432 Uncertainty=0.807
         Limit>1.5 CCD=7 Hardness=4.600 Uncertainty=0.757

    The diagnostic plots created by the filtering tasks should also be examined as they provide an indication of the quality of the data. These have names like mos1S003-hist.qdp, mos2S004-hist.qdp, and pnS005-hist.qdp and cam be plotted using the command, for example, qdp mos1S003-hist.qdp. Figure 1 shows the light-curve screening for the MOS1 instrument.

    Figure 1: MOS1 light curve for the Abell 1795 observation. Notice the residual variation of the nominally good data which indicates the possible existence of residual soft proton contamination.

  3. Run source detection and make point-source masks. Note that in this and many tasks below, the exposure ID must be explicitly noted. This is done by the prefixm, prefixp, and prefix parameters.

     cheese prefixm='1S003 2S004' prefixp=S005 scale=0.5 rate=1.0 dist=40.0 \
       clobber=0 elow=400 ehigh=10000

  4. Examine the MOS soft band images to confirm that the indicated CCDs are operating in an anomalous state, or that any other CCDs are doing so, and verify that the point-source (cheese) masks look reasonable.

    Figure 2: MOS1 (left) and MOS2 (right) images in the soft (0.2-1.0 keV) band. Note the excess counts in the upper left CCD in the MOS1 image most noticeable in the unexposed (to the sky) upper left corner. Displayed by the command: ds9 *soft* &.

    Figure 3: MOS1 (left), MOS2 (middle), and pn (right) cheese masks. Displayed by the command: ds9 *cheese* &

  5. It is quite often useful to create and image of the entire field of view to get a perspective on the observation. In this section a simplified version of the ESAS Imaging thread is provided with minimal comments. For a more extensive discussion see the ESAS Imaging thread and the XMM-ESAS Cookbook.

    Use mos-spectra and pn-spectra to create the required intermediate spectra.

     mos-spectra prefix=1S003 caldb=/PATH/esascaldb region=reg.txt mask=1 \
       elow=400 ehigh=7200 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos-spectra prefix=2S004 caldb=/PATH/esascaldb region=reg.txt mask=1 \
       elow=400 ehigh=7200 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=pn-reg.txt \
       mask=1 elow=400 ehigh=7200 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1

    Use mos_back and pn_back to create the quiescent particle background (QPB) spectra and images (in detector coordinates).

     mos_back prefix=1S003 caldb=/PATH/esascaldb diag=0 elow=400 ehigh=7200 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/PATH/esascaldb diag=0 elow=400 ehigh=7200 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     pn_back prefix=S005 caldb=/PATH/esascaldb diag=0 elow=400 ehigh=7200 \
      quad1=1 quad2=1 quad3=1 quad4=1

    Use rot-im-det-sky to transform the QPB images in detector coordinates into sky coordinates.

     rot-im-det-sky prefix=1S003 mask=0 elow=400 ehigh=7200 mode=1
     rot-im-det-sky prefix=2S004 mask=0 elow=400 ehigh=7200 mode=1
     rot-im-det-sky prefix=S005 mask=0 elow=400 ehigh=7200 mode=1

    comb combines the MOS1, MOS2, and pn images, as well as images from multiple exposures. Since the spectra and images were created with point sources removed, comb must be run using the cheese masking. Note that the soft proton background is not removed for this image.

     comb caldb=/software/XMM/CCF/esas withpartcontrol=1 withsoftcontrol=0 withswcxcontrol=0 \
      nbands=1 elowlist=400 ehighlist=7200 mask=1 ndata=3 prefixlist="1S003 2S004 S005"

    adapt_900 adaptively smooths the image.

     adapt_900 smoothingcounts=50 thresholdmasking=0.02 detector=0 binning=2 elow=400 ehigh=7200 \
      withpartcontrol=yes withsoftcontrol=yes withswcxcontrol=0

    Display the smoothed image.

     ds9 adapt-400-7200.fits &

    Figure 4 shows the final adaptively smoothed, QPB subtracted, and exposure corrected image.

    Figure 4: Abell 1795 adaptively smoothed, background subtracted, and exposure corrected image in the 0.4-7.2 keV band.

    Get rid of some intermediate files before starting on the spectral analysis part of the thread.

     rm -f *-obj-im.fits *obj-im-sp-det.fits *exp-im.fits *mask-im.fits *obj.pi \
      *obj-oot.pi *rmf *arf

  6. Use xmmselect to find the detector coordinates of the center of the cluster for each of the three detectors. Use these coordinates to create the region selection expressions for the annuli. This is a bit tedious. For this thread the annuli are: 0-30", 30-60", 60-120", 120-180", 180-240", 240-300", 300-420", 420-540", 540-720", and 720-840". The region files used below have names starting with reg1, reg2, and reg3 for the MOS1, MOS2, and pn detectors, respectively. Note that the units are 0.05 arc seconds.

     reg1-0-30.txt: &&((DETX,DETY) IN circle(201.0,-219.0,600))
     reg1-30-60.txt: &&((DETX,DETY) IN circle(201.0,-219.0,1200))&&!((DETX,DETY) IN circle(201.0,-219.0,600))
     reg1-60-120.txt: &&((DETX,DETY) IN circle(201.0,-219.0,2400))&&!((DETX,DETY) IN circle(201.0,-219.0,1200))
     reg1-120-180.txt: &&((DETX,DETY) IN circle(201.0,-219.0,3600))&&!((DETX,DETY) IN circle(201.0,-219.0,2400))
     reg1-180-240.txt: &&((DETX,DETY) IN circle(201.0,-219.0,4800))&&!((DETX,DETY) IN circle(201.0,-219.0,3600))
     reg1-240-300.txt: &&((DETX,DETY) IN circle(201.0,-219.0,6000))&&!((DETX,DETY) IN circle(201.0,-219.0,4800))
     reg1-300-420.txt: &&((DETX,DETY) IN circle(201.0,-219.0,8400))&&!((DETX,DETY) IN circle(201.0,-219.0,6000))
     reg1-420-540.txt: &&((DETX,DETY) IN circle(201.0,-219.0,10800))&&!((DETX,DETY) IN circle(201.0,-219.0,8400))
     reg1-540-720.txt: &&((DETX,DETY) IN circle(201.0,-219.0,14400))&&!((DETX,DETY) IN circle(201.0,-219.0,10800))
     reg1-720-840.txt: &&((DETX,DETY) IN circle(201.0,-219.0,16800))&&!((DETX,DETY) IN circle(201.0,-219.0,14400))
     reg2-0-30.txt: &&((DETX,DETY) IN circle(39.5,-26.0,600))
     reg2-120-180.txt: &&((DETX,DETY) IN circle(39.5,-26.0,3600))&&!((DETX,DETY) IN circle(39.5,-26.0,2400))
     reg2-180-240.txt: &&((DETX,DETY) IN circle(39.5,-26.0,4800))&&!((DETX,DETY) IN circle(39.5,-26.0,3600))
     reg2-240-300.txt: &&((DETX,DETY) IN circle(39.5,-26.0,6000))&&!((DETX,DETY) IN circle(39.5,-26.0,4800))
     reg2-300-420.txt: &&((DETX,DETY) IN circle(39.5,-26.0,8400))&&!((DETX,DETY) IN circle(39.5,-26.0,6000))
     reg2-30-60.txt: &&((DETX,DETY) IN circle(39.5,-26.0,1200))&&!((DETX,DETY) IN circle(39.5,-26.0,600))
     reg2-420-540.txt: &&((DETX,DETY) IN circle(39.5,-26.0,10800))&&!((DETX,DETY) IN circle(39.5,-26.0,8400))
     reg2-540-720.txt: &&((DETX,DETY) IN circle(39.5,-26.0,14400))&&!((DETX,DETY) IN circle(39.5,-26.0,10800))
     reg2-60-120.txt: &&((DETX,DETY) IN circle(39.5,-26.0,2400))&&!((DETX,DETY) IN circle(39.5,-26.0,1200))
     reg2-720-840.txt: &&((DETX,DETY) IN circle(39.5,-26.0,16800))&&!((DETX,DETY) IN circle(39.5,-26.0,14400))
     reg3-0-30.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,600))
     reg3-120-180.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,3600))&&!((DETX,DETY) IN circle( 59.0, -147.5,2400))
     reg3-180-240.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,4800))&&!((DETX,DETY) IN circle( 59.0, -147.5,3600))
     reg3-240-300.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,6000))&&!((DETX,DETY) IN circle( 59.0, -147.5,4800))
     reg3-300-420.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,8400))&&!((DETX,DETY) IN circle( 59.0, -147.5,6000))
     reg3-30-60.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,1200))&&!((DETX,DETY) IN circle( 59.0, -147.5,600))
     reg3-420-540.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,10800))&&!((DETX,DETY) IN circle( 59.0, -147.5,8400))
     reg3-540-720.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,14400))&&!((DETX,DETY) IN circle( 59.0, -147.5,10800))
     reg3-60-120.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,2400))&&!((DETX,DETY) IN circle( 59.0, -147.5,1200))
     reg3-720-840.txt: &&((DETX,DETY) IN circle( 59.0, -147.5,16800))&&!((DETX,DETY) IN circle( 59.0, -147.5,14400))

  7. Use mos-spectra and pn-spectra to create the intermediate spectra that are used by mos_back and pn_back to create the QPB spectra for the regions of interest and their associated RMF and ARF files. MOS CCDs affected by anomalous states should be deselected (this is done by the ccd# parameters). MOS1 CCD6 should be deselected if the observation took place after the meteorite damage (See Section 3.3 of the XMM-Newton Users Handbook). Currently the tasks require the path to the additional CalDB files required for XMM-ESAS (in this case, /PATH/esascaldb). As images are not of interest, the energy selection parameters are set to 0. In order to avoid overwriting important files, the processing of an individual region should be done together and the files renamed.

    Skip to the end of the this section.

     mos-spectra prefix=1S003 caldb=/PATH/esascaldb region=reg1-0-30.txt mask=1 \
       elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/PATH/esascaldb diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-0-30.pi
     mv mos1S003-back.pi mos1S003-back-0-30.pi
     mv mos1S003.rmf mos1S003-0-30.rmf
     mv mos1S003.arf mos1S003-0-30.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-0-30.fits
     specgroup spectrumset=mos1S003-obj-0-30.pi mincounts=100 rmfset=mos1S003-0-30.rmf \
       arfset=mos1S003-0-30.arf backgndset=mos1S003-back-0-30.pi \
       groupedset=mos1S003-obj-0-30-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-0-30.txt\
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-0-30.pi
     mv mos2S004-back.pi mos2S004-back-0-30.pi
     mv mos2S004.rmf mos2S004-0-30.rmf
     mv mos2S004.arf mos2S004-0-30.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-0-30.fits
     specgroup spectrumset=mos2S004-obj-0-30.pi mincounts=100 rmfset=mos2S004-0-30.rmf \
       arfset=mos2S004-0-30.arf backgndset=mos2S004-back-0-30.pi \
       groupedset=mos2S004-obj-0-30-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-0-30.txt mask=1 \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-0-30.pi
     mv pnS005-obj.pi pnS005-obj-0-30.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-0-30.pi
     mv pnS005-back.pi pnS005-back-0-30.pi
     mv pnS005.rmf pnS005-0-30.rmf
     mv pnS005.arf pnS005-0-30.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-0-30.fits
     specgroup spectrumset=pnS005-obj-0-30.pi mincounts=100 rmfset=pnS005-0-30.rmf \
       arfset=pnS005-0-30.arf backgndset=pnS005-back-0-30.pi \
       groupedset=pnS005-obj-0-30-grp.pi

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-30-60.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-30-60.pi
     mv mos1S003-back.pi mos1S003-back-30-60.pi
     mv mos1S003.rmf mos1S003-30-60.rmf
     mv mos1S003.arf mos1S003-30-60.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-30-60.fits
     specgroup spectrumset=mos1S003-obj-30-60.pi mincounts=100 rmfset=mos1S003-30-60.rmf \
       arfset=mos1S003-30-60.arf backgndset=mos1S003-back-30-60.pi \
       groupedset=mos1S003-obj-30-60-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-30-60.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-30-60.pi
     mv mos2S004-back.pi mos2S004-back-30-60.pi
     mv mos2S004.rmf mos2S004-30-60.rmf
     mv mos2S004.arf mos2S004-30-60.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-30-60.fits
     specgroup spectrumset=mos2S004-obj-30-60.pi mincounts=100 rmfset=mos2S004-30-60.rmf \
       arfset=mos2S004-30-60.arf backgndset=mos2S004-back-30-60.pi \
       groupedset=mos2S004-obj-30-60-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-30-60.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-30-60.pi
     mv pnS005-obj.pi pnS005-obj-30-60.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-30-60.pi
     mv pnS005-back.pi pnS005-back-30-60.pi
     mv pnS005.rmf pnS005-30-60.rmf
     mv pnS005.arf pnS005-30-60.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-30-60.fits
     specgroup spectrumset=pnS005-obj-30-60.pi mincounts=100 rmfset=pnS005-30-60.rmf \
       arfset=pnS005-30-60.arf backgndset=pnS005-back-30-60.pi \
       groupedset=pnS005-obj-30-60-grp.pi

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-60-120.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-60-120.pi
     mv mos1S003-back.pi mos1S003-back-60-120.pi
     mv mos1S003.rmf mos1S003-60-120.rmf
     mv mos1S003.arf mos1S003-60-120.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-60-120.fits
     specgroup spectrumset=mos1S003-obj-60-120.pi mincounts=100 rmfset=mos1S003-60-120.rmf \
       arfset=mos1S003-60-120.arf backgndset=mos1S003-back-60-120.pi \
       groupedset=mos1S003-obj-60-120-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-60-120.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 \
       ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-60-120.pi
     mv mos2S004-back.pi mos2S004-back-60-120.pi
     mv mos2S004.rmf mos2S004-60-120.rmf
     mv mos2S004.arf mos2S004-60-120.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-60-120.fits
     specgroup spectrumset=mos2S004-obj-60-120.pi mincounts=100 rmfset=mos2S004-60-120.rmf \
       arfset=mos2S004-60-120.arf backgndset=mos2S004-back-60-120.pi \
       groupedset=mos2S004-obj-60-120-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-60-120.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-60-120.pi
     mv pnS005-obj.pi pnS005-obj-60-120.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-60-120.pi
     mv pnS005-back.pi pnS005-back-60-120.pi
     mv pnS005.rmf pnS005-60-120.rmf
     mv pnS005.arf pnS005-60-120.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-60-120.fits
     specgroup spectrumset=pnS005-obj-60-120.pi mincounts=100 rmfset=pnS005-60-120.rmf \
       arfset=pnS005-60-120.arf backgndset=pnS005-back-60-120.pi \
       groupedset=pnS005-obj-60-120-grp.pi

    Skip to the end of the this section.

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-120-180.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-120-180.pi
     mv mos1S003-back.pi mos1S003-back-120-180.pi
     mv mos1S003.rmf mos1S003-120-180.rmf
     mv mos1S003.arf mos1S003-120-180.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-120-180.fits
     specgroup spectrumset=mos1S003-obj-120-180.pi mincounts=100 rmfset=mos1S003-120-180.rmf \
       arfset=mos1S003-120-180.arf backgndset=mos1S003-back-120-180.pi \
       groupedset=mos1S003-obj-120-180-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-120-180.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-120-180.pi
     mv mos2S004-back.pi mos2S004-back-120-180.pi
     mv mos2S004.rmf mos2S004-120-180.rmf
     mv mos2S004.arf mos2S004-120-180.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-120-180.fits
     specgroup spectrumset=mos2S004-obj-120-180.pi mincounts=100 rmfset=mos2S004-120-180.rmf \
       arfset=mos2S004-120-180.arf backgndset=mos2S004-back-120-180.pi \
       groupedset=mos2S004-obj-120-180-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-120-180.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-120-180.pi
     mv pnS005-obj.pi pnS005-obj-120-180.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-120-180.pi
     mv pnS005-back.pi pnS005-back-120-180.pi
     mv pnS005.rmf pnS005-120-180.rmf
     mv pnS005.arf pnS005-120-180.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-120-180.fits
     specgroup spectrumset=pnS005-obj-120-180.pi mincounts=100 rmfset=pnS005-120-180.rmf \
       arfset=pnS005-120-180.arf backgndset=pnS005-back-120-180.pi \
       groupedset=pnS005-obj-120-180-grp.pi

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-180-240.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-180-240.pi
     mv mos1S003-back.pi mos1S003-back-180-240.pi
     mv mos1S003.rmf mos1S003-180-240.rmf
     mv mos1S003.arf mos1S003-180-240.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-180-240.fits
     specgroup spectrumset=mos1S003-obj-180-240.pi mincounts=100 rmfset=mos1S003-180-240.rmf \
       arfset=mos1S003-180-240.arf backgndset=mos1S003-back-180-240.pi \
       groupedset=mos1S003-obj-180-240-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-180-240.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-180-240.pi
     mv mos2S004-back.pi mos2S004-back-180-240.pi
     mv mos2S004.rmf mos2S004-180-240.rmf
     mv mos2S004.arf mos2S004-180-240.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-180-240.fits
     specgroup spectrumset=mos2S004-obj-180-240.pi mincounts=100 rmfset=mos2S004-180-240.rmf \
       arfset=mos2S004-180-240.arf backgndset=mos2S004-back-180-240.pi \
       groupedset=mos2S004-obj-180-240-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-180-240.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-180-240.pi
     mv pnS005-obj.pi pnS005-obj-180-240.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-180-240.pi
     mv pnS005-back.pi pnS005-back-180-240.pi
     mv pnS005.rmf pnS005-180-240.rmf
     mv pnS005.arf pnS005-180-240.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-180-240.fits
     specgroup spectrumset=pnS005-obj-180-240.pi mincounts=100 rmfset=pnS005-180-240.rmf \
       arfset=pnS005-180-240.arf backgndset=pnS005-back-180-240.pi \
       groupedset=pnS005-obj-180-240-grp.pi

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-240-300.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-240-300.pi
     mv mos1S003-back.pi mos1S003-back-240-300.pi
     mv mos1S003.rmf mos1S003-240-300.rmf
     mv mos1S003.arf mos1S003-240-300.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-240-300.fits
     specgroup spectrumset=mos1S003-obj-240-300.pi mincounts=100 rmfset=mos1S003-240-300.rmf \
       arfset=mos1S003-240-300.arf backgndset=mos1S003-back-240-300.pi \
       groupedset=mos1S003-obj-240-300-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-240-300.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-240-300.pi
     mv mos2S004-back.pi mos2S004-back-240-300.pi
     mv mos2S004.rmf mos2S004-240-300.rmf
     mv mos2S004.arf mos2S004-240-300.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-240-300.fits
     specgroup spectrumset=mos2S004-obj-240-300.pi mincounts=100 rmfset=mos2S004-240-300.rmf \
       arfset=mos2S004-240-300.arf backgndset=mos2S004-back-240-300.pi \
       groupedset=mos2S004-obj-240-300-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-240-300.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-240-300.pi
     mv pnS005-obj.pi pnS005-obj-240-300.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-240-300.pi
     mv pnS005-back.pi pnS005-back-240-300.pi
     mv pnS005.rmf pnS005-240-300.rmf
     mv pnS005.arf pnS005-240-300.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-240-300.fits
     specgroup spectrumset=pnS005-obj-240-400.pi mincounts=100 rmfset=pnS005-240-300.rmf \
       arfset=pnS005-240-300.arf backgndset=pnS005-back-240-300.pi \
       groupedset=pnS005-obj-240-300-grp.pi

    Skip to the end of the this section.

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-300-420.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-300-420.pi
     mv mos1S003-back.pi mos1S003-back-300-420.pi
     mv mos1S003.rmf mos1S003-300-420.rmf
     mv mos1S003.arf mos1S003-300-420.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-300-420.fits
     specgroup spectrumset=mos1S003-obj-300-420.pi mincounts=100 rmfset=mos1S003-300-420.rmf \
       arfset=mos1S003-300-420.arf backgndset=mos1S003-back-300-420.pi \
       groupedset=mos1S003-obj-300-420-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-300-420.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-300-420.pi
     mv mos2S004-back.pi mos2S004-back-300-420.pi
     mv mos2S004.rmf mos2S004-300-420.rmf
     mv mos2S004.arf mos2S004-300-420.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-300-420.fits
     specgroup spectrumset=mos2S004-obj-300-420.pi mincounts=100 rmfset=mos2S004-300-420.rmf \
       arfset=mos2S004-300-420.arf backgndset=mos2S004-back-300-420.pi \
       groupedset=mos2S004-obj-300-420-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-300-420.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-300-420.pi
     mv pnS005-obj.pi pnS005-obj-300-420.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-300-420.pi
     mv pnS005-back.pi pnS005-back-300-420.pi
     mv pnS005.rmf pnS005-300-420.rmf
     mv pnS005.arf pnS005-300-420.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-300-420.fits
     specgroup spectrumset=pnS005-obj-300-420.pi mincounts=100 rmfset=pnS005-300-420.rmf \
       arfset=pnS005-300-420.arf backgndset=pnS005-back-300-420.pi \
       groupedset=pnS005-obj-300-420-grp.pi

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-420-540.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-420-540.pi
     mv mos1S003-back.pi mos1S003-back-420-540.pi
     mv mos1S003.rmf mos1S003-420-540.rmf
     mv mos1S003.arf mos1S003-420-540.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-420-540.fits
     specgroup spectrumset=mos1S003-obj-420-540.pi mincounts=100 rmfset=mos1S003-420-540.rmf \
       arfset=mos1S003-420-540.arf backgndset=mos1S003-back-420-540.pi \
       groupedset=mos1S003-obj-420-540-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-420-540.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-420-540.pi
     mv mos2S004-back.pi mos2S004-back-420-540.pi
     mv mos2S004.rmf mos2S004-420-540.rmf
     mv mos2S004.arf mos2S004-420-540.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-420-540.fits
     specgroup spectrumset=mos2S004-obj-420-540.pi mincounts=100 rmfset=mos2S004-420-540.rmf \
       arfset=mos2S004-420-540.arf backgndset=mos2S004-back-420-540.pi \
       groupedset=mos2S004-obj-420-540-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-420-540.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-420-540.pi
     mv pnS005-obj.pi pnS005-obj-420-540.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-420-540.pi
     mv pnS005-back.pi pnS005-back-420-540.pi
     mv pnS005.rmf pnS005-420-540.rmf
     mv pnS005.arf pnS005-420-540.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-420-540.fits
     specgroup spectrumset=pnS005-obj-420-540.pi mincounts=100 rmfset=pnS005-420-540.rmf \
       arfset=pnS005-420-540.arf backgndset=pnS005-back-420-540.pi \
       groupedset=pnS005-obj-420-540-grp.pi

     mos-spectra prefix=1S003 caldb=/software/XMM/CCF/esas region=reg1-540-720.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-540-720.pi
     mv mos1S003-back.pi mos1S003-back-540-720.pi
     mv mos1S003.rmf mos1S003-540-720.rmf
     mv mos1S003.arf mos1S003-540-720.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-540-720.fits
     specgroup spectrumset=mos1S003-obj-540-720.pi mincounts=100 rmfset=mos1S003-540-720.rmf \
       arfset=mos1S003-540-720.arf backgndset=mos1S003-back-540-720.pi \
       groupedset=mos1S003-obj-540-720-grp.pi

     mos-spectra prefix=2S004 caldb=/software/XMM/CCF/esas region=reg2-540-720.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-540-720.pi
     mv mos2S004-back.pi mos2S004-back-540-720.pi
     mv mos2S004.rmf mos2S004-540-720.rmf
     mv mos2S004.arf mos2S004-540-720.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-540-720.fits
     specgroup spectrumset=mos2S004-obj-540-720.pi mincounts=100 rmfset=mos2S004-540-720.rmf \
       arfset=mos2S004-540-720.arf backgndset=mos2S004-back-540-720.pi \
       groupedset=mos2S004-obj-540-720-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-540-720.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-540-720.pi
     mv pnS005-obj.pi pnS005-obj-540-720.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-540-720.pi
     mv pnS005-back.pi pnS005-back-540-720.pi
     mv pnS005.rmf pnS005-540-720.rmf
     mv pnS005.arf pnS005-540-720.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-540-720.fits
     specgroup spectrumset=pnS005-obj-540-270.pi mincounts=100 rmfset=pnS005-540-720.rmf \
       arfset=pnS005-540-720.arf backgndset=pnS005-back-540-720.pi \
       groupedset=pnS005-obj-540-720-grp.pi

     mos-spectra prefix="1S003" caldb=/software/XMM/CCF/esas region=reg1-720-840.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mos_back prefix=1S003 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=0 ccd6=1 ccd7=1
     mv mos1S003-obj.pi mos1S003-obj-720-840.pi
     mv mos1S003-back.pi mos1S003-back-720-840.pi
     mv mos1S003.rmf mos1S003-720-840.rmf
     mv mos1S003.arf mos1S003-720-840.arf
     mv mos1S003-obj-im-sp-det.fits mos1S003-sp-720-840.fits
     specgroup spectrumset=mos1S003-obj-720-840.pi mincounts=100 rmfset=mos1S003-720-840.rmf \
       arfset=mos1S003-720-840.arf backgndset=mos1S003-back-720-840.pi \
       groupedset=mos1S003-obj-720-840-grp.pi

     mos-spectra prefix="2S004" caldb=/software/XMM/CCF/esas region=reg2-720-840.txt \
       mask=1 elow=0 ehigh=0 ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mos_back prefix=2S004 caldb=/software/XMM/CCF/esas diag=0 elow=0 ehigh=0 \
       ccd1=1 ccd2=1 ccd3=1 ccd4=1 ccd5=1 ccd6=1 ccd7=1
     mv mos2S004-obj.pi mos2S004-obj-720-840.pi
     mv mos2S004-back.pi mos2S004-back-720-840.pi
     mv mos2S004.rmf mos2S004-720-840.rmf
     mv mos2S004.arf mos2S004-720-840.arf
     mv mos2S004-obj-im-sp-det.fits mos2S004-sp-720-840.fits
     specgroup spectrumset=mos2S004-obj-720-840.pi mincounts=100 rmfset=mos2S004-720-840.rmf \
       arfset=mos2S004-720-840.arf backgndset=mos2S004-back-720-840.pi \
       groupedset=mos2S004-obj-720-840-grp.pi

     pn-spectra prefix=S005 caldb=/software/XMM/CCF/esas region=reg3-720-840.txt \
       mask=1 elow=0 ehigh=0 pattern=4 quad1=1 quad2=1 quad3=1 quad4=1
     pn_back prefix=S005 caldb=/XMM/sas/CCF/esas diag=0 elow=0 ehigh=0 \
       quad1=1 quad2=1 quad3=1 quad4=1
     mv pnS005-obj-os.pi pnS005-obj-os-720-840.pi
     mv pnS005-obj.pi pnS005-obj-720-840.pi
     mv pnS005-obj-oot.pi pnS005-obj-oot-720-840.pi
     mv pnS005-back.pi pnS005-back-720-840.pi
     mv pnS005.rmf pnS005-720-840.rmf
     mv pnS005.arf pnS005-720-840.arf
     mv pnS005-obj-im-sp-det.fits pnS005-sp-720-840.fits
     specgroup spectrumset=pnS005-obj-720-840.pi mincounts=100 rmfset=pnS005-720-840.rmf \
       arfset=pnS005-720-840.arf backgndset=pnS005-back-720-840.pi \
       groupedset=pnS005-obj-720-840-grp.pi

  8. Because of XMM's finite PSF, if there are strong spectral gradients the spectrum in one region can be significantly affected by the spectrum in a neighboring region. The creation and use of cross arfs provides a means of addressing this cross-talk. Use arfgen and evselect to create cross arfs for the spectral analysis. evselect is first run to create an instrument map. Some care must be taken in choosing the pixel size of the image as arfgen computation time scale as the number of pixels but if the pixels are too coarse the ARF will be underestimated. Generally, the ARF will reach its asymptotic limit about the time it reached infinite computation time. After the detector map has been created, arfgen can be run for the appropriate ARFs - the smaller the region is, the finer the detector map should be. Also, not that arfgen in the cross-arf mode can be run for the appropriate ARF for a given region on its own.

    Skip to the end of the this section.

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=50 ycolumn=DETY yimagebinsize=50 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-0-30-0-30.arf spectrumset=mos1S003-obj-0-30.pi \
       crossreg_spectrumset=mos1S003-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-0-30-30-60.arf spectrumset=mos1S003-obj-30-60.pi \
       crossreg_spectrumset=mos1S003-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-0-30-60-120.arf spectrumset=mos1S003-obj-60-120.pi \
       crossreg_spectrumset=mos1S003-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-30-60-30-60.arf spectrumset=mos1S003-obj-30-60.pi \
       crossreg_spectrumset=mos1S003-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-30-60-0-30.arf spectrumset=mos1S003-obj-0-30.pi \
       crossreg_spectrumset=mos1S003-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-30-60-60-120.arf spectrumset=mos1S003-obj-60-120.pi \
       crossreg_spectrumset=mos1S003-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=50 ycolumn=DETY yimagebinsize=50 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-0-30-0-30.arf spectrumset=mos2S004-obj-0-30.pi \
       crossreg_spectrumset=mos2S004-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-0-30-30-60.arf spectrumset=mos2S004-obj-30-60.pi \
       crossreg_spectrumset=mos2S004-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-0-30-60-120.arf spectrumset=mos2S004-obj-60-120.pi \
       crossreg_spectrumset=mos2S004-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-30-60-30-60.arf spectrumset=mos2S004-obj-30-60.pi \
       crossreg_spectrumset=mos2S004-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-30-60-0-30.arf spectrumset=mos2S004-obj-0-30.pi \
       crossreg_spectrumset=mos2S004-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-30-60-60-120.arf spectrumset=mos2S004-obj-60-120.pi \
       crossreg_spectrumset=mos2S004-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=50 ycolumn=DETY yimagebinsize=50 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-0-30-0-30.arf spectrumset=pnS005-obj-0-30.pi \
       crossreg_spectrumset=pnS005-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-0-30-30-60.arf spectrumset=pnS005-obj-30-60.pi \
       crossreg_spectrumset=pnS005-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-0-30-60-120.arf spectrumset=pnS005-obj-60-120.pi \
       crossreg_spectrumset=pnS005-obj-0-30.pi crossregionarf=yes withrmfset=yes r\
       mfset=pnS005-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-30-60-30-60.arf spectrumset=pnS005-obj-30-60.pi \
       crossreg_spectrumset=pnS005-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-30-60-0-30.arf spectrumset=pnS005-obj-0-30.pi \
       crossreg_spectrumset=pnS005-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-30-60-60-120.arf spectrumset=pnS005-obj-60-120.pi \
       crossreg_spectrumset=pnS005-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

    Skip to the end of the this section.

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       cxcolumn=DETX ximagebinsize=70 ycolumn=DETY yimagebinsize=70 squarepixels=yes \
       cimagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-0-30-120-180.arf spectrumset=mos1S003-obj-120-180.pi \
       ccrossreg_spectrumset=mos1S003-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       crmfset=mos1S003-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       cdetmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       cmodelootcorr=no
     arfgen arfset=mos1S003-30-60-120-180.arf spectrumset=mos1S003-obj-120-180.pi \
       ccrossreg_spectrumset=mos1S003-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       crmfset=mos1S003-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       cdetmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       cmodelootcorr=no
     arfgen arfset=mos1S003-60-120-60-120.arf spectrumset=mos1S003-obj-60-120.pi \
       ccrossreg_spectrumset=mos1S003-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       crmfset=mos1S003-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       cdetmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       cmodelootcorr=no
     arfgen arfset=mos1S003-60-120-0-30.arf spectrumset=mos1S003-obj-0-30.pi \
       ccrossreg_spectrumset=mos1S003-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       crmfset=mos1S003-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       cdetmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       cmodelootcorr=no
     arfgen arfset=mos1S003-60-120-30-60.arf spectrumset=mos1S003-obj-30-60.pi \
       crossreg_spectrumset=mos1S003-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-60-120-120-180.arf spectrumset=mos1S003-obj-120-180.pi \
       crossreg_spectrumset=mos1S003-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=70 ycolumn=DETY yimagebinsize=70 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-0-30-120-180.arf spectrumset=mos2S004-obj-120-180.pi \
       crossreg_spectrumset=mos2S004-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-30-60-120-180.arf spectrumset=mos2S004-obj-120-180.pi \
       crossreg_spectrumset=mos2S004-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-60-120-60-120.arf spectrumset=mos2S004-obj-60-120.pi \
       crossreg_spectrumset=mos2S004-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-60-120-0-30.arf spectrumset=mos2S004-obj-0-30.pi \
       crossreg_spectrumset=mos2S004-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-60-120-30-60.arf spectrumset=mos2S004-obj-30-60.pi \
       crossreg_spectrumset=mos2S004-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-60-120-120-180.arf spectrumset=mos2S004-obj-120-180.pi \
       crossreg_spectrumset=mos2S004-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       mxcolumn=DETX ximagebinsize=70 ycolumn=DETY yimagebinsize=70 squarepixels=yes \
       mimagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-0-30-120-180.arf spectrumset=pnS005-obj-120-180.pi \
       mcrossreg_spectrumset=pnS005-obj-0-30.pi crossregionarf=yes withrmfset=yes \
       mrmfset=pnS005-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       mdetmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       mmodelootcorr=no
     arfgen arfset=pnS005-30-60-120-180.arf spectrumset=pnS005-obj-120-180.pi \
       mcrossreg_spectrumset=pnS005-obj-30-60.pi crossregionarf=yes withrmfset=yes \
       mrmfset=pnS005-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       mdetmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       mmodelootcorr=no
     arfgen arfset=pnS005-60-120-60-120.arf spectrumset=pnS005-obj-60-120.pi \
       mcrossreg_spectrumset=pnS005-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       mrmfset=pnS005-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       mdetmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       mmodelootcorr=no
     arfgen arfset=pnS005-60-120-0-30.arf spectrumset=pnS005-obj-0-30.pi \
       mcrossreg_spectrumset=pnS005-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       mrmfset=pnS005-0-30.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       mdetmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       mmodelootcorr=no
     arfgen arfset=pnS005-60-120-30-60.arf spectrumset=pnS005-obj-30-60.pi \
       mcrossreg_spectrumset=pnS005-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       mrmfset=pnS005-30-60.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       mdetmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       mmodelootcorr=no
     arfgen arfset=pnS005-60-120-120-180.arf spectrumset=pnS005-obj-120-180.pi \
       mcrossreg_spectrumset=pnS005-obj-60-120.pi crossregionarf=yes withrmfset=yes \
       mrmfset=pnS005-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       mdetmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       mmodelootcorr=no

    Skip to the end of the this section.

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=100 ycolumn=DETY yimagebinsize=100 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-120-180-120-180.arf spectrumset=mos1S003-obj-120-180.pi \
       crossreg_spectrumset=mos1S003-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
    arfgen arfset=mos1S003-120-180-60-120.arf spectrumset=mos1S003-obj-60-120.pi \
       crossreg_spectrumset=mos1S003-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-120-180-180-240.arf spectrumset=mos1S003-obj-180-240.pi \
       crossreg_spectrumset=mos1S003-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=100 ycolumn=DETY yimagebinsize=100 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-120-180-120-180.arf spectrumset=mos2S004-obj-120-180.pi \
       crossreg_spectrumset=mos2S004-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-120-180-60-120.arf spectrumset=mos2S004-obj-60-120.pi \
       crossreg_spectrumset=mos2S004-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-120-180-180-240.arf spectrumset=mos2S004-obj-180-240.pi \
       crossreg_spectrumset=mos2S004-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=100 ycolumn=DETY yimagebinsize=100 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-120-180-120-180.arf spectrumset=pnS005-obj-120-180.pi \
       crossreg_spectrumset=pnS005-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-120-180-60-120.arf spectrumset=pnS005-obj-60-120.pi \
       crossreg_spectrumset=pnS005-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-60-120.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-120-180-180-240.arf spectrumset=pnS005-obj-180-240.pi \
       crossreg_spectrumset=pnS005-obj-120-180.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=150 ycolumn=DETY yimagebinsize=150 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-180-240-180-240.arf spectrumset=mos1S003-obj-180-240.pi \
       crossreg_spectrumset=mos1S003-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-180-240-120-180.arf spectrumset=mos1S003-obj-120-180.pi \
       crossreg_spectrumset=mos1S003-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-180-240-240-300.arf spectrumset=mos1S003-obj-240-300.pi \
       crossreg_spectrumset=mos1S003-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=150 ycolumn=DETY yimagebinsize=150 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-180-240-180-240.arf spectrumset=mos2S004-obj-180-240.pi \
       crossreg_spectrumset=mos2S004-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-180-240-120-180.arf spectrumset=mos2S004-obj-120-180.pi \
       crossreg_spectrumset=mos2S004-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-180-240-240-300.arf spectrumset=mos2S004-obj-240-300.pi \
       crossreg_spectrumset=mos2S004-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=150 ycolumn=DETY yimagebinsize=150 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-180-240-180-240.arf spectrumset=pnS005-obj-180-240.pi \
       crossreg_spectrumset=pnS005-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-180-240-120-180.arf spectrumset=pnS005-obj-120-180.pi \
       crossreg_spectrumset=pnS005-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-120-180.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-180-240-240-300.arf spectrumset=pnS005-obj-240-300.pi \
       crossreg_spectrumset=pnS005-obj-180-240.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

    Skip to the end of the this section.

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=200 ycolumn=DETY yimagebinsize=200 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-240-300-240-300.arf spectrumset=mos1S003-obj-240-300.pi \
       crossreg_spectrumset=mos1S003-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-240-300-180-240.arf spectrumset=mos1S003-obj-180-240.pi \
       crossreg_spectrumset=mos1S003-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-240-300-300-420.arf spectrumset=mos1S003-obj-300-420.pi \
       crossreg_spectrumset=mos1S003-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=200 ycolumn=DETY yimagebinsize=200 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-240-300-240-300.arf spectrumset=mos2S004-obj-240-300.pi \
       crossreg_spectrumset=mos2S004-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
    arfgen arfset=mos2S004-240-300-180-240.arf spectrumset=mos2S004-obj-180-240.pi \
       crossreg_spectrumset=mos2S004-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-240-300-300-420.arf spectrumset=mos2S004-obj-300-420.pi \
       crossreg_spectrumset=mos2S004-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=200 ycolumn=DETY yimagebinsize=200 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-240-300-240-300.arf spectrumset=pnS005-obj-240-300.pi \
       crossreg_spectrumset=pnS005-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-240-300-180-240.arf spectrumset=pnS005-obj-180-240.pi \
       crossreg_spectrumset=pnS005-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-180-240.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-240-300-300-420.arf spectrumset=pnS005-obj-300-420.pi \
       crossreg_spectrumset=pnS005-obj-240-300.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=250 ycolumn=DETY yimagebinsize=250 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-300-420-300-420.arf spectrumset=mos1S003-obj-300-420.pi \
       crossreg_spectrumset=mos1S003-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-300-420-240-300.arf spectrumset=mos1S003-obj-240-300.pi \
       crossreg_spectrumset=mos1S003-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-300-420-420-540.arf spectrumset=mos1S003-obj-420-540.pi \
       crossreg_spectrumset=mos1S003-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=250 ycolumn=DETY yimagebinsize=250 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-300-420-300-420.arf spectrumset=mos2S004-obj-300-420.pi \
       crossreg_spectrumset=mos2S004-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-300-420-240-300.arf spectrumset=mos2S004-obj-240-300.pi \
       crossreg_spectrumset=mos2S004-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-300-420-420-540.arf spectrumset=mos2S004-obj-420-540.pi \
       crossreg_spectrumset=mos2S004-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=250 ycolumn=DETY yimagebinsize=250 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-300-420-300-420.arf spectrumset=pnS005-obj-300-420.pi \
       crossreg_spectrumset=pnS005-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-300-420-240-300.arf spectrumset=pnS005-obj-240-300.pi \
       crossreg_spectrumset=pnS005-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-240-300.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-300-420-420-540.arf spectrumset=pnS005-obj-420-540.pi \
       crossreg_spectrumset=pnS005-obj-300-420.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

    Skip to the end of the this section.

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=300 ycolumn=DETY yimagebinsize=300 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-420-540-420-540.arf spectrumset=mos1S003-obj-420-540.pi \
       crossreg_spectrumset=mos1S003-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits\
       modelootcorr=no
     arfgen arfset=mos1S003-420-540-300-420.arf spectrumset=mos1S003-obj-300-420.pi \
       crossreg_spectrumset=mos1S003-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-420-540-540-720.arf spectrumset=mos1S003-obj-540-720.pi \
       crossreg_spectrumset=mos1S003-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=300 ycolumn=DETY yimagebinsize=300 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-420-540-420-540.arf spectrumset=mos2S004-obj-420-540.pi \
       crossreg_spectrumset=mos2S004-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-420-540-300-420.arf spectrumset=mos2S004-obj-300-420.pi \
       crossreg_spectrumset=mos2S004-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-420-540-540-720.arf spectrumset=mos2S004-obj-540-720.pi \
       crossreg_spectrumset=mos2S004-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=300 ycolumn=DETY yimagebinsize=300 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-420-540-420-540.arf spectrumset=pnS005-obj-420-540.pi \
       crossreg_spectrumset=pnS005-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-420-540-300-420.arf spectrumset=pnS005-obj-300-420.pi \
       crossreg_spectrumset=pnS005-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-300-420.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-420-540-540-720.arf spectrumset=pnS005-obj-540-720.pi \
       crossreg_spectrumset=pnS005-obj-420-540.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=350 ycolumn=DETY yimagebinsize=350 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-540-720-540-720.arf spectrumset=mos1S003-obj-540-720.pi \
       crossreg_spectrumset=mos1S003-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-540-720-420-540.arf spectrumset=mos1S003-obj-420-540.pi \
       crossreg_spectrumset=mos1S003-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-540-720-720-840.arf spectrumset=mos1S003-obj-720-840.pi \
       crossreg_spectrumset=mos1S003-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-720-840.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=350 ycolumn=DETY yimagebinsize=350 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-540-720-540-720.arf spectrumset=mos2S004-obj-540-720.pi \
       crossreg_spectrumset=mos2S004-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-540-720-420-540.arf spectrumset=mos2S004-obj-420-540.pi \
       crossreg_spectrumset=mos2S004-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-540-720-720-840.arf spectrumset=mos2S004-obj-720-840.pi \
       crossreg_spectrumset=mos2S004-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-720-840.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=350 ycolumn=DETY yimagebinsize=350 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-540-720-540-720.arf spectrumset=pnS005-obj-540-720.pi \
       crossreg_spectrumset=pnS005-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-540-720-420-540.arf spectrumset=pnS005-obj-420-540.pi \
       crossreg_spectrumset=pnS005-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-420-540.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-540-720-720-840.arf spectrumset=pnS005-obj-720-840.pi \
       crossreg_spectrumset=pnS005-obj-540-720.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-720-840.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

    Skip to the end of the this section.

     evselect table=mos1S003-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=400 ycolumn=DETY yimagebinsize=400 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos1S003-720-840-720-840.arf spectrumset=mos1S003-obj-720-840.pi \
       crossreg_spectrumset=mos1S003-obj-720-840.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-720-840.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no
     arfgen arfset=mos1S003-720-840-540-720.arf spectrumset=mos1S003-obj-540-720.pi \
       crossreg_spectrumset=mos1S003-obj-720-840.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos1S003-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos1S003-clean.fits \
       modelootcorr=no

     evselect table=mos2S004-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=400 ycolumn=DETY yimagebinsize=400 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=mos2S004-720-840-720-840.arf spectrumset=mos2S004-obj-720-840.pi \
       crossreg_spectrumset=mos2S004-obj-720-840.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-720-840.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no
     arfgen arfset=mos2S004-720-840-540-720.arf spectrumset=mos2S004-obj-540-720.pi \
       crossreg_spectrumset=mos2S004-obj-720-840.pi crossregionarf=yes withrmfset=yes \
       rmfset=mos2S004-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=mos2S004-clean.fits \
       modelootcorr=no

     evselect table=pnS005-clean.fits:EVENTS ignorelegallimits=yes imageset=detmap.ds \
       xcolumn=DETX ximagebinsize=400 ycolumn=DETY yimagebinsize=400 squarepixels=yes \
       imagebinning=binSize withimageset=yes

     arfgen arfset=pnS005-720-840-720-840.arf spectrumset=pnS005-obj-720-840.pi \
       crossreg_spectrumset=pnS005-obj-720-840.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-720-840.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no
     arfgen arfset=pnS005-720-840-540-720.arf spectrumset=pnS005-obj-540-720.pi \
       crossreg_spectrumset=pnS005-obj-720-840.pi crossregionarf=yes withrmfset=yes \
       rmfset=pnS005-540-720.rmf extendedsource=yes modelee=no withbadpixcorr=yes \
       detmaptype=dataset detmaparray=detmap.ds badpixlocation=pnS005-clean.fits \
       modelootcorr=no

  9. At this point all of the XMM spectra (source and background) and calibration files (RMF and ARF) files have been created for the spectral fitting process. However, it must be remembered that there are other background components besides the QPB which still contribute to the source spectrum: the cosmic X-ray background and possibly solar wind charge exchange (typically only the OVII and OVIII lines are critical) and soft proton events. Each of these three background components can be addressed as part of the fitting process. For example (in pseudo Xspec terms):

     Model = cosmic + swcx + sp + source
       = (apec + wabs*(apec + apec + pow) + (gauss + gauss) + (pow) + (source)
       -> con

    Unfortunately, the fitting is more complicated than this. In Xspec V12.0 and higher, the soft proton component must be added as a separate model as it requires a separate diagonal response file. The cosmic background component needs to be constrained, which can be accomplished by simultaneously fitting a ROSAT All-Sky Survey (RASS) spectrum of the off-source background. The RASS spectrum and associated response file can be obtained using the HEASARC X-ray Background Tool. As extended emission is being studied (or you wouldn't be looking at this script), the various spectra need to be normalized to the same solid angle, per square arc minute, which is useful as the RASS spectra come in those units. A second constant scale factor can be added to the model and allowed to float to account for residual calibration discrepancies. Finally, additional Gaussian lines must be added to the model to account for the instrumental fluorescent lines associated with the QPD but are not in the model QPB spectra (the lines are so bright that the residual gain variation of an observation will cause significant residuals).

     Model = cosmic + swcx + sp + source
       = (apec + wabs*(apec + apec + pow) + (gauss + gauss) + (pow) + (source)
       -> (gauss + gauss + gauss + gauss + gauss + gauss + gauss) +
          con*con(gauss + gauss + apec + wabs*(apec + apec + pow) + source)

    To further complicate the model setup, all the spectra should be fit simultaneously to maximize the available constraints, e.g., the cosmic and SWCX background are the same for all XMM-Newton spectra and some parameters of the fluorescent line and SP background can be linked. This is all discussed to a much greater extent in the XMM-ESAS Cookbook.

    A sample Xspec XCM file for these data can be found here. Figure 5 shows the results of the fit.

    Figure 5: Spectral fit of the radial profile spectra of Abell 1795.


Last Updated: 19 May 2011


Caveats

For information about known issues, refer to the ESAS Warnings and Watchouts pages.


send this page to a friend
   Copyright 2013© European Space Agency. All rights reserved.
e-mail the webmaster    
This page was last updated on 9 May, 2013.