[WIP]: Allowing for general theory covariance matrices to be included in fits

Created by: RosalynLP

This addresses #628 (closed).

The idea is to allow for generic theory covariance matrices to be loaded from file and added together for use in the fit. Currently I am considering the general case of a higher twist covariance matrix but in general we could have several functions in theorycovariance/construction.py which loads different covariance matrices, and then they can all be added together in config.py.

To do:

• Load in theory covariance matrix from file in theorycovariance/construction.py

• Apply cuts to the relevant datasets in this theory covariance matrix

• Find the cuts of the other datasets in the fit and expand the theory covariance matrix with zeros in all the other dataset entries

• Import the expanded theory covariance matrix into config.py intheory_covmat

• Combine the various theory covariance matrices together

• Deal with the fact that you may have some subset of theory covariance matrices but not all and that vp-setupfit should work in all cases (in particular that you can load other covmats when there is no MHOU covmat, but that you need not specify all the different theories - this leads to a runcard mess)

• Take care of the fact that different namespaces can lead to different output file names for the theory covmat which in turn affects them being loaded in in nnfit.cc

N.B. The JLAB data exists in here because it is necessary for the particular case of higher twist errors.