Skip to content

Commit

Permalink
modify input kwargs to allow for arguments to be passed to Lnu_of_nu …
Browse files Browse the repository at this point in the history
…in the MQ24 code
  • Loading branch information
noahfranz13 committed Aug 2, 2024
1 parent c0bfa7d commit c8e2544
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 7 deletions.
2 changes: 1 addition & 1 deletion src/syncfit/_version.py
Original file line number Diff line number Diff line change
@@ -1 +1 @@
__version__ = "0.3.3"
__version__ = "0.3.4"
8 changes: 2 additions & 6 deletions src/syncfit/models/mq_model.py
Original file line number Diff line number Diff line change
Expand Up @@ -37,14 +37,10 @@ def __init__(self, prior=None, p=None, t=None):
self.t = t

def SED(self, nu, p, log_bG_sh, log_Mdot, log_epsilon_T, log_epsilon_e, log_epsilon_B,
t, lum_dist, **kwargs):
t, lum_dist, ell_dec = 1.0, f = 3.0/16.0, **kwargs):

# set microphysical and geometric parameters
# log_epsilon_e = -1
# log_epsilon_B = log_epsilon_e # assume equipartition
delta = 10**log_epsilon_e/10**log_epsilon_T
f = 3.0/16.0
ell_dec = 1.0

t = (t*u.day).to(u.s).value

Expand All @@ -53,7 +49,7 @@ def SED(self, nu, p, log_bG_sh, log_Mdot, log_epsilon_T, log_epsilon_e, log_epsi
Lnu = Lnu_of_nu(
10**log_bG_sh, Mdot_over_vw, nu, t, p=p,
epsilon_T=10**log_epsilon_T, epsilon_B=10**log_epsilon_B, epsilon_e=10**log_epsilon_e,
f=f,ell_dec=ell_dec,radius_insteadof_time=False
f=f,ell_dec=ell_dec, **kwargs
) * u.erg / (u.s * u.Hz)

lum_dist_cm = lum_dist*u.cm # give it units so the conversion works well
Expand Down

0 comments on commit c8e2544

Please sign in to comment.