models_spatial_src.rst 11.5 KB

Edit Raw Blame History


Spatial source model components
The following sections present the spatial model components that are available
in ctools for gamma-ray sources.

Note
Except of the DiffuseMapCube model, all spatial models are normalized
so that when integrated over the sphere the result is unity. For clarify,
the spatial normalisation is omitted in the formulae below.

Point source

The PointSource model describes a point source
<source name="Crab" type="PointSource">
  <spatialModel type="PointSource">
    <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="1"/>
    <parameter name="DEC" scale="1.0" value="22.0145" min="-90"  max="90"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

RA is the Right Ascension (degrees)
DEC is the Declination (degrees)

The RA and DEC parameters can be remplaced by GLON and GLAT
for Galactic coordinates
<source name="Crab" type="PointSource">
  <spatialModel type="PointSource">
    <parameter name="GLON" scale="1.0" value="184.5575" min="-360" max="360" free="1"/>
    <parameter name="GLAT" scale="1.0" value="-5.7844"  min="-90"  max="90"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

GLON is the Galactic longitude (degrees)
GLAT is the Galactic latitude (degrees)


Note
For compatibility with the Fermi/LAT ScienceTools the model type
PointSource can be replaced by SkyDirFunction.


Radial source

The RadialDisk model describes a uniform intensity distribution within
a given radius
<source name="Crab" type="ExtendedSource">
  <spatialModel type="RadialDisk">
    <parameter name="RA"     scale="1.0" value="83.6331" min="-360" max="360" free="1"/>
    <parameter name="DEC"    scale="1.0" value="22.0145" min="-90"  max="90"  free="1"/>
    <parameter name="Radius" scale="1.0" value="0.20"    min="0.01" max="10"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

RA is the Right Ascension (degrees)
DEC is the Declination (degrees)
Radius is the disk radius (degrees)

The RadialGaussian model describes a Gaussian intensity distribution
<source name="Crab" type="ExtendedSource">
  <spatialModel type="RadialGaussian">
    <parameter name="RA"    scale="1.0" value="83.6331" min="-360" max="360" free="1"/>
    <parameter name="DEC"   scale="1.0" value="22.0145" min="-90"  max="90"  free="1"/>
    <parameter name="Sigma" scale="1.0" value="0.20"    min="0.01" max="10"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

and implements
M_{\rm spatial}(\theta) = \frac{1}{2 \pi \sigma^2} \exp
               \left(-\frac{1}{2}\frac{\theta^2}{\sigma^2} \right),

where

RA is the Right Ascension (degrees)
DEC is the Declination (degrees)
\sigma
 == Sigma (degrees)

The RadialShell model describes a spherical shell projected on the sky
<source name="Crab" type="ExtendedSource">
  <spatialModel type="RadialShell">
    <parameter name="RA"     scale="1.0" value="83.6331" min="-360" max="360" free="1"/>
    <parameter name="DEC"    scale="1.0" value="22.0145" min="-90"  max="90"  free="1"/>
    <parameter name="Radius" scale="1.0" value="0.30"    min="0.01" max="10"  free="1"/>
    <parameter name="Width"  scale="1.0" value="0.10"    min="0.01" max="10"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

