bug fixed

pcigale/sed_modules/fritz2006.py

@@ -142,9 +142,10 @@ class Fritz2006(SedModule):
         lumin_agn_new = self.fritz2006.lumin_agn * ext_fac
         lumin_scatt_new = self.fritz2006.lumin_scatt * ext_fac
         # Calculate the total extincted luminosity averaged over all directions
+        # note that self.opening_angle is different from Fritz's definiting!
         l_ext = np.trapz(self.fritz2006.lumin_intrin_agn * (1-ext_fac),
                          x=self.fritz2006.wave) * \
-                        (1 - np.sin( np.deg2rad(self.opening_angle/2.) ))
+                        (1 - np.cos( np.deg2rad(self.opening_angle) ))
         # Casey (2012) modified black body model
         conv = self.c / (self.fritz2006.wave * self.fritz2006.wave)
         lumin_blackbody = (conv * (1. - np.exp(-(lambda_0 / self.fritz2006.wave)

pcigale/sed_modules/xray.py

@@ -76,7 +76,7 @@ class Xray(SedModule):
         # also have two components:
         #   2.1 high-mass X-ray binaries (HMXB)
         #   2.2 low-mass X-ray binaries (LMXB)
-        # Assuming gamma=2. for both components (Wu et al. 2008)
+        # Assuming gamma=2. for both components (Wu & Gu 2008)
         # normalized such that L(2-10 keV) = 1.
         self.lumin_xrb = self.wave**-1 / np.log(5.)