@davidthomas5412 Let's start with the simplest possible inference: simple Monte Carlo. We'll need the ForegroundCatalog to be able to set_prior and then draw_halo_masses, and then we'll need to compute_likelihood (using the the function @sweverett already implemented). This will give us N_samples * (N_halos + 1) numbers (if we keep the concentrations fixed for now): this is going to get into the Mb quite quickly, so we might want to think about how to store the samples (we need to keep them all). A quick and easy solution is to write out simple plain text halo catalogs that can be read back in with the existing ForegroundCatalog method, but this won't scale well: let's start thinking about database table design.
@davidthomas5412 Let's start with the simplest possible inference: simple Monte Carlo. We'll need the ForegroundCatalog to be able to
set_priorand thendraw_halo_masses, and then we'll need tocompute_likelihood(using the the function @sweverett already implemented). This will give usN_samples * (N_halos + 1)numbers (if we keep the concentrations fixed for now): this is going to get into the Mb quite quickly, so we might want to think about how to store the samples (we need to keep them all). A quick and easy solution is to write out simple plain text halo catalogs that can be read back in with the existingForegroundCatalogmethod, but this won't scale well: let's start thinking about database table design.