Make compliant to latest GammaLib interface change (model XML format)

Jürgen Knödlseder
1 parent 673710bd
Showing 18 changed files with 185 additions and 144 deletions
doc/source/users/user_manual/getting_started/beyond_combine.rst
doc/source/users/user_manual/getting_started/beyond_irf.rst
doc/source/users/user_manual/getting_started/beyond_model.rst
doc/source/users/user_manual/getting_started/beyond_stacked.rst
doc/source/users/user_manual/getting_started/energy_dispersion.rst
doc/source/users/user_manual/getting_started/fitting.rst
doc/source/users/user_manual/getting_started/models.rst
doc/source/users/user_manual/getting_started/python.rst
doc/source/users/user_manual/getting_started/simulating.rst
doc/source/users/user_manual/getting_started/unbinned.rst
doc/source/users/user_manual/iact_analysis/analysis.rst
doc/source/users/user_manual/iact_analysis/high_level.rst
doc/source/users/user_manual/iact_analysis/simulation.rst
doc/source/users/user_manual/iact_analysis/tips.rst
doc/source/users/user_manual/introduction.rst
test/data/sciver/crab_const.xml
test/data/sciver/crab_plaw2.xml
test/science_verification.py
@@ -194,7 +194,7 @@ the log file:
   2016-06-29T19:14:38:  Instrument scale factors ..: unity
   2016-06-29T19:14:38:  Observation identifiers ...: all
   2016-06-29T19:14:38:  Model type ................: PointSource
-  2016-06-29T19:14:38:  Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+  2016-06-29T19:14:38:  Model components ..........: "PointSource" * "PowerLaw" * "Constant"
   2016-06-29T19:14:38:  Number of parameters ......: 6
   2016-06-29T19:14:38:  Number of spatial par's ...: 2
   2016-06-29T19:14:38:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -176,7 +176,7 @@ both observations.
   2016-06-29T20:10:30:  Instrument scale factors ..: unity
   2016-06-29T20:10:30:  Observation identifiers ...: all
   2016-06-29T20:10:30:  Model type ................: PointSource
-  2016-06-29T20:10:30:  Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+  2016-06-29T20:10:30:  Model components ..........: "PointSource" * "PowerLaw" * "Constant"
   2016-06-29T20:10:30:  Number of parameters ......: 6
   2016-06-29T20:10:30:  Number of spatial par's ...: 2
   2016-06-29T20:10:30:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -42,27 +42,27 @@ tied to the ``events2.fits`` file.
   <source_library title="source library">
     <source name="Crab" type="PointSource">
       <spectrum type="PowerLaw">
-         <parameter name="Prefactor" scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
-         <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-         <parameter name="Scale"     scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+         <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+         <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+         <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
       </spectrum>
-      <spatialModel type="SkyDirFunction">
+      <spatialModel type="PointSource">
         <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="0"/>
         <parameter name="DEC" scale="1.0" value="22.0145" min="-90"  max="90"  free="0"/>
       </spatialModel>
     </source>
     <source name="Background_00001" type="CTAIrfBackground" instrument="CTA" id="00001">
       <spectrum type="PowerLaw">
-        <parameter name="Prefactor" scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>
-        <parameter name="Index"     scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
-        <parameter name="Scale"     scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
+        <parameter name="Prefactor"   scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>
+        <parameter name="Index"       scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
+        <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
       </spectrum>
     </source>
     <source name="Background_00002" type="CTAIrfBackground" instrument="CTA" id="00002">
       <spectrum type="PowerLaw">
-        <parameter name="Prefactor" scale="1.0"  value="10.0"  min="1e-3" max="1e+3"   free="1"/>
-        <parameter name="Index"     scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
-        <parameter name="Scale"     scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
+        <parameter name="Prefactor"   scale="1.0"  value="10.0"  min="1e-3" max="1e+3"   free="1"/>
+        <parameter name="Index"       scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
+        <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
       </spectrum>
     </source>
   </source_library>
@@ -104,27 +104,27 @@ This is illustrated in the model definition file below where component
   <source_library title="source library">
     <source name="Crab" type="PointSource">
       <spectrum type="PowerLaw">
-         <parameter name="Prefactor" scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
-         <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-         <parameter name="Scale"     scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+         <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+         <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+         <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
       </spectrum>
-      <spatialModel type="SkyDirFunction">
+      <spatialModel type="PointSource">
         <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="0"/>
         <parameter name="DEC" scale="1.0" value="22.0145" min="-90"  max="90"  free="0"/>
       </spatialModel>
     </source>
     <source name="Background_00001" type="CTAIrfBackground" instrument="CTA" id="00001,00003,00004">
       <spectrum type="PowerLaw">
