Add y, Q2, x, t resolution tests

Relative difference fitted by Gaussian Add sim~rec 2d correlation plot for those variables

benchmarks/dvmp/analysis/vm_invar.cxx

@@ -135,12 +135,12 @@ int vm_invar(const std::string& config_name)
   double hist_range_h[4] = {1., 15., 0.1, 0.};
   
   TString VarName[4] = {"y", "Q2", "x", "t"};
-  TH1D* h_sim[4];
+  TH1D h_sim[4];
   {
   auto h_tmp = d_im.Histo1D({"h_Q2_sim", ";Q^{2};#", 50, hist_range_l[0], hist_range_h[0]}, "Q2_sim");
   //auto& htmp = *h_tmp;
   
-  h_sim[0] = *h_tmp;
+  &h_sim[0] = *h_tmp;
   }
   
   //TH1D h_sim[4];