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'])