Open matthew-w-lundy opened 1 year ago
Currently this is what you have:
SIMPLE = T / conforms to FITS standard
BITPIX = 8 / array data type
NAXIS = 0 / number of array dimensions
EXTEND = T
TELESCOP= 'VERITAS ' / Telescope
LICENSE = ''
COMMENT FITS (Flexible Image Transport System) for COMMENT mat is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; b COMMENT ibcode: 2001A&A...376..359H
END
XTENSION= 'BINTABLE' / binary table extension
BITPIX = 8 / array data type
NAXIS = 2 / number of array dimensions
NAXIS1 = 52 / length of dimension 1
NAXIS2 = 4979 / length of dimension 2
PCOUNT = 0 / number of group parameters
GCOUNT = 1 / number of groups
TFIELDS = 9 / number of table fields
TTYPE1 = 'EVENT_ID'
TFORM1 = '1K '
TTYPE2 = 'TIME '
TFORM2 = '1D '
TUNIT2 = 's '
TTYPE3 = 'RA '
TFORM3 = '1E '
TUNIT3 = 'deg '
TTYPE4 = 'DEC '
TFORM4 = '1E '
TUNIT4 = 'deg '
TTYPE5 = 'ENERGY '
TFORM5 = '1E '
TUNIT5 = 'TeV '
TTYPE6 = 'ALT '
TFORM6 = '1E '
TUNIT6 = 'deg '
TTYPE7 = 'AZ '
TFORM7 = '1E '
TUNIT7 = 'deg '
TTYPE8 = 'MSW '
TFORM8 = '1D '
TTYPE9 = 'MSL '
TFORM9 = '1D '
EXTNAME = 'EVENTS ' / extension name
HDUCLASS= 'GADF ' / This FITS file follows the GADF data format
HDUDOC = 'https://gamma-astro-data-formats.readthedocs.io'
HDUVERS = '0.3 ' / DL3 specification version
HDUCLAS1= 'EVENTS ' / Primary extension class
RADECSYS= 'FK5 ' / equatorial system type
EQUINOX = 2000.0 / base equinox
CREATOR = 'pyV2DL3 v0.5::VEGAS'
ORIGIN = 'VERITAS Collaboration' / Data from VERITAS
TELESCOP= 'VERITAS '
INSTRUME= 'VERITAS '
OBS_ID = 103110 / Run Number
DATE-OBS= '2022-09-25T10:49:18.000' / start date (UTC) of obs yy-mm-dd hh:mm:ss DATE-AVG= '2022-09-25T11:04:18.000' / average date (UTC) of obs
DATE-END= '2022-09-25T11:19:18.000' / end date (UTC) of obs yy-mm-dd hh:mm:ss
TSTART = 556886963.0 / mission time of start of obs [s]
TSTOP = 556888763.0 / mission time of end of obs [s]
MJDREFI = 53402.0 / int part of reference MJD [days]
MJDREFF = 0.0 / fractional part of reference MJD [days]
TIMEUNIT= 's ' / time unit is seconds since MET start
TIMESYS = 'utc ' / time scale is UTC
TIMEREF = 'topocenter' / location from where the observation was made
ONTIME = 1800.0 / sum of good time intervals [s]
LIVETIME= 1571.7390439 / (ontime * deadtime time correction) [s]
DEADC = 0.8731883577222223 / Average dead time correction (LIVETIME/ONTIME) OBJECT = 'Crab ' / observed object
RA_OBJ = 83.63333046211 / right ascension of object [deg]
DEC_OBJ = 22.01447150632731 / declination of object [deg]
RA_PNT = 83.63488735564972 / pointing right ascension [deg]
DEC_PNT = 21.51889080394067 / pointing declination [deg]
ALT_PNT = 65.85945928266763 / average pointing altitude [deg]
AZ_PNT = 108.949695184459 / average pointing azimuth [deg]
TELLIST = 'T1,T2,T3,T4' / comma-separated list of tel IDs
N_TELS = 4 / number of telescopes in event list
EUNIT = 'TeV ' / energy unit
GEOLON = -110.952408 / longitude of array center [deg]
GEOLAT = 31.674951 / latitude of array center [deg]
ALTITUDE= 1277.110044 / altitude of array center [m]
END
Would be nice if the header of the fits file produced contained the following information.