s2min = max(spin2_bounds[0], (mtot*chis - mass1*spin1_bounds[1])/mass2)