Skip to content

Commit

Permalink
update
Browse files Browse the repository at this point in the history
  • Loading branch information
daniel-lyon committed Oct 27, 2023
1 parent 9caaef5 commit 82a9bfa
Show file tree
Hide file tree
Showing 13 changed files with 21 additions and 40 deletions.
27 changes: 10 additions & 17 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -43,12 +43,20 @@ gleam_0856.template()
gleam_0856.fft()
```

Citing zfinder
----------
If you used zfinder in your paper, please be sure to cite it: (link when published)

Contributing
----------
If you have found a bug, please raise an issue. Alternatively, if you know how to fix it, you can be become a contributor by submitting a pull request.

Physical Methodology
----------

A three dimensional `.fits` data cube of right ascension (ra), declination (dec), and frequency is required to find the redshift. Given a target ra and dec, the flux at the target location is calculated across all frequencies. To find the redshift of the source, zfinder searches for emission lines of an element or molecule. Carbon Monoxide (115.2712 GHz, ~2.6mm) is a tracer of cold gas - it is theorised that Active Galactic Nuclei (AGN) at high redshifts have lots of cold gas circling around the black hole. Emission lines appear at integer multiples of the corresponding fundamental frequency, meaning for CO emission lines, they appear at 230.5424 GHz, 345.8136 GHz, etc (theoretically up to infinity). As the redshift of a source increases, the new fundamental frequency becomes equal to 1/(1+z), where z is the redshift. Gaussians at integer multiples of the fundamental frequency are overlayed onto the flux data of the source and the chi-squared is calculated by incrementing z by small changes. The minimum chi-squared corresponds to the most likely redshift of the source. At high redshifts (z > ~5), there should be two emission lines present in the flux data. The presence of two or more evenly spaced emission lines indicates a single real redshift solution; there can only be one redshift at which the gaussians will overlay all lines. Having only one line is a clear indicator of lower redshift (z < ~5) but the redshift solution is also indeterminant as there can be several redshifts where only one emission line is present.

Example Source: J085614 + 022400
Example Source: SPT 0345-47
----------

### Template Fitting Method
Expand All @@ -57,12 +65,6 @@ Example Source: J085614 + 022400
<img src="./figures/template_flux.png">
</h1>

### Template Fitting Animation
![GIF Description](./figures/Animations/flux_animation.gif)
![GIF Description](./figures/Animations/redshift_animation.gif)
![GIF Description](./figures/Animations/chi2_animation.gif)


### FFT Fitting Method
<h1 align="left">
<img src="./figures/fft_chi2.png">
Expand All @@ -73,13 +75,4 @@ Example Source: J085614 + 022400
<h1 align="center">
<img src="./figures/fft_per_pixel.png">
<img src="./figures/zoomed_fft_per_pixel.png">
</h1>


Citing zfinder
----------
If you used zfinder in your paper, please be sure to cite it: (link when published)

Contributing
----------
If you have found a bug, please raise an issue. Alternatively, if you know how to fix it, you can be become a contributor by submitting a pull request.
</h1>
4 changes: 0 additions & 4 deletions docs/conf.py
Original file line number Diff line number Diff line change
Expand Up @@ -21,10 +21,6 @@
'sslf'
]

# from unittest.mock import MagicMock
# for mod in autodoc_mock_imports:
# sys.modules[mod] = MagicMock()s

sys.path.insert(0, os.path.abspath('..'))

import sphinx_rtd_theme
Expand Down
Binary file modified docs/requirements.txt
Binary file not shown.
Binary file modified figures/fft_chi2.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified figures/fft_flux.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified figures/fft_per_pixel.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified figures/template_chi2.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified figures/template_flux.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified figures/zoomed_fft_per_pixel.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
6 changes: 1 addition & 5 deletions runs/GLEAM_0856/GLEAM_0856.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -4,9 +4,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"# Running ZFINDER on GLEAM_0856\n",
"\n",
"Source fits files hosted at: "
"# Running ZFINDER on GLEAM_0856"
]
},
{
Expand All @@ -20,8 +18,6 @@
"fitsfile = '0856_cube_c0.4_nat_80MHz_taper3.fits'\n",
"ra = '08:56:14.8'\n",
"dec = '02:24:00.6'\n",
"# ra = '08:56:14.76'\n",
"# dec = '02:23:59.6'\n",
"aperture_radius = 3\n",
"transition = 115.271\n",
"\n",
Expand Down
4 changes: 1 addition & 3 deletions runs/GLEAM_0856_blob/GLEAM_0856_blob.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -4,9 +4,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"# Running ZFINDER on GLEAM_0856\n",
"\n",
"Source fits files hosted at: "
"# Running ZFINDER on GLEAM_0856"
]
},
{
Expand Down
16 changes: 7 additions & 9 deletions runs/SPT_0345-47/SPT_0345-47.ipynb

Large diffs are not rendered by default.

4 changes: 2 additions & 2 deletions setup.py
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
from setuptools import setup, find_packages

VERSION = '1.0.0'
VERSION = '1.0.2'
DESCRIPTION = 'Redshift finding algorithm'
LONG_DESCRIPTION = 'Find the redshift of high redshift radio galaxies'

Expand All @@ -13,7 +13,7 @@
long_description_content_type="text/markdown",
long_description=LONG_DESCRIPTION,
packages=find_packages(),
install_requires=['numpy', 'matplotlib', 'scipy', 'astropy', 'photutils', 'PyAstronomy', 'sslf', 'radio_beam', 'sphinx', 'sphinx_rtd_theme'],
install_requires=['numpy', 'matplotlib', 'scipy', 'astropy', 'photutils', 'PyAstronomy', 'sslf', 'radio_beam', 'tqdm'],
keywords=['python', 'visualisation', 'redshift', 'galaxy-evolution', 'cosmology',
'epoch-of-reionisation', 'radio-galaxies', 'active-galactic-nuclei',
'high-redshift', 'galaxies', 'black-holes'],
Expand Down

0 comments on commit 82a9bfa

Please sign in to comment.