Emlineqa

descqa/EmlineRatioTest.py

@@ -55,6 +55,9 @@ class EmlineRatioTest(BaseValidationTest):
     sim_drawnum: int, optional, (default: 30000)
         The number of galaxies to draw from the simulated data to perform the comparison.
         The default number is chosen to (hopefully) not make the 2-D KS test too stringent.
+    truncate_cat_name: Bool, optional, (default: False)
+        Specifies whether the catalog name displayed in the summary figure should be 
+        shortened.
     """
     def __init__(self, **kwargs):
 
@@ -88,6 +91,8 @@ def __init__(self, **kwargs):
         self.figlist = []
         self.runcat_name = []
 
+        self.truncate_cat_name = kwargs.get('truncate_cat_name', False)
+
 
     def run_on_single_catalog(self, catalog_instance, catalog_name, output_dir):
 
@@ -142,15 +147,15 @@ def run_on_single_catalog(self, catalog_instance, catalog_name, output_dir):
                                                 'mag_y_lsst'], filters=(uband_maglim | gband_maglim | rband_maglim | iband_maglim | zband_maglim | yband_maglim))
         sz = data['redshift']
         galaxyID = data['galaxyID']
-        Halpha = data['emissionLines/totalLineLuminosity:balmerAlpha6563']* 4.4659e13*u.W/u.Hz
-        Hbeta = data['emissionLines/totalLineLuminosity:balmerBeta4861']* 4.4659e13*u.W/u.Hz
-        NII6584 = data['emissionLines/totalLineLuminosity:nitrogenII6584']* 4.4659e13*u.W/u.Hz
-        OIII5007 = data['emissionLines/totalLineLuminosity:oxygenIII5007']* 4.4659e13*u.W/u.Hz
-        OIII4959 = data['emissionLines/totalLineLuminosity:oxygenIII4959']* 4.4659e13*u.W/u.Hz
-        OII3726 = data['emissionLines/totalLineLuminosity:oxygenII3726']* 4.4659e13*u.W/u.Hz
-        OII3729 = data['emissionLines/totalLineLuminosity:oxygenII3729']* 4.4659e13*u.W/u.Hz
-        SII6716 = data['emissionLines/totalLineLuminosity:sulfurII6716']* 4.4659e13*u.W/u.Hz
-        SII6731 = data['emissionLines/totalLineLuminosity:sulfurII6731']* 4.4659e13*u.W/u.Hz
+        Halpha = data['emissionLines/totalLineLuminosity:balmerAlpha6563'] * 3.839e26*u.W
+        Hbeta = data['emissionLines/totalLineLuminosity:balmerBeta4861'] * 3.839e26*u.W
+        NII6584 = data['emissionLines/totalLineLuminosity:nitrogenII6584'] * 3.839e26*u.W
+        OIII5007 = data['emissionLines/totalLineLuminosity:oxygenIII5007'] * 3.839e26*u.W
+        OIII4959 = data['emissionLines/totalLineLuminosity:oxygenIII4959'] * 3.839e26*u.W
+        OII3726 = data['emissionLines/totalLineLuminosity:oxygenII3726'] * 3.839e26*u.W
+        OII3729 = data['emissionLines/totalLineLuminosity:oxygenII3729'] * 3.839e26*u.W
+        SII6716 = data['emissionLines/totalLineLuminosity:sulfurII6716'] * 3.839e26*u.W
+        SII6731 = data['emissionLines/totalLineLuminosity:sulfurII6731'] * 3.839e26*u.W
         SIItot = SII6716 + SII6731
         OIIItot = OIII5007 + OIII4959
         OIItot = OII3726 + OII3729
@@ -172,7 +177,7 @@ def run_on_single_catalog(self, catalog_instance, catalog_name, output_dir):
 
         for x in range(len(property_list)):
 
-            property_list[x] = (property_list[x][indices]/(4*np.pi*lumdist_small**2)).to('erg/s/cm**2/Hz').value
+            property_list[x] = (property_list[x][indices]/(4*np.pi*lumdist_small**2)).to('erg/s/cm**2').value
 
         Halpha_small, Hbeta_small, NII6584_small, OIII5007_small, OIII4959_small, OII3726_small, OII3729_small, SII6716_small, SII6731_small, SIItot_small, OIIItot_small, OIItot_small = property_list
 
@@ -199,8 +204,11 @@ def run_on_single_catalog(self, catalog_instance, catalog_name, output_dir):
         # Perform the Test and Return Results
         #=========================================
 
+        if self.truncate_cat_name:
+            thisfig, pvalue, medianshift = self.makeplot(catalog_name.split('_')[0])
+        else:
+            thisfig, pvalue, medianshift = self.makeplot(catalog_name)
 
-        thisfig, pvalue, medianshift = self.makeplot(catalog_name)
         self.figlist.append(thisfig)
         self.runcat_name.append(catalog_name)
 

descqa/configs/EmlineLuminosityRatios.yaml

@@ -10,4 +10,5 @@ mag_r_cut: 27.7
 mag_i_cut: 27.0
 mag_z_cut: 26.2
 mag_y_cut: 24.9
+truncate_cat_name: False
 description: 'Plots ratios of emission line luminosities from SDSS data and DESC simulations and compares the two 2-D distributions.'