This is an IDL package to analysis QSO variability and do photometric reverberation mapping (photo-RM)
The core code are learnt from the JAVELIN code (please see https://bitbucket.org/nye17/javelin) and the Python emcee package (http://dan.iel.fm/emcee/current/). There is a modified version of JAVELIN code, including some features such as, adding structure function and changing the transfer functions by Haowen (please see https://github.com/pkuzhw387/Modified_JAVELIN)
Qian Yang
July 30, 2017
Damped random walk (DRW) model QSO continuum variability (probability calculated in lnlike_from_CL.pro)
Continuum + emission line with lag model (covariance matrix in get_S.pro)
Maximum the probability, using zoom in boxes (run drw_zoomin.pro)
MCMC using Metropolis-Hastings method (run drw_mcmc.pro)
MCMC using method learn from emcee (run drw_emcee.pro)
MCMC (emcee.pro debug)
Structure function (compare with previous works)
Continuum lag (no transfer function)
DRW in mag, while scale in flux
Consider host galaxy contribution (not change), which is different from emission line contribution (has lag and transfer function)
Predict light curves (simulate LSST light curve)
Try different format of transfer functions
DRW using different bands
Two emission lines lag
Multi-band analysis
Define conditions (S/N, time baseline, time interval, line ratio)
Any question please feel free to email Qian ([email protected]).