fit_ramps_casertano_no_dq

romanisim.ramp.fit_ramps_casertano_no_dq(resultants, read_noise, read_pattern)

Fit ramps following Casertano+2022, only using full ramps.

This is a simpler implementation of fit_ramps_casertano, which doesn’t address the case of partial ramps broken by CRs. This case is easier and can be done reasonably efficiently in pure python; results can be compared with fit_ramps_casertano in for the case of unbroken ramps.

Parameters:

resultantsnp.ndarry[nresultants, npixel]

the resultants in electrons

read noise: float

the read noise in electrons

read_patternlist[list[int]]

list of lists giving indices of reads entering each resultant

Returns:

parnp.ndarray[nx, ny, 2] (float)

the best fit pedestal and slope for each pixel

varnp.ndarray[nx, ny, 3, 2, 2] (float)

the covariance matrix of par, for each of three noise terms: the read noise, Poisson source noise, and total noise.