-        <parameter name="Prefactor" scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>
-        <parameter name="Index"     scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
-        <parameter name="Scale"     scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
+        <parameter name="Prefactor"   scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>
+        <parameter name="Index"       scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
+        <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
       </spectrum>
     </source>
     <source name="Background_00002" type="CTAIrfBackground" instrument="CTA" id="00002">
       <spectrum type="PowerLaw">
-        <parameter name="Prefactor" scale="1.0"  value="10.0"  min="1e-3" max="1e+3"   free="1"/>
-        <parameter name="Index"     scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
-        <parameter name="Scale"     scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
+        <parameter name="Prefactor"   scale="1.0"  value="10.0"  min="1e-3" max="1e+3"   free="1"/>
+        <parameter name="Index"       scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
+        <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
       </spectrum>
     </source>
   </source_library>
@@ -145,7 +145,7 @@ modelling works:
       <spectrum type="PowerLaw">
         <parameter name="Prefactor" value="5.7" error="0" scale="1e-16" min="1e-07" max="1000" free="1" />
         <parameter name="Index" value="2.48" error="0" scale="-1" min="0" max="5" free="1" />
-        <parameter name="Scale" value="0.3" scale="1e+06" min="0.01" max="1000" free="0" />
+        <parameter name="PivotEnergy" value="0.3" scale="1e+06" min="0.01" max="1000" free="0" />
       </spectrum>
       <spatialModel type="SkyDirFunction">
         <parameter name="RA" value="83.6331" scale="1" min="-360" max="360" free="0" />
@@ -156,7 +156,7 @@ modelling works:
       <spectrum type="PowerLaw">
         <parameter name="Prefactor" value="1" error="0" scale="1" min="0" free="1" />
         <parameter name="Index" value="0" error="0" scale="1" min="-10" max="10" free="1" />
-        <parameter name="Scale" value="1" scale="1e+06" free="0" />
+        <parameter name="PivotEnergy" value="1" scale="1e+06" free="0" />
       </spectrum>
     </source>
   </source_library>
@@ -256,7 +256,7 @@ below.
   2016-06-29T19:54:17:  Instrument scale factors ..: unity
   2016-06-29T19:54:17:  Observation identifiers ...: all
   2016-06-29T19:54:17:  Model type ................: PointSource
-  2016-06-29T19:54:17:  Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+  2016-06-29T19:54:17:  Model components ..........: "PointSource" * "PowerLaw" * "Constant"
   2016-06-29T19:54:17:  Number of parameters ......: 6
   2016-06-29T19:54:17:  Number of spatial par's ...: 2
   2016-06-29T19:54:17:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -102,7 +102,7 @@ Below you see the output from the ``ctlike.log`` file.
   2016-06-29T20:33:50:  Instrument scale factors ..: unity
   2016-06-29T20:33:50:  Observation identifiers ...: all
   2016-06-29T20:33:50:  Model type ................: PointSource
-  2016-06-29T20:33:50:  Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+  2016-06-29T20:33:50:  Model components ..........: "PointSource" * "PowerLaw" * "Constant"
   2016-06-29T20:33:50:  Number of parameters ......: 6
   2016-06-29T20:33:50:  Number of spatial par's ...: 2
   2016-06-29T20:33:50:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -42,7 +42,7 @@ Below we show the XML result file that has been produced by the run:
       <spectrum type="PowerLaw">
         <parameter name="Prefactor" value="5.79064" error="0.101828" scale="1e-16" min="1e-07" max="1000" free="1" />
         <parameter name="Index" value="2.47177" error="0.0154935" scale="-1" min="0" max="5" free="1" />
-        <parameter name="Scale" value="0.3" scale="1e+06" min="0.01" max="1000" free="0" />
+        <parameter name="PivotEnergy" value="0.3" scale="1e+06" min="0.01" max="1000" free="0" />
       </spectrum>
       <spatialModel type="SkyDirFunction">
         <parameter name="RA" value="83.6331" scale="1" min="-360" max="360" free="0" />
@@ -53,7 +53,7 @@ Below we show the XML result file that has been produced by the run:
       <spectrum type="PowerLaw">
         <parameter name="Prefactor" value="1.01642" error="0.0162226" scale="1" min="0.001" max="1000" free="1" />
         <parameter name="Index" value="0.00462945" error="0.00964737" scale="1" min="-5" max="5" free="1" />
-        <parameter name="Scale" value="1" scale="1e+06" min="0.01" max="1000" free="0" />
+        <parameter name="PivotEnergy" value="1" scale="1e+06" min="0.01" max="1000" free="0" />
       </spectrum>
     </source>
   </source_library>
