moved fit routine to their file, added it to the juobs module, fix in fit_routine with combined fits

Improved pvalue functions,now they use less memory.  pvalue function now are in their own file juobs_pvalue.jl

small changes in get_chi and fit_routine, now xdata is typeless in get_chi and AbstractVector in fit_routine

Added methods for read_meson and read_meson_correnction that take a Vector{Tuple{String,String}} as Gamma structure. This method allow to read in one sweep all correlator for any specific gamma structure combination. It is, in a way, complementary to what was already available when one of the two gamma structure was not specified

added keyword parameter 'replica' to corr_obs(cdata::CData,...) for consinstency with corr_obs(cdata::Array{CData,1}...)