ebachelet
commented
6 months ago I am not sure to understand the issue, blend model is slightly negative but it is not the only event with such
Open rachel3834 opened 6 months ago
ebachelet
commented
6 months ago I am not sure to understand the issue, blend model is slightly negative but it is not the only event with such
Example: Gaia23dpi
MOP's fitting process has attempted to fit PSPL models with and without blending to this event, but has not managed to find a model with a reasonable blend. This leads to the target receiving a lower priority than it deserves.
Output from the fit_model_PSPL process:
`Fitting single event: Gaia23dpi Fitting event: Gaia23dpi FIT: Found 1 datasets and a total of 217 datapoints to model for event Gaia23dpi FITTOOLS: established event FITTOOLS: appended 1 telescopes check_event : Everything looks fine... FITTOOLS: Set model 1, static PSPL FITTOOLS: model 1 fit boundaries: t0: [2456900.19677, 2460307.72111] tE: [1.0, 3000.0] u0: [0.0, 2.0] Selecting Tel_0 to estimate u0, tE and fs initial_guess : Initial parameters guess SUCCESS Using guess: [2460293.6819949998, 0.7813635361213648, 133.84285499993712, 4006.182711794462, 656.9324883149255] fit : Trust Region Reflective fit SUCCESS best_model: [2460307.7211099993, 1.0108556059014966, 104.82522903736293, 7244.35960074632, -2574.11220284964] soft_l1 110.41577879330302 /mop/mop/toolbox/fittools.py:28: RuntimeWarning:
invalid value encountered in log10
FITTOOLS: model 1 fitted parameters {'t0': 2460307.72111, 't0_error': 45.1554, 'u0': 1.01086, 'u0_error': 2.98641, 'tE': 104.82523, 'tE_error': 171.70832, 'fsource_Tel_0': 7244.3596, 'fsource_Tel_0_error': 45716.01043, 'fblend_Tel_0': -2574.1122, 'fblend_Tel_0_error': 45714.97093, 'chi2': 146.203, 'piEN': 0.0, 'piEN_error': 0.0, 'piEE': 0.0, 'piEE_error': 0.0, 'red_chi2': 0.69, 'Source_magnitude': 17.75, 'Blend_magnitude': nan, 'Source_mag_error': 6.852, 'Blend_mag_error': -19.282, 'Baseline_magnitude': 17.75, 'Baseline_mag_error': 6.852, 'Fit_covariance': array([[ 2.03901060e+03, 1.29517928e+02, -7.35596159e+03, 1.98604506e+06, -1.98601400e+06], [ 1.29517928e+02, 8.91866866e+00, -5.12255455e+02, 1.36523352e+05, -1.36520219e+05], [-7.35596159e+03, -5.12255455e+02, 2.94837463e+04, -7.83992484e+06, 7.83972952e+06], [ 1.98604506e+06, 1.36523352e+05, -7.83992484e+06, 2.08995361e+09, -2.08990606e+09], [-1.98601400e+06, -1.36520219e+05, 7.83972952e+06, -2.08990606e+09, 2.08985857e+09]]), 'fit_parameters': OrderedDict([('t0', [0, [2456900.19677, 2460307.72111]]), ('u0', [1, [0.0, 2.0]]), ('tE', [2, [1.0, 3000.0]]), ('fsource_Tel_0', [3, (0.0, 7244.359600749891)]), ('fblend_Tel_0', [4, (-7244.359600749891, 7244.359600749891)])]), 'SW_test': 0.986, 'AD_test': 1.313, 'KS_test': 0.112, 'chi2_dof': 0.6896381809507844} FITTOOLS: model 1 evaluated parameters {'t0': 2460307.72111, 't0_error': 45.1554, 'u0': 1.01086, 'u0_error': 2.98641, 'tE': 104.82523, 'tE_error': 171.70832, 'fsource_Tel_0': 7244.3596, 'fsource_Tel_0_error': 45716.01043, 'fblend_Tel_0': -2574.1122, 'fblend_Tel_0_error': 45714.97093, 'chi2': 146.203, 'piEN': 0.0, 'piEN_error': 0.0, 'piEE': 0.0, 'piEE_error': 0.0, 'red_chi2': 0.69, 'Source_magnitude': 17.75, 'Blend_magnitude': nan, 'Source_mag_error': 6.852, 'Blend_mag_error': -19.282, 'Baseline_magnitude': 17.75, 'Baseline_mag_error': 6.852, 'Fit_covariance': array([[ 2.03901060e+03, 1.29517928e+02, -7.35596159e+03, 1.98604506e+06, -1.98601400e+06], [ 1.29517928e+02, 8.91866866e+00, -5.12255455e+02, 1.36523352e+05, -1.36520219e+05], [-7.35596159e+03, -5.12255455e+02, 2.94837463e+04, -7.83992484e+06, 7.83972952e+06], [ 1.98604506e+06, 1.36523352e+05, -7.83992484e+06, 2.08995361e+09, -2.08990606e+09], [-1.98601400e+06, -1.36520219e+05, 7.83972952e+06, -2.08990606e+09, 2.08985857e+09]]), 'fit_parameters': OrderedDict([('t0', [0, [2456900.19677, 2460307.72111]]), ('u0', [1, [0.0, 2.0]]), ('tE', [2, [1.0, 3000.0]]), ('fsource_Tel_0', [3, (0.0, 7244.359600749891)]), ('fblend_Tel_0', [4, (-7244.359600749891, 7244.359600749891)])]), 'SW_test': 0.986, 'AD_test': 1.313, 'KS_test': 0.112, 'chi2_dof': 0.6896381809507844} FITTOOLS: fit no-blend model? True FITTOOLS: Set model 2, static PSPL without blending FITTOOLS: model 2 fit boundaries: t0: [2456900.19677, 2460307.72111] tE: [1.0, 3000.0] u0: [0.0, 2.0] Selecting Tel_0 to estimate u0, tE and fs initial_guess : Initial parameters guess SUCCESS Using guess: [2460293.6819949998, 0.7813635361213648, 133.84285499993712, 4655.591446403771] fit : Trust Region Reflective fit SUCCESS best_model: [2460307.7211099993, 0.8096792066584301, 119.9006507575897, 4669.936050374011] soft_l1 110.56417009397416 FITTOOLS: model 2 fitted parameters {'t0': 2460307.72111, 't0_error': 12.16647, 'u0': 0.80968, 'u0_error': 0.01969, 'tE': 119.90065, 'tE_error': 9.84971, 'fsource_Tel_0': 4669.93605, 'fsource_Tel_0_error': 7.43666, 'chi2': 146.581, 'piEN': 0.0, 'piEN_error': 0.0, 'piEE': 0.0, 'piEE_error': 0.0, 'red_chi2': 0.688, 'Source_magnitude': 18.227, 'Blend_magnitude': nan, 'Source_mag_error': 0.002, 'Blend_mag_error': nan, 'Baseline_magnitude': 18.227, 'Baseline_mag_error': 0.002, 'Fit_covariance': array([[ 1.48023036e+02, -1.97707057e-01, 1.05722956e+02, -1.39474615e+01], [-1.97707057e-01, 3.87542864e-04, -1.30487519e-01, 2.95672126e-02], [ 1.05722956e+02, -1.30487519e-01, 9.70167897e+01, -2.00879906e+01], [-1.39474615e+01, 2.95672126e-02, -2.00879906e+01, 5.53038645e+01]]), 'fit_parameters': OrderedDict([('t0', [0, [2456900.19677, 2460307.72111]]), ('u0', [1, [0.0, 2.0]]), ('tE', [2, [1.0, 3000.0]]), ('fsource_Tel_0', [3, (0.0, 7244.359600749891)])]), 'SW_test': 0.986, 'AD_test': 1.299, 'KS_test': 0.112, 'chi2_dof': 0.6914205644337361} FITTOOLS: model 2 evaluated parameters {'t0': 2460307.72111, 't0_error': 12.16647, 'u0': 0.80968, 'u0_error': 0.01969, 'tE': 119.90065, 'tE_error': 9.84971, 'fsource_Tel_0': 4669.93605, 'fsource_Tel_0_error': 7.43666, 'chi2': 146.581, 'piEN': 0.0, 'piEN_error': 0.0, 'piEE': 0.0, 'piEE_error': 0.0, 'red_chi2': 0.688, 'Source_magnitude': 18.227, 'Blend_magnitude': nan, 'Source_mag_error': 0.002, 'Blend_mag_error': nan, 'Baseline_magnitude': 18.227, 'Baseline_mag_error': 0.002, 'Fit_covariance': array([[ 1.48023036e+02, -1.97707057e-01, 1.05722956e+02, -1.39474615e+01], [-1.97707057e-01, 3.87542864e-04, -1.30487519e-01, 2.95672126e-02], [ 1.05722956e+02, -1.30487519e-01, 9.70167897e+01, -2.00879906e+01], [-1.39474615e+01, 2.95672126e-02, -2.00879906e+01, 5.53038645e+01]]), 'fit_parameters': OrderedDict([('t0', [0, [2456900.19677, 2460307.72111]]), ('u0', [1, [0.0, 2.0]]), ('tE', [2, [1.0, 3000.0]]), ('fsource_Tel_0', [3, (0.0, 7244.359600749891)])]), 'SW_test': 0.986, 'AD_test': 1.299, 'KS_test': 0.112, 'chi2_dof': 0.6914205644337361} FITTOOLS: generated model lightcurve FIT: completed modeling process for Gaia23dpi /usr/local/lib/python3.10/site-packages/django/db/models/fields/init.py:1595: RuntimeWarning:
DateTimeField ReducedDatum.timestamp received a naive datetime (2018-06-29 08:15:27.243860) while time zone support is active.
FIT: Searched for existing models <QuerySet [<ReducedDatum: ReducedDatum object (15507710)>]> FIT: Stored model lightcurve for event Gaia23dpi Fitted parameters for Gaia23dpi: {'Alive': True, 'Last_fit': 2460313.2297376906, 't0': 2460307.72111, 't0_error': 45.1554, 'u0': 1.01086, 'u0_error': 2.98641, 'tE': 104.82523, 'tE_error': 171.70832, 'piEN': 0.0, 'piEN_error': 0.0, 'piEE': 0.0, 'piEE_error': 0.0, 'Source_magnitude': 17.75, 'Source_mag_error': 6.852, 'Blend_magnitude': nan, 'Blend_mag_error': -19.282, 'Baseline_magnitude': 17.75, 'Baseline_mag_error': 6.852, 'Fit_covariance': array([[ 2.03901060e+03, 1.29517928e+02, -7.35596159e+03, 1.98604506e+06, -1.98601400e+06], [ 1.29517928e+02, 8.91866866e+00, -5.12255455e+02, 1.36523352e+05, -1.36520219e+05], [-7.35596159e+03, -5.12255455e+02, 2.94837463e+04, -7.83992484e+06, 7.83972952e+06], [ 1.98604506e+06, 1.36523352e+05, -7.83992484e+06, 2.08995361e+09, -2.08990606e+09], [-1.98601400e+06, -1.36520219e+05, 7.83972952e+06, -2.08990606e+09, 2.08985857e+09]]), 'chi2': 146.203, 'red_chi2': 0.69, 'KS_test': 0.112, 'AD_test': 1.313, 'SW_test': 0.986} Target post save hook: Gaia23dpi created: False FIT: Stored model parameters for event Gaia23dpi`