How do we make catalogues including a column for secondary solution in redshift?

[eclake]

Subject says it all :)

[jacopo-chevallard]

Hi Emma, I'm not sure what you mean, but I don't think there is any way to account for multiple solutions in the current version of pyp-beagle, although I do have scripts allowing you to do so...

[eclake]

Thanks Jacopo. I was hoping to put secondary redshift solutions into the catalogues I'm putting together for Kit's photometry for the HST master catalogue - might you be able to send me scripts if it's not possible with pyp-beagle?

[jacopo-chevallard]

So, there is a pyp_beagle module beagle_multinest_catalogue.py which allows you to load some information (e.g. posterior mean) for different solutions (or modes) from the *_MNstats.dat files created by MultiNest.

A similar approach is also implemented in a script I used for JWST, to creat "Value Added Catalogues" of physical properties from Beagle runs make_value_added_catalogue.py. In this script I also implemented a function get_mode_rows which does a more appropriate thing: it reads the MultiNest *_MNpost_separate.dat files to identify which rows in a Beagle output file pertain to which solution. Once you have matched the rows to the solutions, you can then use the Beagle output in a "standard" way and perform any calculation on the model parameters and derived quantities. My initial idea was to integrate such functionality into pyp_beagle (in the beagle_summary_catalogue.py module), but in the end I never had the time to do it..

Let me know if you need any other detail!

[eclake]

Thank you!