Skip to content
Snippets Groups Projects

corrected missing sigma factor in posterior sampler

  1. Rapidpe samples are generated with Gaussians of standard deviation=grid_size and amplitude=Marginlaized_Likelihood. During the implementation, we call numpy.random.normal function which has an amplitude
    1/2πσ1/\sqrt{2\pi\sigma}
    due to normalization. So, wherever we modify the script to have the amplitude=Marginlaized_Likelihood, need a correction factor of
    2πσ\sqrt{2\pi\sigma}
    .
  1. Rapidpe posterior samples are generated in mchirp-q (or mchirp-mtotal) with a prior uniform in mass1-mass2 as a function of mchirp-q(or mchirp-mtotal). To get the mass1, and mass2 samples we need to convert the mchirp-q (or mchirp-mtotal) sample to mass1-mass2 and assign the prior of the corresponding mchirp-q(or mchirp-mtotal) as a weight for each samples.
  • corrected the weight of m1-m2 samples in the posterior plot
  1. We usually save weighted samples. Additionally save posterior samples after resampling from weighted samples.
  • save posterior samples in addition to weighted posterior samples

Review wiki with these corrections applied: https://git.ligo.org/rapidpe-rift/rapidpe_rift_review_o4/-/wikis/AMR-Review/Review-Version-2

Review wiki before these corrections: https://git.ligo.org/rapidpe-rift/rapidpe_rift_review_o4/-/wikis/AMR-Review

Edited by VINAYA VALSAN

Merge request reports

Merge request pipeline #591228 passed

Merge request pipeline passed for e815d5e5

Test coverage 17.00% (0.33%) from 1 job
Test summary results are being parsed

Merged by VINAYA VALSANVINAYA VALSAN 1 year ago (Jan 21, 2024 5:06pm UTC)

Merge details

  • Changes merged into main with 64122a87 (commits were squashed).
  • Deleted the source branch.

Pipeline #595871 passed

Pipeline passed for 64122a87 on main

Test coverage 16.33% (0.33%) from 1 job

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
Please register or sign in to reply
Loading