From 16355319571052e5dfd593e56d8136711751850d Mon Sep 17 00:00:00 2001 From: Fabian Jankowski Date: Thu, 27 Jun 2024 23:25:36 +0200 Subject: [PATCH] Explain how to estimate the scattering-corrected DM reliably by iterating. --- docs/FAQ.md | 10 ++++++++++ 1 file changed, 10 insertions(+) diff --git a/docs/FAQ.md b/docs/FAQ.md index 0961ba3..bcc9299 100644 --- a/docs/FAQ.md +++ b/docs/FAQ.md @@ -32,3 +32,13 @@ before installing `scatfit`. Then `$ pip3 install scatfit` should work fine. Maybe replace `pip3` by `pip` depending on your Python installation. + +## How do I estimate the scattering-corrected DM reliably? ## + +There are several ways to estimate a pulse or burst's true dispersion measure (DM). For instance, one can optimise the burst profile for the total band-averaged S/N (S/N-maximising DM) or profile structure (structure-maximising DM). Scattering introduces a typically single-sided exponential tail that interferes with the classical DM-estimation techniques by shifting the profile centres in a frequency and DM-dependent way. Hence, one needs to account for scattering by estimating a scattering-corrected DM, which usually differs from the above DM values, sometimes quite significantly. + +`scatfit` blindly dedisperses the data at the fiducial DM you provide in its command line call. This behaviour is advantageous, as it gives the user the flexibility to measure the optimum DM using external software, such as `DM_phase`. `scatfit` then performs its sub-banded scattering fits, after which it analyses the fit residuals to determine a scattering-corrected DM. + +To measure the scattering-corrected DM reliably, one must run `scatfit` iteratively several times by inputting the last determined scattering-corrected DM in the new call until the output DM converges to a stable value. This occurs when the output DM changes less than its fit uncertainty or an externally pre-defined value. The leading edge of the burst profile should be as sharp as possible (in the absence of DM-smearing) and the tail well-described by the scattering model. + +Three to four iterations should be enough to reach a stable scattering-corrected DM if the initial DM was already close to the corrected value. Something must be wrong if the DM does not stabilise after 5-6 iterations.