A robust, flexible and automatic framework for joint estimation of magnitude of completeness, b-value, magnitude binning and their time variation