diff --git a/docs/.buildinfo b/docs/.buildinfo index 446e98f9..30fb32af 100644 --- a/docs/.buildinfo +++ b/docs/.buildinfo @@ -1,4 +1,4 @@ # Sphinx build info version 1 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done. -config: 2acbe60213cec1ff5d56a56e4cbac976 +config: d66d4b627a3641eee5ee7004022669b0 tags: 645f666f9bcd5a90fca523b33c5a78b7 diff --git a/docs/.doctrees/environment.pickle b/docs/.doctrees/environment.pickle index d8c2760e..ee04eef4 100644 Binary files a/docs/.doctrees/environment.pickle and b/docs/.doctrees/environment.pickle differ diff --git a/docs/MulensModel.binarylens.html b/docs/MulensModel.binarylens.html index 84989563..7dc3d146 100644 --- a/docs/MulensModel.binarylens.html +++ b/docs/MulensModel.binarylens.html @@ -6,7 +6,7 @@ - MulensModel.binarylens module — MulensModel 2.22.0 documentation + MulensModel.binarylens module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.binarylensimports.html b/docs/MulensModel.binarylensimports.html index 0df20d99..819bb42e 100644 --- a/docs/MulensModel.binarylensimports.html +++ b/docs/MulensModel.binarylensimports.html @@ -6,7 +6,7 @@ - MulensModel.binarylensimports module — MulensModel 2.22.0 documentation + MulensModel.binarylensimports module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.binarylenswithshear.html b/docs/MulensModel.binarylenswithshear.html index 820887bd..edfd4e91 100644 --- a/docs/MulensModel.binarylenswithshear.html +++ b/docs/MulensModel.binarylenswithshear.html @@ -6,7 +6,7 @@ - MulensModel.binarylenswithshear module — MulensModel 2.22.0 documentation + MulensModel.binarylenswithshear module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.caustics.html b/docs/MulensModel.caustics.html index 73a42188..7a5adb2a 100644 --- a/docs/MulensModel.caustics.html +++ b/docs/MulensModel.caustics.html @@ -6,7 +6,7 @@ - MulensModel.caustics module — MulensModel 2.22.0 documentation + MulensModel.caustics module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.causticspointwithshear.html b/docs/MulensModel.causticspointwithshear.html index 45bb2ed8..4978e29d 100644 --- a/docs/MulensModel.causticspointwithshear.html +++ b/docs/MulensModel.causticspointwithshear.html @@ -6,7 +6,7 @@ - MulensModel.causticspointwithshear module — MulensModel 2.22.0 documentation + MulensModel.causticspointwithshear module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.causticswithshear.html b/docs/MulensModel.causticswithshear.html index 35021ceb..94652c65 100644 --- a/docs/MulensModel.causticswithshear.html +++ b/docs/MulensModel.causticswithshear.html @@ -6,7 +6,7 @@ - MulensModel.causticswithshear module — MulensModel 2.22.0 documentation + MulensModel.causticswithshear module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.coordinates.html b/docs/MulensModel.coordinates.html index f70edac8..5938648c 100644 --- a/docs/MulensModel.coordinates.html +++ b/docs/MulensModel.coordinates.html @@ -6,7 +6,7 @@ - MulensModel.coordinates module — MulensModel 2.22.0 documentation + MulensModel.coordinates module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.event.html b/docs/MulensModel.event.html index 62679111..90cd0848 100644 --- a/docs/MulensModel.event.html +++ b/docs/MulensModel.event.html @@ -6,7 +6,7 @@ - MulensModel.event module — MulensModel 2.22.0 documentation + MulensModel.event module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.fitdata.html b/docs/MulensModel.fitdata.html index ef2eec86..dca5d7ff 100644 --- a/docs/MulensModel.fitdata.html +++ b/docs/MulensModel.fitdata.html @@ -6,7 +6,7 @@ - MulensModel.fitdata module — MulensModel 2.22.0 documentation + MulensModel.fitdata module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.horizons.html b/docs/MulensModel.horizons.html index 4c64df24..2a9d357b 100644 --- a/docs/MulensModel.horizons.html +++ b/docs/MulensModel.horizons.html @@ -6,7 +6,7 @@ - MulensModel.horizons module — MulensModel 2.22.0 documentation + MulensModel.horizons module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.html b/docs/MulensModel.html index d864ca0f..fccbe71b 100644 --- a/docs/MulensModel.html +++ b/docs/MulensModel.html @@ -6,7 +6,7 @@ - MulensModel package — MulensModel 2.22.0 documentation + MulensModel package — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.limbdarkeningcoeffs.html b/docs/MulensModel.limbdarkeningcoeffs.html index 9016a559..ea1adbe2 100644 --- a/docs/MulensModel.limbdarkeningcoeffs.html +++ b/docs/MulensModel.limbdarkeningcoeffs.html @@ -6,7 +6,7 @@ - MulensModel.limbdarkeningcoeffs module — MulensModel 2.22.0 documentation + MulensModel.limbdarkeningcoeffs module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.magnificationcurve.html b/docs/MulensModel.magnificationcurve.html index 5a3aee14..512dd8fe 100644 --- a/docs/MulensModel.magnificationcurve.html +++ b/docs/MulensModel.magnificationcurve.html @@ -6,7 +6,7 @@ - MulensModel.magnificationcurve module — MulensModel 2.22.0 documentation + MulensModel.magnificationcurve module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.model.html b/docs/MulensModel.model.html index 7a8d8a73..08d5d78c 100644 --- a/docs/MulensModel.model.html +++ b/docs/MulensModel.model.html @@ -6,7 +6,7 @@ - MulensModel.model module — MulensModel 2.22.0 documentation + MulensModel.model module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.modelparameters.html b/docs/MulensModel.modelparameters.html index bd698f94..20ec36c5 100644 --- a/docs/MulensModel.modelparameters.html +++ b/docs/MulensModel.modelparameters.html @@ -6,7 +6,7 @@ - MulensModel.modelparameters module — MulensModel 2.22.0 documentation + MulensModel.modelparameters module — MulensModel 2.23.0 documentation @@ -334,7 +334,7 @@
property xi_omega_periapsis

float

-

The argument of periapsis of the xallrap orbit, i.e., the angle +

The argument of periapsis of the xallarap orbit, i.e., the angle between the ascending node and periapsis measured in the direction of motion. The units are degrees.

@@ -348,6 +348,26 @@ If not provided, then it defaults to t_0.

+
+
+property q_source
+

float

+

The mass ratio of the second and the first source. +This is value must be positive and can be > 1. +Defined only for xallarap binary-source models because it does not +affect the magnification for binary-source models without xallarap.

+
+ +
+
+property xallarap_reference_position
+

np.ndarray of shape (2, 1)

+

The position of the first source at t_0_xi relative to +the source center of mass. It is a 2D vector that is subtracted from +the source position along the orbit in order to calculate the shift +caused by xallarap.

