-
Notifications
You must be signed in to change notification settings - Fork 4
Craco-motivated update to scattering #72
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Conversation
… path script directory
…sing zDM and PATH
…some maximum width
…idth distribution
…Next stop - getting actual data values, and converging on them
… and made minor changes to reduce noise in the program
… minor cosmetic things
profxj
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Holy new code!
Great stuff. My main additional request is to add a few
tests related to the new width code.
…ptical.py to allow separation of these modules from standard zdm functionality
profxj
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
great stuff!
JordanHoffmann3
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think it makes sense to me after a lot of staring
This update is to begin modelling CRACO.
I have added code from Andy to calculate CRACO beamshapes.
However, the ability of CRACO to probe wider FRBs - and evidence from observations - suggests the need to update scattering. We now have the ability to set intrinsic width ~ (1+z) and scattering to go as (1+z)^-3
See new plot in scripts/Scattering that shows the effect of different scattering treatments,.
I have NOT updated the default values however. Nor have I determined if the new upper-lognormal (lognormal on lower-side, log-constant on upper side) function is reasonable. Next step is to work out how to fit for this!