SciPy optimise to maximise function for binned data

Question

I have generated some data, using box integration, based on whether or not it lies beneath a straight line. The code to generate this data is given below:

cl