Skip to content

Commit

Permalink
add latex names for params, small changes to model
Browse files Browse the repository at this point in the history
  • Loading branch information
calumgordon committed Jan 17, 2025
1 parent 761b64e commit 706e4a4
Show file tree
Hide file tree
Showing 7 changed files with 2,162 additions and 4 deletions.
10 changes: 6 additions & 4 deletions vega/correlation_func.py
Original file line number Diff line number Diff line change
Expand Up @@ -374,7 +374,8 @@ def _prep_delta_gamma_model(self):

def _compute_gamma_extension(self, g_rf_wave, params):

L1,L2 = 1205.1, 1213
#L1,L2 = 1205.1, 1213
L1,L2 = 1207.6, 1213
p1, p2 = params['Ca1'], params['Ca2']

#this may be a lot slower
Expand All @@ -383,7 +384,7 @@ def _compute_gamma_extension(self, g_rf_wave, params):

if max(g_rf_wave)<L1:
mean_gamma += p1 / (g_rf_wave - L1)**2
mean_gamma += p2 / (g_rf_wave - L2)**2
#mean_gamma += p2 / (g_rf_wave - L2)**2

return mean_gamma

Expand All @@ -398,7 +399,8 @@ def _compute_gamma_extension(self, g_rf_wave, params):

def _compute_gamma_extension_QSO(self, g_rf_wave, params):

L1,L2 = 1205.1, 1213
#L1,L2 = 1205.1, 1213
L1,L2 = 1207.6, 1213
p1, p2 = params['Ca1'], params['Ca2']

#this may be a lot slower
Expand All @@ -407,7 +409,7 @@ def _compute_gamma_extension_QSO(self, g_rf_wave, params):

if max(g_rf_wave)<L1:
mean_gamma += p1 / (g_rf_wave - L1)**2
mean_gamma += p2 / (g_rf_wave - L2)**2
#mean_gamma += p2 / (g_rf_wave - L2)**2

return mean_gamma

Expand Down
Loading

0 comments on commit 706e4a4

Please sign in to comment.