def plot(self, mass=None, sfr=None, z=None, FQdist=None, sfms_prop=None, **mkwargs): ''' Plot the quiescent fraction as a function of stellar mass. The quiescent fraction is calculated based on an evolving sSFR(M*,z) cut. Parameters ---------- FQdist : list List that specifies the mass bins and the quiescent fraction (eq. [massbin, quiescentfraction]) ''' if FQdist is None: if sfr is None or z is None: raise ValueError fq_obj = Fq() masses, fq = fq_obj.Calculate(mass=mass, sfr=sfr, z=z, sfms_prop=sfms_prop) else: masses, fq = FQdist self.masses = masses self.z = z if self.kwargs == {}: kwargs = mkwargs.copy() else: kwargs = (self.kwargs).copy() kwargs.update(mkwargs) if 'label' in kwargs: fq_label = kwargs['label'] else: fq_label = None if 'line_color' in kwargs: line_color = kwargs['line_color'] else: line_color = 'black' if 'line_style' in kwargs: line_style = kwargs['line_style'] else: line_style = '-' if 'lw' in kwargs: line_width = kwargs['lw'] else: line_width = 4 self.subs.plot(masses, fq, color = line_color, lw = line_width, ls = line_style, label = fq_label) return None
def Fq(self, **fq_kwargs): ''' Calculate quiescent fraction of CenQue class object ''' if self.zsnap is None: raise ValueError if self.mass is None: raise ValueError if self.sfr is None: raise ValueError if self.sfr_prop is None: raise ValueError qfrac = Fq() # Star-forming or Quiescent return qfrac.Calculate(mass=self.mass, sfr=self.sfr, z=self.zsnap, sfms_prop=self.sfr_prop['sfms'])