ctpsfcube

Generate point spread function cube for a stacked analysis.

Synopsis

This tool generates a point spread function cube for a stacked analysis. A point spread function cube is a 4-dimensional cube spanned by Right Ascension or Galactic longitude, Declination or Galactic latitude, energy, and offset angle between true and measured arrival direction of a photon. The energy binning of the cube may be either linear, logarithmic, or custom defined using an input file.

ctpsfcube requires on input the event list or observation definition XML file that has been used in the generation of the counts cube using ctbin.

It is not recommended to use the counts cube for the point spread function cube definition, although this is formally possible by specifying the counts cube as incube parameter. This leads however to a large FITS file on output since the number of bins in the counts cube will be multiplied by the number of offset angle bins (typically 200). Since the point spread function varies only slowly over the field of view of the camera it is recommended to use a rather coarse spatial binning to keep the file size manageable (with a typical value of binsz=1.0).

ctpsfcube generates a point spread function cube FITS file comprising three extensions. The primary extension contains a 4-dimensional image that contains the point spread function values. The values are in units of sr**(-1). The next extension named ENERGIES contains a binary table that defines the energies of the point spread function cube. The last extension named DELTAS contains a binary table that defines the offset angles between true and measured arrival direction of the photon.

General parameters

inobs [file]

Input event list or observation definition XML file.

incube [file]

Input counts cube file to extract point spread function cube definition.

caldb [string]

Calibration database.

irf [string]

Instrument response function.

outcube [file]

Output point spread function cube file.

ebinalg <FILE|LIN|LOG|POW> [string]

Algorithm for defining energy bins. For FILE, the energy bins are defined in a FITS file that is specified by the ebinfile parameter, for LIN LOG and POW there will be enumbins energy bins spaced linearly, logarithmically, or following a power law between emin and emax, respectively. For POW, the parameter ebingamma specifies the slope of the power law.

emin [real]

Lower energy value for first energy bin (in TeV) if LIN or LOG energy binning algorithms are used.

emax [real]

Upper energy value for last energy bin (in TeV) if LIN or LOG energy binning algorithms are used.

enumbins [integer]

Number of energy bins if LIN or LOG energy binning algorithms are used.

ebinfile [file]

Name of the file containing the energy binning definition if ebinalg=FILE. You may use csebins to generate a file with appropriate energy binning.

ebingamma [real]

Exponent of the power law for POW energy binning. An exponent of 1.0 corresponds to a logarithmic energy binning.

(addbounds = no) [boolean]

Add energies to the point spread function cube at the observation energy boundaries?

(usepnt = no) [boolean]

Use CTA pointing direction for cube centre instead of xref/yref parameters?

nxpix [integer]

Number of cube bins in Right Ascension or Galactic longitude.

nypix [integer]

Number of cube bins in Declination or Galactic latitude.

binsz [real]

Cube bin size (in degrees/pixel).

coordsys <CEL|GAL> [string]

Coordinate system (CEL - celestial, GAL - galactic).

proj <AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN> [string]

Projection method.

xref [real]

Right Ascension / Galactic longitude of cube centre (J2000, in degrees).

yref [real]

Declination / Galactic latitude of cube centre (J2000, in degrees).

(amax = 0.3) [real]

Upper bound of angular separation between true and measued photon direction (in degrees).

(anumbins = 200) [integer]

Number of angular separation bins.

Standard parameters

(chatter = 2) [integer]
Verbosity of the executable:

chatter = 0: no information will be logged

chatter = 1: only errors will be logged

chatter = 2: errors and actions will be logged

chatter = 3: report about the task execution

chatter = 4: detailed report about the task execution

(clobber = yes) [boolean]

Specifies whether an existing output point spread function cube file should be overwritten.

(debug = no) [boolean]

Enables debug mode. In debug mode the executable will dump any log file output to the console.

(mode = ql) [string]

Mode of automatic parameters (default is ql, i.e. “query and learn”).

(logfile = ctpsfcube.log) [string]

Name of log file.