@@ -112,7 +112,7 @@ this run:
   2016-06-29T16:11:08:  Instrument scale factors ..: unity
   2016-06-29T16:11:08:  Observation identifiers ...: all
   2016-06-29T16:11:08:  Model type ................: PointSource
-  2016-06-29T16:11:08:  Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+  2016-06-29T16:11:08:  Model components ..........: "PointSource" * "PowerLaw" * "Constant"
   2016-06-29T16:11:08:  Number of parameters ......: 6
   2016-06-29T16:11:08:  Number of spatial par's ...: 2
   2016-06-29T16:11:08:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -67,11 +67,11 @@ function that is tabulated in an ASCII file.
   <source_library title="source library">
     <source name="Crab" type="PointSource">
       <spectrum type="PowerLaw">
-         <parameter name="Prefactor" scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
-         <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-         <parameter name="Scale"     scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+         <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+         <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+         <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
       </spectrum>
-      <spatialModel type="SkyDirFunction">
+      <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>
@@ -144,7 +144,7 @@ Point source
   .. code-block:: xml
     <source name="Crab" type="PointSource">
-      <spatialModel type="SkyDirFunction">
+      <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>
@@ -153,6 +153,11 @@ Point source
       </spectrum>
     </source>
+  .. note::
+
+    For compatibility with the Fermi/LAT ScienceTools the model type ``PointSource``
+    can be replaced by ``SkyDirFunction``.
+
 Radial source
 ^^^^^^^^^^^^^
@@ -160,7 +165,7 @@ Radial source
   .. code-block:: xml
     <source name="Crab" type="ExtendedSource">
-      <spatialModel type="DiskFunction">
+      <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"/>
@@ -173,7 +178,7 @@ Radial source
   .. code-block:: xml
     <source name="Crab" type="ExtendedSource">
-      <spatialModel type="GaussFunction">
+      <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"/>
@@ -186,7 +191,7 @@ Radial source
   .. code-block:: xml
     <source name="Crab" type="ExtendedSource">
-      <spatialModel type="ShellFunction">
+      <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"/>
@@ -219,7 +224,7 @@ Elliptical source
   .. code-block:: xml
     <source name="Crab" type="ExtendedSource">
-      <spatialModel type="EllipticalGauss">
+      <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"/>
@@ -238,7 +243,7 @@ Diffuse source
   .. code-block:: xml
     <source name="Crab" type="DiffuseSource">
-      <spatialModel type="ConstantValue">
+      <spatialModel type="DiffuseIsotropic">
          <parameter name="Value" scale="1" value="1" min="1"  max="1" free="0"/>
       </spatialModel>
       <spectrum type="...">
@@ -246,21 +251,32 @@ Diffuse source
       </spectrum>
     </source>
+  .. note::
+
+    For compatibility with the Fermi/LAT ScienceTools the model type
+    ``DiffuseIsotropic`` can be replaced by ``ConstantValue``.
+
   .. code-block:: xml
     <source name="Crab" type="DiffuseSource">
-      <spatialModel type="SpatialMap" file="map.fits">
-         <parameter name="Prefactor" scale="1" value="1" min="0.001" max="1000.0" free="0"/>
+      <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>
+  .. 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``.
+
   .. code-block:: xml
     <source name="Crab" type="DiffuseSource">
-      <spatialModel type="MapCubeFunction" file="map_cube.fits">
+      <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="...">
@@ -268,6 +284,12 @@ Diffuse source
       </spectrum>
     </source>
+  .. 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``.
+
 CTA radial background
 ^^^^^^^^^^^^^^^^^^^^^
@@ -358,10 +380,10 @@ in ctools.
    Exceptions to this rule exist for the following spatial models:
-   - ``ConstantValue``: intensity units are given per steradian, e.g.
+   - ``DiffuseIsotropic``: intensity units are given per steradian, e.g.
      :math:`{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}\,{\rm sr}^{-1}`
-   - ``MapCubeFunction``: intensities are unitless, the spectral model 
+   - ``DiffuseMapCube``: intensities are unitless, the spectral model
      presents a relative scaling of the diffuse model
    If spectral models are used for a background model component, intensity 
@@ -395,7 +417,7 @@ Constant
   where
   * :math:`N_0` = ``Normalization``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1})`
@@ -405,9 +427,9 @@ Power law
   .. code-block:: xml
    <spectrum type="PowerLaw">
-     <parameter name="Prefactor" scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
-     <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-     <parameter name="Scale"     scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+     <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+     <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+     <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
    </spectrum>
   The ``PowerLaw`` model component implements the power law function
@@ -418,28 +440,33 @@ Power law
   where
   * :math:`k_0` = ``Prefactor``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1})`
   * :math:`\gamma` = ``Index``
