-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathLMC_81_mirisim.py
146 lines (114 loc) · 6.44 KB
/
LMC_81_mirisim.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
# script to generate a mirisim scene from the planetary nebula SMP LMC 81
import numpy as np
import matplotlib.pyplot as plt
import pdb
from mirisim.skysim import Background, sed, Point, Skycube, ExternalSed
from mirisim.skysim import wrap_pysynphot as wS
from mirisim.config_parser import SimConfig, SimulatorConfig, SceneConfig
from mirisim import MiriSimulation
# this script will create a scene file for the LMC planetary nebula SMP LMC 81, which is my (placeholder) target for LRS wave calibration. I will use the Spitzer IRS spectrum as external SED for the simulation
# With the ETC I've calculated that I can get a SNR at 12 micron of 120 fro Ngroups = 15, Nint = 1. I will use a standard 2-point nod
# initialise the point source:
pn_sed = ExternalSed(sedfile='LMC_81_Mirisim_spec.txt')
pn = Point(Cen=[0,0])
pn.set_SED(pn_sed)
# now write out to a FITS file
fov = np.array([[-5, 5], [-5,5]])
spat_samp = 0.05
# now we can add a background if requested:
bg = Background(level='low', gradient=5., pa=15.)
scene = pn + bg
# write some output
print(scene)
#scene.writecube(cubefits='lmc_81_scene.fits', FOV=fov, spatsampling=spat_samp, wrange=[5., 15.], wsampling=0.05, overwrite=True, time=0.0)
# now we also want to create the ini file
targetlist = [pn]
## export to ini file -- THIS DOESN'T WORK AS DESCRIBED IN THE JUPYTER NOTEBOOK
scene_config = SceneConfig.makeScene(loglevel=0,
background=bg,
targets = targetlist)
#scene_config.write('lmc_81_scene.ini')
# FIRST SIMULATION: NO DITHERING
#sim_config_nodither = SimConfig.makeSim(
# name = 'LMC-81_lrsslit_nodither', # name given to simulation
# scene = 'lmc_81_scene.ini', # name of scene file to input
# rel_obsdate = 1.0, # relative observation date (0 = launch, 1 = end of 5 yrs)
# POP = 'IMA', # Component on which to center (Imager or MRS)
# ConfigPath = 'LRS_SLIT', # Configure the Optical path (MRS sub-band)
# Dither = False, # Don't Dither
# StartInd = 1, # start index for dither pattern [NOT USED HERE]
# NDither = 2, # number of dither positions [NOT USED HERE]
# DitherPat = 'lrs_recommended_dither.dat', # dither pattern to use [NOT USED HERE]
# disperser = 'SHORT', # [NOT USED HERE]
# detector = 'SW', # [NOT USED HERE]
# mrs_mode = 'SLOW', # [NOT USED HERE]
# mrs_exposures = 10, # [NOT USED HERE]
# mrs_integrations = 3, # [NOT USED HERE]
# mrs_frames = 5, # [NOT USED HERE]
# ima_exposures = 1, # number of exposures
# ima_integrations = 1, # number of integrations
# ima_frames = 15, # number of groups (for MIRI, # Groups = # Frames)
# ima_mode = 'FAST', # Imager read mode (default is FAST ~ 2.3 s)
# filter = 'P750L', # Imager Filter to use
# readDetect = 'FULL' # Portion of detector to read out
#)
#sim_config_nodither.write('LMC-81_simulation_slit_nodither.ini')
#simulator_config = SimulatorConfig.from_default()
#lmc81_sim = MiriSimulation(sim_config_nodither, scene_config, simulator_config)
#lmc81_sim.run()
# SECOND SIMULATION: ALONG SLIT NOD
sim_config_dither = SimConfig.makeSim(
name = 'LMC-81_lrsslit_dither', # name given to simulation
scene = 'lmc_81_scene.ini', # name of scene file to input
rel_obsdate = 1.0, # relative observation date (0 = launch, 1 = end of 5 yrs)
POP = 'IMA', # Component on which to center (Imager or MRS)
ConfigPath = 'LRS_SLIT', # Configure the Optical path (MRS sub-band)
Dither = True, # Nod along the slit
StartInd = 1, # start index for dither pattern [NOT USED HERE]
NDither = 2, # number of dither positions [NOT USED HERE]
DitherPat = 'lrs_recommended_dither.dat', # dither pattern to use [NOT USED HERE]
disperser = 'SHORT', # [NOT USED HERE]
detector = 'SW', # [NOT USED HERE]
mrs_mode = 'SLOW', # [NOT USED HERE]
mrs_exposures = 10, # [NOT USED HERE]
mrs_integrations = 3, # [NOT USED HERE]
mrs_frames = 5, # [NOT USED HERE]
ima_exposures = 1, # number of exposures
ima_integrations = 1, # number of integrations
ima_frames = 15, # number of groups (for MIRI, # Groups = # Frames)
ima_mode = 'FAST', # Imager read mode (default is FAST ~ 2.3 s)
filter = 'P750L', # Imager Filter to use
readDetect = 'FULL' # Portion of detector to read out
)
sim_config_dither.write('LMC-81_simulation_slitdither.ini')
simulator_config = SimulatorConfig.makeSimulator(take_webbPsf=True)
lmc81_sim2 = MiriSimulation(sim_config_dither, scene_config, simulator_config)
lmc81_sim2.run()
# THIRD SIMULATION: SLITLESS
#sim_config_slitless = SimConfig.makeSim(
# name = 'LMC-81_lrsslitless', # name given to simulation
# scene = 'lmc_81_scene.ini', # name of scene file to input
# rel_obsdate = 1.0, # relative observation date (0 = launch, 1 = end of 5 yrs)
# POP = 'IMA', # Component on which to center (Imager or MRS)
# ConfigPath = 'LRS_SLITLESS', # Configure the Optical path (MRS sub-band)
# Dither = False, # Nod along the slit
# StartInd = 1, # start index for dither pattern [NOT USED HERE]
# NDither = 2, # number of dither positions [NOT USED HERE]
# DitherPat = 'lrs_recommended_dither.dat', # dither pattern to use [NOT USED HERE]
# disperser = 'SHORT', # [NOT USED HERE]
# detector = 'SW', # [NOT USED HERE]
# mrs_mode = 'SLOW', # [NOT USED HERE]
# mrs_exposures = 10, # [NOT USED HERE]
# mrs_integrations = 3, # [NOT USED HERE]
# mrs_frames = 5, # [NOT USED HERE]
# ima_exposures = 1, # number of exposures
# ima_integrations = 5, # number of integrations
# ima_frames = 60, # number of groups (for MIRI, # Groups = # Frames)
# ima_mode = 'FAST', # Imager read mode (default is FAST ~ 2.3 s)
# filter = 'P750L', # Imager Filter to use
# readDetect = 'SLITLESSPRISM' # Portion of detector to read out
#)
#sim_config_slitless.write('LMC-81_simulation_slitless.ini')
#simulator_config = SimulatorConfig.from_default()
#lmc81_sim3 = MiriSimulation(sim_config_slitless, scene_config, simulator_config)
#lmc81_sim3.run()