diff --git a/Metallicity_Stack_Commons/analysis/attenuation.py b/Metallicity_Stack_Commons/analysis/attenuation.py
index 414fd6e..08fe2fb 100644
--- a/Metallicity_Stack_Commons/analysis/attenuation.py
+++ b/Metallicity_Stack_Commons/analysis/attenuation.py
@@ -38,7 +38,7 @@ def compute_EBV(fitspath, use_revised=False):
     HDELTA = combine_asc[HD+'_Flux_Observed'].data
 
     HgHb = HGAMMA / HBETA
-    HdHb = HGAMMA / HBETA
+    HdHb = HDELTA / HBETA
 
     EBV = -2.5 * np.log10(HgHb/HgHb_CaseB)/(k_HGAMMA - k_HBETA)