-  * :math:`E_0` = ``Scale``
-    :math:`[{\rm MeV}]`
+  * :math:`E_0` = ``PivotEnergy``
+    :math:`({\rm MeV})`
   .. warning::
-    The ``Scale`` parameter is not intended to be fitted.
+    The ``PivotEnergy`` parameter is not intended to be fitted.
+
+  .. note::
+
+    For compatibility with the Fermi/LAT ScienceTools the parameter
+    ``PivotEnergy`` can be replaced by ``Scale``.
-  An alternative power law function that uses the integral flux as parameter
+  An alternative power law function that uses the integral photon flux as parameter
   rather than the Prefactor is specified by
   .. code-block:: xml
-   <spectrum type="PowerLaw2">
-     <parameter scale="1e-07" name="Integral"   min="1e-07" max="1000.0"    value="1.0" free="1"/>
+   <spectrum type="PowerLawPhotonFlux">
+     <parameter scale="1e-07" name="PhotonFlux" min="1e-07" max="1000.0"    value="1.0" free="1"/>
      <parameter scale="1.0"   name="Index"      min="-5.0"  max="+5.0"      value="-2.0" free="1"/>
      <parameter scale="1.0"   name="LowerLimit" min="10.0"  max="1000000.0" value="100.0" free="0"/>
      <parameter scale="1.0"   name="UpperLimit" min="10.0"  max="1000000.0" value="500000.0" free="0"/>
    </spectrum>
-  The ``PowerLaw2`` model component implements the power law function
+  The ``PowerLawPhotonFlux`` model component implements the power law function
   .. math::
     \frac{dN}{dE} = \frac{N(\gamma+1)E^{\gamma}}
@@ -447,22 +474,28 @@ Power law
   where
-  * :math:`N` = ``Integral``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}]`
+  * :math:`N` = ``PhotonFlux``
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1})`
   * :math:`\gamma` = ``Index``
   * :math:`E_{\rm min}` = ``LowerLimit``
-    :math:`[{\rm MeV}]`
+    :math:`({\rm MeV})`
   * :math:`E_{\rm max}` = ``UpperLimit``
-    :math:`[{\rm MeV}]`
+    :math:`({\rm MeV})`
   .. warning::
     The ``LowerLimit`` and ``UpperLimit`` parameters are always treated as fixed and,
-    the flux given by the ``Integral`` parameter is computed over the 
+    the flux given by the ``PhotonFlux`` parameter is computed over the
     range set by these two parameters.
     Use of this model allows the errors on the integral flux to be evaluated directly
     by :ref:`ctlike`.
+  .. note::
+
+    For compatibility with the Fermi/LAT ScienceTools the model type
+    ``PowerLawPhotonFlux`` can be replaced by ``PowerLaw2`` and the parameter
+    ``PhotonFlux`` can be replaced by ``Integral``.
+
 Exponentially cut-off power law
 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
@@ -470,10 +503,10 @@ Exponentially cut-off power law
   .. code-block:: xml
    <spectrum type="ExpCutoff">
-     <parameter name="Prefactor" scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
-     <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-     <parameter name="Cutoff"    scale="1e6"   value="1.0"  min="0.01"  max="1000.0" free="1"/>
-     <parameter name="Scale"     scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+     <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+     <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+     <parameter name="Cutoff"      scale="1e6"   value="1.0"  min="0.01"  max="1000.0" free="1"/>
+     <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
    </spectrum>
   The ``ExpCutoff`` model component implements the exponentially cut-off power law function
@@ -485,17 +518,21 @@ Exponentially cut-off power law
   where
   * :math:`k_0` = ``Prefactor``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1})`
   * :math:`\gamma` = ``Index``
-  * :math:`E_0` = ``Scale``
-    :math:`[{\rm MeV}]`
+  * :math:`E_0` = ``PivotEnergy``
+    :math:`({\rm MeV})`
   * :math:`E_{\rm cut}` = ``Cutoff``
-    :math:`[{\rm MeV}]`
+    :math:`({\rm MeV})`
   .. warning::
-    The ``Scale`` parameter is not intended to be fitted.
+    The ``PivotEnergy`` parameter is not intended to be fitted.
+
+  .. note::
+    For compatibility with the Fermi/LAT ScienceTools the parameter
+    ``PivotEnergy`` can be replaced by ``Scale``.
 Broken power law
@@ -524,11 +561,11 @@ Broken power law
   where
   * :math:`k_0` = ``Prefactor``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1})`
   * :math:`\gamma_1` = ``Index1``
   * :math:`\gamma_2` = ``Index2``
   * :math:`E_b` = ``BreakValue``
