def Pf_func_gamma(x, samp, d):
"""The equation to calculate :math:`P_x(f)` from the cantilever
parameters and detector noise floor.
f
Frequency, the independent variable.
P_x0
The zero frequency power spectral density of position fluctuations.
f_c
The cantilever resonance frequency.
Q
The quality factor of the cantilever.
P_detector
The detector noise floor.
"""
fc = d['mu_fc'] + samp['dfc']
Q = samp['Q']
kc = samp['kc']
return Pf(x,
calc_P_x0(fc*u.Hz, Q, kc * u('N/m'), d['T']*u.K).to('nm^2/Hz').magnitude,
fc, Q, samp['Pdet'] * d['scale']
)
def Pfi(d, samp, i):
fc = d['mu_fc'] + samp['dfc'][i]
Q = samp['Q'][i]
kc = samp['kc'][i]
return Pf(d['f'],
calc_P_x0(fc*u.Hz, Q, kc*u('N/m'), d['T']*u.K).to('nm^2/Hz').magnitude,
fc, Q, samp['Pdet'][i]*d['scale']
)