| Both sides previous revision
Previous revision
Next revision
|
Previous revision
|
lc_generation [2022/02/15 12:47]
theoastro |
lc_generation [2023/06/02 17:05] (current)
theoastro |
| Note that the parameters vary from model to model - an overview can be found [[https://github.com/nuclear-multimessenger-astronomy/nmma/blob/c9f55663f27625ca7c3aa1b8d74dcb231ca060e7/nmma/em/model.py#L16|here]]. Some models such as ''Bu2019lm'' require to embed a SVD grid for light curve computation which can be found [[https://github.com/nuclear-multimessenger-astronomy/nmma/tree/main/svdmodels|here]]. With these, a light curve object can be instantiated as follows: | Note that the parameters vary from model to model - an overview can be found [[https://github.com/nuclear-multimessenger-astronomy/nmma/blob/c9f55663f27625ca7c3aa1b8d74dcb231ca060e7/nmma/em/model.py#L16|here]]. Some models such as ''Bu2019lm'' require to embed a SVD grid for light curve computation which can be found [[https://github.com/nuclear-multimessenger-astronomy/nmma/tree/main/svdmodels|here]]. With these, a light curve object can be instantiated as follows: |
| |
| t = np.arange(tmin=0.1, tmax=20.0 ,deltat=0.1) | tmin=0.1 |
| | tmax=20.0 |
| | deltat=0.1 |
| | t = np.arange(tmin, tmax, deltat) |
| lc_model = nmma.em.model.SVDLightCurveModel(model='Bu2019lm', sample_times = t, svd_path = "nmma/svdmodels/", parameter_conversion=None, mag_ncoeff=None, lbol_ncoeff=None) | lc_model = nmma.em.model.SVDLightCurveModel(model='Bu2019lm', sample_times = t, svd_path = "nmma/svdmodels/", parameter_conversion=None, mag_ncoeff=None, lbol_ncoeff=None) |
| |
| 'inclination_EM': 0, | 'inclination_EM': 0, |
| 'luminosity_distance': 40} | 'luminosity_distance': 40} |
| lbol, mag= lc_model.generate_lightcurve(t, params) | lbol, mag= lc_model.generate_lightcurve(t, params_range) |
| |
| Through the filter keys such as u,g,r,i,z,y,J,H,K, the light curves in different photometric bands can be obtained, e.g., using ''mag['u']''. An example plot of kilonova light curves in different photometric bands is shown below for AT2017gfo. | Through the filter keys such as u,g,r,i,z,y,J,H,K, the light curves in different photometric bands can be obtained, e.g., using ''mag['u']''. An example plot of kilonova light curves in different photometric bands is shown below for AT2017gfo. |
| grb_model = nmma.em.model.GRBLightCurveModel(t_day, resolution=12, jetType=0) | grb_model = nmma.em.model.GRBLightCurveModel(t_day, resolution=12, jetType=0) |
| params_range = { | params_range = { |
| | |
| # ranges for source parameters | |
| 'inclination_EM': 0, | 'inclination_EM': 0, |
| 'log10_E0': 50., | 'log10_E0': 50., |
| 'log10_epsilon_e':-1.6, | 'log10_epsilon_e':-1.6, |
| 'log10_epsilon_B':-2., | 'log10_epsilon_B':-2., |
| # luminosity distance remains fixed [Mpc] | |
| 'luminosity_distance': 40,} | 'luminosity_distance': 40,} |
| |