-    :math:`[{\rm MeV}]`
+    :math:`({\rm MeV})`
   .. warning::
@@ -543,10 +580,10 @@ Log parabola
   .. code-block:: xml
    <spectrum type="LogParabola">
-     <parameter name="Prefactor" scale="1e-17" value="5.878"   min="1e-07" max="1000.0" free="1"/>
-     <parameter name="Index"     scale="-1"    value="2.32473" min="0.0"   max="+5.0"   free="1"/>
-     <parameter name="Curvature" scale="-1"    value="0.074"   min="-5.0"  max="+5.0"   free="1"/>
-     <parameter name="Scale"     scale="1e6"   value="1.0"     min="0.01"  max="1000.0" free="0"/>
+     <parameter name="Prefactor"   scale="1e-17" value="5.878"   min="1e-07" max="1000.0" free="1"/>
+     <parameter name="Index"       scale="-1"    value="2.32473" min="0.0"   max="+5.0"   free="1"/>
+     <parameter name="Curvature"   scale="-1"    value="0.074"   min="-5.0"  max="+5.0"   free="1"/>
+     <parameter name="PivotEnergy" scale="1e6"   value="1.0"     min="0.01"  max="1000.0" free="0"/>
    </spectrum>
   The ``LogParabola`` model component implements the log parabola function
@@ -557,25 +594,26 @@ Log parabola
   where
   * :math:`k_0` = ``Prefactor``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1})`
   * :math:`\gamma` = ``Index``
   * :math:`\eta` = ``Curvature``
-  * :math:`E_0` = ``Scale``
-    :math:`[{\rm MeV}]`
+  * :math:`E_0` = ``PivotEnergy``
+    :math:`({\rm MeV})`
   .. warning::
-    The ``Scale`` parameter is not intended to be fitted.
+    The ``PivotEnergy`` parameter is not intended to be fitted.
-  An alternative XML format is supported for compatibility with the Fermi/LAT XML format:
+  An alternative XML format is supported for compatibility with the Fermi/LAT
+  XML format:
   .. code-block:: xml
    <spectrum type="LogParabola">
-     <parameter name="Prefactor" scale="1e-17" value="5.878"   min="1e-07" max="1000.0" free="1"/>
-     <parameter name="alpha"     scale="1"     value="2.32473" min="0.0"   max="+5.0"   free="1"/>
-     <parameter name="beta"      scale="1"     value="0.074"   min="-5.0"  max="+5.0"   free="1"/>
-     <parameter name="Scale"     scale="1e6"   value="1.0"     min="0.01"  max="1000.0" free="0"/>
+     <parameter name="norm"  scale="1e-17" value="5.878"   min="1e-07" max="1000.0" free="1"/>
+     <parameter name="alpha" scale="1"     value="2.32473" min="0.0"   max="+5.0"   free="1"/>
+     <parameter name="beta"  scale="1"     value="0.074"   min="-5.0"  max="+5.0"   free="1"/>
+     <parameter name="Eb"    scale="1e6"   value="1.0"     min="0.01"  max="1000.0" free="0"/>
    </spectrum>
   where
@@ -604,11 +642,11 @@ Gaussian
   where
   * :math:`N_0` = ``Normalization``
-    :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}]`
+    :math:`({\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1})`
   * :math:`\bar{E}` = ``Mean``
-    :math:`[{\rm MeV}]`
+    :math:`({\rm MeV})`
   * :math:`\sigma` = ``Sigma``
-    :math:`[{\rm MeV}]`
+    :math:`({\rm MeV})`
 File function
@@ -625,9 +663,9 @@ File function
   The energy and intensity values are defined using an ASCII file with
   columns of energy and differential flux values.
   Energies are given in units of
-  :math:`[{\rm MeV}]`,
+  :math:`{\rm MeV}`,
   intensities are given in units of