+
+
property t_0_1
diff --git a/docs/MulensModel.mulensdata.html b/docs/MulensModel.mulensdata.html index b170038e..dc00c5be 100644 --- a/docs/MulensModel.mulensdata.html +++ b/docs/MulensModel.mulensdata.html @@ -6,7 +6,7 @@ - MulensModel.mulensdata module — MulensModel 2.22.0 documentation + MulensModel.mulensdata module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.mulensobjects.html b/docs/MulensModel.mulensobjects.html index 51dc65e8..ded2a1b0 100644 --- a/docs/MulensModel.mulensobjects.html +++ b/docs/MulensModel.mulensobjects.html @@ -6,7 +6,7 @@ - MulensModel.mulensobjects package — MulensModel 2.22.0 documentation + MulensModel.mulensobjects package — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.mulensobjects.lens.html b/docs/MulensModel.mulensobjects.lens.html index 7d896e44..d6ac72c7 100644 --- a/docs/MulensModel.mulensobjects.lens.html +++ b/docs/MulensModel.mulensobjects.lens.html @@ -6,7 +6,7 @@ - MulensModel.mulensobjects.lens module — MulensModel 2.22.0 documentation + MulensModel.mulensobjects.lens module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.mulensobjects.mulenssystem.html b/docs/MulensModel.mulensobjects.mulenssystem.html index f406df44..e3cc7ac2 100644 --- a/docs/MulensModel.mulensobjects.mulenssystem.html +++ b/docs/MulensModel.mulensobjects.mulenssystem.html @@ -6,7 +6,7 @@ - MulensModel.mulensobjects.mulenssystem module — MulensModel 2.22.0 documentation + MulensModel.mulensobjects.mulenssystem module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.mulensobjects.source.html b/docs/MulensModel.mulensobjects.source.html index 37418951..8e8daa78 100644 --- a/docs/MulensModel.mulensobjects.source.html +++ b/docs/MulensModel.mulensobjects.source.html @@ -6,7 +6,7 @@ - MulensModel.mulensobjects.source module — MulensModel 2.22.0 documentation + MulensModel.mulensobjects.source module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.pointlens.html b/docs/MulensModel.pointlens.html index 0f246db0..288408ab 100644 --- a/docs/MulensModel.pointlens.html +++ b/docs/MulensModel.pointlens.html @@ -6,7 +6,7 @@ - MulensModel.pointlens module — MulensModel 2.22.0 documentation + MulensModel.pointlens module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.pointlenswithshear.html b/docs/MulensModel.pointlenswithshear.html index e1edf389..8589d71e 100644 --- a/docs/MulensModel.pointlenswithshear.html +++ b/docs/MulensModel.pointlenswithshear.html @@ -6,7 +6,7 @@ - MulensModel.pointlenswithshear module — MulensModel 2.22.0 documentation + MulensModel.pointlenswithshear module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.satelliteskycoord.html b/docs/MulensModel.satelliteskycoord.html index 13635e5e..af842d39 100644 --- a/docs/MulensModel.satelliteskycoord.html +++ b/docs/MulensModel.satelliteskycoord.html @@ -6,7 +6,7 @@ - MulensModel.satelliteskycoord module — MulensModel 2.22.0 documentation + MulensModel.satelliteskycoord module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.trajectory.html b/docs/MulensModel.trajectory.html index 7740d762..558882ad 100644 --- a/docs/MulensModel.trajectory.html +++ b/docs/MulensModel.trajectory.html @@ -6,7 +6,7 @@ - MulensModel.trajectory module — MulensModel 2.22.0 documentation + MulensModel.trajectory module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.uniformcausticsampling.html b/docs/MulensModel.uniformcausticsampling.html index 41ef1945..6fd94515 100644 --- a/docs/MulensModel.uniformcausticsampling.html +++ b/docs/MulensModel.uniformcausticsampling.html @@ -6,7 +6,7 @@ - MulensModel.uniformcausticsampling module — MulensModel 2.22.0 documentation + MulensModel.uniformcausticsampling module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.utils.html b/docs/MulensModel.utils.html index 410d7ac1..2434e306 100644 --- a/docs/MulensModel.utils.html +++ b/docs/MulensModel.utils.html @@ -6,7 +6,7 @@ - MulensModel.utils module — MulensModel 2.22.0 documentation + MulensModel.utils module — MulensModel 2.23.0 documentation diff --git a/docs/MulensModel.version.html b/docs/MulensModel.version.html index 971e39c4..d3edc284 100644 --- a/docs/MulensModel.version.html +++ b/docs/MulensModel.version.html @@ -6,7 +6,7 @@ - MulensModel.version module — MulensModel 2.22.0 documentation + MulensModel.version module — MulensModel 2.23.0 documentation diff --git a/docs/_static/documentation_options.js b/docs/_static/documentation_options.js index e030fb03..5cdaab5b 100644 --- a/docs/_static/documentation_options.js +++ b/docs/_static/documentation_options.js @@ -1,6 +1,6 @@ var DOCUMENTATION_OPTIONS = { URL_ROOT: document.getElementById("documentation_options").getAttribute('data-url_root'), - VERSION: '2.22.0', + VERSION: '2.23.0', LANGUAGE: 'None', COLLAPSE_INDEX: false, BUILDER: 'html', diff --git a/docs/genindex.html b/docs/genindex.html index 8e93d247..eae67e2d 100644 --- a/docs/genindex.html +++ b/docs/genindex.html @@ -5,7 +5,7 @@ - Index — MulensModel 2.22.0 documentation + Index — MulensModel 2.23.0 documentation @@ -1035,6 +1035,10 @@

Q

+

R

@@ -1285,6 +1289,8 @@

X

  • x_caustic_in (MulensModel.modelparameters.ModelParameters property)
  • x_caustic_out (MulensModel.modelparameters.ModelParameters property) +
    • +
  • xallarap_reference_position (MulensModel.modelparameters.ModelParameters property)
  • xi_argument_of_latitude_reference (MulensModel.modelparameters.ModelParameters property)
    • diff --git a/docs/index.html b/docs/index.html index 792e69aa..2367d513 100644 --- a/docs/index.html +++ b/docs/index.html @@ -6,7 +6,7 @@ - Welcome to MulensModel’s documentation! — MulensModel 2.22.0 documentation + Welcome to MulensModel’s documentation! — MulensModel 2.23.0 documentation diff --git a/docs/install.html b/docs/install.html index bcae2ffe..0c53e680 100644 --- a/docs/install.html +++ b/docs/install.html @@ -6,7 +6,7 @@ - How to install? — MulensModel 2.22.0 documentation + How to install? — MulensModel 2.23.0 documentation diff --git a/docs/modules.html b/docs/modules.html index 684fa0f7..73189318 100644 --- a/docs/modules.html +++ b/docs/modules.html @@ -6,7 +6,7 @@ - MulensModel — MulensModel 2.22.0 documentation + MulensModel — MulensModel 2.23.0 documentation diff --git a/docs/objects.inv b/docs/objects.inv index 8a4b7b4d..f1f8d776 100644 Binary files a/docs/objects.inv and b/docs/objects.inv differ diff --git a/docs/py-modindex.html b/docs/py-modindex.html index 19a2232a..904999c1 100644 --- a/docs/py-modindex.html +++ b/docs/py-modindex.html @@ -5,7 +5,7 @@ - Python Module Index — MulensModel 2.22.0 documentation + Python Module Index — MulensModel 2.23.0 documentation diff --git a/docs/search.html b/docs/search.html index fbc5bd8f..c5fae9af 100644 --- a/docs/search.html +++ b/docs/search.html @@ -5,7 +5,7 @@ - Search — MulensModel 2.22.0 documentation + Search — MulensModel 2.23.0 documentation diff --git a/docs/searchindex.js b/docs/searchindex.js index b1b2d2e8..3c3cec70 100644 --- a/docs/searchindex.js +++ b/docs/searchindex.js @@ -1 +1 @@ -Search.setIndex({docnames:["MulensModel","MulensModel.binarylens","MulensModel.binarylensimports","MulensModel.binarylenswithshear","MulensModel.caustics","MulensModel.causticspointwithshear","MulensModel.causticswithshear","MulensModel.coordinates","MulensModel.event","MulensModel.fitdata","MulensModel.horizons","MulensModel.limbdarkeningcoeffs","MulensModel.magnificationcurve","MulensModel.model","MulensModel.modelparameters","MulensModel.mulensdata","MulensModel.mulensobjects","MulensModel.mulensobjects.lens","MulensModel.mulensobjects.mulenssystem","MulensModel.mulensobjects.source","MulensModel.pointlens","MulensModel.pointlenswithshear","MulensModel.satelliteskycoord","MulensModel.trajectory","MulensModel.uniformcausticsampling","MulensModel.utils","MulensModel.version","index","install","modules","tutorial","tutorial_fit_pi_E","tutorial_fit_pspl"],envversion:{"sphinx.domains.c":2,"sphinx.domains.changeset":1,"sphinx.domains.citation":1,"sphinx.domains.cpp":5,"sphinx.domains.index":1,"sphinx.domains.javascript":2,"sphinx.domains.math":2,"sphinx.domains.python":3,"sphinx.domains.rst":2,"sphinx.domains.std":2,"sphinx.ext.intersphinx":1,sphinx:56},filenames:["MulensModel.rst","MulensModel.binarylens.rst","MulensModel.binarylensimports.rst","MulensModel.binarylenswithshear.rst","MulensModel.caustics.rst","MulensModel.causticspointwithshear.rst","MulensModel.causticswithshear.rst","MulensModel.coordinates.rst","MulensModel.event.rst","MulensModel.fitdata.rst","MulensModel.horizons.rst","MulensModel.limbdarkeningcoeffs.rst","MulensModel.magnificationcurve.rst","MulensModel.model.rst","MulensModel.modelparameters.rst","MulensModel.mulensdata.rst","MulensModel.mulensobjects.rst","MulensModel.mulensobjects.lens.rst","MulensModel.mulensobjects.mulenssystem.rst","MulensModel.mulensobjects.source.rst","MulensModel.pointlens.rst","MulensModel.pointlenswithshear.rst","MulensModel.satelliteskycoord.rst","MulensModel.trajectory.rst","MulensModel.uniformcausticsampling.rst","MulensModel.utils.rst","MulensModel.version.rst","index.rst","install.rst","modules.rst","tutorial.rst","tutorial_fit_pi_E.rst","tutorial_fit_pspl.rst"],objects:{"":[[0,0,0,"-","MulensModel"]],"MulensModel.BinaryLens":[[0,2,1,"","adaptive_contouring_magnification"],[0,2,1,"","hexadecapole_magnification"],[0,2,1,"","point_source_magnification"],[0,2,1,"","vbbl_magnification"]],"MulensModel.binarylens":[[1,1,1,"","BinaryLens"]],"MulensModel.binarylens.BinaryLens":[[1,2,1,"","adaptive_contouring_magnification"],[1,2,1,"","hexadecapole_magnification"],[1,2,1,"","point_source_magnification"],[1,2,1,"","vbbl_magnification"]],"MulensModel.binarylenswithshear":[[3,1,1,"","BinaryLensWithShear"]],"MulensModel.binarylenswithshear.BinaryLensWithShear":[[3,2,1,"","point_source_magnification"]],"MulensModel.caustics":[[4,1,1,"","Caustics"]],"MulensModel.caustics.Caustics":[[4,1,1,"","CriticalCurve"],[4,3,1,"","critical_curve"],[4,2,1,"","get_caustics"],[4,2,1,"","plot"]],"MulensModel.causticspointwithshear":[[5,1,1,"","CausticsPointWithShear"]],"MulensModel.causticswithshear":[[6,1,1,"","CausticsWithShear"]],"MulensModel.coordinates":[[7,1,1,"","Coordinates"]],"MulensModel.coordinates.Coordinates":[[7,3,1,"","east_projected"],[7,3,1,"","ecliptic_lat"],[7,3,1,"","ecliptic_lon"],[7,3,1,"","galactic_b"],[7,3,1,"","galactic_l"],[7,3,1,"","north_projected"],[7,2,1,"","v_Earth_projected"]],"MulensModel.event":[[8,1,1,"","Event"]],"MulensModel.event.Event":[[8,3,1,"","best_chi2"],[8,3,1,"","best_chi2_parameters"],[8,3,1,"","blend_fluxes"],[8,2,1,"","calculate_chi2_gradient"],[8,3,1,"","chi2"],[8,3,1,"","chi2_gradient"],[8,3,1,"","coords"],[8,3,1,"","data_ref"],[8,3,1,"","datasets"],[8,2,1,"","fit_fluxes"],[8,3,1,"","fits"],[8,3,1,"","fluxes"],[8,2,1,"","get_chi2"],[8,2,1,"","get_chi2_for_dataset"],[8,2,1,"","get_chi2_gradient"],[8,2,1,"","get_chi2_per_point"],[8,2,1,"","get_flux_for_dataset"],[8,2,1,"","get_ref_fluxes"],[8,3,1,"","model"],[8,2,1,"","plot"],[8,2,1,"","plot_data"],[8,2,1,"","plot_model"],[8,2,1,"","plot_residuals"],[8,2,1,"","plot_source_for_datasets"],[8,2,1,"","plot_trajectory"],[8,2,1,"","reset_best_chi2"],[8,3,1,"","source_fluxes"],[8,3,1,"","sum_function"]],"MulensModel.fitdata":[[9,1,1,"","FitData"]],"MulensModel.fitdata.FitData":[[9,1,1,"","FSPL_Derivatives"],[9,3,1,"","blend_flux"],[9,2,1,"","calculate_chi2_gradient"],[9,3,1,"","chi2"],[9,3,1,"","chi2_gradient"],[9,3,1,"","chi2_per_point"],[9,3,1,"","data_magnification"],[9,3,1,"","dataset"],[9,2,1,"","fit_fluxes"],[9,3,1,"","gamma"],[9,2,1,"","get_chi2_gradient"],[9,2,1,"","get_d_A_d_params_for_point_lens_model"],[9,2,1,"","get_d_A_d_rho"],[9,2,1,"","get_d_A_d_u_for_PSPL_model"],[9,2,1,"","get_d_A_d_u_for_point_lens_model"],[9,2,1,"","get_data_magnification"],[9,2,1,"","get_dataset_trajectory"],[9,2,1,"","get_model_fluxes"],[9,2,1,"","get_model_magnitudes"],[9,2,1,"","get_residuals"],[9,3,1,"","model"],[9,2,1,"","scale_fluxes"],[9,3,1,"","source_flux"],[9,3,1,"","source_flux_ratio"],[9,3,1,"","source_fluxes"],[9,2,1,"","update"]],"MulensModel.fitdata.FitData.FSPL_Derivatives":[[9,2,1,"","get_d_A_d_rho"],[9,2,1,"","get_gradient"]],"MulensModel.horizons":[[10,1,1,"","Horizons"]],"MulensModel.horizons.Horizons":[[10,3,1,"","time"],[10,3,1,"","xyz"]],"MulensModel.limbdarkeningcoeffs":[[11,1,1,"","LimbDarkeningCoeffs"]],"MulensModel.limbdarkeningcoeffs.LimbDarkeningCoeffs":[[11,2,1,"","get_limb_coeff_gamma"],[11,2,1,"","get_limb_coeff_u"],[11,2,1,"","get_weighted_limb_coeff_gamma"],[11,2,1,"","set_limb_coeff_gamma"],[11,2,1,"","set_limb_coeff_u"]],"MulensModel.magnificationcurve":[[12,1,1,"","MagnificationCurve"]],"MulensModel.magnificationcurve.MagnificationCurve":[[12,2,1,"","get_binary_lens_magnification"],[12,2,1,"","get_magnification"],[12,2,1,"","get_point_lens_magnification"],[12,3,1,"","methods_for_epochs"],[12,2,1,"","set_magnification_methods"],[12,2,1,"","set_magnification_methods_parameters"]],"MulensModel.model":[[13,1,1,"","Model"]],"MulensModel.model.Model":[[13,3,1,"","bandpasses"],[13,3,1,"","caustics"],[13,3,1,"","coords"],[13,3,1,"","data_magnification"],[13,3,1,"","datasets"],[13,3,1,"","default_magnification_method"],[13,2,1,"","get_data_magnification"],[13,2,1,"","get_lc"],[13,2,1,"","get_limb_coeff_gamma"],[13,2,1,"","get_limb_coeff_u"],[13,2,1,"","get_magnification"],[13,2,1,"","get_magnification_methods"],[13,2,1,"","get_magnification_methods_parameters"],[13,2,1,"","get_parallax"],[13,2,1,"","get_ref_fluxes"],[13,2,1,"","get_residuals"],[13,2,1,"","get_satellite_coords"],[13,2,1,"","get_trajectory"],[13,2,1,"","is_static"],[13,2,1,"","magnification"],[13,3,1,"","methods"],[13,3,1,"","n_lenses"],[13,3,1,"","n_sources"],[13,2,1,"","parallax"],[13,3,1,"","parameters"],[13,2,1,"","plot_caustics"],[13,2,1,"","plot_data"],[13,2,1,"","plot_lc"],[13,2,1,"","plot_magnification"],[13,2,1,"","plot_residuals"],[13,2,1,"","plot_source"],[13,2,1,"","plot_source_for_datasets"],[13,2,1,"","plot_trajectory"],[13,2,1,"","reset_plot_properties"],[13,2,1,"","set_datasets"],[13,2,1,"","set_default_magnification_method"],[13,2,1,"","set_limb_coeff_gamma"],[13,2,1,"","set_limb_coeff_u"],[13,2,1,"","set_magnification_methods"],[13,2,1,"","set_magnification_methods_parameters"],[13,2,1,"","set_source_flux_ratio"],[13,2,1,"","set_source_flux_ratio_for_band"],[13,2,1,"","set_times"],[13,2,1,"","update_caustics"]],"MulensModel.modelparameters":[[14,1,1,"","ModelParameters"],[14,4,1,"","which_parameters"]],"MulensModel.modelparameters.ModelParameters":[[14,3,1,"","alpha"],[14,2,1,"","as_dict"],[14,3,1,"","convergence_K"],[14,3,1,"","dalpha_dt"],[14,3,1,"","ds_dt"],[14,3,1,"","gamma"],[14,3,1,"","gamma_parallel"],[14,3,1,"","gamma_perp"],[14,2,1,"","get_alpha"],[14,2,1,"","get_s"],[14,3,1,"","is_external_mass_sheet"],[14,3,1,"","is_external_mass_sheet_with_shear"],[14,2,1,"","is_finite_source"],[14,2,1,"","is_static"],[14,3,1,"","is_xallarap"],[14,3,1,"","n_lenses"],[14,3,1,"","n_sources"],[14,3,1,"","pi_E"],[14,3,1,"","pi_E_E"],[14,3,1,"","pi_E_N"],[14,3,1,"","pi_E_mag"],[14,3,1,"","q"],[14,3,1,"","rho"],[14,3,1,"","rho_1"],[14,3,1,"","rho_2"],[14,3,1,"","s"],[14,3,1,"","shear_G"],[14,3,1,"","source_1_parameters"],[14,3,1,"","source_2_parameters"],[14,3,1,"","t_0"],[14,3,1,"","t_0_1"],[14,3,1,"","t_0_2"],[14,3,1,"","t_0_kep"],[14,3,1,"","t_0_par"],[14,3,1,"","t_0_xi"],[14,3,1,"","t_E"],[14,3,1,"","t_caustic_in"],[14,3,1,"","t_caustic_out"],[14,3,1,"","t_eff"],[14,3,1,"","t_star"],[14,3,1,"","t_star_1"],[14,3,1,"","t_star_2"],[14,3,1,"","u_0"],[14,3,1,"","u_0_1"],[14,3,1,"","u_0_2"],[14,3,1,"","uniform_caustic_sampling"],[14,3,1,"","x_caustic_in"],[14,3,1,"","x_caustic_out"],[14,3,1,"","xi_Omega_node"],[14,3,1,"","xi_argument_of_latitude_reference"],[14,3,1,"","xi_eccentricity"],[14,3,1,"","xi_inclination"],[14,3,1,"","xi_omega_periapsis"],[14,3,1,"","xi_period"],[14,3,1,"","xi_semimajor_axis"]],"MulensModel.mulensdata":[[15,1,1,"","MulensData"]],"MulensModel.mulensdata.MulensData":[[15,3,1,"","bad"],[15,3,1,"","bandpass"],[15,3,1,"","chi2_fmt"],[15,3,1,"","coords"],[15,2,1,"","copy"],[15,2,1,"","data_and_err_in_chi2_fmt"],[15,2,1,"","data_and_err_in_input_fmt"],[15,3,1,"","ephemerides_file"],[15,3,1,"","err_flux"],[15,3,1,"","err_mag"],[15,3,1,"","errorbars_scale_factors"],[15,3,1,"","errorbars_scaling_equation"],[15,3,1,"","flux"],[15,3,1,"","good"],[15,3,1,"","input_fmt"],[15,3,1,"","mag"],[15,3,1,"","n_epochs"],[15,2,1,"","plot"],[15,3,1,"","plot_properties"],[15,3,1,"","satellite_skycoord"],[15,2,1,"","scale_errorbars"],[15,2,1,"","set_limb_darkening_weights"],[15,3,1,"","time"]],"MulensModel.mulensobjects":[[16,1,1,"","Lens"],[16,1,1,"","MulensSystem"],[16,1,1,"","Source"],[17,0,0,"-","lens"],[18,0,0,"-","mulenssystem"],[19,0,0,"-","source"]],"MulensModel.mulensobjects.Lens":[[16,3,1,"","a_proj"],[16,3,1,"","caustics"],[16,3,1,"","distance"],[16,3,1,"","epsilon"],[16,3,1,"","mass"],[16,3,1,"","mass_1"],[16,3,1,"","mass_2"],[16,3,1,"","mass_3"],[16,3,1,"","n_masses"],[16,3,1,"","pi_L"],[16,2,1,"","plot_caustics"],[16,3,1,"","q"],[16,3,1,"","s"],[16,3,1,"","total_mass"]],"MulensModel.mulensobjects.MulensSystem":[[16,3,1,"","lens"],[16,3,1,"","mu_rel"],[16,3,1,"","pi_E"],[16,3,1,"","pi_rel"],[16,2,1,"","plot_caustics"],[16,2,1,"","plot_magnification"],[16,3,1,"","r_E"],[16,3,1,"","r_E_tilde"],[16,3,1,"","source"],[16,3,1,"","t_E"],[16,3,1,"","theta_E"]],"MulensModel.mulensobjects.Source":[[16,3,1,"","angular_radius"],[16,3,1,"","distance"],[16,3,1,"","pi_S"]],"MulensModel.mulensobjects.lens":[[17,1,1,"","Lens"]],"MulensModel.mulensobjects.lens.Lens":[[17,3,1,"","a_proj"],[17,3,1,"","caustics"],[17,3,1,"","distance"],[17,3,1,"","epsilon"],[17,3,1,"","mass"],[17,3,1,"","mass_1"],[17,3,1,"","mass_2"],[17,3,1,"","mass_3"],[17,3,1,"","n_masses"],[17,3,1,"","pi_L"],[17,2,1,"","plot_caustics"],[17,3,1,"","q"],[17,3,1,"","s"],[17,3,1,"","total_mass"]],"MulensModel.mulensobjects.mulenssystem":[[18,1,1,"","MulensSystem"]],"MulensModel.mulensobjects.mulenssystem.MulensSystem":[[18,3,1,"","lens"],[18,3,1,"","mu_rel"],[18,3,1,"","pi_E"],[18,3,1,"","pi_rel"],[18,2,1,"","plot_caustics"],[18,2,1,"","plot_magnification"],[18,3,1,"","r_E"],[18,3,1,"","r_E_tilde"],[18,3,1,"","source"],[18,3,1,"","t_E"],[18,3,1,"","theta_E"]],"MulensModel.mulensobjects.source":[[19,1,1,"","Source"]],"MulensModel.mulensobjects.source.Source":[[19,3,1,"","angular_radius"],[19,3,1,"","distance"],[19,3,1,"","pi_S"]],"MulensModel.pointlens":[[20,1,1,"","PointLens"],[20,4,1,"","get_pspl_magnification"]],"MulensModel.pointlens.PointLens":[[20,2,1,"","get_point_lens_LD_integrated_magnification"],[20,2,1,"","get_point_lens_finite_source_magnification"],[20,2,1,"","get_point_lens_large_LD_integrated_magnification"],[20,2,1,"","get_point_lens_large_finite_source_magnification"],[20,2,1,"","get_point_lens_limb_darkening_magnification"],[20,2,1,"","get_point_lens_uniform_integrated_magnification"]],"MulensModel.pointlenswithshear":[[21,1,1,"","PointLensWithShear"]],"MulensModel.pointlenswithshear.PointLensWithShear":[[21,2,1,"","get_point_lens_LD_integrated_magnification"],[21,2,1,"","get_point_lens_finite_source_magnification"],[21,2,1,"","get_point_lens_large_LD_integrated_magnification"],[21,2,1,"","get_point_lens_large_finite_source_magnification"],[21,2,1,"","get_point_lens_limb_darkening_magnification"],[21,2,1,"","get_point_lens_uniform_integrated_magnification"],[21,2,1,"","get_point_source_magnification"]],"MulensModel.satelliteskycoord":[[22,1,1,"","SatelliteSkyCoord"]],"MulensModel.satelliteskycoord.SatelliteSkyCoord":[[22,2,1,"","get_satellite_coords"]],"MulensModel.trajectory":[[23,1,1,"","Trajectory"]],"MulensModel.trajectory.Trajectory":[[23,3,1,"","d_perp"],[23,2,1,"","get_xy"],[23,3,1,"","parallax_delta_N_E"],[23,3,1,"","x"],[23,3,1,"","y"]],"MulensModel.uniformcausticsampling":[[24,1,1,"","UniformCausticSampling"]],"MulensModel.uniformcausticsampling.UniformCausticSampling":[[24,2,1,"","caustic_point"],[24,2,1,"","check_valid_trajectory"],[24,2,1,"","get_standard_parameters"],[24,2,1,"","get_uniform_sampling"],[24,2,1,"","get_x_in_x_out"],[24,2,1,"","jacobian"],[24,3,1,"","n_caustics"],[24,3,1,"","q"],[24,3,1,"","s"],[24,2,1,"","which_caustic"]],"MulensModel.utils":[[25,1,1,"","PlotUtils"],[25,1,1,"","Utils"]],"MulensModel.utils.PlotUtils":[[25,2,1,"","find_subtract"],[25,2,1,"","find_subtract_xlabel"],[25,2,1,"","get_color_differences"],[25,2,1,"","get_y_value_y_err"]],"MulensModel.utils.Utils":[[25,2,1,"","astropy_version_check"],[25,2,1,"","complex_fsum"],[25,2,1,"","date_change"],[25,2,1,"","dot"],[25,2,1,"","gamma_to_u"],[25,2,1,"","get_flux_and_err_from_mag"],[25,2,1,"","get_flux_from_mag"],[25,2,1,"","get_mag_and_err_from_flux"],[25,2,1,"","get_mag_from_flux"],[25,2,1,"","get_n_caustics"],[25,2,1,"","u_to_gamma"],[25,2,1,"","vector_product_normalized"],[25,2,1,"","velocity_of_Earth"]],MulensModel:[[0,1,1,"","BinaryLens"],[1,0,0,"-","binarylens"],[2,0,0,"-","binarylensimports"],[3,0,0,"-","binarylenswithshear"],[4,0,0,"-","caustics"],[5,0,0,"-","causticspointwithshear"],[6,0,0,"-","causticswithshear"],[7,0,0,"-","coordinates"],[8,0,0,"-","event"],[9,0,0,"-","fitdata"],[10,0,0,"-","horizons"],[11,0,0,"-","limbdarkeningcoeffs"],[12,0,0,"-","magnificationcurve"],[13,0,0,"-","model"],[14,0,0,"-","modelparameters"],[15,0,0,"-","mulensdata"],[16,0,0,"-","mulensobjects"],[20,0,0,"-","pointlens"],[21,0,0,"-","pointlenswithshear"],[22,0,0,"-","satelliteskycoord"],[23,0,0,"-","trajectory"],[24,0,0,"-","uniformcausticsampling"],[25,0,0,"-","utils"],[26,0,0,"-","version"]]},objnames:{"0":["py","module","Python module"],"1":["py","class","Python class"],"2":["py","method","Python method"],"3":["py","property","Python property"],"4":["py","function","Python function"]},objtypes:{"0":"py:module","1":"py:class","2":"py:method","3":"py:property","4":"py:function"},terms:{"0":[0,1,3,6,8,9,12,13,14,16,17,20,23,24,25,28,30,31,32],"00":[7,8,25],"000":7,"00096":30,"001":[0,1,31],"0039":30,"005":[12,13],"01":31,"04":31,"06":30,"086":31,"0j":6,"1":[0,1,3,4,6,9,11,12,13,14,15,16,17,20,23,24,25,28,30,31,32],"10":[8,12,24,25,31],"100":[12,31],"1000":13,"10000":24,"102":31,"1025":31,"11":31,"110":31,"1103":31,"12":25,"120":30,"120755":31,"133":30,"139":[12,20],"15":[12,13,31],"150":31,"1593":[0,1,12,25],"16":[30,31],"1679":[0,1,12],"17":28,"18":[7,8,31],"180":[7,14,23],"18h00m00":7,"19":30,"1979":21,"1994":[12,20],"2":[4,6,8,9,11,13,14,15,16,17,24,25,30,31,32],"20":32,"200":[12,20,24],"2002":11,"2003":[20,30],"2004":[12,20,30],"2005":31,"2007":[0,1,12],"2008":[0,1,12,14,24,25],"2009":[12,20,24],"2010":[0,1,12,24],"2011":[14,23],"2012":24,"2015":[25,31],"2018":[12,20],"2022":3,"21":31,"2188":[0,1,12],"218982":31,"22":[13,25],"2228":24,"223":30,"23":24,"235":30,"24":3,"245":25,"2450000":[8,13,14,15,25,31],"2452810":30,"2452833":30,"2452845":30,"2452848":30,"2452890":30,"2453200":31,"2453628":31,"2453630":31,"2453950":31,"2455746":[12,13],"2455747":[12,13],"2455748":[12,13],"2456789":13,"246":25,"2460000":[8,13,15,25],"25":30,"26":31,"262":20,"27":24,"270":7,"2719":31,"2810":30,"282":21,"2890":30,"29062":31,"2972":31,"2f":32,"3":[9,11,13,14,22,24,25,31],"30":[7,8,25,31],"30d00m00":7,"3291":31,"3400":31,"35":14,"35507":31,"360":7,"37":31,"37263":31,"377":[0,1,12],"3800":31,"387105":31,"395":24,"3d":25,"3f":32,"4":[0,1,12,24],"40":31,"400":31,"408":[0,1,12],"41":31,"415677":31,"421l":[12,20],"426":24,"430":[12,20],"44":31,"45":31,"479":[12,20],"488817":31,"491":24,"4f":32,"5":[11,13,15,20,28,30,31,32],"50":31,"500":31,"5000":[4,13,16,17,18],"505":[12,20],"515":24,"5157":[12,20],"52":24,"521":11,"53":30,"5380":32,"561":21,"572":11,"587":24,"59":31,"5f":[31,32],"5th":[0,1],"6":[12,13,24,30,31],"603":[12,20],"61":30,"611301":31,"6700":31,"67778":31,"681":[0,1,12,25],"695":[12,20],"7":[12,13,24,31],"71":[12,20,31],"738":14,"787":24,"8":[30,32],"84":31,"87":14,"927":3,"93":31,"947":31,"949":31,"95":31,"9th":3,"boolean":[0,1,8,9,13,14,15,20,23],"case":[7,11,13,14,15,16,17,24,31,32],"class":[0,1,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,27,30,32],"default":[0,1,4,7,8,9,12,13,14,15,23,24,25,30],"do":[0,1,3,8,13,14,24,30,31],"export":28,"final":[0,1,32],"float":[0,1,3,4,5,6,7,8,9,11,12,13,14,15,16,17,18,19,20,22,23,24,25,32],"function":[0,1,8,9,12,13,14,15,16,17,23,24,25,27,31,32],"import":[7,8,10,13,24,30],"int":[4,8,13,14,15,16,17,18,24,25],"long":[15,31],"public":20,"return":[0,1,3,4,7,8,9,10,11,12,13,14,15,20,21,22,24,25,27,31,32],"static":[8,14,25],"true":[0,1,3,8,9,12,13,14,15,23,24,30,31],"try":[31,32],"while":15,A:[7,9,11,13,14,15,16,17,18,20,23,24,25,32],And:[30,31],As:[0,1,31],At:31,Be:20,But:31,For:[0,1,8,9,10,12,13,14,20,23,24,27,30],IN:[],If:[0,1,3,7,8,9,12,13,14,15,16,17,18,19,20,25,28,30,31,32],In:[7,8,13,14,24,28,31,32],Is:25,It:[0,1,8,11,13,14,16,18,23,30,31],Its:[],NOT:30,No:9,Not:[13,16,17,21],ON:13,ONE:13,Or:30,That:27,The:[0,1,3,4,8,9,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,30,31,32],Then:[13,14,30,31,32],There:8,These:8,To:[4,9,13,14,24,28,32],Will:[8,15],__future__:30,__init__:15,_lens_:12,_not_:20,_point_:12,a_proj:[16,17],abil:32,about:[13,31,32],abov:[8,9,24,28,32],acc:[0,1],accept:[4,8,9,13,15],access:[8,14],accord:[14,24],account:[13,23],accur:[8,25],accuraci:[0,1,12,13],across:15,actual:[16,17],ad:[8,15,16,17,30],adapt:[0,1],adaptive_contour:12,adaptive_contouring_magnif:[0,1,12],adaptivecontour:[0,1,12,28],add:[0,1,15,28],add_2450000:[15,31],add_2460000:15,addit:[12,13,15],addition:31,adequ:13,advanc:[0,1,24],advis:30,advoc:14,affect:[0,1,13,23],after:[9,16,17],al:[3,11,12,14,20,23,24,25,30,31],algorithm:[0,1,24,25,31],all:[8,9,13,15,16,17,23,24,27,28,30,32],all_approxim:[0,1],allow:[8,9,12,13,14,32],almost:15,along:[4,13],alpha:[13,14,15,16,18,23,24,30,31],alpha_list:13,alreadi:8,also:[0,1,8,9,13,15,16,18,20,24,30,31],altern:[13,24],alwai:[4,6,14,15,16,18],amplitud:32,an:[7,8,9,10,11,13,14,15,16,17,22,24,25,27,31,32],analysi:24,analyz:13,angl:[7,13,14,24],angular:[16,18,19],angular_radiu:[16,19],ani:[13,25],annual:31,annuli:[12,20],anomali:14,anyth:9,apj:[0,1,3,11,12,14,20,25],appendix:23,appli:[9,13,24,32],applic:[0,1,8,9,15],approach:[20,31],appropri:[8,13,15],approxim:[0,1,12,20,32],ar:[0,1,3,4,7,8,9,12,13,14,15,16,17,18,19,20,23,24,25,27,28,30,31,32],arbitrari:[9,27],archiv:[28,30],area:[0,1],aren:15,arg:[7,13,14,21,31,32],argument:[8,9,10,12,13,14,16,17,18,20,21,23,24,25,32],around:[31,32],arrai:[7,8,9,16,17,20],arrow:13,arrow_kwarg:13,as_dict:14,ascend:14,associ:32,assum:[8,12,13,14,15,16,17,18,19,20,23],astropi:[7,8,10,12,13,14,15,16,17,18,19,22,23,25,28],astropy_version_check:25,attach:13,attr:[8,9],attribut:[0,1,3,4,5,6,8,12,13,14,16,19,20,22,23,32],au:[16,17,18,23],author:[15,28],avail:[12,30],axi:[0,1,3,4,8,13,14,15,24,25,31,32],b:[],back:31,background:[16,19],bad:[8,9,13,15],band:[8,11,13,15],bandpass:[9,11,13,15],barycentr:25,base:[0,1,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,32],baselin:32,bash:28,bashrc:28,basic:[8,14],bayesian:24,becaus:[8,15,24,31],been:8,befor:[0,1,3,16,17,28],begin:[14,25,30],behavior:[8,13],below:[24,31],best:[9,13,14,27,30,31,32],best_chi2:8,best_chi2_paramet:8,better:[31,32],between:[0,1,3,4,6,9,12,13,14,16,17,18,24,25,31],beyond:22,bfg:32,binari:[0,1,3,6,8,9,12,13,14,23,24,25],binarylen:[0,3,29],binarylensimport:[0,29],binarylenswithshear:[0,29],bjd:15,bjd_tdb:15,black:[13,25,30,31],blend:[8,9,13],blend_flux:[8,9,13],blend_flux_0:9,blg234:31,blg:[20,30,31],bodi:[0,1,3,4,6,14,16,17,18,30],bond:30,bool:[8,9,13,14,24,25],both:[0,1,3,4,9,11,13,24,25,31],bottom:24,bozza:[0,1,12,20],bright:[15,31],brighter:32,burn:31,button:28,c:[20,24,28,32],caclul:8,calcul:[0,1,3,4,7,8,9,12,13,14,15,20,21,22,23,24,25,31,32],calculate_chi2_gradi:[8,9],call:[8,9,13,14,15,24,32],can:[0,1,7,8,9,11,12,13,14,15,20,24,27,28,30,31,32],cannot:14,cartesian:25,cartesianrepresent:[10,25],cassan:[14,24],caustic:[0,5,6,8,13,14,16,17,18,24,25,29],caustic_point:24,causticspointwithshear:[0,29],causticswithshear:[0,29],center:[0,1,3,4,11,14,23],central:24,cg:32,chain:31,chang:[5,8,9,13,14,15,21,24,25,32],check:[14,15,24,25,28,31],check_valid_trajectori:24,chi2:[8,9,15,27,30,31],chi2_fmt:[8,15],chi2_for_model:32,chi2_gradi:[8,9],chi2_initi:32,chi2_min:31,chi2_min_theta:31,chi2_per_point:9,chi:[8,9,15,31,32],choic:32,choos:31,chosen:[0,1],circl:13,circular:[12,13],cite:[3,24],clariti:31,clockwis:24,close:24,cobyla:32,code:[12,15,24,25,28,30,31,32],coeff:13,coeffici:[0,1,9,11,12,13,15,16,19,20,25],color:[13,15,24,25,30,31],color_differ:25,color_list:[13,25],colorbar:24,column:15,combin:[8,13,14,16,17,31,32],come:31,command:[13,28],comment:30,common:31,companion:[16,17],compar:[12,13],complex128:24,complex:[0,1,3,5,6,14,24,25],complex_fsum:25,complic:24,compon:[7,9,14,16,17,23,24],compos:21,comput:[9,12,20,31],condit:24,configur:24,congratul:[28,32],conniv:7,consid:[12,21,24],consist:[15,16,18],constraint:[13,32],contact:[15,28],contain:[8,15,32],content:29,contour:[0,1,12],contrari:[8,13],contribut:[8,9,11],control:[8,13],conveni:31,convent:[0,1,11,12,13,14,15,23,24,25],converg:[3,5,6,12,14,21],convergence_k:[3,5,6,14],convert:32,coord:[8,12,13,15,23,31],coordin:[0,1,3,4,8,12,13,14,15,16,17,22,23,24,29,31],copi:15,correspond:[4,9,13,16,17,24,31,32],could:[15,30],counter:24,counterclockwis:14,cover:15,cpp:28,creat:9,critic:4,critical_curv:4,criticalcurv:4,cross:[14,16,18,24],cshrc:28,cubic:22,cunstructor:32,current:[8,9,13,20,28,31],curv:[4,8,12,13,16,18,20,27,31],curvelinear:[14,24],d:[8,9],d_a_d_rho:9,d_perp:23,da:9,da_dparam:9,da_du:9,dai:[14,16,18,32],dalpha_dt:14,darken:[0,1,8,9,11,12,13,15,16,19,20,25],dash:31,dat:[22,31,32],data:[8,9,13,15,22,25,27,31],data_and_err_in_chi2_fmt:15,data_and_err_in_input_fmt:15,data_err:15,data_list:15,data_magnif:[9,13],data_num:8,data_path:[30,31,32],data_ref:[8,13],datapoint:9,dataset:[8,9,13,15,30,31,32],date:[12,13,15,25],date_chang:25,dec:[7,10,13,15],decid:12,def:[31,32],default_magnification_method:13,default_method:12,defin:[0,1,3,4,7,8,9,13,14,15,16,17,23,24,28,32],definit:[16,17,23],deg:[7,14,23],degener:31,degre:[7,13,14,16,18,30],depend:[8,10,32],deprec:[8,9,13,15],deriv:[9,13],descript:[12,24,27],desir:9,detail:[8,13],determin:[0,1,9,13,30],dfm:31,dict:[8,9,11,12,13,14,15,23,24,31],dictionari:[11,12,13,14,15,23,24],didn:28,differ:[8,13,15,23,24,25,31,32],dimensionless:[16,18,23],direct:[13,14,20],directli:[12,13,14],directori:28,dispers:31,distanc:[10,16,17,19],distribut:24,doctr:[8,13],doe:[0,1,8,9,24,30,31],doesn:12,dominik:[0,1,12],don:[12,24,30,31],done:[15,31,32],dot:25,dot_product:25,doubl:[13,14,25],download:[28,30,31,32],dozen:32,draw:24,drawn:24,ds_dt:14,dt:13,du:9,due:13,dure:[0,1],e9967a:25,e:[4,5,7,8,9,10,11,12,13,14,15,20,21,23,24,25,28,32],each:[8,9,11,12,13,15,16,17,20,21,23,24],earth:[7,13,23,25,31],earth_coord:23,earth_orbit:[12,13,23],easi:[20,31],easili:32,east:[7,14,23],east_project:7,eccentr:14,eclipt:7,ecliptic_lat:7,ecliptic_lon:7,effect:[12,13,14,15,20,23,30,31],einstein:[0,1,3,4,6,14,16,17,18,20,24],either:[0,1,8,13,14,15,16,19,28],element:[9,13],ellips:13,els:32,emce:31,en:25,enclos:[0,1],encourag:32,enough:31,ensemblesampl:31,ensur:9,entranc:[14,24],enumer:32,ephemeri:15,ephemerid:[10,12,13,15,22,23],ephemerides_fil:[9,13,15,22],ephemeris_fil:22,epheri:22,epoch:[7,8,12,13,14,15,22,23,24,25,31],epsilon:[16,17],eq:24,equal:[4,13,14,16,18,24],equat:[0,1,3,12,15,20],err:15,err_flux:[9,15,25],err_mag:[15,25,32],error:[0,1,12,15],errorbar:[8,9,13,15,32],errorbars_scale_factor:15,errorbars_scaling_equ:15,es:[8,9,13],estim:[0,1,32],et:[3,11,12,14,20,23,24,25,30,31],etc:[0,1,3,8,14,15],evalu:[9,15,24,32],even:9,evenli:13,event:[0,1,7,12,13,15,23,24,27,29,30,31,32],event_0:31,event_1:31,event_2:31,everi:[11,13,28,31],everyth:30,ew:31,exactli:13,exampl:[8,13,14,15,22,24,30,31,32],except:[13,23,24],exclud:15,execut:[0,1,9],exist:[8,24],exit:[14,24],exoplanet:30,expand:30,expect:13,explan:[8,13],explicitli:[13,14,30],extend:15,extern:[3,14],extract:28,extrapol:22,f_blend:13,f_sourc:13,factor:15,faint:15,fals:[0,1,8,9,12,13,15,20,23,24,25],far:12,fast:20,favourit:32,featur:[13,30],few:32,fig:24,figur:31,file:[10,12,13,15,22,23,25,28,31,32],file_nam:[10,15,30,31,32],fill:13,filter:11,find:[13,25,27,31],find_subtract:25,find_subtract_xlabel:25,fine:[24,31,32],finer:13,finit:[0,1,9,12,14,20,30],finite_source_ld_lee09:12,finite_source_ld_wittmao94:12,finite_source_ld_yoo04:12,finite_source_ld_yoo04_direct:12,finite_source_uniform_gould94:12,finite_source_uniform_gould94_direct:12,finite_source_uniform_lee09:12,finite_source_uniform_wittmao94:12,first:[8,11,13,15,21,25,30,31,32],fit:[8,9,13,14,15,24,27,30],fit_blend:[8,13],fit_flux:[8,9],fitdata:[0,8,29],fix:[8,9,11,13,14],fix_blend:8,fix_blend_flux:[8,9],fix_source_flux:[8,9],fix_source_flux_ratio:[8,9],flag:[8,15],floor:15,flux:[8,9,11,13,15,25,30,32],flux_err:25,flux_ratio_constraint:13,flux_source_1:13,flux_source_2:13,flux_source_band_1:13,flux_source_band_2:13,fmt:[15,32],fold:24,follow:[23,25,30],form:[12,13,31],formal:20,format:[8,15,25,31,32],found:32,four:24,fpic:28,fraction:[0,1,3,4,6,14,16,17,18,20],frame:[10,25],free:9,frequent:[11,24],from:[7,8,9,10,12,13,14,15,20,22,23,24,25,28,30,31,32],fsbl:14,fspl:[9,14],fspl_deriv:9,fsplderiv:9,fsum:8,full:[25,31],full_bjd:[7,25],fulli:28,fun:32,func:9,further:31,g:[7,8,9,11,12,13,14,15,20,23,25,28,32],galact:7,galactic_b:7,galactic_l:7,gamma:[0,1,9,11,12,13,14,15,20,25],gamma_parallel:14,gamma_perp:14,gamma_to_u:25,gaussian:15,gave:8,gca:[4,13,32],gcc:28,gener:[12,23,25,27],get:[8,11,13,14,15,24,25,30,32],get_:14,get_alpha:14,get_binary_lens_magnif:12,get_caust:4,get_chi2:[8,30,31,32],get_chi2_for_dataset:8,get_chi2_gradi:[8,9],get_chi2_per_point:8,get_color_differ:25,get_d_a_d_params_for_point_lens_model:9,get_d_a_d_rho:9,get_d_a_d_u_for_point_lens_model:9,get_d_a_d_u_for_pspl_model:9,get_data_magnif:[9,13],get_dataset_trajectori:9,get_flux_and_err_from_mag:25,get_flux_for_dataset:8,get_flux_from_mag:25,get_gradi:9,get_lc:13,get_limb_coeff_gamma:[9,11,13],get_limb_coeff_u:[11,13],get_mag_and_err_from_flux:25,get_mag_from_flux:25,get_magnif:[12,13],get_magnification_method:13,get_magnification_methods_paramet:13,get_model_flux:9,get_model_magnitud:9,get_n_caust:25,get_parallax:13,get_point_lens_finite_source_magnif:[20,21],get_point_lens_large_finite_source_magnif:[20,21],get_point_lens_large_ld_integrated_magnif:[20,21],get_point_lens_ld_integrated_magnif:[20,21],get_point_lens_limb_darkening_magnif:[20,21],get_point_lens_magnif:12,get_point_lens_uniform_integrated_magnif:[20,21],get_point_source_magnif:21,get_pspl_magnif:20,get_ref_flux:[8,13],get_residu:[9,13],get_satellite_coord:[13,22],get_standard_paramet:24,get_trajectori:13,get_uniform_sampl:24,get_weighted_limb_coeff_gamma:11,get_x_in_x_out:24,get_xi:23,get_y_value_y_err:25,github:28,give:[4,8,9,11,13,14,15,22,32],given:[0,1,3,4,7,8,9,11,12,13,14,15,16,17,19,22,23,24,25,27,32],go:[13,24,28],goe:31,good:[9,15,30,32],got:31,gould:[0,1,12,20],gradient:[8,9],green:[0,1,28],ground:13,guarante:30,guess:32,h:20,ha:[8,9,11,13,14,16,18,23,24,31,32],had:30,hand:[0,1,3,24],handl:14,happen:13,have:[8,11,13,15,28,31,32],haven:[8,31],he:31,henc:[12,13,31,32],here:[13,15,20,24,31,32],hexadecapol:[0,1,12,13],hexadecapole_magnif:[0,1,12],high:13,higher:[0,1,3,23,31],hjd:[15,31],hold:10,home:28,hope:31,horizon:[0,22,29],hour:13,hourangl:7,how:[0,1,15,23,30,31,32],howev:24,http:[25,31],human:[8,13,15],i:[0,1,4,5,8,9,10,11,12,13,14,15,20,21,24,31,32],i_caust:24,id:31,ignor:[0,1,13],imag:[0,1,3,24],impact:[16,18,24,31,32],implement:[0,1,4,13,16,17,20,21],importantli:25,inclin:14,includ:[8,9,12,13,14,15,20,23,28,31,32],index:[8,13,27,31],index_dataset:8,indic:24,individu:9,inf:31,info:10,inform:[8,13,14,15,31,32],inherit:7,initi:32,initial_guess:32,input:[10,14,15,20,23,24,25,27],input_fmt:15,inspect:32,instal:[25,27],instanc:[7,8,11,13,14,15,23,32],instantan:[14,20],instead:[0,1,8,13,14,20,24,31],instruct:[10,15],integr:[0,1,12,20,24],intermedi:24,intern:[4,13,14,23,24],interpol:[12,15,20,22],interv:13,introduc:[11,15,25],invert:32,invert_yaxi:32,invok:28,io:31,is_external_mass_sheet:14,is_external_mass_sheet_with_shear:14,is_finite_sourc:14,is_stat:[13,14],is_xallarap:14,isfinit:31,isinst:32,isnan:31,iter:[12,25],its:[0,1,8,11,15],j:20,jacobian:[8,9,14,24],jd:32,join:[30,31,32],jpl:[10,22],julian:[12,13,25],jupyt:30,just:[8,13,22,28,32],k2_ephemeris_01:22,k:[8,9],kain:24,keep:[8,24],kei:[8,9,11,12,13,15,23,24,32],keyword:[4,7,8,9,13,15,16,17,18,20,21,22,24],km:[7,25],know:[0,1,3,30],kpc:[16,17,19],kwarg:[4,7,8,9,12,13,15,16,17,18,21],l:7,label:[15,31,32],label_list:13,larg:[8,12,20],larger:[4,13,15,20,24],last:[28,32],later:[20,25,31],latitud:[7,14],latter:[8,15,31],ld_accuraci:[0,1],learn:[31,32],least:[9,14],leav:30,lee:[12,20],left:[0,1,3,4],legend:[8,31,32],len:[0,1,3,4,5,6,8,9,12,13,14,16,18,20,21,23,24,25,30,31,32],length:[4,13,15,24,25],lens:[4,8,13,20,27],less:32,let:[31,32],librari:[0,1],light:[8,12,13,20,27,31],like:[13,24],likelihood:[27,31],limb:[0,1,8,9,11,12,13,15,16,19,20,25],limb_darken:[16,19],limbdarkeningcoeff:[0,13,15,16,19,20,29],limit:13,line:[0,1],linear:[0,1,9,11],linestyl:30,linewidth:13,link:[8,13],linspac:24,list:[4,8,9,12,13,14,15,16,17,22,23,24,25,31,32],lm:28,ln_like:31,ln_like_:31,ln_prior:31,ln_prior_:31,ln_prob:31,load:15,loadtxt:15,loc:[31,32],logarithm:31,longer:[8,13],longitud:[7,14],look:[13,31,32],lower:[0,1,3,24],ls:31,lstdc:28,lumin:[13,14],lw:31,m:[0,1],m_1:[16,17],m_i:[16,17],ma:[16,18,19],macho:20,mag:[8,9,13,15,25,32],mag_zeropoint:13,magnif:[0,1,3,9,12,13,16,18,20,21,32],magnification_vector:13,magnificationcurv:[0,13,29],magnitud:[8,9,13,14,15,25,32],mai:[0,1,3,4,7,8,13,14,15,16,17,19,31,32],main:[14,27,30],major:14,make:[0,1,3,8,13,15,28,31],mani:14,mao:[12,20],mark:[8,13,15],marker:[8,13,15],marker_list:13,markers:15,martin:[0,1,12],mass:[0,1,3,4,5,6,14,16,17,18,20,21,23,24,25],mass_1:[0,1,3,16,17],mass_2:[0,1,3,16,17],mass_3:[16,17],materi:32,math:8,matplotlib:[4,8,13,15,24,25,30,31,32],matter:15,mayb:13,mead:32,mean:[0,1,3,8,13,24,25],measur:[14,15],mention:[8,9],messag:32,method:[0,1,12,13,20,30,32],method_paramet:13,methods_for_epoch:12,methods_paramet:[12,13],microarcsecond:[16,19],microlens:[8,11,12,13,14,15,16,18,20,23,24,27,30],might:30,millarcsecond:[16,19],milliarcsecond:[16,17,18,19],mind:24,minim:32,minimum:[14,15,24,25],minut:22,miss:[16,17],mitig:[4,13],mm:[30,31,32],mnra:[0,1,12,20,24],moa:30,moa_data:30,mode:31,model:[0,8,9,12,14,15,16,18,20,21,23,24,27,29,31],model_0:31,model_1:31,model_2:31,model_flux:9,model_mag:9,model_paramet:14,modelparamet:[0,12,13,16,18,20,21,23,24,27,29],modifi:31,modul:[27,29,31,32],month:25,more:[8,9,11,13,15,16,17,18,24,31,32],most:[7,13,14,15,25,30],motion:[13,14,16,18,20,31],mu:27,mu_rel:[16,18],much:32,mulens_data:15,mulensdata:[0,8,9,12,13,27,29,30,31,32],mulensmodel:[28,30,31,32],mulensobject:[0,29],mulenssystem:[0,16],multi:14,multinest:31,multipl:[12,13,15,20],multipli:15,must:[8,9,13,14,15,16,18,24],my_1s2l_model:30,my_data:[8,31,32],my_ev:[30,31,32],my_friends_data:32,my_model:31,my_pspl_model:30,n:[23,24],n_burn:31,n_caustic:[24,25],n_dim:31,n_epoch:[13,15],n_lens:[13,14],n_mass:[16,17],n_min_for_caust:24,n_point:[4,13,16,17,18,24],n_sourc:[9,13,14],n_step:31,n_walker:31,nad:31,name:[9,10,11,12,13,22,31],nan:8,nasa:30,natur:21,nbest:31,ndarrai:[8,9,10,12,13,14,15,16,17,20,21,22,23,24,25,32],ndim:32,ne:14,necessari:[15,28],need:[8,9,12,13,14,15,25,28,31,32],neg:[0,1,3,8,23,25,31],neither:[0,1],nelder:32,net:23,next:32,nfev:32,nice:[13,31,32],nicer:15,node:14,non:[14,31],none:[0,1,3,8,9,12,13,15,16,17,18,19,20,21,22,23,24,25],normal:[13,25],north:[7,14,23],north_project:7,notat:14,note:[0,1,3,4,7,8,11,12,13,14,15,16,17,18,24,30,31,32],notebook:30,noth:13,now:[30,31,32],np:[7,8,9,10,12,13,14,15,16,17,20,21,22,23,24,25,31,32],nu:14,number:[4,8,13,14,15,16,17,18,24,25,31],numer:20,numpi:[8,13,15,16,17,24,31,32],o:28,ob03235:30,ob03235_moa:30,ob03235_ogl:30,ob05086:31,ob08092:32,object:[0,1,3,4,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25],observ:[16,18,31],observatori:13,obtain:12,obviou:24,obvious:31,oct:25,off:[13,31],offset:23,ogl:[20,30,31],ogle_data:30,ok:[31,32],omega:14,onc:13,one:[8,9,13,14,15,16,18,31,32],ones:[8,25],onli:[4,8,9,13,14,15,16,17,18,20,24,31,32],onto:[16,18],op:32,open:28,opposit:[14,15],optim:[8,24,30,32],option:[0,1,4,8,9,12,13,15,22,23],orbit:[13,14],order:[0,1,3,8,13,15,24,28,32],org:25,origin:[0,1,3,4,23],os:[30,31,32],other:[0,1,3,4,8,9,11,12,13,15,24,25,31,32],otherwis:[8,9,25],our:32,out:[25,31],out_text:25,outlin:4,output:[10,23,25,32],outsid:12,over:[8,13,15],overal:9,overestim:[0,1],p:[24,31],packag:[20,25,27,28,29,31,32],paczynski:[12,20],page:[25,27,28],pair:[7,24],panel:24,paper:[0,1,24],parallax:[12,13,14,15,16,17,18,19,23,27],parallax_delta_n_:23,parallel:14,param:[14,31],paramet:[0,1,3,4,7,8,9,11,12,13,14,15,16,17,18,20,21,22,23,24,25,27,30,31],parameter:14,parameteris:24,parameters_to_fit:[31,32],pariti:[0,1,3],part:[9,25,28],partial:9,particular:[8,14],pass:[8,12,13,15,16,17,18,28],passag:14,path:[15,22,28,30,31,32],pc:[16,19],peak:32,peirson:3,per:24,percentil:31,perform:9,periapsi:14,perihelion:14,period:[14,15],perpendicular:[14,25],phot_fmt:[8,9,13,15,25,30],phot_ob08092_o4:32,photometr:[9,15],photometri:[13,15],photometry_fil:[30,31,32],physic:[16,18,19],pi_:[16,19],pi_e:[14,16,18,31],pi_e_:[8,9,14,31],pi_e_mag:14,pi_e_n:[8,9,14,31],pi_l:[16,17],pi_rel:[16,18],piec:[31,32],pip:28,place:25,plane:[4,7,14,16,18,23,25],planet:30,planetari:24,pleas:[3,28,30],plot:[4,8,9,13,15,16,17,18,24,25,30],plot_caust:[13,16,17,18],plot_data:[8,13,15,30,31,32],plot_lc:[8,13,15],plot_magnif:[13,16,18,30],plot_model:[8,30,31,32],plot_param:31,plot_properti:15,plot_residu:[8,13,15],plot_sourc:[8,13],plot_source_for_dataset:[8,13],plot_trajectori:[8,13],plotutil:25,plt:[4,13,24,30,31,32],plu:21,point:[0,1,3,4,8,9,12,13,14,15,16,17,18,20,21,24,25,30,31,32],point_sourc:12,point_source_magnif:[0,1,3],point_source_point_len:12,pointlen:[0,21,29],pointlenswithshear:[0,29],pointlik:20,polynomi:[0,1,3],posit:[0,1,3,8,10,12,13,15,20,23,24,31],possibl:[0,1,3,4,13,15],powel:32,pre:12,preced:8,precis:20,prepar:10,prescript:12,present:[0,1,20,32],previou:[8,13],previous:[8,9,13],primari:[0,1,3,14,15,16,17],print:[7,8,11,13,14,15,30,31,32],print_funct:30,prior:31,probabl:31,problem:[15,28],procedur:32,produc:[5,9,24],product:25,profil:[0,1,28],project:[7,14,16,17,18,23,25],proper:[14,16,18,20],properti:[4,7,8,9,10,12,13,14,15,16,17,18,19,23,24,31],propos:20,provid:[0,1,8,13,14,15,20,22,24,31,32],psbl:14,pspl:[9,14,20],pspl_magnif:[20,21],pspl_model:32,put:24,py:[8,9,13,28],pyplot:[4,13,15,16,17,18,24,25,30,31,32],python2:30,python3:30,python:[28,31],pythonpath:28,q:[0,1,3,4,6,14,16,17,23,24,25,30],quadrupol:[0,1,12,13],qualiti:13,quantiti:[14,16,17,18,19,25],quit:24,r:[28,31],r_e:[16,18],r_e_tild:[16,18],ra:[7,10,13,15],rad:14,radian:30,radiu:[0,1,13,14,16,18,19,20],rais:24,randn:31,random:31,randomli:24,rang:[8,12,13,14,15,22,24,30,31],rapper:28,rare:13,rate:14,ratio:[4,6,8,9,13,14,16,17,24,25,32],re:[0,1,3,9,28,30,31],read:[4,10,16,17,18],readi:31,real:24,recalcul:[8,9],recent:28,recommend:[0,1,3],rectilinear:31,red:[13,30,32],refer:[8,9,10,13,14,31],refit:9,refsdal:[5,21],regard:20,regardless:14,regress:13,reject:31,rel:[0,1,9,11,12,14,16,17,18,20,24,25,31],relat:[0,1,25,32],releas:28,relev:8,rememb:[8,11,13,31],remov:13,repeat:[24,28],report:[9,31],repositori:28,repr:[],repres:[7,25],represent:13,request:[0,1,7,8,9,13,25],requir:[8,9,15,23,24,28],rescal:9,reset:[8,13],reset_best_chi2:8,reset_plot_properti:13,reshap:31,residu:[8,9,13,15,31],reson:24,respect:[8,9,13,15],result:[0,1,8,30,31,32],result_fun:32,retriev:9,rho:[0,1,9,12,13,14,20,30],rho_1:14,rho_2:14,right:[0,1,3,24,28],ring:[0,1,3,4,6,14,16,17,18,24],rough:32,round:31,routin:[24,32],run:[8,9,12,24,28],run_mcmc:31,s:[0,1,3,4,6,7,13,14,16,17,18,20,24,25,28,30,31,32],sai:31,same:[8,12,13,15,20,24,25,32],sampl:[24,31],sample_file_01:15,satellit:[10,12,13,15,22,23],satellite_skycoord:[12,13,15,22,23],satelliteskycoord:[0,13,29],save:[9,15],scale:[4,8,9,13,15,16,17,25,30],scale_errorbar:15,scale_flux:9,scatter:[4,8,13,15,16,17,18,24],scheme:24,scipi:32,search:[27,31],sec:[0,1],second:[15,25,31,32],secondari:[0,1,3,14,16,17],see:[0,1,8,9,10,12,13,14,15,16,17,18,20,22,23,24,27,30,31,32],seem:32,select:24,self:[8,13,15],semi:14,separ:[0,1,3,4,6,8,9,12,13,14,16,17,20,23,24,25,31],sequenc:[0,1,32],set:[0,1,7,8,9,11,12,13,14,15,16,17,18,19,23,25,27,30,32],set_aspect:[4,13],set_dataset:13,set_default_magnification_method:13,set_limb_coeff_gamma:[11,13],set_limb_coeff_u:[11,13],set_limb_darkening_weight:15,set_magnification_method:[12,13,30],set_magnification_methods_paramet:[12,13],set_source_flux_ratio:13,set_source_flux_ratio_for_band:13,set_tim:13,setattr:[14,31,32],setenv:28,setup:28,share:28,shear:[3,5,6,12,14,21],shear_g:[3,5,6,14],sheat:[3,21],sheet:[3,14],shift:[14,23],should:[0,1,8,9,12,13,15,16,19,23,24,25,28,32],shouldn:15,show:[8,13,15,24,30,31,32],show_bad:[8,13,15],show_data:13,show_errorbar:[8,13,15],shown:[8,13],sigma:31,sign:14,signific:13,significantli:[0,1],similar:[15,24,25,28,32],simpl:30,simpli:14,simplifi:20,simpson:[0,1],sinc:30,singl:[5,8,9,13,16,17,21,25,31,32],size:[0,1,13,14,15,16,17,18,20,24,25],size_list:13,skip:28,skowron:[14,23],sky:[7,12,13,14,15,23,25],sky_coordin:7,skycoord:[7,8,12,13,15,22,23],slightli:[8,13,31,32],slow:[12,20],slower:[8,31],small:[0,1,20,25,31],smaller:[12,13,24],smallest:[8,31,32],smooth:24,so:[0,1,8,13,23,24,28,30,32],solar:[16,17],solmass:[16,17,18],solut:[0,1,3,31],some:[0,1,3,8,24,27,30,31,32],someth:[8,9],somewhat:30,sonam:28,sourc:[0,1,3,8,9,12,13,14,16,18,20,21,23,24,28,30,31,32],source_1_paramet:14,source_2_paramet:14,source_flux:[8,9,13],source_flux_0:9,source_flux_1:[9,13],source_flux_2:13,source_flux_ratio:[9,13],source_i:[0,1,3],source_x:[0,1,3],space:[8,9,11,13,15,24,25],special:15,specif:[13,14,24],specifi:[7,8,9,12,13,15,16,17,18,19,20,21,23,24,30,31,32],spheric:13,sphinx:[],spitzer_ephemeris_01:22,split:31,squar:[0,1,4,9,13,15],stabl:20,stamp:23,standard:[10,12,14,16,17,24],star:[16,18,19,20,31],starblg234:31,start:[0,1,3,8,13,24,30,31,32],start_1:31,startup:28,state:[0,1],step:[28,31,32],still:32,stoke:12,stop:13,store:[4,13,16,17,32],str:[7,8,9,10,11,12,13,14,15,22,23,25,32],string:[8,15,25],structur:[4,5,6,13,16,18],studi:11,subdirectori:28,submodul:29,subpackag:[29,32],subtract:[8,13,15,25],subtract_2450000:[8,13,15,25,30,31],subtract_2460000:[8,13,15,25],success:[13,32],suitabl:[0,1,13],sum:[0,1,8,9,13,14,16,17,23,25],sum_funct:8,sun:[13,31],suppli:15,support:14,suppos:[25,30],suppress:15,sure:[0,1,3,8,13,15,28],symbol:[13,24],symmetr:24,system:[0,1,3,5,6,9,13,14,15,16,17,18,23],systemat:24,t:[8,12,15,24,28,30,31],t_0:[8,9,13,14,24,30,31,32],t_0_1:[13,14],t_0_2:[13,14],t_0_kep:[13,14],t_0_par:[7,14,31],t_0_xi:14,t_1:31,t_2:31,t_caustic_in:[14,24],t_caustic_out:[14,24],t_e:[8,9,13,14,16,18,24,30,31,32],t_e_1:14,t_e_2:14,t_eff:[8,9,14],t_rang:[8,13,30],t_star:14,t_star_1:14,t_star_2:14,t_start:[13,31],t_stop:[13,31],tabl:12,take:[8,9,24,31],taken:[8,13,15,23],tbd:[15,16,17],tbl:30,tcsh:28,tdb:10,termin:28,tertiari:[16,17],test:[0,1,28],text:25,th:[8,9],than:[0,1,3,9,11,12,13,15,16,17,18,24,31,32],thei:[8,9,12,13,15,24,31],them:[0,1,8,15,31,32],theorem:[0,1,12],theta:[31,32],theta_:[14,16,18,25],thi:[0,1,3,4,7,8,9,11,12,13,14,15,16,17,20,21,23,24,25,28,30,31,32],thing:[30,32],third:32,those:[8,13,30],three:[0,1,15,31],thu:24,time:[0,1,8,10,12,13,14,15,16,18,22,23,24,25,28,30,31,32],timescal:[14,31],titl:[8,31],tnc:32,togeth:[30,32],tolist:32,topocentr:[12,13,23],topolog:24,total:[0,1,3,9,11,15,16,17],total_mass:[16,17],toward:7,trajectori:[0,8,9,12,13,14,15,20,21,24,29],transform:[13,25],translat:15,trend:31,tupl:[8,13],turn:[13,31,32],tutori:27,two:[0,1,3,4,9,13,14,15,24,30,31,32],txt:[28,30],type:[9,11,13,14,15,20,23,25,32],typic:15,u:[0,1,7,11,12,13,14,16,19,20,25],u_0:[8,9,14,16,18,24,30,31,32],u_0_1:14,u_0_2:14,u_limb_darken:[0,1],u_to_gamma:25,uncertainti:[0,1,9,15,25],under:24,underli:12,uniform:[12,20,24],uniform_caustic_sampl:14,uniformcausticsampl:[0,14,29],uniformli:24,unit:[7,13,14,15,16,17,18,19,28],unless:[13,15],unnorm:31,unpack:28,up:[22,30],updat:[8,9,13,28],update_caust:13,upper:24,us:[0,1,3,4,8,9,11,12,13,14,15,16,17,18,20,22,23,24,25,28,30,31,32],usag:[14,15],user:13,user_nam:28,util:[0,13,29],v:[11,15],v_earth_perp_:7,v_earth_perp_n:7,v_earth_project:[7,25],val:[],valerio:[0,1,12],valid:14,valu:[0,1,7,8,9,11,12,13,14,15,16,17,19,20,23,24,25,27,31,32],valueerror:24,variat:[0,1],vastli:[0,1],vbbinarylens:[12,20],vbbinarylensinglibrari:28,vbbinarylensinglibrary_wrapp:28,vbbl:[0,1,12,13,28,30],vbbl_magnif:[0,1,12],vbbl_on:3,ve:28,vector:[4,12,13,14,15,23,25,31],vector_1:25,vector_2:25,vector_product_norm:25,veloc:[7,14,25],velocity_of_earth:25,veri:[0,1,12,13,20,24,25,31,32],version:[0,25,28,29],via:[13,32],visibl:15,visual:24,vs:15,wa:[8,9,13,20,21,31,32],wai:[13,24],walker:31,want:[8,9,13,15,16,17,25,30,31,32],we:[8,13,24,25,30,31,32],websit:[0,1,31],weight:[11,15],well:[9,12,13,20,31],were:[15,30],what:[0,1,3,12,23],when:[8,11,12,13,15,24],where:[7,13,24,32],whether:[8,9,13,14,15],which:[7,8,9,12,13,14,15,16,17,22,23,24,31,32],which_caust:24,which_paramet:14,wide:24,width:13,wiki:25,wikipedia:25,window:28,without:[9,12,13,15,16,17,31],witt:[12,20],wl:28,won:15,work:[12,15,20,24,28,30],would:[9,32],written:30,wrong:31,wyrzykowski:31,x0:32,x:[0,1,3,4,8,10,13,15,16,17,23,24,30,32],x_caustic:24,x_caustic_in:[14,24],x_caustic_out:[14,24],x_caustic_point:24,x_caustic_sep:24,xallarap:[13,14],xallrap:14,xi_argument_of_latitude_refer:14,xi_eccentr:14,xi_inclin:14,xi_omega_nod:14,xi_omega_periapsi:14,xi_period:14,xi_semimajor_axi:14,xlabel:25,xlim:[30,31],xy:23,xyz:10,y:[0,1,3,4,10,16,17,23,24,25],year:[14,31],yerr:32,yet:[8,13,24,31],ylim:30,yoo:[12,20],you:[0,1,3,7,8,9,11,12,13,14,15,16,17,24,25,28,30,31,32],your:[24,28,30,32],yr:[14,16,18],z:10,zero:[8,14,24,25],zeropoint:25,zip:31,zoom:[31,32],zorder:15,zorder_list:13},titles:["MulensModel package","MulensModel.binarylens module","MulensModel.binarylensimports module","MulensModel.binarylenswithshear module","MulensModel.caustics module","MulensModel.causticspointwithshear module","MulensModel.causticswithshear module","MulensModel.coordinates module","MulensModel.event module","MulensModel.fitdata module","MulensModel.horizons module","MulensModel.limbdarkeningcoeffs module","MulensModel.magnificationcurve module","MulensModel.model module","MulensModel.modelparameters module","MulensModel.mulensdata module","MulensModel.mulensobjects package","MulensModel.mulensobjects.lens module","MulensModel.mulensobjects.mulenssystem module","MulensModel.mulensobjects.source module","MulensModel.pointlens module","MulensModel.pointlenswithshear module","MulensModel.satelliteskycoord module","MulensModel.trajectory module","MulensModel.uniformcausticsampling module","MulensModel.utils module","MulensModel.version module","Welcome to MulensModel\u2019s documentation!","How to install?","MulensModel","Very Basic Tutorial","Microlensing Parallax Fitting Tutorial","Basic Fitting Tutorial"],titleterms:{"import":[31,32],basic:[30,32],binarylen:1,binarylensimport:2,binarylenswithshear:3,caustic:4,causticspointwithshear:5,causticswithshear:6,combin:30,content:[0,16,27],coordin:7,data:[30,32],defin:30,document:27,event:8,exercis:[31,32],fit:[31,32],fitdata:9,horizon:10,how:28,indic:27,instal:28,introduc:30,len:17,limbdarkeningcoeff:11,magnificationcurv:12,microlens:31,model:[13,30,32],modelparamet:14,modul:[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26],mulensdata:15,mulensmodel:[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,29],mulensobject:[16,17,18,19],mulenssystem:18,old:28,packag:[0,16],parallax:31,paramet:32,plot:[31,32],pointlen:20,pointlenswithshear:21,prepar:32,run:31,s:27,sampler:31,satelliteskycoord:22,set:31,sourc:19,submodul:[0,16],subpackag:0,tabl:27,trajectori:23,tutori:[30,31,32],uniformcausticsampl:24,util:25,veri:30,version:26,welcom:27}}) \ No newline at end of file +Search.setIndex({docnames:["MulensModel","MulensModel.binarylens","MulensModel.binarylensimports","MulensModel.binarylenswithshear","MulensModel.caustics","MulensModel.causticspointwithshear","MulensModel.causticswithshear","MulensModel.coordinates","MulensModel.event","MulensModel.fitdata","MulensModel.horizons","MulensModel.limbdarkeningcoeffs","MulensModel.magnificationcurve","MulensModel.model","MulensModel.modelparameters","MulensModel.mulensdata","MulensModel.mulensobjects","MulensModel.mulensobjects.lens","MulensModel.mulensobjects.mulenssystem","MulensModel.mulensobjects.source","MulensModel.pointlens","MulensModel.pointlenswithshear","MulensModel.satelliteskycoord","MulensModel.trajectory","MulensModel.uniformcausticsampling","MulensModel.utils","MulensModel.version","index","install","modules","tutorial","tutorial_fit_pi_E","tutorial_fit_pspl"],envversion:{"sphinx.domains.c":2,"sphinx.domains.changeset":1,"sphinx.domains.citation":1,"sphinx.domains.cpp":5,"sphinx.domains.index":1,"sphinx.domains.javascript":2,"sphinx.domains.math":2,"sphinx.domains.python":3,"sphinx.domains.rst":2,"sphinx.domains.std":2,"sphinx.ext.intersphinx":1,sphinx:56},filenames:["MulensModel.rst","MulensModel.binarylens.rst","MulensModel.binarylensimports.rst","MulensModel.binarylenswithshear.rst","MulensModel.caustics.rst","MulensModel.causticspointwithshear.rst","MulensModel.causticswithshear.rst","MulensModel.coordinates.rst","MulensModel.event.rst","MulensModel.fitdata.rst","MulensModel.horizons.rst","MulensModel.limbdarkeningcoeffs.rst","MulensModel.magnificationcurve.rst","MulensModel.model.rst","MulensModel.modelparameters.rst","MulensModel.mulensdata.rst","MulensModel.mulensobjects.rst","MulensModel.mulensobjects.lens.rst","MulensModel.mulensobjects.mulenssystem.rst","MulensModel.mulensobjects.source.rst","MulensModel.pointlens.rst","MulensModel.pointlenswithshear.rst","MulensModel.satelliteskycoord.rst","MulensModel.trajectory.rst","MulensModel.uniformcausticsampling.rst","MulensModel.utils.rst","MulensModel.version.rst","index.rst","install.rst","modules.rst","tutorial.rst","tutorial_fit_pi_E.rst","tutorial_fit_pspl.rst"],objects:{"":[[0,0,0,"-","MulensModel"]],"MulensModel.BinaryLens":[[0,2,1,"","adaptive_contouring_magnification"],[0,2,1,"","hexadecapole_magnification"],[0,2,1,"","point_source_magnification"],[0,2,1,"","vbbl_magnification"]],"MulensModel.binarylens":[[1,1,1,"","BinaryLens"]],"MulensModel.binarylens.BinaryLens":[[1,2,1,"","adaptive_contouring_magnification"],[1,2,1,"","hexadecapole_magnification"],[1,2,1,"","point_source_magnification"],[1,2,1,"","vbbl_magnification"]],"MulensModel.binarylenswithshear":[[3,1,1,"","BinaryLensWithShear"]],"MulensModel.binarylenswithshear.BinaryLensWithShear":[[3,2,1,"","point_source_magnification"]],"MulensModel.caustics":[[4,1,1,"","Caustics"]],"MulensModel.caustics.Caustics":[[4,1,1,"","CriticalCurve"],[4,3,1,"","critical_curve"],[4,2,1,"","get_caustics"],[4,2,1,"","plot"]],"MulensModel.causticspointwithshear":[[5,1,1,"","CausticsPointWithShear"]],"MulensModel.causticswithshear":[[6,1,1,"","CausticsWithShear"]],"MulensModel.coordinates":[[7,1,1,"","Coordinates"]],"MulensModel.coordinates.Coordinates":[[7,3,1,"","east_projected"],[7,3,1,"","ecliptic_lat"],[7,3,1,"","ecliptic_lon"],[7,3,1,"","galactic_b"],[7,3,1,"","galactic_l"],[7,3,1,"","north_projected"],[7,2,1,"","v_Earth_projected"]],"MulensModel.event":[[8,1,1,"","Event"]],"MulensModel.event.Event":[[8,3,1,"","best_chi2"],[8,3,1,"","best_chi2_parameters"],[8,3,1,"","blend_fluxes"],[8,2,1,"","calculate_chi2_gradient"],[8,3,1,"","chi2"],[8,3,1,"","chi2_gradient"],[8,3,1,"","coords"],[8,3,1,"","data_ref"],[8,3,1,"","datasets"],[8,2,1,"","fit_fluxes"],[8,3,1,"","fits"],[8,3,1,"","fluxes"],[8,2,1,"","get_chi2"],[8,2,1,"","get_chi2_for_dataset"],[8,2,1,"","get_chi2_gradient"],[8,2,1,"","get_chi2_per_point"],[8,2,1,"","get_flux_for_dataset"],[8,2,1,"","get_ref_fluxes"],[8,3,1,"","model"],[8,2,1,"","plot"],[8,2,1,"","plot_data"],[8,2,1,"","plot_model"],[8,2,1,"","plot_residuals"],[8,2,1,"","plot_source_for_datasets"],[8,2,1,"","plot_trajectory"],[8,2,1,"","reset_best_chi2"],[8,3,1,"","source_fluxes"],[8,3,1,"","sum_function"]],"MulensModel.fitdata":[[9,1,1,"","FitData"]],"MulensModel.fitdata.FitData":[[9,1,1,"","FSPL_Derivatives"],[9,3,1,"","blend_flux"],[9,2,1,"","calculate_chi2_gradient"],[9,3,1,"","chi2"],[9,3,1,"","chi2_gradient"],[9,3,1,"","chi2_per_point"],[9,3,1,"","data_magnification"],[9,3,1,"","dataset"],[9,2,1,"","fit_fluxes"],[9,3,1,"","gamma"],[9,2,1,"","get_chi2_gradient"],[9,2,1,"","get_d_A_d_params_for_point_lens_model"],[9,2,1,"","get_d_A_d_rho"],[9,2,1,"","get_d_A_d_u_for_PSPL_model"],[9,2,1,"","get_d_A_d_u_for_point_lens_model"],[9,2,1,"","get_data_magnification"],[9,2,1,"","get_dataset_trajectory"],[9,2,1,"","get_model_fluxes"],[9,2,1,"","get_model_magnitudes"],[9,2,1,"","get_residuals"],[9,3,1,"","model"],[9,2,1,"","scale_fluxes"],[9,3,1,"","source_flux"],[9,3,1,"","source_flux_ratio"],[9,3,1,"","source_fluxes"],[9,2,1,"","update"]],"MulensModel.fitdata.FitData.FSPL_Derivatives":[[9,2,1,"","get_d_A_d_rho"],[9,2,1,"","get_gradient"]],"MulensModel.horizons":[[10,1,1,"","Horizons"]],"MulensModel.horizons.Horizons":[[10,3,1,"","time"],[10,3,1,"","xyz"]],"MulensModel.limbdarkeningcoeffs":[[11,1,1,"","LimbDarkeningCoeffs"]],"MulensModel.limbdarkeningcoeffs.LimbDarkeningCoeffs":[[11,2,1,"","get_limb_coeff_gamma"],[11,2,1,"","get_limb_coeff_u"],[11,2,1,"","get_weighted_limb_coeff_gamma"],[11,2,1,"","set_limb_coeff_gamma"],[11,2,1,"","set_limb_coeff_u"]],"MulensModel.magnificationcurve":[[12,1,1,"","MagnificationCurve"]],"MulensModel.magnificationcurve.MagnificationCurve":[[12,2,1,"","get_binary_lens_magnification"],[12,2,1,"","get_magnification"],[12,2,1,"","get_point_lens_magnification"],[12,3,1,"","methods_for_epochs"],[12,2,1,"","set_magnification_methods"],[12,2,1,"","set_magnification_methods_parameters"]],"MulensModel.model":[[13,1,1,"","Model"]],"MulensModel.model.Model":[[13,3,1,"","bandpasses"],[13,3,1,"","caustics"],[13,3,1,"","coords"],[13,3,1,"","data_magnification"],[13,3,1,"","datasets"],[13,3,1,"","default_magnification_method"],[13,2,1,"","get_data_magnification"],[13,2,1,"","get_lc"],[13,2,1,"","get_limb_coeff_gamma"],[13,2,1,"","get_limb_coeff_u"],[13,2,1,"","get_magnification"],[13,2,1,"","get_magnification_methods"],[13,2,1,"","get_magnification_methods_parameters"],[13,2,1,"","get_parallax"],[13,2,1,"","get_ref_fluxes"],[13,2,1,"","get_residuals"],[13,2,1,"","get_satellite_coords"],[13,2,1,"","get_trajectory"],[13,2,1,"","is_static"],[13,2,1,"","magnification"],[13,3,1,"","methods"],[13,3,1,"","n_lenses"],[13,3,1,"","n_sources"],[13,2,1,"","parallax"],[13,3,1,"","parameters"],[13,2,1,"","plot_caustics"],[13,2,1,"","plot_data"],[13,2,1,"","plot_lc"],[13,2,1,"","plot_magnification"],[13,2,1,"","plot_residuals"],[13,2,1,"","plot_source"],[13,2,1,"","plot_source_for_datasets"],[13,2,1,"","plot_trajectory"],[13,2,1,"","reset_plot_properties"],[13,2,1,"","set_datasets"],[13,2,1,"","set_default_magnification_method"],[13,2,1,"","set_limb_coeff_gamma"],[13,2,1,"","set_limb_coeff_u"],[13,2,1,"","set_magnification_methods"],[13,2,1,"","set_magnification_methods_parameters"],[13,2,1,"","set_source_flux_ratio"],[13,2,1,"","set_source_flux_ratio_for_band"],[13,2,1,"","set_times"],[13,2,1,"","update_caustics"]],"MulensModel.modelparameters":[[14,1,1,"","ModelParameters"],[14,4,1,"","which_parameters"]],"MulensModel.modelparameters.ModelParameters":[[14,3,1,"","alpha"],[14,2,1,"","as_dict"],[14,3,1,"","convergence_K"],[14,3,1,"","dalpha_dt"],[14,3,1,"","ds_dt"],[14,3,1,"","gamma"],[14,3,1,"","gamma_parallel"],[14,3,1,"","gamma_perp"],[14,2,1,"","get_alpha"],[14,2,1,"","get_s"],[14,3,1,"","is_external_mass_sheet"],[14,3,1,"","is_external_mass_sheet_with_shear"],[14,2,1,"","is_finite_source"],[14,2,1,"","is_static"],[14,3,1,"","is_xallarap"],[14,3,1,"","n_lenses"],[14,3,1,"","n_sources"],[14,3,1,"","pi_E"],[14,3,1,"","pi_E_E"],[14,3,1,"","pi_E_N"],[14,3,1,"","pi_E_mag"],[14,3,1,"","q"],[14,3,1,"","q_source"],[14,3,1,"","rho"],[14,3,1,"","rho_1"],[14,3,1,"","rho_2"],[14,3,1,"","s"],[14,3,1,"","shear_G"],[14,3,1,"","source_1_parameters"],[14,3,1,"","source_2_parameters"],[14,3,1,"","t_0"],[14,3,1,"","t_0_1"],[14,3,1,"","t_0_2"],[14,3,1,"","t_0_kep"],[14,3,1,"","t_0_par"],[14,3,1,"","t_0_xi"],[14,3,1,"","t_E"],[14,3,1,"","t_caustic_in"],[14,3,1,"","t_caustic_out"],[14,3,1,"","t_eff"],[14,3,1,"","t_star"],[14,3,1,"","t_star_1"],[14,3,1,"","t_star_2"],[14,3,1,"","u_0"],[14,3,1,"","u_0_1"],[14,3,1,"","u_0_2"],[14,3,1,"","uniform_caustic_sampling"],[14,3,1,"","x_caustic_in"],[14,3,1,"","x_caustic_out"],[14,3,1,"","xallarap_reference_position"],[14,3,1,"","xi_Omega_node"],[14,3,1,"","xi_argument_of_latitude_reference"],[14,3,1,"","xi_eccentricity"],[14,3,1,"","xi_inclination"],[14,3,1,"","xi_omega_periapsis"],[14,3,1,"","xi_period"],[14,3,1,"","xi_semimajor_axis"]],"MulensModel.mulensdata":[[15,1,1,"","MulensData"]],"MulensModel.mulensdata.MulensData":[[15,3,1,"","bad"],[15,3,1,"","bandpass"],[15,3,1,"","chi2_fmt"],[15,3,1,"","coords"],[15,2,1,"","copy"],[15,2,1,"","data_and_err_in_chi2_fmt"],[15,2,1,"","data_and_err_in_input_fmt"],[15,3,1,"","ephemerides_file"],[15,3,1,"","err_flux"],[15,3,1,"","err_mag"],[15,3,1,"","errorbars_scale_factors"],[15,3,1,"","errorbars_scaling_equation"],[15,3,1,"","flux"],[15,3,1,"","good"],[15,3,1,"","input_fmt"],[15,3,1,"","mag"],[15,3,1,"","n_epochs"],[15,2,1,"","plot"],[15,3,1,"","plot_properties"],[15,3,1,"","satellite_skycoord"],[15,2,1,"","scale_errorbars"],[15,2,1,"","set_limb_darkening_weights"],[15,3,1,"","time"]],"MulensModel.mulensobjects":[[16,1,1,"","Lens"],[16,1,1,"","MulensSystem"],[16,1,1,"","Source"],[17,0,0,"-","lens"],[18,0,0,"-","mulenssystem"],[19,0,0,"-","source"]],"MulensModel.mulensobjects.Lens":[[16,3,1,"","a_proj"],[16,3,1,"","caustics"],[16,3,1,"","distance"],[16,3,1,"","epsilon"],[16,3,1,"","mass"],[16,3,1,"","mass_1"],[16,3,1,"","mass_2"],[16,3,1,"","mass_3"],[16,3,1,"","n_masses"],[16,3,1,"","pi_L"],[16,2,1,"","plot_caustics"],[16,3,1,"","q"],[16,3,1,"","s"],[16,3,1,"","total_mass"]],"MulensModel.mulensobjects.MulensSystem":[[16,3,1,"","lens"],[16,3,1,"","mu_rel"],[16,3,1,"","pi_E"],[16,3,1,"","pi_rel"],[16,2,1,"","plot_caustics"],[16,2,1,"","plot_magnification"],[16,3,1,"","r_E"],[16,3,1,"","r_E_tilde"],[16,3,1,"","source"],[16,3,1,"","t_E"],[16,3,1,"","theta_E"]],"MulensModel.mulensobjects.Source":[[16,3,1,"","angular_radius"],[16,3,1,"","distance"],[16,3,1,"","pi_S"]],"MulensModel.mulensobjects.lens":[[17,1,1,"","Lens"]],"MulensModel.mulensobjects.lens.Lens":[[17,3,1,"","a_proj"],[17,3,1,"","caustics"],[17,3,1,"","distance"],[17,3,1,"","epsilon"],[17,3,1,"","mass"],[17,3,1,"","mass_1"],[17,3,1,"","mass_2"],[17,3,1,"","mass_3"],[17,3,1,"","n_masses"],[17,3,1,"","pi_L"],[17,2,1,"","plot_caustics"],[17,3,1,"","q"],[17,3,1,"","s"],[17,3,1,"","total_mass"]],"MulensModel.mulensobjects.mulenssystem":[[18,1,1,"","MulensSystem"]],"MulensModel.mulensobjects.mulenssystem.MulensSystem":[[18,3,1,"","lens"],[18,3,1,"","mu_rel"],[18,3,1,"","pi_E"],[18,3,1,"","pi_rel"],[18,2,1,"","plot_caustics"],[18,2,1,"","plot_magnification"],[18,3,1,"","r_E"],[18,3,1,"","r_E_tilde"],[18,3,1,"","source"],[18,3,1,"","t_E"],[18,3,1,"","theta_E"]],"MulensModel.mulensobjects.source":[[19,1,1,"","Source"]],"MulensModel.mulensobjects.source.Source":[[19,3,1,"","angular_radius"],[19,3,1,"","distance"],[19,3,1,"","pi_S"]],"MulensModel.pointlens":[[20,1,1,"","PointLens"],[20,4,1,"","get_pspl_magnification"]],"MulensModel.pointlens.PointLens":[[20,2,1,"","get_point_lens_LD_integrated_magnification"],[20,2,1,"","get_point_lens_finite_source_magnification"],[20,2,1,"","get_point_lens_large_LD_integrated_magnification"],[20,2,1,"","get_point_lens_large_finite_source_magnification"],[20,2,1,"","get_point_lens_limb_darkening_magnification"],[20,2,1,"","get_point_lens_uniform_integrated_magnification"]],"MulensModel.pointlenswithshear":[[21,1,1,"","PointLensWithShear"]],"MulensModel.pointlenswithshear.PointLensWithShear":[[21,2,1,"","get_point_lens_LD_integrated_magnification"],[21,2,1,"","get_point_lens_finite_source_magnification"],[21,2,1,"","get_point_lens_large_LD_integrated_magnification"],[21,2,1,"","get_point_lens_large_finite_source_magnification"],[21,2,1,"","get_point_lens_limb_darkening_magnification"],[21,2,1,"","get_point_lens_uniform_integrated_magnification"],[21,2,1,"","get_point_source_magnification"]],"MulensModel.satelliteskycoord":[[22,1,1,"","SatelliteSkyCoord"]],"MulensModel.satelliteskycoord.SatelliteSkyCoord":[[22,2,1,"","get_satellite_coords"]],"MulensModel.trajectory":[[23,1,1,"","Trajectory"]],"MulensModel.trajectory.Trajectory":[[23,3,1,"","d_perp"],[23,2,1,"","get_xy"],[23,3,1,"","parallax_delta_N_E"],[23,3,1,"","x"],[23,3,1,"","y"]],"MulensModel.uniformcausticsampling":[[24,1,1,"","UniformCausticSampling"]],"MulensModel.uniformcausticsampling.UniformCausticSampling":[[24,2,1,"","caustic_point"],[24,2,1,"","check_valid_trajectory"],[24,2,1,"","get_standard_parameters"],[24,2,1,"","get_uniform_sampling"],[24,2,1,"","get_x_in_x_out"],[24,2,1,"","jacobian"],[24,3,1,"","n_caustics"],[24,3,1,"","q"],[24,3,1,"","s"],[24,2,1,"","which_caustic"]],"MulensModel.utils":[[25,1,1,"","PlotUtils"],[25,1,1,"","Utils"]],"MulensModel.utils.PlotUtils":[[25,2,1,"","find_subtract"],[25,2,1,"","find_subtract_xlabel"],[25,2,1,"","get_color_differences"],[25,2,1,"","get_y_value_y_err"]],"MulensModel.utils.Utils":[[25,2,1,"","astropy_version_check"],[25,2,1,"","complex_fsum"],[25,2,1,"","date_change"],[25,2,1,"","dot"],[25,2,1,"","gamma_to_u"],[25,2,1,"","get_flux_and_err_from_mag"],[25,2,1,"","get_flux_from_mag"],[25,2,1,"","get_mag_and_err_from_flux"],[25,2,1,"","get_mag_from_flux"],[25,2,1,"","get_n_caustics"],[25,2,1,"","u_to_gamma"],[25,2,1,"","vector_product_normalized"],[25,2,1,"","velocity_of_Earth"]],MulensModel:[[0,1,1,"","BinaryLens"],[1,0,0,"-","binarylens"],[2,0,0,"-","binarylensimports"],[3,0,0,"-","binarylenswithshear"],[4,0,0,"-","caustics"],[5,0,0,"-","causticspointwithshear"],[6,0,0,"-","causticswithshear"],[7,0,0,"-","coordinates"],[8,0,0,"-","event"],[9,0,0,"-","fitdata"],[10,0,0,"-","horizons"],[11,0,0,"-","limbdarkeningcoeffs"],[12,0,0,"-","magnificationcurve"],[13,0,0,"-","model"],[14,0,0,"-","modelparameters"],[15,0,0,"-","mulensdata"],[16,0,0,"-","mulensobjects"],[20,0,0,"-","pointlens"],[21,0,0,"-","pointlenswithshear"],[22,0,0,"-","satelliteskycoord"],[23,0,0,"-","trajectory"],[24,0,0,"-","uniformcausticsampling"],[25,0,0,"-","utils"],[26,0,0,"-","version"]]},objnames:{"0":["py","module","Python module"],"1":["py","class","Python class"],"2":["py","method","Python method"],"3":["py","property","Python property"],"4":["py","function","Python function"]},objtypes:{"0":"py:module","1":"py:class","2":"py:method","3":"py:property","4":"py:function"},terms:{"0":[0,1,3,6,8,9,12,13,14,16,17,20,23,24,25,28,30,31,32],"00":[7,8,25],"000":7,"00096":30,"001":[0,1,31],"0039":30,"005":[12,13],"01":31,"04":31,"06":30,"086":31,"0j":6,"1":[0,1,3,4,6,9,11,12,13,14,15,16,17,20,23,24,25,28,30,31,32],"10":[8,12,24,25,31],"100":[12,31],"1000":13,"10000":24,"102":31,"1025":31,"11":31,"110":31,"1103":31,"12":25,"120":30,"120755":31,"133":30,"139":[12,20],"15":[12,13,31],"150":31,"1593":[0,1,12,25],"16":[30,31],"1679":[0,1,12],"17":28,"18":[7,8,31],"180":[7,14,23],"18h00m00":7,"19":30,"1979":21,"1994":[12,20],"2":[4,6,8,9,11,13,14,15,16,17,24,25,30,31,32],"20":32,"200":[12,20,24],"2002":11,"2003":[20,30],"2004":[12,20,30],"2005":31,"2007":[0,1,12],"2008":[0,1,12,14,24,25],"2009":[12,20,24],"2010":[0,1,12,24],"2011":[14,23],"2012":24,"2015":[25,31],"2018":[12,20],"2022":3,"21":31,"2188":[0,1,12],"218982":31,"22":[13,25],"2228":24,"223":30,"23":24,"235":30,"24":3,"245":25,"2450000":[8,13,14,15,25,31],"2452810":30,"2452833":30,"2452845":30,"2452848":30,"2452890":30,"2453200":31,"2453628":31,"2453630":31,"2453950":31,"2455746":[12,13],"2455747":[12,13],"2455748":[12,13],"2456789":13,"246":25,"2460000":[8,13,15,25],"25":30,"26":31,"262":20,"27":24,"270":7,"2719":31,"2810":30,"282":21,"2890":30,"29062":31,"2972":31,"2d":14,"2f":32,"3":[9,11,13,14,22,24,25,31],"30":[7,8,25,31],"30d00m00":7,"3291":31,"3400":31,"35":14,"35507":31,"360":7,"37":31,"37263":31,"377":[0,1,12],"3800":31,"387105":31,"395":24,"3d":25,"3f":32,"4":[0,1,12,24],"40":31,"400":31,"408":[0,1,12],"41":31,"415677":31,"421l":[12,20],"426":24,"430":[12,20],"44":31,"45":31,"479":[12,20],"488817":31,"491":24,"4f":32,"5":[11,13,15,20,28,30,31,32],"50":31,"500":31,"5000":[4,13,16,17,18],"505":[12,20],"515":24,"5157":[12,20],"52":24,"521":11,"53":30,"5380":32,"561":21,"572":11,"587":24,"59":31,"5f":[31,32],"5th":[0,1],"6":[12,13,24,30,31],"603":[12,20],"61":30,"611301":31,"6700":31,"67778":31,"681":[0,1,12,25],"695":[12,20],"7":[12,13,24,31],"71":[12,20,31],"738":14,"787":24,"8":[30,32],"84":31,"87":14,"927":3,"93":31,"947":31,"949":31,"95":31,"9th":3,"boolean":[0,1,8,9,13,14,15,20,23],"case":[7,11,13,14,15,16,17,24,31,32],"class":[0,1,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,27,30,32],"default":[0,1,4,7,8,9,12,13,14,15,23,24,25,30],"do":[0,1,3,8,13,14,24,30,31],"export":28,"final":[0,1,32],"float":[0,1,3,4,5,6,7,8,9,11,12,13,14,15,16,17,18,19,20,22,23,24,25,32],"function":[0,1,8,9,12,13,14,15,16,17,23,24,25,27,31,32],"import":[7,8,10,13,24,30],"int":[4,8,13,14,15,16,17,18,24,25],"long":[15,31],"public":20,"return":[0,1,3,4,7,8,9,10,11,12,13,14,15,20,21,22,24,25,27,31,32],"static":[8,14,25],"true":[0,1,3,8,9,12,13,14,15,23,24,30,31],"try":[31,32],"while":15,A:[7,9,11,13,14,15,16,17,18,20,23,24,25,32],And:[30,31],As:[0,1,31],At:31,Be:20,But:31,For:[0,1,8,9,10,12,13,14,20,23,24,27,30],If:[0,1,3,7,8,9,12,13,14,15,16,17,18,19,20,25,28,30,31,32],In:[7,8,13,14,24,28,31,32],Is:25,It:[0,1,8,11,13,14,16,18,23,30,31],NOT:30,No:9,Not:[13,16,17,21],ON:13,ONE:13,Or:30,That:27,The:[0,1,3,4,8,9,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,30,31,32],Then:[13,14,30,31,32],There:8,These:8,To:[4,9,13,14,24,28,32],Will:[8,15],__future__:30,__init__:15,_lens_:12,_not_:20,_point_:12,a_proj:[16,17],abil:32,about:[13,31,32],abov:[8,9,24,28,32],acc:[0,1],accept:[4,8,9,13,15],access:[8,14],accord:[14,24],account:[13,23],accur:[8,25],accuraci:[0,1,12,13],across:15,actual:[16,17],ad:[8,15,16,17,30],adapt:[0,1],adaptive_contour:12,adaptive_contouring_magnif:[0,1,12],adaptivecontour:[0,1,12,28],add:[0,1,15,28],add_2450000:[15,31],add_2460000:15,addit:[12,13,15],addition:31,adequ:13,advanc:[0,1,24],advis:30,advoc:14,affect:[0,1,13,14,23],after:[9,16,17],al:[3,11,12,14,20,23,24,25,30,31],algorithm:[0,1,24,25,31],all:[8,9,13,15,16,17,23,24,27,28,30,32],all_approxim:[0,1],allow:[8,9,12,13,14,32],almost:15,along:[4,13,14],alpha:[13,14,15,16,18,23,24,30,31],alpha_list:13,alreadi:8,also:[0,1,8,9,13,15,16,18,20,24,30,31],altern:[13,24],alwai:[4,6,14,15,16,18],amplitud:32,an:[7,8,9,10,11,13,14,15,16,17,22,24,25,27,31,32],analysi:24,analyz:13,angl:[7,13,14,24],angular:[16,18,19],angular_radiu:[16,19],ani:[13,25],annual:31,annuli:[12,20],anomali:14,anyth:9,apj:[0,1,3,11,12,14,20,25],appendix:23,appli:[9,13,24,32],applic:[0,1,8,9,15],approach:[20,31],appropri:[8,13,15],approxim:[0,1,12,20,32],ar:[0,1,3,4,7,8,9,12,13,14,15,16,17,18,19,20,23,24,25,27,28,30,31,32],arbitrari:[9,27],archiv:[28,30],area:[0,1],aren:15,arg:[7,13,14,21,31,32],argument:[8,9,10,12,13,14,16,17,18,20,21,23,24,25,32],around:[31,32],arrai:[7,8,9,16,17,20],arrow:13,arrow_kwarg:13,as_dict:14,ascend:14,associ:32,assum:[8,12,13,14,15,16,17,18,19,20,23],astropi:[7,8,10,12,13,14,15,16,17,18,19,22,23,25,28],astropy_version_check:25,attach:13,attr:[8,9],attribut:[0,1,3,4,5,6,8,12,13,14,16,19,20,22,23,32],au:[16,17,18,23],author:[15,28],avail:[12,30],axi:[0,1,3,4,8,13,14,15,24,25,31,32],back:31,background:[16,19],bad:[8,9,13,15],band:[8,11,13,15],bandpass:[9,11,13,15],barycentr:25,base:[0,1,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,32],baselin:32,bash:28,bashrc:28,basic:[8,14],bayesian:24,becaus:[8,14,15,24,31],been:8,befor:[0,1,3,16,17,28],begin:[14,25,30],behavior:[8,13],below:[24,31],best:[9,13,14,27,30,31,32],best_chi2:8,best_chi2_paramet:8,better:[31,32],between:[0,1,3,4,6,9,12,13,14,16,17,18,24,25,31],beyond:22,bfg:32,binari:[0,1,3,6,8,9,12,13,14,23,24,25],binarylen:[0,3,29],binarylensimport:[0,29],binarylenswithshear:[0,29],bjd:15,bjd_tdb:15,black:[13,25,30,31],blend:[8,9,13],blend_flux:[8,9,13],blend_flux_0:9,blg234:31,blg:[20,30,31],bodi:[0,1,3,4,6,14,16,17,18,30],bond:30,bool:[8,9,13,14,24,25],both:[0,1,3,4,9,11,13,24,25,31],bottom:24,bozza:[0,1,12,20],bright:[15,31],brighter:32,burn:31,button:28,c:[20,24,28,32],caclul:8,calcul:[0,1,3,4,7,8,9,12,13,14,15,20,21,22,23,24,25,31,32],calculate_chi2_gradi:[8,9],call:[8,9,13,14,15,24,32],can:[0,1,7,8,9,11,12,13,14,15,20,24,27,28,30,31,32],cannot:14,cartesian:25,cartesianrepresent:[10,25],cassan:[14,24],caus:14,caustic:[0,5,6,8,13,14,16,17,18,24,25,29],caustic_point:24,causticspointwithshear:[0,29],causticswithshear:[0,29],center:[0,1,3,4,11,14,23],central:24,cg:32,chain:31,chang:[5,8,9,13,14,15,21,24,25,32],check:[14,15,24,25,28,31],check_valid_trajectori:24,chi2:[8,9,15,27,30,31],chi2_fmt:[8,15],chi2_for_model:32,chi2_gradi:[8,9],chi2_initi:32,chi2_min:31,chi2_min_theta:31,chi2_per_point:9,chi:[8,9,15,31,32],choic:32,choos:31,chosen:[0,1],circl:13,circular:[12,13],cite:[3,24],clariti:31,clockwis:24,close:24,cobyla:32,code:[12,15,24,25,28,30,31,32],coeff:13,coeffici:[0,1,9,11,12,13,15,16,19,20,25],color:[13,15,24,25,30,31],color_differ:25,color_list:[13,25],colorbar:24,column:15,combin:[8,13,14,16,17,31,32],come:31,command:[13,28],comment:30,common:31,companion:[16,17],compar:[12,13],complex128:24,complex:[0,1,3,5,6,14,24,25],complex_fsum:25,complic:24,compon:[7,9,14,16,17,23,24],compos:21,comput:[9,12,20,31],condit:24,configur:24,congratul:[28,32],conniv:7,consid:[12,21,24],consist:[15,16,18],constraint:[13,32],contact:[15,28],contain:[8,15,32],content:29,contour:[0,1,12],contrari:[8,13],contribut:[8,9,11],control:[8,13],conveni:31,convent:[0,1,11,12,13,14,15,23,24,25],converg:[3,5,6,12,14,21],convergence_k:[3,5,6,14],convert:32,coord:[8,12,13,15,23,31],coordin:[0,1,3,4,8,12,13,14,15,16,17,22,23,24,29,31],copi:15,correspond:[4,9,13,16,17,24,31,32],could:[15,30],counter:24,counterclockwis:14,cover:15,cpp:28,creat:9,critic:4,critical_curv:4,criticalcurv:4,cross:[14,16,18,24],cshrc:28,cubic:22,cunstructor:32,current:[8,9,13,20,28,31],curv:[4,8,12,13,16,18,20,27,31],curvelinear:[14,24],d:[8,9],d_a_d_rho:9,d_perp:23,da:9,da_dparam:9,da_du:9,dai:[14,16,18,32],dalpha_dt:14,darken:[0,1,8,9,11,12,13,15,16,19,20,25],dash:31,dat:[22,31,32],data:[8,9,13,15,22,25,27,31],data_and_err_in_chi2_fmt:15,data_and_err_in_input_fmt:15,data_err:15,data_list:15,data_magnif:[9,13],data_num:8,data_path:[30,31,32],data_ref:[8,13],datapoint:9,dataset:[8,9,13,15,30,31,32],date:[12,13,15,25],date_chang:25,dec:[7,10,13,15],decid:12,def:[31,32],default_magnification_method:13,default_method:12,defin:[0,1,3,4,7,8,9,13,14,15,16,17,23,24,28,32],definit:[16,17,23],deg:[7,14,23],degener:31,degre:[7,13,14,16,18,30],depend:[8,10,32],deprec:[8,9,13,15],deriv:[9,13],descript:[12,24,27],desir:9,detail:[8,13],determin:[0,1,9,13,30],dfm:31,dict:[8,9,11,12,13,14,15,23,24,31],dictionari:[11,12,13,14,15,23,24],didn:28,differ:[8,13,15,23,24,25,31,32],dimensionless:[16,18,23],direct:[13,14,20],directli:[12,13,14],directori:28,dispers:31,distanc:[10,16,17,19],distribut:24,doctr:[8,13],doe:[0,1,8,9,14,24,30,31],doesn:12,dominik:[0,1,12],don:[12,24,30,31],done:[15,31,32],dot:25,dot_product:25,doubl:[13,14,25],download:[28,30,31,32],dozen:32,draw:24,drawn:24,ds_dt:14,dt:13,du:9,due:13,dure:[0,1],e9967a:25,e:[4,5,7,8,9,10,11,12,13,14,15,20,21,23,24,25,28,32],each:[8,9,11,12,13,15,16,17,20,21,23,24],earth:[7,13,23,25,31],earth_coord:23,earth_orbit:[12,13,23],easi:[20,31],easili:32,east:[7,14,23],east_project:7,eccentr:14,eclipt:7,ecliptic_lat:7,ecliptic_lon:7,effect:[12,13,14,15,20,23,30,31],einstein:[0,1,3,4,6,14,16,17,18,20,24],either:[0,1,8,13,14,15,16,19,28],element:[9,13],ellips:13,els:32,emce:31,en:25,enclos:[0,1],encourag:32,enough:31,ensemblesampl:31,ensur:9,entranc:[14,24],enumer:32,ephemeri:15,ephemerid:[10,12,13,15,22,23],ephemerides_fil:[9,13,15,22],ephemeris_fil:22,epheri:22,epoch:[7,8,12,13,14,15,22,23,24,25,31],epsilon:[16,17],eq:24,equal:[4,13,14,16,18,24],equat:[0,1,3,12,15,20],err:15,err_flux:[9,15,25],err_mag:[15,25,32],error:[0,1,12,15],errorbar:[8,9,13,15,32],errorbars_scale_factor:15,errorbars_scaling_equ:15,es:[8,9,13],estim:[0,1,32],et:[3,11,12,14,20,23,24,25,30,31],etc:[0,1,3,8,14,15],evalu:[9,15,24,32],even:9,evenli:13,event:[0,1,7,12,13,15,23,24,27,29,30,31,32],event_0:31,event_1:31,event_2:31,everi:[11,13,28,31],everyth:30,ew:31,exactli:13,exampl:[8,13,14,15,22,24,30,31,32],except:[13,23,24],exclud:15,execut:[0,1,9],exist:[8,24],exit:[14,24],exoplanet:30,expand:30,expect:13,explan:[8,13],explicitli:[13,14,30],extend:15,extern:[3,14],extract:28,extrapol:22,f_blend:13,f_sourc:13,factor:15,faint:15,fals:[0,1,8,9,12,13,15,20,23,24,25],far:12,fast:20,favourit:32,featur:[13,30],few:32,fig:24,figur:31,file:[10,12,13,15,22,23,25,28,31,32],file_nam:[10,15,30,31,32],fill:13,filter:11,find:[13,25,27,31],find_subtract:25,find_subtract_xlabel:25,fine:[24,31,32],finer:13,finit:[0,1,9,12,14,20,30],finite_source_ld_lee09:12,finite_source_ld_wittmao94:12,finite_source_ld_yoo04:12,finite_source_ld_yoo04_direct:12,finite_source_uniform_gould94:12,finite_source_uniform_gould94_direct:12,finite_source_uniform_lee09:12,finite_source_uniform_wittmao94:12,first:[8,11,13,14,15,21,25,30,31,32],fit:[8,9,13,14,15,24,27,30],fit_blend:[8,13],fit_flux:[8,9],fitdata:[0,8,29],fix:[8,9,11,13,14],fix_blend:8,fix_blend_flux:[8,9],fix_source_flux:[8,9],fix_source_flux_ratio:[8,9],flag:[8,15],floor:15,flux:[8,9,11,13,15,25,30,32],flux_err:25,flux_ratio_constraint:13,flux_source_1:13,flux_source_2:13,flux_source_band_1:13,flux_source_band_2:13,fmt:[15,32],fold:24,follow:[23,25,30],form:[12,13,31],formal:20,format:[8,15,25,31,32],found:32,four:24,fpic:28,fraction:[0,1,3,4,6,14,16,17,18,20],frame:[10,25],free:9,frequent:[11,24],from:[7,8,9,10,12,13,14,15,20,22,23,24,25,28,30,31,32],fsbl:14,fspl:[9,14],fspl_deriv:9,fsplderiv:9,fsum:8,full:[25,31],full_bjd:[7,25],fulli:28,fun:32,func:9,further:31,g:[7,8,9,11,12,13,14,15,20,23,25,28,32],galact:7,galactic_b:7,galactic_l:7,gamma:[0,1,9,11,12,13,14,15,20,25],gamma_parallel:14,gamma_perp:14,gamma_to_u:25,gaussian:15,gave:8,gca:[4,13,32],gcc:28,gener:[12,23,25,27],get:[8,11,13,14,15,24,25,30,32],get_:14,get_alpha:14,get_binary_lens_magnif:12,get_caust:4,get_chi2:[8,30,31,32],get_chi2_for_dataset:8,get_chi2_gradi:[8,9],get_chi2_per_point:8,get_color_differ:25,get_d_a_d_params_for_point_lens_model:9,get_d_a_d_rho:9,get_d_a_d_u_for_point_lens_model:9,get_d_a_d_u_for_pspl_model:9,get_data_magnif:[9,13],get_dataset_trajectori:9,get_flux_and_err_from_mag:25,get_flux_for_dataset:8,get_flux_from_mag:25,get_gradi:9,get_lc:13,get_limb_coeff_gamma:[9,11,13],get_limb_coeff_u:[11,13],get_mag_and_err_from_flux:25,get_mag_from_flux:25,get_magnif:[12,13],get_magnification_method:13,get_magnification_methods_paramet:13,get_model_flux:9,get_model_magnitud:9,get_n_caust:25,get_parallax:13,get_point_lens_finite_source_magnif:[20,21],get_point_lens_large_finite_source_magnif:[20,21],get_point_lens_large_ld_integrated_magnif:[20,21],get_point_lens_ld_integrated_magnif:[20,21],get_point_lens_limb_darkening_magnif:[20,21],get_point_lens_magnif:12,get_point_lens_uniform_integrated_magnif:[20,21],get_point_source_magnif:21,get_pspl_magnif:20,get_ref_flux:[8,13],get_residu:[9,13],get_satellite_coord:[13,22],get_standard_paramet:24,get_trajectori:13,get_uniform_sampl:24,get_weighted_limb_coeff_gamma:11,get_x_in_x_out:24,get_xi:23,get_y_value_y_err:25,github:28,give:[4,8,9,11,13,14,15,22,32],given:[0,1,3,4,7,8,9,11,12,13,14,15,16,17,19,22,23,24,25,27,32],go:[13,24,28],goe:31,good:[9,15,30,32],got:31,gould:[0,1,12,20],gradient:[8,9],green:[0,1,28],ground:13,guarante:30,guess:32,h:20,ha:[8,9,11,13,14,16,18,23,24,31,32],had:30,hand:[0,1,3,24],handl:14,happen:13,have:[8,11,13,15,28,31,32],haven:[8,31],he:31,henc:[12,13,31,32],here:[13,15,20,24,31,32],hexadecapol:[0,1,12,13],hexadecapole_magnif:[0,1,12],high:13,higher:[0,1,3,23,31],hjd:[15,31],hold:10,home:28,hope:31,horizon:[0,22,29],hour:13,hourangl:7,how:[0,1,15,23,30,31,32],howev:24,http:[25,31],human:[8,13,15],i:[0,1,4,5,8,9,10,11,12,13,14,15,20,21,24,31,32],i_caust:24,id:31,ignor:[0,1,13],imag:[0,1,3,24],impact:[16,18,24,31,32],implement:[0,1,4,13,16,17,20,21],importantli:25,inclin:14,includ:[8,9,12,13,14,15,20,23,28,31,32],index:[8,13,27,31],index_dataset:8,indic:24,individu:9,inf:31,info:10,inform:[8,13,14,15,31,32],inherit:7,initi:32,initial_guess:32,input:[10,14,15,20,23,24,25,27],input_fmt:15,inspect:32,instal:[25,27],instanc:[7,8,11,13,14,15,23,32],instantan:[14,20],instead:[0,1,8,13,14,20,24,31],instruct:[10,15],integr:[0,1,12,20,24],intermedi:24,intern:[4,13,14,23,24],interpol:[12,15,20,22],interv:13,introduc:[11,15,25],invert:32,invert_yaxi:32,invok:28,io:31,is_external_mass_sheet:14,is_external_mass_sheet_with_shear:14,is_finite_sourc:14,is_stat:[13,14],is_xallarap:14,isfinit:31,isinst:32,isnan:31,iter:[12,25],its:[0,1,8,11,15],j:20,jacobian:[8,9,14,24],jd:32,join:[30,31,32],jpl:[10,22],julian:[12,13,25],jupyt:30,just:[8,13,22,28,32],k2_ephemeris_01:22,k:[8,9],kain:24,keep:[8,24],kei:[8,9,11,12,13,15,23,24,32],keyword:[4,7,8,9,13,15,16,17,18,20,21,22,24],km:[7,25],know:[0,1,3,30],kpc:[16,17,19],kwarg:[4,7,8,9,12,13,15,16,17,18,21],l:7,label:[15,31,32],label_list:13,larg:[8,12,20],larger:[4,13,15,20,24],last:[28,32],later:[20,25,31],latitud:[7,14],latter:[8,15,31],ld_accuraci:[0,1],learn:[31,32],least:[9,14],leav:30,lee:[12,20],left:[0,1,3,4],legend:[8,31,32],len:[0,1,3,4,5,6,8,9,12,13,14,16,18,20,21,23,24,25,30,31,32],length:[4,13,15,24,25],lens:[4,8,13,20,27],less:32,let:[31,32],librari:[0,1],light:[8,12,13,20,27,31],like:[13,24],likelihood:[27,31],limb:[0,1,8,9,11,12,13,15,16,19,20,25],limb_darken:[16,19],limbdarkeningcoeff:[0,13,15,16,19,20,29],limit:13,line:[0,1],linear:[0,1,9,11],linestyl:30,linewidth:13,link:[8,13],linspac:24,list:[4,8,9,12,13,14,15,16,17,22,23,24,25,31,32],lm:28,ln_like:31,ln_like_:31,ln_prior:31,ln_prior_:31,ln_prob:31,load:15,loadtxt:15,loc:[31,32],logarithm:31,longer:[8,13],longitud:[7,14],look:[13,31,32],lower:[0,1,3,24],ls:31,lstdc:28,lumin:[13,14],lw:31,m:[0,1],m_1:[16,17],m_i:[16,17],ma:[16,18,19],macho:20,mag:[8,9,13,15,25,32],mag_zeropoint:13,magnif:[0,1,3,9,12,13,14,16,18,20,21,32],magnification_vector:13,magnificationcurv:[0,13,29],magnitud:[8,9,13,14,15,25,32],mai:[0,1,3,4,7,8,13,14,15,16,17,19,31,32],main:[14,27,30],major:14,make:[0,1,3,8,13,15,28,31],mani:14,mao:[12,20],mark:[8,13,15],marker:[8,13,15],marker_list:13,markers:15,martin:[0,1,12],mass:[0,1,3,4,5,6,14,16,17,18,20,21,23,24,25],mass_1:[0,1,3,16,17],mass_2:[0,1,3,16,17],mass_3:[16,17],materi:32,math:8,matplotlib:[4,8,13,15,24,25,30,31,32],matter:15,mayb:13,mead:32,mean:[0,1,3,8,13,24,25],measur:[14,15],mention:[8,9],messag:32,method:[0,1,12,13,20,30,32],method_paramet:13,methods_for_epoch:12,methods_paramet:[12,13],microarcsecond:[16,19],microlens:[8,11,12,13,14,15,16,18,20,23,24,27,30],might:30,millarcsecond:[16,19],milliarcsecond:[16,17,18,19],mind:24,minim:32,minimum:[14,15,24,25],minut:22,miss:[16,17],mitig:[4,13],mm:[30,31,32],mnra:[0,1,12,20,24],moa:30,moa_data:30,mode:31,model:[0,8,9,12,14,15,16,18,20,21,23,24,27,29,31],model_0:31,model_1:31,model_2:31,model_flux:9,model_mag:9,model_paramet:14,modelparamet:[0,12,13,16,18,20,21,23,24,27,29],modifi:31,modul:[27,29,31,32],month:25,more:[8,9,11,13,15,16,17,18,24,31,32],most:[7,13,14,15,25,30],motion:[13,14,16,18,20,31],mu:27,mu_rel:[16,18],much:32,mulens_data:15,mulensdata:[0,8,9,12,13,27,29,30,31,32],mulensmodel:[28,30,31,32],mulensobject:[0,29],mulenssystem:[0,16],multi:14,multinest:31,multipl:[12,13,15,20],multipli:15,must:[8,9,13,14,15,16,18,24],my_1s2l_model:30,my_data:[8,31,32],my_ev:[30,31,32],my_friends_data:32,my_model:31,my_pspl_model:30,n:[23,24],n_burn:31,n_caustic:[24,25],n_dim:31,n_epoch:[13,15],n_lens:[13,14],n_mass:[16,17],n_min_for_caust:24,n_point:[4,13,16,17,18,24],n_sourc:[9,13,14],n_step:31,n_walker:31,nad:31,name:[9,10,11,12,13,22,31],nan:8,nasa:30,natur:21,nbest:31,ndarrai:[8,9,10,12,13,14,15,16,17,20,21,22,23,24,25,32],ndim:32,ne:14,necessari:[15,28],need:[8,9,12,13,14,15,25,28,31,32],neg:[0,1,3,8,23,25,31],neither:[0,1],nelder:32,net:23,next:32,nfev:32,nice:[13,31,32],nicer:15,node:14,non:[14,31],none:[0,1,3,8,9,12,13,15,16,17,18,19,20,21,22,23,24,25],normal:[13,25],north:[7,14,23],north_project:7,notat:14,note:[0,1,3,4,7,8,11,12,13,14,15,16,17,18,24,30,31,32],notebook:30,noth:13,now:[30,31,32],np:[7,8,9,10,12,13,14,15,16,17,20,21,22,23,24,25,31,32],nu:14,number:[4,8,13,14,15,16,17,18,24,25,31],numer:20,numpi:[8,13,15,16,17,24,31,32],o:28,ob03235:30,ob03235_moa:30,ob03235_ogl:30,ob05086:31,ob08092:32,object:[0,1,3,4,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25],observ:[16,18,31],observatori:13,obtain:12,obviou:24,obvious:31,oct:25,off:[13,31],offset:23,ogl:[20,30,31],ogle_data:30,ok:[31,32],omega:14,onc:13,one:[8,9,13,14,15,16,18,31,32],ones:[8,25],onli:[4,8,9,13,14,15,16,17,18,20,24,31,32],onto:[16,18],op:32,open:28,opposit:[14,15],optim:[8,24,30,32],option:[0,1,4,8,9,12,13,15,22,23],orbit:[13,14],order:[0,1,3,8,13,14,15,24,28,32],org:25,origin:[0,1,3,4,23],os:[30,31,32],other:[0,1,3,4,8,9,11,12,13,15,24,25,31,32],otherwis:[8,9,25],our:32,out:[25,31],out_text:25,outlin:4,output:[10,23,25,32],outsid:12,over:[8,13,15],overal:9,overestim:[0,1],p:[24,31],packag:[20,25,27,28,29,31,32],paczynski:[12,20],page:[25,27,28],pair:[7,24],panel:24,paper:[0,1,24],parallax:[12,13,14,15,16,17,18,19,23,27],parallax_delta_n_:23,parallel:14,param:[14,31],paramet:[0,1,3,4,7,8,9,11,12,13,14,15,16,17,18,20,21,22,23,24,25,27,30,31],parameter:14,parameteris:24,parameters_to_fit:[31,32],pariti:[0,1,3],part:[9,25,28],partial:9,particular:[8,14],pass:[8,12,13,15,16,17,18,28],passag:14,path:[15,22,28,30,31,32],pc:[16,19],peak:32,peirson:3,per:24,percentil:31,perform:9,periapsi:14,perihelion:14,period:[14,15],perpendicular:[14,25],phot_fmt:[8,9,13,15,25,30],phot_ob08092_o4:32,photometr:[9,15],photometri:[13,15],photometry_fil:[30,31,32],physic:[16,18,19],pi_:[16,19],pi_e:[14,16,18,31],pi_e_:[8,9,14,31],pi_e_mag:14,pi_e_n:[8,9,14,31],pi_l:[16,17],pi_rel:[16,18],piec:[31,32],pip:28,place:25,plane:[4,7,14,16,18,23,25],planet:30,planetari:24,pleas:[3,28,30],plot:[4,8,9,13,15,16,17,18,24,25,30],plot_caust:[13,16,17,18],plot_data:[8,13,15,30,31,32],plot_lc:[8,13,15],plot_magnif:[13,16,18,30],plot_model:[8,30,31,32],plot_param:31,plot_properti:15,plot_residu:[8,13,15],plot_sourc:[8,13],plot_source_for_dataset:[8,13],plot_trajectori:[8,13],plotutil:25,plt:[4,13,24,30,31,32],plu:21,point:[0,1,3,4,8,9,12,13,14,15,16,17,18,20,21,24,25,30,31,32],point_sourc:12,point_source_magnif:[0,1,3],point_source_point_len:12,pointlen:[0,21,29],pointlenswithshear:[0,29],pointlik:20,polynomi:[0,1,3],posit:[0,1,3,8,10,12,13,14,15,20,23,24,31],possibl:[0,1,3,4,13,15],powel:32,pre:12,preced:8,precis:20,prepar:10,prescript:12,present:[0,1,20,32],previou:[8,13],previous:[8,9,13],primari:[0,1,3,14,15,16,17],print:[7,8,11,13,14,15,30,31,32],print_funct:30,prior:31,probabl:31,problem:[15,28],procedur:32,produc:[5,9,24],product:25,profil:[0,1,28],project:[7,14,16,17,18,23,25],proper:[14,16,18,20],properti:[4,7,8,9,10,12,13,14,15,16,17,18,19,23,24,31],propos:20,provid:[0,1,8,13,14,15,20,22,24,31,32],psbl:14,pspl:[9,14,20],pspl_magnif:[20,21],pspl_model:32,put:24,py:[8,9,13,28],pyplot:[4,13,15,16,17,18,24,25,30,31,32],python2:30,python3:30,python:[28,31],pythonpath:28,q:[0,1,3,4,6,14,16,17,23,24,25,30],q_sourc:14,quadrupol:[0,1,12,13],qualiti:13,quantiti:[14,16,17,18,19,25],quit:24,r:[28,31],r_e:[16,18],r_e_tild:[16,18],ra:[7,10,13,15],rad:14,radian:30,radiu:[0,1,13,14,16,18,19,20],rais:24,randn:31,random:31,randomli:24,rang:[8,12,13,14,15,22,24,30,31],rapper:28,rare:13,rate:14,ratio:[4,6,8,9,13,14,16,17,24,25,32],re:[0,1,3,9,28,30,31],read:[4,10,16,17,18],readi:31,real:24,recalcul:[8,9],recent:28,recommend:[0,1,3],rectilinear:31,red:[13,30,32],refer:[8,9,10,13,14,31],refit:9,refsdal:[5,21],regard:20,regardless:14,regress:13,reject:31,rel:[0,1,9,11,12,14,16,17,18,20,24,25,31],relat:[0,1,25,32],releas:28,relev:8,rememb:[8,11,13,31],remov:13,repeat:[24,28],report:[9,31],repositori:28,repres:[7,25],represent:13,request:[0,1,7,8,9,13,25],requir:[8,9,15,23,24,28],rescal:9,reset:[8,13],reset_best_chi2:8,reset_plot_properti:13,reshap:31,residu:[8,9,13,15,31],reson:24,respect:[8,9,13,15],result:[0,1,8,30,31,32],result_fun:32,retriev:9,rho:[0,1,9,12,13,14,20,30],rho_1:14,rho_2:14,right:[0,1,3,24,28],ring:[0,1,3,4,6,14,16,17,18,24],rough:32,round:31,routin:[24,32],run:[8,9,12,24,28],run_mcmc:31,s:[0,1,3,4,6,7,13,14,16,17,18,20,24,25,28,30,31,32],sai:31,same:[8,12,13,15,20,24,25,32],sampl:[24,31],sample_file_01:15,satellit:[10,12,13,15,22,23],satellite_skycoord:[12,13,15,22,23],satelliteskycoord:[0,13,29],save:[9,15],scale:[4,8,9,13,15,16,17,25,30],scale_errorbar:15,scale_flux:9,scatter:[4,8,13,15,16,17,18,24],scheme:24,scipi:32,search:[27,31],sec:[0,1],second:[14,15,25,31,32],secondari:[0,1,3,14,16,17],see:[0,1,8,9,10,12,13,14,15,16,17,18,20,22,23,24,27,30,31,32],seem:32,select:24,self:[8,13,15],semi:14,separ:[0,1,3,4,6,8,9,12,13,14,16,17,20,23,24,25,31],sequenc:[0,1,32],set:[0,1,7,8,9,11,12,13,14,15,16,17,18,19,23,25,27,30,32],set_aspect:[4,13],set_dataset:13,set_default_magnification_method:13,set_limb_coeff_gamma:[11,13],set_limb_coeff_u:[11,13],set_limb_darkening_weight:15,set_magnification_method:[12,13,30],set_magnification_methods_paramet:[12,13],set_source_flux_ratio:13,set_source_flux_ratio_for_band:13,set_tim:13,setattr:[14,31,32],setenv:28,setup:28,shape:14,share:28,shear:[3,5,6,12,14,21],shear_g:[3,5,6,14],sheat:[3,21],sheet:[3,14],shift:[14,23],should:[0,1,8,9,12,13,15,16,19,23,24,25,28,32],shouldn:15,show:[8,13,15,24,30,31,32],show_bad:[8,13,15],show_data:13,show_errorbar:[8,13,15],shown:[8,13],sigma:31,sign:14,signific:13,significantli:[0,1],similar:[15,24,25,28,32],simpl:30,simpli:14,simplifi:20,simpson:[0,1],sinc:30,singl:[5,8,9,13,16,17,21,25,31,32],size:[0,1,13,14,15,16,17,18,20,24,25],size_list:13,skip:28,skowron:[14,23],sky:[7,12,13,14,15,23,25],sky_coordin:7,skycoord:[7,8,12,13,15,22,23],slightli:[8,13,31,32],slow:[12,20],slower:[8,31],small:[0,1,20,25,31],smaller:[12,13,24],smallest:[8,31,32],smooth:24,so:[0,1,8,13,23,24,28,30,32],solar:[16,17],solmass:[16,17,18],solut:[0,1,3,31],some:[0,1,3,8,24,27,30,31,32],someth:[8,9],somewhat:30,sonam:28,sourc:[0,1,3,8,9,12,13,14,16,18,20,21,23,24,28,30,31,32],source_1_paramet:14,source_2_paramet:14,source_flux:[8,9,13],source_flux_0:9,source_flux_1:[9,13],source_flux_2:13,source_flux_ratio:[9,13],source_i:[0,1,3],source_x:[0,1,3],space:[8,9,11,13,15,24,25],special:15,specif:[13,14,24],specifi:[7,8,9,12,13,15,16,17,18,19,20,21,23,24,30,31,32],spheric:13,spitzer_ephemeris_01:22,split:31,squar:[0,1,4,9,13,15],stabl:20,stamp:23,standard:[10,12,14,16,17,24],star:[16,18,19,20,31],starblg234:31,start:[0,1,3,8,13,24,30,31,32],start_1:31,startup:28,state:[0,1],step:[28,31,32],still:32,stoke:12,stop:13,store:[4,13,16,17,32],str:[7,8,9,10,11,12,13,14,15,22,23,25,32],string:[8,15,25],structur:[4,5,6,13,16,18],studi:11,subdirectori:28,submodul:29,subpackag:[29,32],subtract:[8,13,14,15,25],subtract_2450000:[8,13,15,25,30,31],subtract_2460000:[8,13,15,25],success:[13,32],suitabl:[0,1,13],sum:[0,1,8,9,13,14,16,17,23,25],sum_funct:8,sun:[13,31],suppli:15,support:14,suppos:[25,30],suppress:15,sure:[0,1,3,8,13,15,28],symbol:[13,24],symmetr:24,system:[0,1,3,5,6,9,13,14,15,16,17,18,23],systemat:24,t:[8,12,15,24,28,30,31],t_0:[8,9,13,14,24,30,31,32],t_0_1:[13,14],t_0_2:[13,14],t_0_kep:[13,14],t_0_par:[7,14,31],t_0_xi:14,t_1:31,t_2:31,t_caustic_in:[14,24],t_caustic_out:[14,24],t_e:[8,9,13,14,16,18,24,30,31,32],t_e_1:14,t_e_2:14,t_eff:[8,9,14],t_rang:[8,13,30],t_star:14,t_star_1:14,t_star_2:14,t_start:[13,31],t_stop:[13,31],tabl:12,take:[8,9,24,31],taken:[8,13,15,23],tbd:[15,16,17],tbl:30,tcsh:28,tdb:10,termin:28,tertiari:[16,17],test:[0,1,28],text:25,th:[8,9],than:[0,1,3,9,11,12,13,15,16,17,18,24,31,32],thei:[8,9,12,13,15,24,31],them:[0,1,8,15,31,32],theorem:[0,1,12],theta:[31,32],theta_:[14,16,18,25],thi:[0,1,3,4,7,8,9,11,12,13,14,15,16,17,20,21,23,24,25,28,30,31,32],thing:[30,32],third:32,those:[8,13,30],three:[0,1,15,31],thu:24,time:[0,1,8,10,12,13,14,15,16,18,22,23,24,25,28,30,31,32],timescal:[14,31],titl:[8,31],tnc:32,togeth:[30,32],tolist:32,topocentr:[12,13,23],topolog:24,total:[0,1,3,9,11,15,16,17],total_mass:[16,17],toward:7,trajectori:[0,8,9,12,13,14,15,20,21,24,29],transform:[13,25],translat:15,trend:31,tupl:[8,13],turn:[13,31,32],tutori:27,two:[0,1,3,4,9,13,14,15,24,30,31,32],txt:[28,30],type:[9,11,13,14,15,20,23,25,32],typic:15,u:[0,1,7,11,12,13,14,16,19,20,25],u_0:[8,9,14,16,18,24,30,31,32],u_0_1:14,u_0_2:14,u_limb_darken:[0,1],u_to_gamma:25,uncertainti:[0,1,9,15,25],under:24,underli:12,uniform:[12,20,24],uniform_caustic_sampl:14,uniformcausticsampl:[0,14,29],uniformli:24,unit:[7,13,14,15,16,17,18,19,28],unless:[13,15],unnorm:31,unpack:28,up:[22,30],updat:[8,9,13,28],update_caust:13,upper:24,us:[0,1,3,4,8,9,11,12,13,14,15,16,17,18,20,22,23,24,25,28,30,31,32],usag:[14,15],user:13,user_nam:28,util:[0,13,29],v:[11,15],v_earth_perp_:7,v_earth_perp_n:7,v_earth_project:[7,25],valerio:[0,1,12],valid:14,valu:[0,1,7,8,9,11,12,13,14,15,16,17,19,20,23,24,25,27,31,32],valueerror:24,variat:[0,1],vastli:[0,1],vbbinarylens:[12,20],vbbinarylensinglibrari:28,vbbinarylensinglibrary_wrapp:28,vbbl:[0,1,12,13,28,30],vbbl_magnif:[0,1,12],vbbl_on:3,ve:28,vector:[4,12,13,14,15,23,25,31],vector_1:25,vector_2:25,vector_product_norm:25,veloc:[7,14,25],velocity_of_earth:25,veri:[0,1,12,13,20,24,25,31,32],version:[0,25,28,29],via:[13,32],visibl:15,visual:24,vs:15,wa:[8,9,13,20,21,31,32],wai:[13,24],walker:31,want:[8,9,13,15,16,17,25,30,31,32],we:[8,13,24,25,30,31,32],websit:[0,1,31],weight:[11,15],well:[9,12,13,20,31],were:[15,30],what:[0,1,3,12,23],when:[8,11,12,13,15,24],where:[7,13,24,32],whether:[8,9,13,14,15],which:[7,8,9,12,13,14,15,16,17,22,23,24,31,32],which_caust:24,which_paramet:14,wide:24,width:13,wiki:25,wikipedia:25,window:28,without:[9,12,13,14,15,16,17,31],witt:[12,20],wl:28,won:15,work:[12,15,20,24,28,30],would:[9,32],written:30,wrong:31,wyrzykowski:31,x0:32,x:[0,1,3,4,8,10,13,15,16,17,23,24,30,32],x_caustic:24,x_caustic_in:[14,24],x_caustic_out:[14,24],x_caustic_point:24,x_caustic_sep:24,xallarap:[13,14],xallarap_reference_posit:14,xi_argument_of_latitude_refer:14,xi_eccentr:14,xi_inclin:14,xi_omega_nod:14,xi_omega_periapsi:14,xi_period:14,xi_semimajor_axi:14,xlabel:25,xlim:[30,31],xy:23,xyz:10,y:[0,1,3,4,10,16,17,23,24,25],year:[14,31],yerr:32,yet:[8,13,24,31],ylim:30,yoo:[12,20],you:[0,1,3,7,8,9,11,12,13,14,15,16,17,24,25,28,30,31,32],your:[24,28,30,32],yr:[14,16,18],z:10,zero:[8,14,24,25],zeropoint:25,zip:31,zoom:[31,32],zorder:15,zorder_list:13},titles:["MulensModel package","MulensModel.binarylens module","MulensModel.binarylensimports module","MulensModel.binarylenswithshear module","MulensModel.caustics module","MulensModel.causticspointwithshear module","MulensModel.causticswithshear module","MulensModel.coordinates module","MulensModel.event module","MulensModel.fitdata module","MulensModel.horizons module","MulensModel.limbdarkeningcoeffs module","MulensModel.magnificationcurve module","MulensModel.model module","MulensModel.modelparameters module","MulensModel.mulensdata module","MulensModel.mulensobjects package","MulensModel.mulensobjects.lens module","MulensModel.mulensobjects.mulenssystem module","MulensModel.mulensobjects.source module","MulensModel.pointlens module","MulensModel.pointlenswithshear module","MulensModel.satelliteskycoord module","MulensModel.trajectory module","MulensModel.uniformcausticsampling module","MulensModel.utils module","MulensModel.version module","Welcome to MulensModel\u2019s documentation!","How to install?","MulensModel","Very Basic Tutorial","Microlensing Parallax Fitting Tutorial","Basic Fitting Tutorial"],titleterms:{"import":[31,32],basic:[30,32],binarylen:1,binarylensimport:2,binarylenswithshear:3,caustic:4,causticspointwithshear:5,causticswithshear:6,combin:30,content:[0,16,27],coordin:7,data:[30,32],defin:30,document:27,event:8,exercis:[31,32],fit:[31,32],fitdata:9,horizon:10,how:28,indic:27,instal:28,introduc:30,len:17,limbdarkeningcoeff:11,magnificationcurv:12,microlens:31,model:[13,30,32],modelparamet:14,modul:[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26],mulensdata:15,mulensmodel:[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,29],mulensobject:[16,17,18,19],mulenssystem:18,old:28,packag:[0,16],parallax:31,paramet:32,plot:[31,32],pointlen:20,pointlenswithshear:21,prepar:32,run:31,s:27,sampler:31,satelliteskycoord:22,set:31,sourc:19,submodul:[0,16],subpackag:0,tabl:27,trajectori:23,tutori:[30,31,32],uniformcausticsampl:24,util:25,veri:30,version:26,welcom:27}}) \ No newline at end of file diff --git a/docs/tutorial.html b/docs/tutorial.html index 30a940b1..a4f1cf07 100644 --- a/docs/tutorial.html +++ b/docs/tutorial.html @@ -6,7 +6,7 @@ - Very Basic Tutorial — MulensModel 2.22.0 documentation + Very Basic Tutorial — MulensModel 2.23.0 documentation diff --git a/docs/tutorial_fit_pi_E.html b/docs/tutorial_fit_pi_E.html index 9138a82a..de35ecc5 100644 --- a/docs/tutorial_fit_pi_E.html +++ b/docs/tutorial_fit_pi_E.html @@ -6,7 +6,7 @@ - Microlensing Parallax Fitting Tutorial — MulensModel 2.22.0 documentation + Microlensing Parallax Fitting Tutorial — MulensModel 2.23.0 documentation diff --git a/docs/tutorial_fit_pspl.html b/docs/tutorial_fit_pspl.html index fff6543f..39cc1b10 100644 --- a/docs/tutorial_fit_pspl.html +++ b/docs/tutorial_fit_pspl.html @@ -6,7 +6,7 @@ - Basic Fitting Tutorial — MulensModel 2.22.0 documentation + Basic Fitting Tutorial — MulensModel 2.23.0 documentation diff --git a/documents/parameter_names.pdf b/documents/parameter_names.pdf index 85c3e26f..fd1d35a5 100644 Binary files a/documents/parameter_names.pdf and b/documents/parameter_names.pdf differ diff --git a/documents/parameter_names.tex b/documents/parameter_names.tex index 824c74ff..70cbb08e 100644 --- a/documents/parameter_names.tex +++ b/documents/parameter_names.tex @@ -24,7 +24,7 @@ {\LARGE Microlensing parameters in \MM}\\ \bigskip Radek Poleski\\ -last update: Jul 2023 +last update: Apr 2024 \end{center} \bigskip @@ -63,6 +63,7 @@ $\xi_u$ & \texttt{xi\_argument\_of\_latitude\_reference} & deg & The argument of latitude at the reference epoch ($t_{0,\chi}$). The argument of latitude is a sum of true anomaly ($\nu$, changes with time) and the argument of periapsis ($\omega$, orbit parameter, i.e., does not change with time): $u = \nu + \omega$.$^c$\\ $\xi_e$ & \texttt{xi\_eccentricity} & & The eccentricity of a xallarap orbit.\\ $\xi_\omega$ & \texttt{xi\_omega\_periapsis} & deg & The argument of periapsis of a xallarap orbit.$^c$\\ +$q_\mathrm{source}$ & \texttt{q\_source} & & Mass ratio of source components: $m_{s,2}/m_{s,1}$. It is valid only for xallarap models.\\ $t_{0,\xi}$ & \texttt{t\_0\_xi} & & The reference epoch for parameters in xallarap models.$^a$\\ \hline \caption{Notes: \newline diff --git a/source/MulensModel/modelparameters.py b/source/MulensModel/modelparameters.py index 7cae8618..bdc67673 100644 --- a/source/MulensModel/modelparameters.py +++ b/source/MulensModel/modelparameters.py @@ -3,6 +3,7 @@ import warnings from MulensModel.uniformcausticsampling import UniformCausticSampling +from MulensModel.orbits.orbit import Orbit # For definition of class ModelParameters see below. @@ -217,60 +218,137 @@ def __init__(self, parameters): "as a parameter\ne.g., ModelParameters({'t_0': " + "2456789.0, 'u_0': 0.123, 't_E': 23.45})") - self._count_sources(parameters.keys()) - self._count_lenses(parameters.keys()) self._set_type(parameters.keys()) + self._count_sources(set(parameters.keys())) + self._count_lenses(parameters.keys()) self._check_types('alpha' in parameters.keys()) if self.n_sources == 1: - self._check_valid_combination_1_source(parameters.keys()) - if self._type['Cassan08']: - self._uniform_caustic = None - self._standard_parameters = None + self._init_1_source(parameters) elif self.n_sources == 2: - self._check_valid_combination_2_sources(parameters.keys()) - if 't_E' not in parameters.keys(): - raise KeyError('Currently, the binary source calculations ' + - 'require t_E to be directly defined, i.e., ' + - 'has to be the same for both sources.') - (params_1, params_2) = self._divide_parameters(parameters) - try: - self._source_1_parameters = ModelParameters(params_1) - except Exception: - print("ERROR IN ITIALIZING SOURCE 1") - raise - try: - self._source_2_parameters = ModelParameters(params_2) - except Exception: - print("ERROR IN ITIALIZING SOURCE 2") - raise - # The block above forces checks from "== 1" block above to be - # run on each source parameters separately. + self._init_2_sources(parameters) else: raise ValueError('wrong number of sources') + self._set_parameters(parameters) + def _init_1_source(self, parameters): + """ + initialize model with 1 source + """ + self._check_valid_combination_1_source(parameters.keys()) + if self._type['Cassan08']: + self._uniform_caustic = None + self._standard_parameters = None + if self.is_xallarap: + delta_1 = self._get_xallarap_position(parameters) + self._xallarap_reference_position = delta_1 + + def _init_2_sources(self, parameters): + """ + initialize model with s sources + """ + self._check_valid_combination_2_sources(parameters.keys()) + if 't_E' not in parameters.keys(): + raise KeyError('Currently, the binary source calculations ' + + 'require t_E to be directly defined, i.e., ' + + 'has to be the same for both sources.') + + (params_1, params_2) = self._divide_parameters(parameters) + + try: + self._source_1_parameters = ModelParameters(params_1) + except Exception: + print("ERROR IN ITIALIZING SOURCE 1") + raise + try: + self._source_2_parameters = ModelParameters(params_2) + except Exception: + print("ERROR IN ITIALIZING SOURCE 2") + raise + # The try/except blocks above force checks from ._init_1_source() + # to be run on each source parameters separately. + + if self.is_xallarap: + self._update_sources_xallarap_reference() + + def _update_sources_xallarap_reference(self): + """ + Update .xallarap_reference_position for each source parameters + + Note: below we're calling private function and set private + properties NOT of self, but self._source_X_parameters, + which both are of the same type as self. + """ + delta_1 = self._source_1_parameters._get_xallarap_position() + self._source_1_parameters._xallarap_reference_position = delta_1 + self._source_2_parameters._xallarap_reference_position = delta_1 + + def _get_xallarap_position(self, parameters=None): + """ + Get position at t_0_xi from xallarap Orbit object. + + Note: this function is called in 2 different ways: + - directly, i.e., self._get_xallarap_orbit(), and + - indirectly, i.e., self._source_1_parameters._get_xallarap_orbit(). + """ + if parameters is None: + parameters = self.parameters + t_0_xi = parameters.get('t_0_xi', parameters['t_0']) + + zip_ = parameters.items() + orbit_parameters = {key[3:]: value + for (key, value) in zip_ if key[:3] == "xi_"} + orbit_parameters['epoch_reference'] = t_0_xi + orbit = Orbit(**orbit_parameters) + return orbit.get_reference_plane_position([t_0_xi]) + def _count_sources(self, keys): - """How many sources there are?""" - binary_params = ['t_0_1', 't_0_2', 'u_0_1', 'u_0_2', 'rho_1', 'rho_2', - 't_star_1', 't_star_2'] - common = set(binary_params).intersection(set(keys)) - if len(common) == 0: + """ + How many luminous sources there are? + We're also checking for ill-defined xallarap with + *_1 and *_2 binary source parameters. + """ + finite_source_params = ['rho_1', 'rho_2', 't_star_1', 't_star_2'] + binary_params = ['t_0_1', 't_0_2', 'u_0_1', 'u_0_2'] + binary_params += finite_source_params + common = set(binary_params).intersection(keys) + finite_params = common.intersection(finite_source_params) + n_finite = len(finite_params) + + if len(common) == 0 and 'q_source' not in keys: self._n_sources = 1 + elif len(common) == 0 and 'q_source' in keys: + self._n_sources = 2 elif len(common) == 1: - raise ValueError('Wrong parameters - the only binary source ' + - 'parameter is {:}'.format(common)) + if self.is_xallarap and n_finite == 1: + self._n_sources = int(list(finite_params)[0][-1]) + else: + raise ValueError('Wrong parameters - the only binary ' + + 'source parameter is {:}'.format(common)) + elif len(common) == 2 and n_finite == 2 and 'q_source' in keys: + if common in [{'rho_1', 't_star_1'}, {'rho_2', 't_star_2'}]: + raise ValueError('source size overdefined: {:}'.format(common)) + else: + self._n_sources = 2 else: - common_no_1_2 = {param[:-2] for param in common} - condition_1 = (len(common_no_1_2) == len(common)) - condition_2 = not ( - 'rho' in common_no_1_2 and 't_star' in common_no_1_2) - if condition_1 and condition_2: - raise ValueError( - 'Given binary source parameters do not allow defining ' + - 'the Model: {:}'.format(common)) + self._check_for_underdefined_source(common) self._n_sources = 2 + def _check_for_underdefined_source(self, common): + """ + Make sure that finite source size in binary source model + is not underdefined. + """ + common_no_1_2 = {param[:-2] for param in common} + condition_1 = (len(common_no_1_2) == len(common)) + condition_2 = not ( + 'rho' in common_no_1_2 and 't_star' in common_no_1_2) + if condition_1 and condition_2: + raise ValueError( + 'Given binary source parameters do not allow defining ' + + 'the Model: {:}'.format(common)) + def _count_lenses(self, keys): """How many lenses there are?""" self._n_lenses = 1 @@ -294,7 +372,8 @@ def _set_type(self, keys): 'lens 2-parameter orbital motion': 'dalpha_dt ds_dt', 'mass sheet': 'convergence_K shear_G', 'xallarap': ('xi_period xi_semimajor_axis xi_inclination ' - 'xi_Omega_node xi_argument_of_latitude_reference')} + 'xi_Omega_node xi_argument_of_latitude_reference ' + 'xi_eccentricity xi_omega_periapsis q_source')} parameter_to_type = dict() for (key, values) in temp.items(): @@ -312,7 +391,6 @@ def _check_types(self, alpha_defined): Check if self._type values make sense """ n_lenses = self._n_lenses - n_sources = self._n_sources # Lens orbital motion requires binary lens: if self._type['lens 2-parameter orbital motion'] and n_lenses == 1: @@ -327,10 +405,6 @@ def _check_types(self, alpha_defined): 'You defined alpha for single lens model ' 'without external mass sheet. This is not allowed.') - if n_sources > 1 and self._type['xallarap']: - raise NotImplementedError('We have not yet implemented xallarap ' - 'and multiple luminous sources') - def _check_valid_combination_1_source(self, keys): """ Check that the user hasn't over-defined the ModelParameters. @@ -343,7 +417,7 @@ def _check_valid_combination_1_source(self, keys): 'x_caustic_in x_caustic_out t_caustic_in t_caustic_out ' 'xi_period xi_semimajor_axis xi_inclination xi_Omega_node ' 'xi_argument_of_latitude_reference xi_eccentricity ' - 'xi_omega_periapsis t_0_xi').split()) + 'xi_omega_periapsis t_0_xi q_source').split()) difference = set(keys) - allowed_keys if len(difference) > 0: derived_1 = ['gamma', 'gamma_perp', 'gamma_parallel'] @@ -388,6 +462,7 @@ def _divide_parameters(self, parameters): """ separate_parameters = ( 't_0_1 t_0_2 u_0_1 u_0_2 rho_1 rho_2 t_star_1 t_star_2'.split()) + skipped_parameters = ['q_source'] parameters_1 = {} parameters_2 = {} for (key, value) in parameters.items(): @@ -398,13 +473,54 @@ def _divide_parameters(self, parameters): parameters_2[key[:-2]] = value else: raise ValueError('unexpected error') + elif key in skipped_parameters: + continue else: parameters_1[key] = value parameters_2[key] = value + + if self.n_sources == 2 and self._type['xallarap']: + self._set_changed_parameters_2nd_source(parameters['q_source'], + parameters_2) + return (parameters_1, parameters_2) + def _set_changed_parameters_2nd_source(self, q_source, parameters_2): + """ + For xallarap model with 2 sources, the orbit of the second source + must have 2 parameters changed. + Functions starts with tests of input + """ + if q_source <= 0.: + raise ValueError('q_source cannot be negative') + + check_keys = ['xi_semimajor_axis', 'xi_argument_of_latitude_reference'] + for key in check_keys: + if key not in parameters_2: + raise KeyError('xallarap model with 2 sources requires ' + key) + + parameters_2['xi_semimajor_axis'] /= q_source + parameters_2['xi_argument_of_latitude_reference'] += 180. + def __repr__(self): """A nice way to represent a ModelParameters object as a string""" + out = self._get_main_parameters_to_print() + + if self.is_xallarap: + fmt = "\nxallarap reference position: ({:.4f}, {:.4f})" + if self.n_sources == 1: + source = self + else: + source = self._source_1_parameters + position = source.xallarap_reference_position + out += fmt.format(position[0, 0], position[1, 0]) + + return out + + def _get_main_parameters_to_print(self): + """ + prepare all the standard parameters to be printed + """ keys = self._get_keys_for_repr() formats = self._get_formats_dict_for_repr() ordered_keys = self._get_ordered_keys_for_repr() @@ -488,6 +604,7 @@ def _get_formats_dict_for_repr(self): 'name': 'xallarap eccentricity'}, 'xi_omega_periapsis': {'width': 11, 'precision': 5, 'unit': 'deg', 'name': 'xallarap omega periapsis'}, + 'q_source': {'width': 12, 'precision': 8}, 't_0_xi': {'width': 13, 'precision': 5, 'unit': 'HJD'}, } # Add binary source parameters with the same settings. @@ -516,7 +633,7 @@ def _get_ordered_keys_for_repr(self): 'x_caustic_in', 'x_caustic_out', 't_caustic_in', 't_caustic_out', 'xi_period', 'xi_semimajor_axis', 'xi_inclination', 'xi_Omega_node', 'xi_argument_of_latitude_reference', - 'xi_eccentricity', 'xi_omega_periapsis', 't_0_xi' + 'xi_eccentricity', 'xi_omega_periapsis', 'q_source', 't_0_xi' ] return ordered_keys @@ -546,15 +663,28 @@ def _get_formats_for_repr(self, form, full_name): def _check_valid_combination_2_sources(self, keys): """ - make sure that there is no conflict between t_0 and t_0_1 etc. + Make sure that there is no conflict between t_0 and t_0_1 etc. + Also make sure that xallarap is not mixed with t_0_1, u_0_1 etc. """ - binary_params = ( - 't_0_1 t_0_2 u_0_1 u_0_2 rho_1 rho_2 t_star_1 t_star_2'.split()) + binary_params_nonFS = 't_0_1 t_0_2 u_0_1 u_0_2'.split() + binary_params_FS_2 = ['rho_2', 't_star_2'] + binary_params_FS = ['rho_1', 't_star_1'] + binary_params_FS_2 + binary_params = binary_params_nonFS + binary_params_FS for parameter in binary_params: if (parameter in keys) and (parameter[:-2] in keys): raise ValueError('You cannot set {:} and {:}'.format( parameter, parameter[:-2])) + common = set(keys).intersection(binary_params_nonFS) + if self.is_xallarap and len(common) > 0: + msg = 'xallarap parameters cannot be mixed with {:}' + raise NotImplementedError(msg.format(common)) + + common = set(keys).intersection(binary_params_FS_2) + if self.is_xallarap and len(common) > 0 and 'q_source' not in keys: + raise KeyError('You cannot define xallarap model without ' + 'q_star but with rho_2 or t_star_2') + def _check_valid_combination_1_source_standard(self, keys): """ Here we check parameters for non-Cassan08 parameterization. @@ -792,9 +922,21 @@ def _update_sources(self, parameter, value): if parameter in self._source_1_parameters.parameters: setattr(self._source_1_parameters, parameter, value) + if parameter == 'xi_semimajor_axis': + value /= self.parameters['q_source'] + elif parameter == 'xi_argument_of_latitude_reference': + value += 180. + if parameter in self._source_2_parameters.parameters: setattr(self._source_2_parameters, parameter, value) + if parameter == 'q_source': + value_ = self.parameters['xi_semimajor_axis'] / value + setattr(self._source_2_parameters, 'xi_semimajor_axis', value_) + + if self.is_xallarap and self.n_sources > 1: + self._update_sources_xallarap_reference() + def _set_time_quantity(self, key, new_time): """ Save a variable with units of time (e.g. t_E, t_star, @@ -1442,6 +1584,7 @@ def xi_period(self, new_value): if new_value < 0.: raise ValueError('Xallarap period cannot be negative') self.parameters['xi_period'] = new_value + self._update_sources('xi_period', new_value) @property def xi_semimajor_axis(self): @@ -1457,6 +1600,7 @@ def xi_semimajor_axis(self, new_value): if new_value < 0.: raise ValueError('Xallarap semimajor axis cannot be negative') self.parameters['xi_semimajor_axis'] = new_value + self._update_sources('xi_semimajor_axis', new_value) @property def xi_Omega_node(self): @@ -1475,6 +1619,7 @@ def xi_Omega_node(self, new_value): self._warn_if_angle_outside_reasonable_range(new_value, 'xi_Omega_node') self.parameters['xi_Omega_node'] = new_value + self._update_sources('xi_Omega_node', new_value) @property def xi_inclination(self): @@ -1492,6 +1637,7 @@ def xi_inclination(self, new_value): self._warn_if_angle_outside_reasonable_range(new_value, 'xi_inclination') self.parameters['xi_inclination'] = new_value + self._update_sources('xi_inclination', new_value) @property def xi_argument_of_latitude_reference(self): @@ -1512,6 +1658,7 @@ def xi_argument_of_latitude_reference(self, new_value): self._warn_if_angle_outside_reasonable_range( new_value, 'xi_argument_of_latitude_reference') self.parameters['xi_argument_of_latitude_reference'] = new_value + self._update_sources('xi_argument_of_latitude_reference', new_value) @property def xi_eccentricity(self): @@ -1527,13 +1674,14 @@ def xi_eccentricity(self, new_value): if new_value < 0. or new_value >= 1.: raise ValueError('xallarap eccentricity has to be between 0 and 1') self.parameters['xi_eccentricity'] = new_value + self._update_sources('xi_eccentricity', new_value) @property def xi_omega_periapsis(self): """ *float* - The argument of periapsis of the xallrap orbit, i.e., the angle + The argument of periapsis of the xallarap orbit, i.e., the angle between the ascending node and periapsis measured in the direction of motion. The units are degrees. @@ -1545,6 +1693,7 @@ def xi_omega_periapsis(self, new_value): self._warn_if_angle_outside_reasonable_range( new_value, 'xi_omega_periapsis') self.parameters['xi_omega_periapsis'] = new_value + self._update_sources('xi_omega_periapsis', new_value) @property def t_0_xi(self): @@ -1562,6 +1711,38 @@ def t_0_xi(self): @t_0_xi.setter def t_0_xi(self, new_value): self.parameters['t_0_xi'] = new_value + self._update_sources('t_0_xi', new_value) + + @property + def q_source(self): + """ + *float* + + The mass ratio of the second and the first source. + This is value must be positive and can be > 1. + Defined only for xallarap binary-source models because it does not + affect the magnification for binary-source models without xallarap. + """ + return self.parameters['q_source'] + + @q_source.setter + def q_source(self, new_value): + if new_value < 0.: + raise ValueError('q_source cannot be negative') + self.parameters['q_source'] = new_value + self._update_sources('q_source', new_value) + + @property + def xallarap_reference_position(self): + """ + *np.ndarray* of shape (2, 1) + + The position of the first source at :py:attr:`~t_0_xi` relative to + the source center of mass. It is a 2D vector that is subtracted from + the source position along the orbit in order to calculate the shift + caused by xallarap. + """ + return self._xallarap_reference_position @property def t_0_1(self): diff --git a/source/MulensModel/orbits/orbit.py b/source/MulensModel/orbits/orbit.py index 1e62a009..3f0d0e5c 100644 --- a/source/MulensModel/orbits/orbit.py +++ b/source/MulensModel/orbits/orbit.py @@ -213,7 +213,7 @@ def _get_periapsis_epoch(self, u_reference, epoch_reference): the argument_of_latitude (u) at given epoch """ time_shift = self._period * u_reference / (2. * np.pi) - return epoch_reference - time_shift + return np.float64(epoch_reference - time_shift) def get_orbital_plane_position(self, time): """ diff --git a/source/MulensModel/tests/test_ModelParameters.py b/source/MulensModel/tests/test_ModelParameters.py index 8d9e9733..670b19a5 100644 --- a/source/MulensModel/tests/test_ModelParameters.py +++ b/source/MulensModel/tests/test_ModelParameters.py @@ -45,6 +45,8 @@ def test_init_for_2_sources(self): with self.assertRaises(KeyError): mm.ModelParameters({'t_01': 1, 'u_0_1': 0.1, 't_eff_1': 10, 't_0_2': 10., 'u_0_2': 0.01, 't_eff_2': 20.}) + with self.assertRaises(Exception): + mm.ModelParameters({'t_0_1': 1, 'u_0_2': 0.1, 't_E': 10}) def test_init_parameters(): @@ -174,12 +176,12 @@ def test_rho_t_e_t_star(): class test(unittest.TestCase): def test_too_much_rho_t_e_t_star(self): + t_0 = 2450000. + u_0 = 0.1 + t_E = 20. * u.day + rho = 0.001 + t_star = t_E * rho with self.assertRaises(KeyError): - t_0 = 2450000. - u_0 = 0.1 - t_E = 20. * u.day - rho = 0.001 - t_star = t_E * rho mm.ModelParameters({ 't_0': t_0, 'u_0': u_0, 't_E': t_E, 'rho': rho, 't_star': t_star}) @@ -458,7 +460,7 @@ def test_failing_single_lens_with_mass_sheet(self): xallarap_parameters = { - 't_0': 0, 't_E': 9., 'u_0': 0.1, 'xi_period': 12.345, + 't_0': 2., 't_E': 9., 'u_0': 0.1, 'xi_period': 12.345, 'xi_semimajor_axis': 0.54321, 'xi_Omega_node': 0.123, 'xi_inclination': 9.8765, 'xi_argument_of_latitude_reference': 24.68, 'xi_eccentricity': 0.5, 'xi_omega_periapsis': 12.3456, 't_0_xi': 1.} @@ -468,7 +470,7 @@ def setup_xallarap(key): """ Setup for xallarap tests. """ - model = mm.ModelParameters(xallarap_parameters) + model = mm.ModelParameters({**xallarap_parameters}) return (model, xallarap_parameters[key]) @@ -510,8 +512,6 @@ def test_missing_xallarap_parameters(self): parameters.pop(parameter) with self.assertRaises(KeyError): mm.ModelParameters(parameters) - print("Test failed (i.e. KeyError was not raised) for ", - parameter) def test_negative_period(self): """ @@ -558,6 +558,89 @@ def test_xallarap_and_Cassan08(self): with self.assertRaises(NotImplementedError): mm.ModelParameters(parameters) + def test_negative_source_mass_ratio_1(self): + """ + q_source must be positive in __init__() + """ + parameters = {**xallarap_parameters, 'q_source': -0.12345} + with self.assertRaises(ValueError): + _ = mm.ModelParameters(parameters) + + def test_negative_source_mass_ratio_2(self): + """ + q_source must be positive + """ + parameters = {**xallarap_parameters, 'q_source': 0.12345} + model = mm.ModelParameters(parameters) + with self.assertRaises(ValueError): + setattr(model, 'q_source', -0.12345) + + def test_overdefined_source_size(self): + """ + overdefine first sourece size + """ + parameters = {**xallarap_parameters, + 'rho_1': 0.1, 't_star_1': 0.1, 'q_source': 1.0} + with self.assertRaises(ValueError): + mm.ModelParameters(parameters) + + def test_mixed_binary_source_and_xallarap(self): + """ + Xallarap cannot be combined with just + one stardard binary source parameter. + """ + parameters = {**xallarap_parameters, 'u_0_2': 0.123} + with self.assertRaises(ValueError): + mm.ModelParameters(parameters) + + def test_no_q_source_but_with_rho_2(self): + """ + xallarap model without q_source cannot have rho_2 + """ + parameters = {**xallarap_parameters, 'rho_2': 0.1} + with self.assertRaises(KeyError): + mm.ModelParameters(parameters) + + def test_PSPL_and_q_source(self): + """ + Make sure one cannot provide only PSPL parameters and q_source. + """ + with self.assertRaises(KeyError): + mm.ModelParameters({'t_0': 1, 'u_0': 2, 't_E': 3, 'q_source': 1}) + + +def test_print_xallarap(): + """ + Test if printing of printing of xallarap model works as expected. + """ + model = mm.ModelParameters(xallarap_parameters) + expected = ( + " t_0 (HJD) u_0 t_E (d) xallarap period (d) xallarap " + "semimajor axis xallarap inclination (deg) xallarap Omega node (deg) " + "xallarap argument of latitude reference (deg) xallarap eccentricity " + "xallarap omega periapsis (deg) t_0_xi (HJD) " + "\n 2.00000 0.100000 9.0000 12.3450 " + " 0.543210 9.87650 0.12300 " + " 24.68000 0.500000 " + " 12.34560 1.00000 " + "\nxallarap reference position: (0.2673, 0.0582)" + ) + assert model.__repr__() == expected + + +def test_print_xallarap_with_q_source(): + """ + Test if printing of printing of xallarap model with q_source works + as expected. Most stuff was tested in test_print_xallarap(), so we + check only the parts that are important here. + """ + parameters = {**xallarap_parameters, 'q_source': 0.12345} + model = mm.ModelParameters(parameters) + lines = model.__repr__().split("\n") + assert lines[0][-27:] == " q_source t_0_xi (HJD) " + assert lines[1][-27:] == " 0.12345000 1.00000 " + assert lines[2] == "xallarap reference position: (0.2673, 0.0582)" + @pytest.mark.parametrize( "parameter", @@ -601,3 +684,102 @@ def test_is_xallarap_2(): parameters = {'t_0_1': 0, 'u_0_1': 1, 't_0_2': 5, 'u_0_2': 0.1, 't_E': 9} model_params = mm.ModelParameters(parameters) assert not model_params.is_xallarap + + +def test_xallarap_n_sources(): + """ + Make sure that number of sources in xallarap models is properly calculated + """ + parameters = {**xallarap_parameters} + model_1S = mm.ModelParameters(parameters) + assert model_1S.n_sources == 1 + + parameters['q_source'] = 1. + model_2S = mm.ModelParameters(parameters) + assert model_2S.n_sources == 2 + + parameters['rho_1'] = 0.1 + parameters['rho_2'] = 0.2 + model_3 = mm.ModelParameters(parameters) + assert model_3.n_sources == 2 + + parameters = {**xallarap_parameters, 't_star_1': 2.} + model_4 = mm.ModelParameters(parameters) + assert model_4.n_sources == 1 + + +def test_2S1L_xallarap_individual_source_parameters(): + """ + Make sure that parameters of both sources are properly set. + Most importantly, xi_u is shifted by 180 deg and xi_a is scaled by + q_source. + """ + q_source = 1.23456 + parameters_1st = {**xallarap_parameters} + + parameters_2nd = {**parameters_1st} + parameters_2nd['xi_semimajor_axis'] /= q_source + parameters_2nd['xi_argument_of_latitude_reference'] += 180. + + parameters = {'q_source': q_source, **parameters_1st} + model = mm.ModelParameters(parameters) + check_1st = model.source_1_parameters.as_dict() + check_1st['t_E'] = check_1st['t_E'].value + check_2nd = model.source_2_parameters.as_dict() + check_2nd['t_E'] = check_2nd['t_E'].value + + assert check_1st == parameters_1st + assert check_2nd == parameters_2nd + + +tested_keys_3 = tested_keys_2 + ['q_source'] + + +@pytest.mark.parametrize("key", tested_keys_3) +def test_changes_of_xallrap_parameters_for_both_sources(key): + """ + Make sure that chainging a xallarap parameter in a binary source event + with binary sources model properly changes parameters of each parameter. + For q_source make sure that it actually it's not passed to the parameters + of each source. + """ + q_source = 1.23456 + factor = 1.1 + parameters = {'q_source': q_source, **xallarap_parameters} + model = mm.ModelParameters(parameters) + old_value = getattr(model, key) + new_value = factor * old_value + setattr(model, key, new_value) + + assert getattr(model, key) == new_value + + if key == 'q_source': + assert 'q_source' not in model.source_1_parameters.parameters + assert 'q_source' not in model.source_2_parameters.parameters + key_a = 'xi_semimajor_axis' + xi_a = xallarap_parameters[key_a] + assert model.source_1_parameters.parameters[key_a] == xi_a + assert model.source_2_parameters.parameters[key_a] == xi_a / new_value + return + + assert getattr(model.source_1_parameters, key) == new_value + + new_value_2 = new_value + if key == 'xi_argument_of_latitude_reference': + new_value_2 += 180. + elif key == 'xi_semimajor_axis': + new_value_2 /= q_source + + assert getattr(model.source_2_parameters, key) == new_value_2 + + +def test_reference_position(): + """ + Make sure that for xallarap model bith sources have the same + reference position. + """ + parameters = {**xallarap_parameters, 'q_source': 0.12345} + model = mm.ModelParameters(parameters) + text_1 = model.source_1_parameters.__repr__().split("\n")[-1] + text_2 = model.source_2_parameters.__repr__().split("\n")[-1] + assert text_1 == text_2 diff --git a/source/MulensModel/trajectory.py b/source/MulensModel/trajectory.py index ddce3d10..718cd106 100644 --- a/source/MulensModel/trajectory.py +++ b/source/MulensModel/trajectory.py @@ -338,6 +338,5 @@ def _get_shifts_xallarap(self): key[3:]: value for (key, value) in zip_ if key[:3] == "xi_"} orbit_parameters['epoch_reference'] = t_0_xi orbit = Orbit(**orbit_parameters) - reference_position = orbit.get_reference_plane_position(t_0_xi) positions = orbit.get_reference_plane_position(self.times) - return positions - reference_position.reshape((2, 1)) + return positions - self.parameters.xallarap_reference_position diff --git a/source/MulensModel/version.py b/source/MulensModel/version.py index 4d8af451..db969b3a 100644 --- a/source/MulensModel/version.py +++ b/source/MulensModel/version.py @@ -1 +1 @@ -__version__ = "2.22.2" +__version__ = "2.23.0"