Homec4science
Diffusion ARRAKIHS (master) Browse

arrakihs/Workspace/CumNumOfSat/satComputeMstarvsMhalo.py
0d50a2737ed2master

Last Change
Raw File

satComputeMstarvsMhalo.py
View Options

#!/usr/bin/env python3
import sys,os,string
import argparse
import satlib
import numpy as np
import Ptools as pt
import pickle
from scipy.interpolate import splrep,splev
####################################################################
# option parser
####################################################################
description=""
epilog =""""""
parser = argparse.ArgumentParser(description=description,epilog=epilog,formatter_class=argparse.RawDescriptionHelpFormatter)
parser.add_argument("--LvvsMh_model",
action="store",
type=str,
dest="LvvsMh_model",
default = None,
help="Lv Mh model")
parser.add_argument("-p",
action="store",
dest="ps",
metavar='FILE NAME',
type=str,
default=None,
help='output file name')
parser.add_argument("--dpi",
action="store",
dest="dpi",
type=float,
default=300,
help="DPI of the saved file",
metavar=" FLOAT")
####################################################################
# main
####################################################################
def MakePlot(opt):
# get Mh vs Lv relation
Mhbins,Lvmins,Lvmaxs = satlib.GetStellarMassHaloMassRelation(model=opt.LvvsMh_model)
ax = pt.gca()
ax.fill_between(Mhbins,Lvmins,Lvmaxs,facecolor='orange',alpha=0.5,label=r"$\textrm{model=%s}$"%opt.LvvsMh_model)
#pt.plot(Mhbins,Lvmins)
#pt.plot(Mhbins,Lvmaxs)
###########################
# finalize
###########################
xmin = 1e8
xmax = 1e11
ymin = 1e2
ymax = 1e10
xlabel = r"$\rm{Halo\,\,Mass\,\,(M_{\rm halo})}\,\,[M_{\odot}]$"
ylabel = r"$\rm{Stellar\,\,Mass\,\,(M_{\star})}\,\,[M_{\odot}]$"
pt.SetAxis(xmin,xmax,ymin,ymax,log="xy")
pt.xlabel(xlabel,fontsize=pt.labelfont)
pt.ylabel(ylabel,fontsize=pt.labelfont)
pt.grid(False)
pt.legend()
pt.title(r"$\textrm{model=%s}$"%opt.LvvsMh_model,fontsize=pt.labelfont)
if __name__ == '__main__':
opt = parser.parse_args()
files = []
pt.InitPlot(files, opt)
pt.labelfont=20
# pt.figure(figsize=(8*2,6*2))
# pt.figure(dpi=10)
#fig = pt.gcf()
# fig.subplots_adjust(left=0.1)
# fig.subplots_adjust(right=1)
# fig.subplots_adjust(bottom=0.12)
# fig.subplots_adjust(top=0.95)
# fig.subplots_adjust(wspace=0.25)
# fig.subplots_adjust(hspace=0.02)
MakePlot(opt)
pt.EndPlot(files, opt)
c4science · Help