Bayesian Optimization of Terahertz Quantum Cascade Lasers