Fix mc output shapes

Author: andreicuceu

bin/run_vega_mc_mpi.py

@@ -45,6 +45,9 @@ def print_func(message):
         raise ValueError('Warning: You called "run_vega_mc_mpi.py" without asking'
                          ' for monte carlo. Add "run_montecarlo = True" to the "[control]" section.')
 
+    # Activate monte carlo mode
+    vega.monte_carlo = True
+
     # Check if we need to run a forecast
     forecast = vega.main_config['control'].getboolean('forecast', False)
     if forecast:

vega/output.py

@@ -422,8 +422,8 @@ def write_monte_carlo(self, cpu_id=None):
         covariances = np.array(self.analysis.mc_covariances)
 
         names = np.array(list(bestfits.keys()))
-        bestfit_table = np.array([bestfits[name][:, 0] for name in names])
-        errors_table = np.array([bestfits[name][:, 1] for name in names])
+        bestfit_table = np.array([bestfits[name][:, 0] for name in names]).T
+        errors_table = np.array([bestfits[name][:, 1] for name in names]).T
         covariances = covariances.reshape(bestfit_table.shape[0]*len(names), len(names))
 
         # Get the data types for the columns