Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision 		Previous revision
comb_analysis [2022/02/22 10:33]
theoastro [EOS analysis] 		comb_analysis [2022/02/27 08:46] (current)
theoastro
Line 5: Line 5:
   *  the Hubble Constant.   *  the Hubble Constant.
  
-In the following, we will take as an example the EOS analysis.+In the following, we will take the EOS analysis as an example.
  
 ==== EOS analysis ==== ==== EOS analysis ====
  
 **Generate a simulation set** **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. 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   mkdir -p ./output
Line 16: Line 17:
  
   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   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:33

Page Tools