Data

+ System Overview
+ Software Download
+ Documentation
+ Cicerone
+ Analysis Threads
+ User Contributions

XML Model Definitions for Likelihood

Example model definitions for Likelihood:

PowerLaw


<source name="PowerLaw_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="PowerLaw">

<parameter free="1" max="1000.0" min="0.001" name="Prefactor" scale="1e-09" value="1"/>

<parameter free="1" max="-1.0" min="-5." name="Index" scale="1.0" value="-2.1"/>

<parameter free="0" max="2000.0" min="30.0" name="Scale" scale="1.0" value="100.0"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

BrokenPowerLaw


<source name="BrokenPowerLaw_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="BrokenPowerLaw">

<parameter free="1" max="1000.0" min="0.001" name="Prefactor" scale="1e-09" value="1"/>

<parameter free="1" max="-1.0" min="-5." name="Index1" scale="1.0" value="-1.8"/>

<parameter free="1" max="2000.0" min="30.0" name="BreakValue" scale="1.0" value="1000.0"/>

<parameter free="1" max="-1.0" min="-5." name="Index2" scale="1.0" value="-2.3"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

PowerLaw2


<source name="PowerLaw2_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="PowerLaw2">

<parameter free="1" max="1000.0" min="1e-05" name="Integral" scale="1e-06" value="1.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index" scale="1.0" value="-2.0"/>

<parameter free="0" max="200000.0" min="20.0" name="LowerLimit" scale="1.0" value="20.0"/>

<parameter free="0" max="200000.0" min="20.0" name="UpperLimit" scale="1.0" value="2e5"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>    </source>

BrokenPowerLaw2


<source name="BrokenPowerLaw2_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="BrokenPowerLaw2">

<parameter free="1" max="1000.0" min="0.001" name="Integral" scale="1e-04" value="1.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index1" scale="1.0" value="-1.8"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index2" scale="1.0" value="-2.3"/>

<parameter free="1" max="10000.0" min="30.0" name="BreakValue" scale="1.0" value="1000.0"/>

<parameter free="0" max="200000.0" min="20.0" name="LowerLimit" scale="1.0" value="20.0"/>

<parameter free="0" max="200000.0" min="20.0" name="UpperLimit" scale="1.0" value="2e5"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

SmoothBrokenPowerLaw


<source name="SmoothBrokenPowerLaw_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="SmoothBrokenPowerLaw">

<parameter free="1" max="1e10" min="0.0" name="Prefactor" scale="1e-06" value="1.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index1" scale="1.0" value="-2.0"/>

<parameter free="0" max="2000.0" min="30.0" name="Scale" scale="1.0" value="100.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index2" scale="1.0" value="-2.0"/>

<parameter free="1" max="5e5" min="20" name="BreakValue" scale="1.0" value="1e3"/>

<parameter free="1" max="10" min="0.01" name="Beta" scale="1.0" value="0.2"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

LogParabola


<source name="LogParabola_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="LogParabola">

<parameter free="1" max="1000.0" min="0.001" name="norm" scale="1e-9" value="1"/>

<parameter free="1" max="10" min="0" name="alpha" scale="1.0" value="1"/>

<parameter free="1" max="1e4" min="20" name="Eb" scale="1" value="300."/>

<parameter free="1" max="10" min="0" name="beta" scale="1.0" value="2"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

ExpCutoff


<source name="ExpCutoff_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="ExpCutoff">

<parameter free="1" max="100000.0" min="0.01" name="Prefactor" scale="1e-09" value="50"/>

<parameter free="1" max="-1.0" min="-5." name="Index" scale="1.0" value="-2.1"/>

<parameter free="0" max="2000.0" min="30.0" name="Scale" scale="1.0" value="100.0"/>

<parameter free="1" max="300.0" min="1.0" name="Ebreak" scale="1.0" value="10.0"/>

<parameter free="1" max="300.0" min="0.1" name="P1" scale="1000.0" value="100."/>

<parameter free="0" max="1.0" min="-1.0" name="P2" scale="1.0" value="0"/>

<parameter free="0" max="1.0" min="-1.0" name="P3" scale="1.0" value="0"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

BPLExpCutoff


<source name="BPLExpCutoff_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="BPLExpCutoff">

<parameter free="1" max="100000.0" min="0.01" name="Prefactor" scale="1e-09" value="1"/>

<parameter free="1" max="-1.001" min="-5." name="Index1" scale="1.0" value="-2.1"/>

<parameter free="1" max="-1.001" min="-5." name="Index2" scale="1.0" value="-2.1"/>

<parameter free="1" max="10000.0" min="1.0" name="BreakValue" scale="1.0" value="1000.0"/>

<parameter free="1" max="300.0" min="1.0" name="Eabs" scale="1.0" value="10.0"/>

<parameter free="1" max="300.0" min="0.1" name="P1" scale="1000.0" value="100."/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

Gaussian


<source name="Gaussian_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="Gaussian">

<parameter free="1" max="1000.0" min="0.001" name="Prefactor" scale="1e-09" value="1"/>

