Multi variable fitting 

Implementation of multi-variable fit
 - [x] extend Minuit to allow additional variables
 - [x] add helper functions to recast fit-variables into the fit range
 - [x] change fitting script
    - [x] add fit_variables parameter
    - [x] load additional fit_variables
    - [x] store templates for additional variables
 - [x] add scripts to run the fit checks
This gives 1000x lumi, as far as I can see on the plots