MPPNPy/iso_plot.py at master · Bigash94/MPPNPy · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
import matplotlib.pyplot as pl
import pandas as pd

factor =0.2
df = pd.read_csv('table_selected_isotopes_for_a_mass_coord.dat', skiprows=[0,1], usecols=[1,2,3], sep='\&')

# data sets
x = df[' specie '].values
#y1 = df[' M3.z2m2 KADoNiS 0.3 '].values
#y2 = df[' M3.z2m2 KADoNiS 1.0 '].values
y1 = df[' M15.z6m3 KADoNiS 0.3 '].values
y2 = df[' M15.z6m3 KADoNiS 1.0 '].values

A=[]
for i in x:
    a=i[3:]
    A.append(float(a))

A_plot=[]
pf1=[]
pf2=[]
for i in range(len(A)):
    if y1[i] and y2[i] >= 3.:
        if y1[i]-y2[i] >= factor*y1[i]:
            A_plot.append(A[i])
            pf1.append(y1[i])
            pf2.append(y2[i])
        if y2[i]-y1[i] >= factor*y2[i]:
            A_plot.append(A[i])
            pf1.append(y1[i])
            pf2.append(y2[i])

fig = pl.figure()            # Figure object
ax = fig.add_subplot(1,1,1)

pl.xticks(size=25)
pl.yticks(size=25)

xminorlocator = pl.MultipleLocator(1)
xmajorlocator = pl.MultipleLocator(10)
ax.xaxis.set_major_locator(xmajorlocator)
ax.xaxis.set_minor_locator(xminorlocator)
yminorlocator = pl.MultipleLocator(1)
ymajorlocator = pl.MultipleLocator(0.5)
ax.yaxis.set_major_locator(ymajorlocator)
ax.yaxis.set_minor_locator(yminorlocator)

#pl.semilogy(A_plot,pf1,'gs',markersize=15.,label='M3.z2m2 KADoNiS 0.3')
#pl.semilogy(A_plot,pf2,'ro',markersize=15.,label='M3.z2m2 KADoNiS 1.0')
pl.plot(A,y2/y1,'ko',markersize=10.,label='kadonis 1.0/0.3 for M3.z2m2')
pl.legend(numpoints=1,loc='best',prop={'size':30})

pl.xlabel('$Mass$ $number$', fontsize=40)
pl.ylabel('$X_{i}/X_{sun}$', fontsize=40)

pl.show()