<parameter free="1" max="1e5" min="1e3" name="Mean" scale="1.0" value="7e4"/>

<parameter free="1" max="30" min="1e4" name="Sigma" scale="1.0" value="1e3"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

ConstantValue


<source name="ConstantValue_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="ConstantValue">

<parameter free="0" max="10" min="0" name="Value" scale="1" value="1"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

FileFunction


<source name="FileFunction_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum file="foo.dat" type="FileFunction">

<parameter free="1" max="1e5" min="1e-5" name="Normalization" scale="1" value="1.0"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

BandFunction


<source name="BandFunction_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="BandFunction">

<parameter free="1" max="1000.0" min="1e-05" name="norm" scale="1e-09" value="1.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="alpha" scale="1.0" value="-1.8"/>

<parameter free="1" max="-1.0" min="-5.0" name="beta" scale="1.0" value="-2.5"/>

<parameter free="1" max="-1.0" min="-5.0" name="Ep" scale="1." value="0.1"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

PLSuperExpCutoff


<source name="PLSuperExpCutoff_source" type="PointSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 -->

<spectrum type="PLSuperExpCutoff">

<parameter free="1" max="1000" min="1e-05" name="Prefactor" scale="1e-07" value="1"/>

<parameter free="1" max="0" min="-5" name="Index1" scale="1" value="-1.7"/>

<parameter free="0" max="1000" min="50" name="Scale" scale="1" value="200"/>

<parameter free="1" max="30000" min="500" name="Cutoff" scale="1" value="3000"/>

<parameter free="1" max="5" min="0" name="Index2" scale="1" value="1.5"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360." min="-360." name="RA" scale="1.0" value="83.45"/>

<parameter free="0" max="90." min="-90." name="DEC" scale="1.0" value="21.72"/>

</spatialModel>

</source>

DMFitFunction


<source name="DM_Example" type="PointSource">

<spectrum file="$(BASE_DIR)/data/Likelihood/gammamc_dif.dat" type="DMFitFunction">

<parameter error="1." free="0" max="1.e+5" min="1.e-5" name="norm" scale="1.e+20" value="5.0" />

<parameter error="1." free="0" max="5000.0" min="0." name="sigmav" scale="1.e-26" value="3.0" />

<parameter error="1." free="0" max="5000.0" min="1." name="mass" scale="1.0" value="10"/>

<parameter error="0.1" free="0" max="1.0" min="0.0" name="bratio" scale="1.0" value="1"/>

<parameter free="0" max="10" min="1" name="channel0" scale="1.0" value="4"/>

<parameter free="0" max="10" min="1" name="channel1" scale="1.0" value="1"/>

</spectrum>

<spatialModel type="SkyDirFunction">

<parameter free="0" max="360" min="-360" name="RA" scale="1.0" value="128.8272"/>

<parameter free="0" max="90" min="-90" name="DEC" scale="1.0" value="-45.1762"/>

</spatialModel>

</source>

Spatial Models

To use an extended spatial model such as SpatialMap or MapCubeFunction with Unbinned Likelihood, one must include the map_based_integral attribute in the spatialModel portion of the XML model definition. This attribute should not be used in a Binned Likelihood analysis. An example is provided below:

SpatialMap


<source name="SpatialMap_source" type="DiffuseSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 sr^-1 -->

<spectrum type="PowerLaw2">

<parameter free="1" max="1000.0" min="1e-05" name="Integral" scale="1e-06" value="1.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index" scale="1.0" value="-2.0"/>

<parameter free="0" max="200000.0" min="20.0" name="LowerLimit" scale="1.0" value="20.0"/>

<parameter free="0" max="200000.0" min="20.0" name="UpperLimit" scale="1.0" value="2e5"/>

</spectrum>

<spatialModel file="$(PATH_TO_FILE)/SpatialMap_source.fits" type="SpatialMap"  map_based_integral="true">

<parameter free="0" max="1000.0" min="0.001" name="Prefactor" scale= "1.0" value="1.0"/>

</spatialModel>

</source>

MapCubeFunction


<source name="MapCubeFunction_source" type="DiffuseSource">

<!-- point source units are cm^-2 s^-1 MeV^-1 sr^-1 -->

<spectrum type="PowerLaw2">

<parameter free="1" max="1000.0" min="1e-05" name="Integral" scale="1e-06" value="1.0"/>

<parameter free="1" max="-1.0" min="-5.0" name="Index" scale="1.0" value="-2.0"/>

<parameter free="0" max="200000.0" min="20.0" name="LowerLimit" scale="1.0" value="20.0"/>

<parameter free="0" max="200000.0" min="20.0" name="UpperLimit" scale="1.0" value="2e5"/>

</spectrum>

<spatialModel file="$(PATH_TO_FILES)/MapCubeFunction_source.fits" type="SpatialMap"  map_based_integral="true">

<parameter free="0" max="1000.0" min="0.001" name="Normalization" scale= "1.0" value="1.0"/>

</spatialModel>

</source>

National Aeronautics and Space Administration

Goddard Space Flight Center