chis_max = (mass1*spin1_bounds[1] + mass2*spin2_bounds[1])/mtot