PSPL Model Classes
PSPL
- class model.PSPL_PhotAstrom_noPar_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_noParallax
,PSPL_PhotAstromParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])noParallax: Position of the observed source position in arcsec.
get_astrometry_unlensed
(t_obs)noParallax: Get the astrometry of the source if the lens didn't exist.
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_lens_astrometry
(t_obs)Equation of motion for just the foreground lens.
Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Get the x, y astrometry for each of the two source images, which we label plus and minus.
calc_piE_ecliptic
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t_obs)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
PSPL Parallax
- class model.PSPL_PhotAstrom_Par_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_PhotAstromParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_PhotAstrom_LumLens_Par_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_PhotAstromParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_PhotAstrom_LumLens_Par_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_PhotAstromParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_PhotAstrom_LumLens_Par_Param4(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_PhotAstromParam4
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
PSPL_parallax2 / PSPL_multiphot_parallax
- class model.PSPL_PhotAstrom_Par_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_PhotAstromParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_PhotAstrom_noPar_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_noParallax
,PSPL_PhotAstromParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])noParallax: Position of the observed source position in arcsec.
get_astrometry_unlensed
(t_obs)noParallax: Get the astrometry of the source if the lens didn't exist.
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_lens_astrometry
(t_obs)Equation of motion for just the foreground lens.
Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Get the x, y astrometry for each of the two source images, which we label plus and minus.
calc_piE_ecliptic
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t_obs)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- class model.PSPL_PhotAstrom_Par_Param3(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_PhotAstromParam3
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_PhotAstrom_Par_Param4(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_PhotAstromParam4
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_PhotAstrom_noPar_Param4(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_PhotAstrom
,PSPL_noParallax
,PSPL_PhotAstromParam4
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])noParallax: Position of the observed source position in arcsec.
get_astrometry_unlensed
(t_obs)noParallax: Get the astrometry of the source if the lens didn't exist.
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_lens_astrometry
(t_obs)Equation of motion for just the foreground lens.
Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Get the x, y astrometry for each of the two source images, which we label plus and minus.
calc_piE_ecliptic
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t_obs)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- class model.PSPL_Astrom_Par_Param4(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_Astrom
,PSPL_Parallax
,PSPL_AstromParam4
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_Astrom_Par_Param3(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_Astrom
,PSPL_Parallax
,PSPL_AstromParam3
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
PSPL_phot
- class model.PSPL_Phot_noPar_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_Phot
,PSPL_noParallax
,PSPL_PhotParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
calc_piE_ecliptic
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- class model.PSPL_Phot_noPar_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_Phot
,PSPL_noParallax
,PSPL_PhotParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
calc_piE_ecliptic
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
PSPL_phot_parallax / PSPL_phot_multiphot_parallax
- class model.PSPL_Phot_Par_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_Phot
,PSPL_Parallax
,PSPL_PhotParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- class model.PSPL_Phot_Par_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_Phot
,PSPL_Parallax
,PSPL_PhotParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Parallax: Get lens astrometry
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
PSPL Phot parallax with GP
- class model.PSPL_Phot_Par_GP_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_Phot
,PSPL_Parallax
,PSPL_GP_PhotParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_Phot_Par_GP_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_Phot
,PSPL_Parallax
,PSPL_GP_PhotParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_Phot_Par_GP_Param1_2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_Phot
,PSPL_Parallax
,PSPL_GP_PhotParam1_2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_Phot_Par_GP_Param2_2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_Phot
,PSPL_Parallax
,PSPL_GP_PhotParam2_2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
PSPL Phot, no parallax with GP
- class model.PSPL_Phot_noPar_GP_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_Phot
,PSPL_noParallax
,PSPL_GP_PhotParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
calc_piE_ecliptic
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_Phot_noPar_GP_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_Phot
,PSPL_noParallax
,PSPL_GP_PhotParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
calc_piE_ecliptic
get_astrometry
get_astrometry_unlensed
get_centroid_shift
get_chi2_photometry
get_lens_astrometry
get_lnL_constant
get_photometry
get_resolved_amplification
get_resolved_astrometry
log_likely_astrometry
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
PSPL PhotAstrom, parallax with GP
- class model.PSPL_PhotAstrom_Par_GP_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_GP_PhotAstromParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_GP_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_GP_PhotAstromParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_GP_Param3(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_GP_PhotAstromParam3
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_GP_Param4(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax
,PSPL_GP_PhotAstromParam4
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_LumLens_GP_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_GP_PhotAstromParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_LumLens_GP_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_GP_PhotAstromParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_LumLens_GP_Param3(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_GP_PhotAstromParam3
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_Par_LumLens_GP_Param4(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_Parallax_LumLens
,PSPL_GP_PhotAstromParam4
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
Parallax: Get piE_ecliptic
Parallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])Parallax: Get astrometry
get_astrometry_unlensed
(t_obs)Get the astrometry of the source if the lens didn't exist.
get_centroid_shift
(t[, ast_filt_idx])Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_geoproj_ast_params
(t0par)get_geoproj_params
(t0par)get_lens_astrometry
(t_obs)Parallax: Get lens astrometry
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
start
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- calc_piE_ecliptic()
Parallax: Get piE_ecliptic
- get_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
Parallax: Get astrometry
- get_astrometry_unlensed(t_obs)
Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t, ast_filt_idx=0)
Parallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Parallax: Get lens astrometry
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Parallax: Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Parallax: Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image.
xS_minus is the vector position of the plus image.
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
PSPL PhotAstrom, no parallax with GP
- class model.PSPL_PhotAstrom_noPar_GP_Param1(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_noParallax
,PSPL_GP_PhotAstromParam1
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])noParallax: Position of the observed source position in arcsec.
get_astrometry_unlensed
(t_obs)noParallax: Get the astrometry of the source if the lens didn't exist.
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_lens_astrometry
(t_obs)Equation of motion for just the foreground lens.
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
calc_piE_ecliptic
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t_obs)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.
- class model.PSPL_PhotAstrom_noPar_GP_Param2(*args, **kwargs)
Bases:
ModelClassABC
,PSPL_GP
,PSPL_PhotAstrom
,PSPL_noParallax
,PSPL_GP_PhotAstromParam2
Helper class that provides a standard way to create an ABC using inheritance.
Methods
animate
(tE, time_steps, frame_time, name, ...)Produces animation of microlensing event.
noParallax: Get the photometric amplification term at a set of times, t.
get_astrometry
(t_obs[, ast_filt_idx])noParallax: Position of the observed source position in arcsec.
get_astrometry_unlensed
(t_obs)noParallax: Get the astrometry of the source if the lens didn't exist.
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
get_lens_astrometry
(t_obs)Equation of motion for just the foreground lens.
get_log_det_covariance
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
get_photometry_with_gp
(t_obs, mag_obs, ...)Returns photometry with GP noise added in.
Get the photometric amplification term at a set of times, t for both the plus and minus images.
get_resolved_astrometry
(t_obs)Get the x, y astrometry for each of the two source images, which we label plus and minus.
log_likely_photometry
(t_obs, mag_obs, ...[, ...])Calculate the log-likelihood for the PSPL + GP model and photometric data.
calc_piE_ecliptic
get_chi2_astrometry
get_chi2_photometry
get_lnL_constant
get_photometry
log_likely_astrometry
log_likely_astrometry_each
log_likely_photometry_each
- animate(tE, time_steps, frame_time, name, size, zoom, astrometry)
Produces animation of microlensing event. This function takes the PSPL and makes an animation, the input variables are as follows
- Parameters
- tE:
- number of einstein crossings times before/after the peak you want the animation to plot
e.g tE = 2 => graph will go from -2 tE to 2 tE
- time_steps:
number of time steps before/after peak, so total number of time steps will be 2 times this value
- frame_time:
times in ms of each frame in the animation
- name: string
the animation will be saved as name.html
- size: list
[horizontal, vertical] cm’s
- zoom:
# of einstein radii plotted in vertical direction
- get_amplification(t)
noParallax: Get the photometric amplification term at a set of times, t.
- Parameters
- t:
Array of times in MJD.DDD
- get_astrometry(t_obs, ast_filt_idx=0)
noParallax: Position of the observed source position in arcsec.
- get_astrometry_unlensed(t_obs)
noParallax: Get the astrometry of the source if the lens didn’t exist.
- Returns
- xS_unlensednumpy array, dtype=float, shape = len(t_obs) x 2
The unlensed positions of the source in arcseconds.
- get_centroid_shift(t)
noParallax: Get the centroid shift (in mas) for a list of observation times (in MJD).
- get_lens_astrometry(t_obs)
Equation of motion for just the foreground lens.
- Parameters
- t_obsarray_like
Time (in MJD).
- get_log_det_covariance(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_photometry_with_gp(t_obs, mag_obs, mag_err_obs, filt_index=0, t_pred=None)
Returns photometry with GP noise added in.
Note
This will throw an error if this is a filter with use_gp_phot[filt_index] = False.
- get_resolved_amplification(t)
Get the photometric amplification term at a set of times, t for both the plus and minus images.
- Parameters
- t:
Array of times in MJD.DDD
- get_resolved_astrometry(t_obs)
Get the x, y astrometry for each of the two source images, which we label plus and minus.
- Returns
- [xS_plus, xS_minus]list of numpy arrays
xS_plus is the vector position of the plus image in arcsec
xS_minus is the vector position of the plus image in arcsec
- log_likely_photometry(t_obs, mag_obs, mag_err_obs, filt_index=0)
Calculate the log-likelihood for the PSPL + GP model and photometric data.
Note
The GP will only be used for filters where use_gp_phot[filt_index] = True.