Some analysis notebooks for IXPE data using the maximum likelihood techique outlined in https://arxiv.org/abs/2204.00140.
IXPE Reduction:
- Reduce-CenX3.ipynb : reduce the data for Cen X-3 including source barycentring, period finding and folding
- Reduce-Generic-IXPE.ipynb : reduce the data for a generic IXPE pulsar observation including Solar-System and source barycentring, period from GBM and folding. Solar-system barycenter correction requires the HEASOFT barycen program. This example is Her X-1 during the main-on.
Julia MCMC:
- MCMC-Julia.jl : Julia script to perform MCMC estimation of confidence regions on RVM parameters (uses photondata.txt produced by the reduce notebooks): photondata contains phase, 2cos(2PA)modf, 2sin(2PA)*modf for each photon (IXPE convention)
- PBS-Julia-Script : Run the Julia script on a PBS server with 32 CPUs
Precession:
- Precession-HerX1-Final.ipynb : Python notebook to derive formulae for a precessing triaxial body and applying them to the polarization angle measurements of Her X-1