Skip to content

Commit

Permalink
pixel noise updates
Browse files Browse the repository at this point in the history
turn off flat field errors by default
fix bug to specify ngroup as an image array
  • Loading branch information
JarronL committed Sep 27, 2019
1 parent 223dae5 commit 9d961ab
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 2 deletions.
3 changes: 2 additions & 1 deletion pynrc/nrc_utils.py
Original file line number Diff line number Diff line change
Expand Up @@ -2801,7 +2801,7 @@ def sat_limit_webbpsf(filter_or_bp, pupil=None, mask=None, module='A', pix_scale

def pix_noise(ngroup=2, nf=1, nd2=0, tf=10.737, rn=15.0, ktc=29.0, p_excess=(0,0),
fsrc=0.0, idark=0.003, fzodi=0, fbg=0, ideal_Poisson=False,
ff_noise = True, **kwargs):
ff_noise=False, **kwargs):
"""Noise per pixel
Theoretical noise calculation of a generalized MULTIACCUM ramp in terms of e-/sec.
Expand Down Expand Up @@ -2839,6 +2839,7 @@ def pix_noise(ngroup=2, nf=1, nd2=0, tf=10.737, rn=15.0, ktc=29.0, p_excess=(0,0
otherwise MULTIACCUM equation is used?
ff_noise : bool
Include flat field errors in calculation? From JWST-CALC-003894.
Default=False.
Notes
-----
Expand Down
5 changes: 4 additions & 1 deletion pynrc/pynrc_core.py
Original file line number Diff line number Diff line change
Expand Up @@ -225,8 +225,11 @@ def pixel_noise(self, fsrc=0.0, fzodi=0.0, fbg=0.0, ng=None, verbose=False, **kw
"""

ma = self.multiaccum
if ng is none:
ng = ma.ngroup

# Pixel noise per ramp (e-/sec/pix)
pn = pix_noise(ma.ngroup, ma.nf, ma.nd2, tf=self.time_frame, \
pn = pix_noise(ng, ma.nf, ma.nd2, tf=self.time_frame, \
rn=self.read_noise, ktc=self.ktc, p_excess=self.p_excess, \
idark=self.dark_current, fsrc=fsrc, fzodi=fzodi, fbg=fbg, **kwargs)

Expand Down

0 comments on commit 9d961ab

Please sign in to comment.