-  :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`.
+  :math:`{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}`.
   The only parameter is a multiplicative normalization:
   .. math::
@@ -665,9 +703,9 @@ Node function
   power law which is defined by a set of energy and intensity values
   (the so called nodes) that are piecewise connected by power laws.
   Energies are given in units of
-  :math:`[{\rm MeV}]`,
+  :math:`{\rm MeV}`,
   intensities are given in units of
-  :math:`[{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}]`.
+  :math:`{\rm ph}\,\,{\rm cm}^{-2}\,{\rm s}^{-1}\,{\rm MeV}^{-1}`.
   .. warning::
@@ -170,7 +170,7 @@ member of the observation container:
     Instrument scale factors ..: unity
     Observation identifiers ...: all
     Model type ................: PointSource
-    Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+    Model components ..........: "PointSource" * "PowerLaw" * "Constant"
     Number of parameters ......: 6
     Number of spatial par's ...: 2
      RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -68,9 +68,9 @@ The source and background model is defined by the XML file
   <source_library title="source library">
     <source name="Crab" type="PointSource">
       <spectrum type="PowerLaw">
-         <parameter name="Prefactor" scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
-         <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-         <parameter name="Scale"     scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+         <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+         <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+         <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
       </spectrum>
       <spatialModel type="SkyDirFunction">
         <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="0"/>
@@ -78,12 +78,12 @@ The source and background model is defined by the XML file
       </spatialModel>
     </source>
     <source name="CTABackgroundModel" type="CTAIrfBackground" instrument="CTA">
-      <spectrum type="PowerLaw">  
-        <parameter name="Prefactor" scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>  
-        <parameter name="Index"     scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>  
-        <parameter name="Scale"     scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>  
+      <spectrum type="PowerLaw">
+        <parameter name="Prefactor"   scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>
+        <parameter name="Index"       scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
+        <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
       </spectrum>
-    </source>     
+    </source>
   </source_library>
 The model consists of a source library that contains two components:
@@ -103,8 +103,8 @@ maximum (``max``) value, and a parameter may be free (``free=&quot;1&quot;``) or fixed
 relevant here for the simulations, but they will be important for the model 
 fitting later.
-The spectral intensity I(E) (in units of photons/cm2/s/MeV) of the power
-law is given by 
+The spectral intensity I(E) (in units of photons/cm2/s/MeV) of the power law is
+given by
 .. math::
     \frac{dN}{dE} = N_0 \left( \frac{E}{E_0} \right)^{\gamma}
@@ -113,7 +113,7 @@ where the parameters in the XML definition have the following mappings:
 * :math:`N_0` = ``Prefactor``
 * :math:`\gamma` = ``Index``
-* :math:`E_0` = ``Scale``
+* :math:`E_0` = ``PivotEnergy``
 .. warning::
@@ -175,7 +175,7 @@ the same events for both analyses.
   2016-06-29T16:24:45:  Instrument scale factors ..: unity
   2016-06-29T16:24:45:  Observation identifiers ...: all
   2016-06-29T16:24:45:  Model type ................: PointSource
-  2016-06-29T16:24:45:  Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
+  2016-06-29T16:24:45:  Model components ..........: "PointSource" * "PowerLaw" * "Constant"
   2016-06-29T16:24:45:  Number of parameters ......: 6
   2016-06-29T16:24:45:  Number of spatial par's ...: 2
   2016-06-29T16:24:45:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
@@ -229,31 +229,31 @@ Model definition XML file
 	<source_library title="source library">
 	  <source name="bkg_11111" type="CTAIrfBackground" instrument="HESS" id="11111">
 	    <spectrum type="ConstantValue">
-	      <parameter name="Value" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
+	      <parameter name="Normalization" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
 	    </spectrum>
 	  </source>
 	  <source name="bkg_11112" type="CTAIrfBackground" instrument="HESS" id="11112">
 	    <spectrum type="ConstantValue">
-	      <parameter name="Value" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
+	      <parameter name="Normalization" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
 	    </spectrum>
 	  </source>
 	  <source name="bkg_11113" type="CTAIrfBackground" instrument="HESS" id="11113">
 	    <spectrum type="ConstantValue">
-	      <parameter name="Value" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
+	      <parameter name="Normalization" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
 	    </spectrum>
 	  </source>
 	  <source name="bkg_11114" type="CTAIrfBackground" instrument="HESS" id="11114">
 	    <spectrum type="ConstantValue">
-	      <parameter name="Value" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
+	      <parameter name="Normalization" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
 	    </spectrum>
 	  </source>
 	    <source name="Crab" type="PointSource">
 	    <spectrum type="PowerLaw">
-	       <parameter name="Prefactor" scale="1e-17" value="3.0"  min="1e-07" max="1000.0" free="1"/>
-	       <parameter name="Index"     scale="-1"    value="2.5" min="0.0"   max="+5.0"   free="1"/>
-	       <parameter name="Scale"     scale="1e6"   value="1.0"  min="0.01"  max="1000.0" free="0"/>
+	       <parameter name="Prefactor"   scale="1e-17" value="3.0"  min="1e-07" max="1000.0" free="1"/>
+	       <parameter name="Index"       scale="-1"    value="2.5" min="0.0"   max="+5.0"   free="1"/>
+	       <parameter name="PivotEnergy" scale="1e6"   value="1.0"  min="0.01"  max="1000.0" free="0"/>
 	    </spectrum>
-	    <spatialModel type="SkyDirFunction">
+	    <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>
@@ -514,18 +514,18 @@ The output model of :ref:`ctbkgcube` looks the following:
 	<source_library title="source library">
 	  <source name="BackgroundModel" type="CTACubeBackground" instrument="CTA,HESS,MAGIC,VERITAS">
 	    <spectrum type="PowerLaw">
-	      <parameter name="Prefactor" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
-	      <parameter name="Index" value="0" error="0" scale="1" min="-5" max="5" free="1" />
-	      <parameter name="Scale" value="1" scale="1e+06" free="0" />
+	      <parameter name="Prefactor"   value="1" error="0" scale="1" min="0.01" max="100" free="1" />
+	      <parameter name="Index"       value="0" error="0" scale="1" min="-5" max="5" free="1" />
+	      <parameter name="PivotEnergy" value="1" scale="1e+06" free="0" />
 	    </spectrum>
 	  </source>
 	  <source name="Crab" type="PointSource">
 	    <spectrum type="PowerLaw">
-	       <parameter name="Prefactor" scale="1e-17" value="3.0"  min="1e-07" max="1000.0" free="1"/>
-	       <parameter name="Index"     scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
-	       <parameter name="Scale"     scale="1e6"   value="1."  min="0.01"  max="1000.0" free="0"/>
+	       <parameter name="Prefactor"   scale="1e-17" value="3.0"  min="1e-07" max="1000.0" free="1"/>
+	       <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+	       <parameter name="PivotEnergy" scale="1e6"   value="1."  min="0.01"  max="1000.0" free="0"/>
 	    </spectrum>
-	    <spatialModel type="SkyDirFunction">
+	    <spatialModel type="PointSource">
 	      <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="0"/>
 	      <parameter name="DEC" scale="1.0" value="22.0145" min="-90"  max="90"  free="0"/>
 	    </spatialModel>
@@ -3,7 +3,8 @@
 High level analysis tools for IACT data
 =======================================
-This section describes and explains how to create high level output. The following tools are explained:
+This section describes and explains how to create high level output. The
+following scripts and tools are explained:
 * :ref:`csobsinfo` (inspect and observation XML files)
 * :ref:`csmodelinfo` (inspect model XML files)
@@ -16,9 +17,10 @@ This section describes and explains how to create high level output. The followi
 Inspecting observation definition files
 ---------------------------------------
-After having created an observation definition file, the user might want to know a summary of its content.
-Important quantities usually are the number of observations, the total livetime, or the mean zenith angle.
-For this purpose, :ref:`csobsinfo` inspects an observation XML container for its content.
+After having created an observation definition file, the user might want to
+know a summary of its content. Important quantities usually are the number of
+observations, the total livetime, or the mean zenith angle. For this purpose,
+:ref:`csobsinfo` inspects an observation XML container for its content.
 .. code-block:: bash
@@ -142,20 +142,20 @@ e.g. for the Crab Nebula:
 	<source_library title="source library">
 	  <source name="Crab" type="PointSource">
 	    <spectrum type="PowerLaw">
-	       <parameter name="Prefactor" scale="1e-17" value="3.5"  min="1e-07" max="1000.0" free="1"/>
-	       <parameter name="Index"     scale="-1"    value="2.5"  min="0.0"   max="+5.0"   free="1"/>
-	       <parameter name="Scale"     scale="1e6"   value="1.0"  min="0.01"  max="1000.0" free="0"/>
+	       <parameter name="Prefactor"   scale="1e-17" value="3.5"  min="1e-07" max="1000.0" free="1"/>
+	       <parameter name="Index"       scale="-1"    value="2.5"  min="0.0"   max="+5.0"   free="1"/>
+	       <parameter name="PivotEnergy" scale="1e6"   value="1.0"  min="0.01"  max="1000.0" free="0"/>
 	    </spectrum>
-	    <spatialModel type="SkyDirFunction">
+	    <spatialModel type="PointSource">
 	      <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="0"/>
 	      <parameter name="DEC" scale="1.0" value="22.0145" min="-90"  max="90"  free="0"/>
 	    </spatialModel>
 	  </source>
 	  <source name="CTABackgroundModel" type="CTAIrfBackground" instrument="CTA">
 	    <spectrum type="PowerLaw">	
-	      <parameter name="Prefactor" scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>	
-	      <parameter name="Index"     scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>	
-	      <parameter name="Scale"     scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>	
+	      <parameter name="Prefactor"   scale="1.0"  value="1.0"  min="1e-3" max="1e+3"   free="1"/>	
+	      <parameter name="Index"       scale="1.0"  value="0.0"  min="-5.0" max="+5.0"   free="1"/>
+	      <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
 	    </spectrum>
 	  </source>	
 	</source_library>
@@ -264,20 +264,20 @@ excess map. In the section about :ref:`high level analysis tools &lt;sec_high_level
 the usage of :ref:`cttsmap` was demonstrated. This tool can also be used to
 generate flux maps. For this, the user has to incorporate a test source in the
 XML model file. This source should consist of a spectral component using a
