scicoder_project/segue.py at master · sukhdeep2/scicoder_project · 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
#!/anaconda/bin/python

import sys
import astropy
from astropy.io import fits

FileName=sys.argv[1]

Data = fits.open(FileName,memmap=True)

print Data[1].header[20]
tbdata = Data[1].data
feh = tbdata.field("FEH_ADOP")
dist = tbdata.field("DIST_ADOP")
v = tbdata.field("RV_ADOP")
spec = tbdata.field("SPECTYPE_HAMMER")
Data.close()

inx = ((feh > -9999) & (dist > -9999) & (v > -9999))
feh = feh[inx]
dist = dist[inx]
v = v[inx]
spec = spec[x]
unique_spectpyes = set(spec)


import matplotlib.pyplot as plt

fig,ax = plt.subplots(3,1,figsize=(10,10))
ax[0].hist(feh,40)
ax[0].set_xlim(-4,1)
ax[0].set_xlabel('[FE/H]',fontsize=20)
ax[1].hist(dist,800)
ax[1].set_xlim(0,50)
ax[1].set_xlabel('Distance (kpc =??)',fontsize=20)
ax[2].hist(v,50)
ax[2].set_xlim(-500,500)
ax[2].set_xlabel('V (km/s)',fontsize=20)
fig.subplots_adjust(hspace=0.3)


import matplotlib.pyplot as plt

fig,ax = plt.subplots(3,1,figsize=(10,10))
ax[0].hist(feh,40)
ax[0].set_xlim(-4,1)
ax[0].set_xlabel('[FE/H]',fontsize=20)
ax[1].hist(dist,800)
ax[1].set_xlim(0,50)
ax[1].set_xlabel('Distance (kpc =??)',fontsize=20)
ax[2].hist(v,50)
ax[2].set_xlim(-500,500)
ax[2].set_xlabel('V (km/s)',fontsize=20)
fig.subplots_adjust(hspace=0.3)