In [2]:

    

import matplotlib.pyplot as plt
from pylab import *

params = {'backend': 'pdf',
            'figure.figsize': [4.3, 3.0],
            'font.family':'serif',
            'font.size':10,
            'font.serif': 'Times Roman',
            'axes.titlesize': 'medium',
            'axes.labelsize': 'medium',
            'legend.fontsize': 8,
            'legend.frameon' : False,
            'text.usetex': True,
            'figure.dpi': 600,
            'lines.markersize': 2,
            'lines.linewidth': 2,
            'lines.antialiased': False,
            'path.simplify': False,
            'legend.handlelength':3,
            'figure.subplot.bottom':0.15,
            'figure.subplot.top':0.9,
            'figure.subplot.left':0.15,
            'figure.subplot.right':0.92}

mpl.rcParams.update(params)


    

In [ ]:

    




    

In [70]:

    

plt.hist(peris_c,bins=np.linspace(0,100,100), alpha = 0.6, label="$3.23M_\\odot$, one pulse")
plt.hist(peris_a,bins=np.linspace(0,100,100), alpha = 0.6, label="Same but in 7 equal pulses",lw=0.5, ec="k", ls="--", histtype='stepfilled')
plt.hist(peris_b,bins=np.linspace(0,100,100), alpha = 0.6, label="7 pulses from MESA",lw=0.5, ec="k", histtype='stepfilled')
plt.hist(peris_d,bins=np.linspace(0,100,100), alpha = 0.6, label="Same+10 km/s kick",lw=0.5, ec="k", ls=":", histtype='stepfilled')
plt.gca().set_xlabel("distance at periastron")
plt.gca().set_ylabel("number")
plt.legend(loc=1, fontsize=6)
plt.title("$39M_\odot+39M_\odot, a=50R_\odot$ at core collapse")
plt.yscale('log', nonposy='clip')
plt.savefig("post_kick_peri.pdf")
plt.show()


    

In [ ]: