Update calc_sfr.py

Updated IMF correction. Values don't change, but the factor is consistent with what's quoted in Madau & Dickinson (0.66 = 1./1.5)
sagasurvey · Dec 21, 2023 · 41f85fb · 41f85fb
1 parent e5feb3b
commit 41f85fb
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/SAGA/objects/calc_sfr.py b/SAGA/objects/calc_sfr.py
@@ -33,8 +33,8 @@ def calc_SFR_NUV(NUV_mag, NUV_mag_err, dist_mpc, internal_ext=0.9):
     l_nuv_msun = log_L_nuv - np.log10(3.826e33)
 
     # CONVVERT TO SFR:   EQ 3, inglesias- paramo 2006
-    # AND ACCOUNT FOR IMF
-    log_SFR_NUV = l_nuv_msun - 9.33 - np.log10(1.5)
+    # Salpeter -> Koupa IMF using Madua & Dickinson 2014 (Figure 4)
+    log_SFR_NUV = l_nuv_msun - 9.33 + np.log10(0.66)
     log_SFR_NUV_err = log_flux_nuv_err
 
     return log_SFR_NUV, log_SFR_NUV_err
@@ -60,9 +60,9 @@ def calc_SFR_Halpha(EW_Halpha, EW_Halpha_err, spec_z, Mr, EWc=2.5, BD=3.25):
     L_Halpha_err = term1_err * term2 * term3
 
     # EQ 3, Bauer et al above
-    # Account for IMF
-    SFR = L_Halpha / (1.27e34 * 1.5)
-    SFR_err = L_Halpha_err / (1.27e34 * 1.5)
+    # Account for IMF, Madua & Dickinson 2014 (Figure $=4)
+    SFR = (L_Halpha * 0.66) / 1.27e34 
+    SFR_err = (L_Halpha_err * 0.66)/ 1.27e34 
 
     log_Ha_SFR = np.log10(SFR)