Fermi Science Support Center
There will be a network outage, Sunday, April 26 from 13:00 EDT - Monday April 27 01:00 EDT. All services, web sites, and systems that require access to the internet will be affected by this work.

LAT Background Models

NOTE: The Galactic diffuse model for Pass 7 is singificantly larger than in previous releases. For some computers this can produce memory-related issues when running tools that use this model, such as gtlike, gtsrcprob, gtdiffrsp, and gtsrcmaps.

Galactic interstellar
emission model
Event Selection/
IRF Name
Isotropic spectral template
(see below for usage notes)

gll_iem_v05_rev1.fit (above)
replaces the previous version,
(details below)
Pass 7 Reprocessed Source (front+back)
Pass 7 Reprocessed Source (front only)
Pass 7 Reprocessed Source (back only)
Pass 7 Reprocessed Clean (front+back)
Pass 7 Reprocessed Clean (front only)
Pass 7 Reprocessed Clean (back only)
gal_2yearp7v6_v0.fits or
Pass 7 (V6) Source (front+back)
Pass 7 (V6) Source (front)
Pass 7 (V6) Source (back)
Pass 7 (V6) Clean (front+back)
Pass 7 (V6) Clean (front)
Pass 7 (V6) Clean (back)
gll_iem_v02_P6_V11_DIFFUSE.fit Pass 6 (V11) Diffuse (front+back)
Pass 6 (V11) Diffuse (front only)
Pass 6 (V11) Diffuse (back only)
gll_iem_v02.fit Pass 6 (V3) Diffuse (front+back)
Pass 6 (V3) Diffuse (front only)
Pass 6 (V3) Diffuse (back only)
Pass 6 (V3) DataClean (front+back)

Many analyses of LAT data require models of Galactic diffuse and isotropic emission (see the likelihood analysis tutorial for an example).

The Galactic Interstellar Emission Model

The Galactic model is a spatial and spectral template. The model for the Galactic diffuse emission was developed using spectral line surveys of HI and CO (as a tracer of H2) to derive the distribution of interstellar gas in Galactocentric rings. Infrared tracers of dust column density were used to correct column densities in directions where the optical depth of HI was either over or under-estimated. The model of the diffuse gamma-ray emission was then constructed by fitting the gamma-ray emissivities of the rings in several energy bands to the LAT observations. The fitting also required a model of the inverse Compton emission calculated using GALPROP and a model for the isotropic diffuse emission. A complete description of the Galactic Diffuse Model for Pass 7 Reprocessed Data is available.

June 3, 2014 Update: A new version of the Interstellar Emission Model for the analysis of P7REP data is available. This new version, gll_iem_v05_rev1, corrects a bug in gll_iem_v05 that was built with anomalously low molecular hydrogen emissivities in the outer Galaxy. More than 10° from the plane or in the longitude range l < 70° and l > 300°, the effect on sources is minimal. In the affected region below 10° latitude and in the longitude range 70° to 300°, it is advised to reanalyze sources with the corrected interstellar emission model. Differences are largest on top of CO clouds, and the scatter on source parameters before/after the correction is larger than the statistical error. A plot showing the ratio of the two models (new/old) in the affected region is available here.

Caveats related to the use of these diffuse models are regularly updated, and are available from the LAT Caveats page.

Isotropic Spectral Templates

An isotropic spectral template provides the spectral form from a fit to the all-sky emission (|b|>30 deg) that is not represented in the Galactic diffuse model and therefore includes both extragalactic diffuse gamma rays and remaining residual (misclassified) cosmic-ray emission. The Pass 7 Source class contains a significant background contribution at high energies. An isotropic model is determined under the assumption of a specified Galactic diffuse model and for one of the LAT event class selections. The resulting spectral template should be used with the same Galactic diffuse model and event class selection that were used to derive it. The three columns in the isotropic template are: the central energy (in MeV) for the band, the differential flux in that band (in ph/cm2/s/MeV/sr), and the 1-sigma uncertainty on the differential flux. The uncertainty is not used for likelihood analysis by the Science Tools.

The Pass 7 Ultraclean class does not have an isotropic spectral template because it is expected to be used to study the type of emission the template would model out of the data. Analysis using the Pass 7 Transient class typically requires case-specific estimates of the background and so no general template is provided.

Updating your XML files

Both the isotropic and Galactic diffuse models are now provided within the Science Tools release. This means your XML model will need to include the proper path to find the model files. You can place them in your local directory, or you can copy them into a single location for easy reference. It is necessary for your XML model to properly represent the source name (as is used in the LAT events files) as well as the proper path and filename for the model files. For this example, the models have been placed in the "/refdata/fermi/diffuseModels" directory under the Fermi Science Tools.

An example of the correct XML for these models is:

<source_library title="source library">

<source name="iso_source_v05" type="DiffuseSource">
<spectrum file="$(FERMI_DIR)/refdata/fermi/diffuseModels/iso_source_v05.txt" type="FileFunction">
<parameter free="1" max="1000" min="1e-05" name="Normalization" scale="1" value="1" />
<spatialModel type="ConstantValue">
<parameter free="0" max="10.0" min="0.0" name="Value" scale="1.0" value="1.0"/>

<source name="gll_iem_v05_rev1" type="DiffuseSource">
<spectrum type="ConstantValue">
<parameter free="1" max="10.0" min="0.0" name="Value" scale="1.0" value= "1.0"/>
<spatialModel file="$(FERMI_DIR)/refdata/fermi/diffuseModels/gll_iem_v05_rev1.fit" type="MapCubeFunction">
<parameter free="0" max="1000.0" min="0.001" name="Normalization" scale= "1.0" value="1.0"/>


The easiest way to generate XML models that use the new diffuse models is to run the user contributed tool "make2FGLxml.py" (available here). If the model files are in your current directory, a sample call that provides the proper output is:

>>> from make2FGLxml import *
This is make2FGLxml version 03.
NOTE: You must have run gtselect on the event file you use as input.
>>> mymodel = srcList('gll_psc_v08.fit','filtered_gti.fits','model.xml')
>>> mymodel.makeModel('gll_iem_v05_rev1.fit', 'gll_iem_v05_rev1', 'iso_source_v05.txt', 'iso_source_v05')
Creating file and adding sources for 2FGL

Using the Earth Limb Template

The 2-year LAT catalog analysis required an additional all-sky component to account for emission from the earth limb that was not completely removed from the data at energies below 200 MeV. These are gamma-rays that are in the broad tails of the PSF and so pass the 100 degree zenith angle cut. The spatial template was derived from residuals in the 50–68 MeV energy range and assumed that the spatial shape is independent of energy. The very soft spectral template was derived by adding the limb template to the model. The template is specific to the 2-year data set because the residual earth limb emission depends on the orientation of the LAT. You can find the Earth limb template FITS file on the 2-year LAT catalog data products page.

The Earth limb template is the only all-sky model component to use files for both the spatial and spectral models. It is intended for use only for analysis of the 2-year dataset used in the second LAT catalog. To include the Earth limb in your analysis, download the model files from the 2-year LAT catalog data products page. Add the following source to your XML model, and ensure the path information to each model component is correct.

<source name="EarthLimb" type="DiffuseSource">
<spectrum file="$(PATH_TO_FILES)/limb_2year_P76_source_v0_smooth.txt" type="FileFunction">
<parameter free="1" max="1000" min="1e-05" name="Normalization" scale="1" value="1" />
<spatialModel file="$(PATH_TO_FILES)/limb_2year_smooth.fits" type="SpatialMap">
<parameter free="0" max="1000.0" min="0.001" name="Normalization" scale= "1.0" value="1.0"/>