Add LMC survey

Jürgen Knödlseder
1 parent e950b6d6
Showing 1 changed file with 63 additions and 4 deletions
examples/make_pointings.py
@@ -240,6 +240,59 @@ def set_gc(lst=True):
     return obsdef
+# ============= #
+# Setup LMC KSP #
+# ============= #
+def set_lmc(hours=250.0, lst=True):
+    """
+    Setup LMC KSP.
+    """
+    # Initialise observation definition
+    obsdef = []
+
+    # Set LMC centre
+    centre = gammalib.GSkyDir()
+    centre.radec_deg(80.0, -69.0)
+
+    # Set offset angle and number of pointings
+    offset = 1.0 # degrees
+    n_pnt  = 12
+
+    # Prepare computations
+    dphi     = 360.0/n_pnt
+    duration = hours*3600.0/float(n_pnt)
+
+    # Loop over pointings
+    for ipnt in range(n_pnt):
+
+        # Compute pointing direction
+        pnt = centre.copy()
+        pnt.rotate_deg(ipnt*dphi, offset)
+        lon = pnt.l_deg()
+        lat = pnt.b_deg()
+
+        # Set positions and duration
+        obs = {'lon': lon, 'lat': lat, 'duration': duration}
+
+        # Add IRF
+        caldb = "aar"
+        if lst:
+            irf = "DESY20140105_50h"
+        else:
+            irf = "DESY20140105_50h_noLST"
+        obs['caldb'] = caldb
+        obs['irf']   = irf
+
+        # Append observation
+        obsdef.append(obs)
+
+    # Dump statistics
+    print("Number of pointings: %d (%.2f sec)" % (n_pnt,duration))
+
+    # Return observation definition
+    return obsdef
+
+
 # ======================================= #
 # Write observation definition dictionary #
 # ======================================= #
@@ -302,13 +355,14 @@ if __name__ == &#39;__main__&#39;:
     if need_help:
         print("Usage: make_pointing.py [OPTIONS]")
         print("     -h       Display this usage message")
-        print("     gps      Galactic plane scan (2 row scheme)")
-        print("     gps3     Galactic plane scan (3 row scheme)")
+        print("     gps      Galactic plane survey (2 row scheme)")
+        print("     gps3     Galactic plane survey (3 row scheme)")
         print("     extgal   Extragalactic survey")
-        print("     gc       Galactic centre")
+        print("     gc       Galactic centre survey")
+        print("     lmc      LMC survey")
         sys.exit()
-    # Galactic plane scan
+    # Galactic plane survey
     if obsname == "gps":
         obsdef = set_gps(lst=True)
         write_obsdef("gps.dat", obsdef)
@@ -326,6 +380,11 @@ if __name__ == &#39;__main__&#39;:
         obsdef = set_gc(lst=True)
         write_obsdef("gc.dat", obsdef)
+    # LMC
+    elif obsname == "lmc":
+        obsdef = set_lmc(lst=True)
+        write_obsdef("lmc.dat", obsdef)
+
     # Invalid pattern
     else:
         print("Unknown option \""+obsname+"\"")