# README_3PC_SEDs.txt
# README version 1, created by david.smith@u-bordeaux.fr on 2023 March 5.
# Update 2023 May 25 (das for mtk)
# Update 2023 July 25 (das)
# Update 2023 August 3 (das) -- 11 of the more recent gamma-ray pulsars
  did not have flux points or limits in the SED plots. Added the 4FGL name to hdu[1] of the .fits file, regenerated those plots, and updated the .tgz and 
  the concatenated .pdf file. The pulsars are
  J1259-8148 J1306-6043  J1402+1306 J1623-6936 J1757-6032 J1803-6707 
  J1824-0621 J1835-3259B J1852-1310 J2029-4239 J2116+1345

This README describes the collection of figures and data showing the
Spectral Energy Distributions (SEDs) of the pulsars provided with the
Fermi LAT 3rd catalog of gamma-ray pulsars ("3PC").

SPECTRAL ENERGY DISTRIBUTIONS:
There is one SED for each gamma-ray pulsar for which there is a TS > 25 point source in 4FGL-DR3.
It derives from the 4FGL-DR3 catalog, https://fermi.gsfc.nasa.gov/ssc/data/access/lat/12yr_catalog/, 
but with the following differences:
i) Only 295 rows are kept (Pulsars, and sources associated with pulsars).
ii) The original PLEC columns have been replaced by PLEC columns with b=2/3 fixed (all pulsars) ; 
and with b left free in the likelihood fit (119 pulsars). See 3PC Section 6 for details.

The SED plots are currently provided in .pdf, .png, and .txt formats.
Plots for each pulsar are provided in the gzipped tar file: 3PC_SEDs_YYYYMMDD.tgz
This tar file unpacks into the subdirectory ./3PC_SEDs/
In addition, the file 3PC_AllSEDs_YYYMMDD.pdf is a single PDF with all the SEDs, one per page.

The data for generating the plots are in the "LAT_Point_Source_Catalog" extension (HDU 3) of 
the main catalog FITS file, 3PC_Catalog+SEDs_YYMMDD.fits.  In the "LAT_Point_Source_Catalog" HDU 
statistical uncertainties for values XXX are designated Unc_XXX.

We provide the script that generated the plots, which reads the
3PC_Catalog+SEDs_YYYYMMDD.fits version of the catalog.

This is run like this:
% python 3PC_SEDPlotter.py J0030+0451
By default, this will generate PDF, PNG and .txt files of the SED for the specified pulsar,
in the directory ./3PC_SEDs/

You can put several pulsar names. If you put NO pulsar name, you generate the output for _all_ gamma-ray pulsars.
The .txt file lists the data points and the PLEC4 model curve.
PLEC is an abbreviation for Power Law with Exponential Cutoff.

REQUIREMENTS -- 3PC_SEDPlotter.py uses these libraries: numpy, astropy, matplotlib