-PowerLaw2 model:
+``PowerLawPhotonFlux`` model:
 .. code-block:: xml
 	<source name="TestSource" type="PointSource">
-	 <spectrum type="PowerLaw2">
-	  <parameter scale="1e-07" name="Integral"   min="1e-07" max="1000.0"    value="1.0" free="1"/>
+	 <spectrum type="PowerLawPhotonFlux">
+	  <parameter scale="1e-07" name="PhotonFlux" min="1e-07" max="1000.0"    value="1.0" free="1"/>
 	  <parameter scale="1.0"   name="Index"      min="-5.0"  max="+5.0"      value="-2.0" free="1"/>
 	  <parameter scale="1.0"   name="LowerLimit" min="10.0"  max="1000000.0" value="100.0" free="0"/>
 	  <parameter scale="1.0"   name="UpperLimit" min="10.0"  max="1000000.0" value="500000.0" free="0"/>
 	 </spectrum>
-	 <spatialModel type="SkyDirFunction">
-	  <parameter free="0" max="360" min="-360" name="RA" scale="1" value="83.6331" />
-	  <parameter free="0" max="90" min="-90" name="DEC" scale="1" value="22.0145" />
+	 <spatialModel type="PointSource">
+	  <parameter free="0" max="360" min="-360" name="RA"  scale="1" value="83.6331" />
+	  <parameter free="0" max="90"  min="-90"  name="DEC" scale="1" value="22.0145" />
 	 </spatialModel>
 	</source>
