TODs¶
The result of a simulation is a TOD
, which encapsulates the generated time-ordered data (and all metadata, like the pointing coordinates).
The overall signal in the
tod.plot()
It has a few useful features. We can see what our data looks like with
tod.plot()
Components¶
The total signal in each detector can be accessed as
tod.signal # returns an array
which is the sum of all of the simulated fields (e.g. noise, atmosphere, CMB) separately, contributing to the incident power. For convenience, we can also access the individual fields as
tod.get_field("atmosphere") # returns an array
We can see all the available fields with tod.fields
.
Units¶
TODs are by default in units of picowatts, but we can convert to any unit that is a combination of an SI prefix and a base unit (one of K_RJ, K_CMB, or W).
tod_in_rj_units = tod.to(units="mK_RJ")
tod_in_cmb_units = tod.to(units="uK_CMB")
Load and save¶
We can output TODs to disk, .. code-block:: python
tod.to_fits(‘filename.fits’)
or if you rather work with hdf files,
tod.to_hdf('filename.hdf5')
you can load fits files back with
tod = TOD.from_fits('filename.fits', format='MUSTANG-2')
note, that you can also load in real MUSTANG-2 data with the same command.