In this talk I will show how to derive analytical gravitational waveforms associated to the coalescence of a “hairy” black hole binary in Einstein-scalar-Gauss-Bonnet gravity (ESGB). I will present the post-Newtonian (PN) Lagrangian and gravitational wave fluxes radiated by such systems. The PN framework relies on reducing the black holes to point particles with scalar field-dependent masses. I will show that this procedure amounts to fixing the Wald entropy of each black hole during their slow, adiabatic inspiral. I will then extend the domain of validity of the PN approximation to the strong field regime near merger, by generalizing the effective-one-body (EOB) formalism to ESGB gravity.