PROTON

Description:

proton uses the spectral fitting results from Xspec and model soft proton detector maps to create model soft proton contamination maps for a given observation.

Requirements:

The current location should be the working directory.

mosspectra and mosback or pnspectra and pnback must have been successfully run for the exposure to be used. The spectra produced by these tasks must be fit to determine the spectral parameters of the soft proton contamination. Currently the spectrummust be modeled as a power law or a broken power law not folded through the instrumental effective areas. In Xspec V11 this is accomplished by including pow/b or bknpow/b in the fitted model. For Xspec V12 this requires adding an additional model and response where the response is a diagonal matrix which is included in the XMM-ESAS calibration files.

The input spectral parameters are simply the fitted values from Xspec (assuming that no additional normalization has been applied to the component).

Calling Parameters:

imagefile – Input object image file in detector coordinates.

specfile – File name of spectrum file used in the spectral fit to determine the residual SP contamination.

withsky – Parameter that does not currently work. If set to `yes', is supposed to rotate the image from detector to sky coordinates.

spmapdet - File name for the output SPF map in detector coordinates.

ccds – These are flags to include or exclude individual CCDs for the MOS detectors (e.g., “T T T T T T T”), or to include or exclude the quadrants of the pn (e.g., “T T T T”).

elow – The low energy limit for the energy band in eV (e.g., elow=400).

ehigh – The high energy limit for the energy band in eV (e.g., ehigh=1250).

speccontrol – 1 for a power law model, 2 for a broken power law.

spectrumcontrol=1

pindex – The fitted power law index.

pnorm – The fitted scale factor for the power law.

spectrumcontrol=2

bindl – The fitted lower power law index.

bbreak – The break energy, typically fixed to the default value of 3.0.

bindh – The fitted upper power law index.

bnorm – The fitted scale factor for the power law.

Output:

prefix-protimdet-elow-ehigh.fits – The model soft proton background image count image in detector coordinates.

Example:

proton imagefile=mos1S001-fovimdet-350-1100.fits specfile=mos1S001-fovt.pi elow=350 ehigh=1100 ccds="T T T F F F T" speccontrol=1 pindex=0.1 pnorm=2.208e-2

proton imagefile=mos1S001-fovimdet-350-1100.fits specfile=mos1S001-fovt.pi elow=350 ehigh=1100 ccds="T T T F F F T" speccontrol=2 bindl=0.85 bbreak=2.0 bindh=1.42 bnorm=0.025