solution output analysis now takes into account foci list

Author: ltetrel

src/openlifu/plan/solution.py

@@ -199,8 +199,8 @@ def analyze(self, transducer: Transducer, options: SolutionOptions = SolutionOpt
             #     distance=options.beamwidth_radius,
             #     options=mask_options)
 
-            pk = np.max(pnp_MPa.data * mainlobe_mask)  #TODO: pnp_MPa supposed to be a list for each focus: pnp_MPa(focus_index)
-            solution_analysis.mainlobe_pnp_MPa = pk
+            pk = np.max(pnp_MPa.data[focus_index] * mainlobe_mask)  #TODO: pnp_MPa supposed to be a list for each focus: pnp_MPa(focus_index)
+            solution_analysis.mainlobe_pnp_MPa += [pk]
 
         #     thresh_m3dB = pk*10**(-3 / 20)
         #     thresh_m6dB = pk*10**(-6 / 20)