Karim Sabri / ctools

Commit 9f4171afedbe56958493132a609720d729e7cd8b

Authored by Jürgen Knödlseder
1 parent d8896b9b

Add site dependent IRF information

Inline Side-by-side
Showing 1 changed file with 22 additions and 9 deletions
examples/make_pointings.py
  View file @9f4171a
... ... @@ -31,7 +31,7 @@ import math
31 31 # ======================= #
32 32 # Setup double row scheme #
33 33 # ======================= #
34   -def set_double_row(lmin=-30.0, lmax=30.0, separation=3.0, hours=100.0):
  34 +def set_double_row(lmin=-30.0, lmax=30.0, separation=3.0, hours=100.0, site=None):
35 35 """
36 36 Setup double row scheme.
37 37 """
... ... @@ -62,6 +62,18 @@ def set_double_row(lmin=-30.0, lmax=30.0, separation=3.0, hours=100.0):
62 62 obsdef.append(obs)
63 63 lon += separation
64 64  
  65 + # Optionally add-in site dependent IRF
  66 + if site != None:
  67 + if site == "North":
  68 + caldb = "tenerife"
  69 + irf = "DESY20140105_50h"
  70 + else:
  71 + caldb = "aar"
  72 + irf = "DESY20140105_50h"
  73 + for obs in obsdef:
  74 + obs['caldb'] = caldb
  75 + obs['irf'] = irf
  76 +
65 77 # Return observation definition
66 78 return obsdef
67 79  
... ... @@ -77,19 +89,19 @@ def set_gps():
77 89 obsdef = []
78 90  
79 91 # Add inner region South
80   - obsdef.extend(set_double_row(lmin=-60.0, lmax=60.0, hours=780))
  92 + obsdef.extend(set_double_row(lmin=-60.0, lmax=60.0, hours=780, site="South"))
81 93  
82 94 # Add Vela & Carina region
83   - obsdef.extend(set_double_row(lmin=240.0, lmax=300.0, hours=180))
  95 + obsdef.extend(set_double_row(lmin=240.0, lmax=300.0, hours=180, site="South"))
84 96  
85 97 # Add 210-240 region
86   - obsdef.extend(set_double_row(lmin=210.0, lmax=240.0, hours=60))
  98 + obsdef.extend(set_double_row(lmin=210.0, lmax=240.0, hours=60, site="South"))
87 99  
88 100 # Add Cygnus, Perseus
89   - obsdef.extend(set_double_row(lmin=60.0, lmax=150.0, hours=450))
  101 + obsdef.extend(set_double_row(lmin=60.0, lmax=150.0, hours=450, site="North"))
90 102  
91 103 # Add Anticentre
92   - obsdef.extend(set_double_row(lmin=150.0, lmax=210.0, hours=150))
  104 + obsdef.extend(set_double_row(lmin=150.0, lmax=210.0, hours=150, site="North"))
93 105  
94 106 # Return observation definition
95 107 return obsdef
... ... @@ -106,7 +118,7 @@ def write_obsdef(filename, obsdef):
106 118 file = open(filename, 'w')
107 119  
108 120 # Write header
109   - file.write("ra,dec,duration\n")
  121 + file.write("ra,dec,duration,caldb,irf\n")
110 122  
111 123 # Loop over pointings
112 124 for obs in obsdef:
... ... @@ -124,7 +136,8 @@ def write_obsdef(filename, obsdef):
124 136 dec = obs["dec"]
125 137  
126 138 # Write information
127   - file.write("%8.4f,%8.4f,%.4f\n" % (ra, dec, obs["duration"]))
  139 + file.write("%8.4f,%8.4f,%.4f,%s,%s\n" % \
  140 + (ra, dec, obs['duration'], obs['caldb'], obs['irf']))
128 141  
129 142 # Close file
130 143 file.close()
... ... @@ -161,4 +174,4 @@ if __name__ == '__main__':
161 174 obsdef = set_gps()
162 175  
163 176 # Write observation definition file
164   - write_obsdef("test.dat", obsdef)
  177 + write_obsdef("gps.dat", obsdef)
... ...