@@ -20,7 +20,8 @@ also be specified as arguments separated by a whitespace to the ctool:
 .. code-block:: bash
-  $ ctobssim ra=83.63 dec=22.01 rad=5.0 tmin=0 tmax=1800 emin=0.1 emax=100.0 caldb=prod2 irf=South_50h inmodel=$CTOOLS/share/models/crab.xml outevents=events.fits
+  $ ctobssim ra=83.63 dec=22.01 rad=5.0 tmin=0 tmax=1800 emin=0.1 emax=100.0 caldb=prod2 irf=South_0.5h \
+             inmodel=$CTOOLS/share/models/crab.xml outevents=events.fits
 If you need help about the usage of a ctool or a cscript, type the name of 
 the ctool or cscript followed by the ``--help`` option:
@@ -7,7 +7,7 @@
       <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
     </spectrum>
     <spatialModel type="DiffuseIsotropic">
-      <parameter name="Normalization" scale="1" value="1" min="0.1" max="10" free="0"/>
+      <parameter name="Value" scale="1" value="1" min="0.1" max="10" free="0"/>
     </spatialModel>
   </source>
   <source name="Background" type="CTAIrfBackground" instrument="CTA">
@@ -2,7 +2,7 @@
 <source_library title="source library">
   <source name="Crab" type="PointSource">
     <spectrum type="PowerLawPhotonFlux">
-      <parameter name="PhotonFlux" scale="1e-5" value="1.617"      min="1e-07" max="1000.0"    free="1"/>
+      <parameter name="PhotonFlux" scale="1e-5" value="1.617"    min="1e-07" max="1000.0"    free="1"/>
       <parameter name="Index"      scale="-1.0" value="2.48"     min="-5.0"  max="+5.0"      free="1"/>
       <parameter name="LowerLimit" scale="1.0"  value="100.0"    min="10.0"  max="1000000.0" free="0"/>
       <parameter name="UpperLimit" scale="1.0"  value="500000.0" min="10.0"  max="1000000.0" free="0"/>
@@ -26,13 +26,13 @@
 # along with this program.  If not, see <http://www.gnu.org/licenses/>.
 #
 # ==========================================================================
-import gammalib
-import ctools
-import cscripts
 import os
 import csv
 import math
 import sys
+import gammalib
+import ctools
+import cscripts
 # ========================== #
@@ -375,7 +375,7 @@ class sciver(gammalib.GPythonTestSuite):
         Test power law 2 model
         """
         self.pull('data/sciver/crab_plaw2')
-        self.test('Pull_Crab_Integral')
+        self.test('Pull_Crab_PhotonFlux')
         self.test('Pull_Crab_Index')
         self.test('Pull_Background_Prefactor')
         self.test('Pull_Background_Index')