def rp(self, lam, unit='kev'):
assert unit in c.ALLOWED_LAM_UNITS
lam_cm = c._lam_cm(lam, unit)
mm1 = self.rho / (2.0 * c.m_p) * c.r_e / (2.0 * np.pi) * np.power(
lam_cm, 2)
return mm1 + 1.0
'''# Returns 1 if the wavelength supplied is too low energy (i.e. inappropriate for applying Drude)
def ip(self, lam, unit='kev'):
lam_cm = c._lam_cm(lam, unit)
return self._interp_helper(lam_cm, self.interps[1])
def rp(self, lam, unit='kev'):
lam_cm = c._lam_cm(lam, unit)
return self._interp_helper(lam_cm, self.interps[0], rp=True)