#ifndef PM_SUBTRACTION_IO_H
#define PM_SUBTRACTION_IO_H

#include <pslib.h>

#include <pmHDU.h>
#include <pmFPA.h>

/// Write subtraction kernels within a readout to a FITS file
bool pmReadoutWriteSubtractionKernels(
    pmReadout *readout,                 ///< Readout for which to write subtraction kernels (in analysis MD)
    psFits *fits                        ///< FITS file to which to write
    );


bool pmReadoutReadSubtractionKernels(
    pmReadout *readout,                 ///< Readout for which to read subtraction kernels (into analysis MD)
    psFits *fits                        ///< FITS file to which to write
    );


#endif
