LAT Data Selection Recommendations (adapted from the Cicerone)

These recommendations apply only to an analysis of LAT data using the P6_V3 instrument response functions.

Limitations for Data Use. When using Fermi data, it is vital to recognize the limitations associated with these analyses. Each instrument team maintains a set of caveats which must be understood by the user. These are listed on the following pages, and will be kept up-to-date as the teams refine their understanding of the data, and develop their science processing. (See Caveats on the Use of LAT Data.)

Event Selection Recommendations (P6_V3)

Analysis Type Minimum energy
(emin)
Maximum energy
(emax)
Max zenith angle
(zmax)
Minimum event class
(evclsmin)
Maximum event class
(evclsmax)
Galactic Point Source Analysis
100 (MeV)
-
105 (degrees)
3
3
Off-plane Point Source Analysis
100 (MeV)
-
105 (degrees)
3
3
Transient Spectral Analysis (<200s)
100 (MeV)
-
105 (degrees)
1
3
Galactic Diffuse Analysis
100 (MeV)
100000 (MeV)
105 (degrees)
3
3

Time Selection Recommendations

Analysis Type ROI-based zenith angle cut
(roicut)
Relational filter expression
(filter)
Galactic Point Source Analysis
yes  
"DATA_QUAL==1"
Off-plane Point Source Analysis
yes  
"DATA_QUAL==1"
Transient Spectral Analysis
yes  
"DATA_QUAL==1"
Galactic Diffuse Analysis
   no**
"DATA_QUAL==1"
** An exposure correction must be made using the zmax option in gtltcube

Prior to performing any analysis on Fermi data, it is necessary to select the set of photons that are best suited for the science being performed. In most cases, the set of data selection criteria (cuts) you make will depend on the analysis being performed. Since the photon data is meant to serve for any analysis type, you will need to process the cuts before you start your analysis. There are several standard analysis types, each requiring their own sets of data cuts. These are:

  • Point Source Spectral analysis
  • Transient Spectral analysis (for timescales < 200 seconds)
  • Timing analysis (pulsar phase, GRB temporal analysis, etc.)
  • Galactic Diffuse analysis
  • Extragalactic Diffuse analysis
    Note: This analysis is not yet supported due to ongoing validation of the LAT instrument.

There are many configurable parameters that can be used in preparing the LAT data for analysis. These parameters are inputs in various tools like gtselect. Certain cuts are strongly recommended for all data analysis processes, while others are flexible, depending on the type of analysis being performed. The tables above outline the recommended cuts to be used for the event selection (gtselect) and time selection (gtmktime) appropriate for the different analysis types. Analysis of more complex sources, or confused regions, will require optimization of some of these parameters beyond the values listed below. If a parameter is not listed, then there is no specific recommended value for its use within that analysis chain.

Caution: Once you have selected a value for a parameter within your analysis, it is strongly recommended that you do not modify that value throughout the analysis chain for that dataset. Changing such values (even for positions) could generate duplicate header keywords, causing the data file to fail processing in later steps.

Preparation of data for analysis

Making cuts using the events file keywords

Applying selection cuts to the LAT data is typically performed using gtselect. Cuts are made to bound the dataset's time range, energy range, position, region of interest radius, and maximum zenith angle. The maximum zenith angle selection is designed to exclude time periods when any portion of the region of interest is too close to the Earth's limb, resulting in elevated background levels. The Earth's limb lies at a zenith angle of 113 degrees, so a suggested value of 105o provides protection against significant contamination by atmospheric gammas.

Note: This cut requires careful handling when calculating livetime. (See Likelihood Livetime and Exposure.)

Additionally, gtselect can be used to include/exclude events that have a lower probability of being photons (event classes 1 or 2). In general, it is recommended that only diffuse class (event class = 3) photons be used for LAT analysis. However, analysis of short duration transients such as prompt emission from gamma-ray bursts is not as sensitive to the presence of moderate background  contamination, and so do not suffer from the presence of lower probability events. The instrument response functions used in the analysis tools must match the event classes selected.

The reconstruction and classification of LAT events is performed by the LAT instrument team as part of the Level 1 science data processing. The version of the event reconstruction algorithms used is referred to as passN. Currently, the LAT data pipeline uses Pass 6 reconstruction algorithms. This processing is inherent in the science data, and cannot be changed. For a given reconstruction, the LAT team generates an accompanying set of parameterized instrument response functions (IRFs) that are designed for analysis of that particular dataset. These are labeled VN (where N is the version of the IRF). There may be several updates of IRFs for a given reconstruction, with one recommended for analysis (currently P6_V3). The IRFs selected later in the analysis chain must match the data reconstruction version (currently Pass 6) and the event class selected at this stage (i.e., P6_V3_DIFFUSE) for a point source analysis. (Also see Event Classes and Instrument Response Function.)

The acceptance cone radius will vary with analysis type and source location. Performing a spectral analysis on a point source in the galactic plane will require a larger initial analysis region than for a source off the plane, to allow for fitting of multiple nearby sources. An acceptance cone radius of 10 deg is appropriate for spectral analysis of point sources off the Galactic plane, while 10-20 deg may be necessary for point sources located near the Galactic Plane. A timing analysis that does not use spectral modeling (for example, pulsar lightcurves) might select a smaller region to reduce background levels. These studies (timing analysis, GRB spectral analysis, extended source analysis, and others) require optimization of the acceptance cone on a case-by-case basis. Such optimization will need to consider the largest PSF of the LAT in the energy range of interest, as well as providing sufficient statistics to account for the surrounding sky structure. An all-sky analysis does not require this selection.

Making cuts using the spacecraft file keywords

The spacecraft files are used to select the good time intervals (GTIs) for data analysis by using gtmktime. Selecting times when the data quality is good (DATA_QUAL==1) will exclude time periods when some spacecraft event has affected the quality of the data. gtmktime can be used to generate and combine GTIs for any logical expression of keywords in the spacecraft file. Contamination from albedo gamma rays form the Earth is reduced by applying an additional timing selection to the maximum zenith angle cut made with gtselect. Excluding time intervals when the Earth's limb intersects the selected region of interest is one method of correcting the exposure for the zenith cut. This is accomplished by using the ROIcut=yes option in gtmktime. For other exposure correction options, see the Likelihood Livetime and Exposure section.

GBM data does not typically require selection cuts prior to use. However, cuts can be applied should they be needed for a tailored analysis chain.


Last updated by: Chuck Patterson 12/06/2010