edges.frequencies

Tools for dealing with frequency arrays.

edges.frequencies.clip_freqs(freq: Annotated[Quantity, PhysicalType('frequency')], low: Annotated[Quantity, PhysicalType('frequency')], high: Annotated[Quantity, PhysicalType('frequency')]) → Annotated[Quantity, Unit('MHz')]

Clip an array of frequencies within a range.

edges.frequencies.edges_freq_mask(low: Annotated[Quantity, PhysicalType('frequency')], high: Annotated[Quantity, PhysicalType('frequency')]) → ndarray[tuple[Any, ...], dtype[bool]]

Create a mask for the raw EDGES spectrum frequencies.

edges.frequencies.edges_raw_freqs(f_low: Quantity, PhysicalType('frequency')]=<Quantity 0. MHz>, f_high: Quantity, PhysicalType('frequency')]=<Quantity inf MHz>) → Annotated[Quantity, Unit('MHz')]

Get the raw frequency array of the EDGES spectrometer.

Parameters:

f_low, f_high – A frequency range to keep.

Returns:

freqs – The raw frequencies of the spectrometer.

Notes

This is correct. The channel width is the important thing. The channel width is given by the FFT. We actually take 32678*2 samples of data at 400 Mega-samples per second. We only use the first half of the samples (since it’s real input). Regardless, the frequency channel width is thus 400 MHz / (32678*2) == 200 MHz / 32678 ~ 6.103 kHz

edges.frequencies.get_mask(freq: ~typing.Annotated[~astropy.units.quantity.Quantity, PhysicalType('frequency')], low: ~typing.Annotated[~astropy.units.quantity.Quantity, PhysicalType('frequency')] = <Quantity 0. MHz>, high: ~typing.Annotated[~astropy.units.quantity.Quantity, PhysicalType('frequency')] = <Quantity inf MHz>) → ndarray[tuple[Any, ...], dtype[bool]]

Get a mask from a frequency array between a given range.