Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision 		Previous revision
comb_analysis [2022/02/22 10:23]
theoastro created 		comb_analysis [2022/02/27 08:46] (current)
theoastro
Line 1: Line 1:
 ===== Combined Analysis ===== ===== Combined Analysis =====
 +
 +NMMA is capable of performing combined analyses to constrain:
 +  *  the neutron star equation of state (EOS), and
 +  *  the Hubble Constant.
 +
 +In the following, we will take the EOS analysis as an example.
 +
 +==== EOS analysis ====
 +
 +**Generate a simulation set**
 +
 +First of all, you need to create an output directory, this output will host all the data that will be used to constrain the EOS.
 +  mkdir -p ./output
 +
 +Running the following command line will generate a json file (injection.json) with the BILBY processing of compact binary merging events. We take here binaries of type BNS (NSBH is also an option). This injection comprises a simulation set of parameters: ''luminosity_distance'', ''log10_mej_wind'', ''KNphi'', ''inclination_EM'', ''KNtimeshift'', ''geocent_time'' for the Bu2019lm model. This creates an injection.json file in the ./output directory.
 +
 +  nmma_create_injection --prior-file ./priors/Bu2019lm.prior --eos-file ./example_files/eos/ALF2.dat --binary-type BNS -f ./output/injection --n-injection 100 --original-parameters --extension json
 +
 +**Light curve posterior**
 +
 +EM data will house the posteriors of the electromagnetic data you will produce: in particular the lc.csv (./example_files/csv_lightcurve/outdir/macroeventID, where macroeventID in range(0, 100)) light curves. We now compute posteriors using NMMA on this simulated set of 100 events, of which we assume a fraction is detectable by Zwicky Transient Facility (ZTF). The result can be find at ./output/EMdata
 +
 +for macroeventID in {0..99}
 +
 +  do
 +    mkdir -p ./output/EMdata/outdir/$macroeventID/
 +    light_curve_analysis --model Bu2019lm --svd-path ./svdmodels --gptype tensorflow --outdir ./output/EMdata/outdir/$macroeventID --label injection_Bu2019lm --prior ./priors/Bu2019lm.prior --tmin 0 --tmax 7 --dt 0.5 --error-budget 1.0 --nlive 256 --Ebv-max 0 --injection ./output/injection.json --injection-num $macroeventID --injection-detection-limit 22,22,22 --injection-outfile ./output/EMdata/outdir/$macroeventID/lc.csv --generation-seed 42 --filters g,r,i --ztf-sampling --ztf-uncertainties --plot --remove-nondetections --optimal-augmentation --optimal-augmentation-filters u,g,r,i,z,y,J,H,K --optimal-augmentation-N 100
 +  done
 +
 +
 +
 +
 +
Last modified: le 2022/02/22 10:23

Page Tools