BGLS calculates the Bayesian Generalized Lomb-Scargle periodogram. It takes as input arrays with a time series, a dataset and errors on those data, and returns arrays with sampled periods and the periodogram values at those periods.
https://ui.adsabs.harvard.edu/abs/2015A%26A...573A.101M