My code for the CTI/deferred charge correction is available at: https://github.com/snyder18/ctisim
The script I use for correcting the SLAC data is: https://github.com/snyder18/ctisim/blob/master/scripts/correct_images.py
I do bias/overscan/gain calibration but no assembly because I am just saving the results to another FITs file with the same multiple HDUs as the original image file.
The parameterization per amplifier are two variables for the electronics component (scale and decay time), a global CTI parameter, and the trap definition. For traps I have been using a look-up table with values taken from the flat field overscan results, but one could also define the traps via a functional form and parameters.