/
update_tracker.py
150 lines (119 loc) · 5.65 KB
/
update_tracker.py
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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
import functions
import mysql_query
from numpy import *
import os
import sys
import string
import pyfits
###################
### Description ###
###################
### Update hatsouth.intersigma
########################
### Start of program ###
########################
start_date = "2015-03-28"
end_date = "2015-04-31"
query_entry = "select SPECtype,SPECobject,SPECmjd,SPEChjd,SPECrv,SPECrv_err,SPECtelescope,SPECresolution,SPECteff,SPECteff_err,SPEClogg,SPEClogg_err,SPECfeh,SPECfeh_err,SPECccfheight,SPECexptime,SPECsn"
query_entry = query_entry + " from SPEC where SPECutdate >= \""+start_date+"\" and SPECutdate <=\""+end_date+"\" and SPECobject like \"HATS%\" and SPECtelescope=\"ANU23\""
#query_entry = query_entry + " from SPEC where SPECobject = \"HATS778-008\" "
exposure_info = mysql_query.query_hsmso(query_entry)
#print exposure_info
if len(exposure_info) > 0:
for entry in exposure_info:
output = ""
if entry[7] == 23000:
inst = "ANU23"
else:
inst = "WiFeS"
if entry[0] == "RV":
#output = output + entry[1] + " " #name
output = output + str(2400000 + entry[3]) + " " #hjd
output = output + str(entry[4]) + " " #RV
output = output + str(entry[5]) + " " #RVerr
#output = output + str(entry[6]) + " " #tel
output = output + inst + " " #tel
output = output + str(entry[7]) + " " #res
output = output + "0 0 0 0 0 0 0 0 " #teff terr logg loggerr feh feherr vrot vroterr
output = output + str(entry[14]) + " " #ccfheight
output = output + str(entry[15]) + " " #exptime
output = output + str(entry[16]) + "\n" #S/N
if entry[0] == "ST":
#output = output + entry[1] + " " #name
output = output + str(2400000 + entry[2]) + " " #mjd
output = output + "-999 -999 " #rv rverr
#output = output + str(entry[6]) + " " #tel
output = output + inst + " " #tel
output = output + str(entry[7]) + " " #res
output = output + str(entry[8]) + " 300 " #teff
#output = output + str(entry[9]) + " "
output = output + str(entry[10]) + " 0.3 " #logg
#output = output + str(entry[11]) + " "
output = output + str(entry[12]) + " 0.5 " #feh
output = output + "0 0 0 " #feh feherr vrot vroterr ccfheight
output = output + str(entry[15]) + " " #exptime
output = output + str(entry[16]) + "\n" #S/N
print entry[1],output
f = open("tracker_temp.txt","w")
f.write(output)
f.close()
candidate = entry[1]
tracker_command = "./updatecandidate.py -c "+candidate
os.system(tracker_command)
####################
### RV Standards ###
####################
query_entry = "select SPECtype,SPECobject,SPECmjd,SPEChjd,SPECrv,SPECrv_err,SPECtelescope,SPECresolution,SPECteff,SPECteff_err,SPEClogg,SPEClogg_err,SPECfeh,SPECfeh_err,SPECccfheight,SPECexptime,SPECsn"
query_entry = query_entry + " from SPEC where SPECutdate >= \""+start_date+"\" and SPECutdate <=\""+end_date+"\" and SPECinstrum=\"echelle\""
query_entry = query_entry + " and (SPECobject=\"HD100623\""
query_entry = query_entry + " or SPECobject=\"HD97343\""
query_entry = query_entry + " or SPECobject=\"HD96700\""
query_entry = query_entry + " or SPECobject=\"HD37213\""
query_entry = query_entry + " or SPECobject=\"HD34721\""
query_entry = query_entry + " or SPECobject=\"HD196761\""
query_entry = query_entry + " or SPECobject=\"HD189625\""
query_entry = query_entry + " or SPECobject=\"HD198802\""
query_entry = query_entry + ")"
exposure_info = mysql_query.query_hsmso(query_entry)
#print exposure_info
if len(exposure_info) > 0:
for entry in exposure_info:
output = ""
if entry[7] == 23000:
inst = "ANU23"
else:
inst = "WiFeS"
if entry[0] == "RV":
#output = output + entry[1] + " " #name
output = output + str(2400000 + entry[3]) + " " #hjd
output = output + str(entry[4]) + " " #RV
output = output + str(entry[5]) + " " #RVerr
#output = output + str(entry[6]) + " " #tel
output = output + inst + " " #tel
output = output + str(entry[7]) + " " #res
output = output + "0 0 0 0 0 0 0 0 " #teff terr logg loggerr feh feherr vrot vroterr
output = output + str(entry[14]) + " " #ccfheight
output = output + str(entry[15]) + " " #exptime
output = output + str(entry[16]) + "\n" #S/N
if entry[0] == "ST":
#output = output + entry[1] + " " #name
output = output + str(2400000 + entry[2]) + " " #mjd
output = output + "-999 -999 " #rv rverr
#output = output + str(entry[6]) + " " #tel
output = output + inst + " " #tel
output = output + str(entry[7]) + " " #res
output = output + str(entry[8]) + " 300 " #teff
#output = output + str(entry[9]) + " "
output = output + str(entry[10]) + " 0.3 " #logg
#output = output + str(entry[11]) + " "
output = output + str(entry[12]) + " 0.5 " #feh
output = output + "0 0 0 " #feh feherr vrot vroterr ccfheight
output = output + str(entry[15]) + " " #exptime
output = output + str(entry[16]) + "\n" #S/N
print entry[1],output
f = open("tracker_temp.txt","w")
f.write(output)
f.close()
candidate = entry[1]
tracker_command = "./updatecandidate.py -c "+candidate
os.system(tracker_command)