and implements
M_{\rm spatial}(\theta) =  n_0 \left \{
\begin{array}{l l}
   \displaystyle
   \sqrt{ \theta_{\rm out}^2 - \theta^2 } - \sqrt{ \theta_{\rm in}^2 - \theta^2 }
   & \mbox{if $\theta \le \theta_{\rm in}$} \\
   \\
  \displaystyle
   \sqrt{ \theta_{\rm out}^2 - \theta^2 }
   & \mbox{if $\theta_{\rm in} < \theta \le \theta_{\rm out}$} \\
   \\
  \displaystyle
  0 & \mbox{if $\theta > \theta_{\rm out}$}
\end{array}
\right .

where

RA is the Right Ascension (degrees)
DEC is the Declination (degrees)
\theta_{\rm out}
 = Radius + Width (degrees)
\theta_{\rm in}
 = Radius (degrees)


Note
For all radial source models the RA and DEC parameters can be
remplaced by GLON and GLAT for Galactic coordinates.


Elliptical source

The EllipticalDisk model describes a uniform intensity distribution within
an elliptical circumference:
<source name="Crab" type="ExtendedSource">
  <spatialModel type="EllipticalDisk">
    <parameter name="RA"          scale="1.0" value="83.6331" min="-360"  max="360" free="1"/>
    <parameter name="DEC"         scale="1.0" value="22.0145" min="-90"   max="90"  free="1"/>
    <parameter name="PA"          scale="1.0" value="45.0"    min="-360"  max="360" free="1"/>
    <parameter name="MinorRadius" scale="1.0" value="0.5"     min="0.001" max="10"  free="1"/>
    <parameter name="MajorRadius" scale="1.0" value="2.0"     min="0.001" max="10"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

RA is the Right Ascension (degrees)
DEC is the Declination (degrees)
PA is the position angle, counted counterclockwise from North (degrees)
MinorRadius is the minor radius of the ellipse (degrees)
MajorRadius is the major radius of the ellipse (degrees)

The EllipticalGaussian model describes a Gaussian intensity distribution
<source name="Crab" type="ExtendedSource">
  <spatialModel type="EllipticalGaussian">
    <parameter name="RA"          scale="1.0" value="83.6331" min="-360"  max="360" free="1"/>
    <parameter name="DEC"         scale="1.0" value="22.0145" min="-90"   max="90"  free="1"/>
    <parameter name="PA"          scale="1.0" value="45.0"    min="-360"  max="360" free="1"/>
    <parameter name="MinorRadius" scale="1.0" value="0.5"     min="0.001" max="10"  free="1"/>
    <parameter name="MajorRadius" scale="1.0" value="2.0"     min="0.001" max="10"  free="1"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

and implements
M_{\rm spatial}(\theta, \phi) = \exp \left( -\frac{\theta^2}{2 r_\mathrm{eff}^2} \right),

with
r_\mathrm{eff} = \frac{ab} {\sqrt{\left( a \sin (\phi - \phi_0) \right)^2 +
                 \sqrt{\left( b \cos (\phi - \phi_0) \right)^2}}}

where

RA is the Right Ascension (degrees)
DEC is the Declination (degrees)
PA is the position angle, counted counterclockwise from North (degrees)
a
 = MinorRadius (degrees)
b
 = MajorRadius (degrees)
\phi_0
 is the position angle of the ellipse, counted counterclockwise
from North
\phi
 is the azimuth angle with respect to North.


Note
For all elliptical source models the RA and DEC parameters can be
remplaced by GLON and GLAT for Galactic coordinates.


Diffuse source

The DiffuseIsotropic model describes an isotropic intensity distribution
<source name="Crab" type="DiffuseSource">
  <spatialModel type="DiffuseIsotropic">
    <parameter name="Value" scale="1" value="1" min="1"  max="1" free="0"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

Value is isotropic intensity


Note
For compatibility with the Fermi/LAT ScienceTools the model type
DiffuseIsotropic can be replaced by ConstantValue.

The DiffuseMap model describes an arbitrary intensity distribution in
form of a sky map
<source name="Crab" type="DiffuseSource">
  <spatialModel type="DiffuseMap" file="map.fits">
    <parameter name="Normalization" scale="1" value="1" min="0.001" max="1000.0" free="0"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

Normalization is a normalization value

and the file attribute specifies a sky map FITS file name. If a file name
without path is specified it is assumed that the FITS file resides in the same
directory as the model definition XML file.

Note
For compatibility with the Fermi/LAT ScienceTools the model type
DiffuseMap can be replaced by SpatialMap and the parameter
Normalization can be replaced by Prefactor.


Note
You may add the normalize="0" to the spatial model tag to avoid the
normalization of the diffuse map to an integral of unity
<source name="Crab" type="DiffuseSource">
  <spatialModel type="DiffuseMap" file="map.fits" normalize="0">
    <parameter name="Normalization" scale="1" value="1" min="0.001" max="1000.0" free="0"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>


The DiffuseMapCube model describes an arbitrary energy-dependent intensity
distribution in form of a map cube
<source name="Crab" type="DiffuseSource">
  <spatialModel type="DiffuseMapCube" file="map_cube.fits">
    <parameter name="Normalization" scale="1" value="1" min="0.001" max="1000.0" free="0"/>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

where

Normalization is a normalization value


Note
For compatibility with the Fermi/LAT ScienceTools the model type
DiffuseMapCube can be replaced by MapCubeFunction and the parameter
Normalization can be replaced by Value.


Composite model

The Composite model implements a composite model that is the sum of an
arbitrary number of spatial models
<source name="Crab" type="CompositeSource">
  <spatialModel type="Composite">
    <spatialModel type="PointSource" component="PointSource">
      <parameter name="RA"    scale="1.0" value="83.6331" min="-360" max="360" free="1"/>
      <parameter name="DEC"   scale="1.0" value="22.0145" min="-90"  max="90"  free="1"/>
    </spatialModel>
    <spatialModel type="RadialGaussian">
      <parameter name="RA"    scale="1.0" value="83.6331" min="-360" max="360" free="1"/>
      <parameter name="DEC"   scale="1.0" value="22.0145" min="-90"  max="90"  free="1"/>
      <parameter name="Sigma" scale="1.0" value="0.20"    min="0.01" max="10"  free="1"/>
    </spatialModel>
  </spatialModel>
  <spectrum type="...">
    ...
  </spectrum>
</source>

which implements
M_{\rm spatial}(x,y|E) = \frac{1}{N} \sum_{i=0}^{N-1} M_{\rm spatial}^{(i)}(x,y|E)

where M_{\rm spatial}^{(i)}(x,y|E)
 is any spatial model component
(including another composite model), and N
 is the number of
model components that are combined.