forked from gnarayan/decat_pointings
-
Notifications
You must be signed in to change notification settings - Fork 0
/
get_silicon.py
53 lines (44 loc) · 2.2 KB
/
get_silicon.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
import numpy as np
import pandas as pd
from copy import copy
c1_corners = pd.read_csv('cornerCoords_SN-C1.dat',delim_whitespace=True)
c1_field_center_ra = 54.2743
c1_field_center_dec = -27.1116
def get_ccd_deg_offsets(df):
ccd_ra_centers = {}
ccd_dec_centers = {}
for ccd in range(1,63):
ww = df['CCD'] == ccd
ccd_ra_centers[ccd] = -1*(np.mean(df['RA'][ww].to_numpy()) - c1_field_center_ra)
ccd_dec_centers[ccd] = -1*(np.mean(df['Dec'][ww].to_numpy()) - c1_field_center_dec)
return ccd_ra_centers,ccd_dec_centers
raoff,decoff = get_ccd_deg_offsets(c1_corners)
def get_field_center_for_target_on_specific_ccd(ras,decs,ccds,extra_offset_in_arcsec=50):
field_cen_ras = []
field_cen_decs = []
for ra,dec,ccd in zip(ras,decs,ccds):
field_cen_ras.append(ra+raoff[int(ccd)]+(np.random.randint(2)-.5)*2*extra_offset_in_arcsec*0.000277778)
field_cen_decs.append(dec+decoff[int(ccd)]+(np.random.randint(2)-.5)*2*extra_offset_in_arcsec*0.000277778)
return np.array(field_cen_ras),np.array(field_cen_decs)
def is_on_silicon(RApoints,DECpoints,RAcands,DECcands):
ccds = []
raminarcsecs = []
decminarcsecs = []
for RApoint,DECpoint,RAcand,DECcand in zip(RApoints,DECpoints,RAcands,DECcands):
found = np.nan
raminarcsec = np.nan
decminarcsec = np.nan
for ccd in range(1,63):
ww = c1_corners['CCD'] == ccd
maxra = np.max(c1_corners['RA'][ww].to_numpy())+RApoint-c1_field_center_ra
minra = np.min(c1_corners['RA'][ww].to_numpy())+RApoint-c1_field_center_ra
maxdec = np.max(c1_corners['Dec'][ww].to_numpy())+DECpoint-c1_field_center_dec
mindec = np.min(c1_corners['Dec'][ww].to_numpy())+DECpoint-c1_field_center_dec
if (RAcand < maxra) & (RAcand > minra) & (DECcand < maxdec) & (DECcand > mindec):
found = copy(ccd)
raminarcsec = min([maxra-RAcand,RAcand-minra])*3600
decminarcsec = min([maxdec-DECcand,DECcand-mindec])*3600
ccds.append(found)
raminarcsecs.append(raminarcsec)
decminarcsecs.append(decminarcsec)
return np.array(ccds),np.array(raminarcsecs),np.array(decminarcsecs)