From 4b6e1a176bfecc3e2b7101250d55cb7913234c45 Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Thu, 2 Oct 2025 17:34:09 -0700 Subject: [PATCH 1/4] Extract spectral laser intensity --- .../staging_injector/templates/analyze_simulation.py | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/simulation_data/staging_injector/templates/analyze_simulation.py b/simulation_data/staging_injector/templates/analyze_simulation.py index 77430712..407cba36 100755 --- a/simulation_data/staging_injector/templates/analyze_simulation.py +++ b/simulation_data/staging_injector/templates/analyze_simulation.py @@ -50,7 +50,6 @@ def analyze_simulation(): nitrogen_density_function = re.findall(r'nitrogen1\.density_function\(x,y,z\) = (.+)', text)[0] plateau_length = eval( re.findall(r'my_constants\.plateau_length = (.+)', text)[0] ) - # Compute red/blue shift: wavelength such that 13.5%/86.5% of the spectrum energy is below S, info = ts.get_laser_spectral_intensity( iteration=ts.iterations[-1], pol=pol) @@ -190,11 +189,11 @@ def visualize_iteration(iteration): # Plot of the laser spectrum fig.add_subplot(gs[3,1]) - S, info = ts.get_spectrum(iteration=iteration, pol=pol) + S, info = ts.get_laser_spectral_intensity(iteration=iteration, pol=pol) lambd = 2*np.pi*c/info.omega[1:] plt.xlabel(r'Wavelength [$\mu m$]') - plt.title('Laser spectrum') - plt.plot( 1.e6*lambd, 1e3*S[1:], color='r' ) + plt.title('Laser spectrum [J/m]') + plt.plot( 1.e6*lambd, S[1:]/lambd**2, color='r' ) plt.xlim(0.5,1.2) plt.subplots_adjust(hspace=0.5) From 2854d22b91d72a8acc930a401aaf450bf2a1209c Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Thu, 2 Oct 2025 17:42:50 -0700 Subject: [PATCH 2/4] Update units --- .../staging_injector/templates/analyze_simulation.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/simulation_data/staging_injector/templates/analyze_simulation.py b/simulation_data/staging_injector/templates/analyze_simulation.py index 407cba36..5b86a20a 100755 --- a/simulation_data/staging_injector/templates/analyze_simulation.py +++ b/simulation_data/staging_injector/templates/analyze_simulation.py @@ -192,8 +192,8 @@ def visualize_iteration(iteration): S, info = ts.get_laser_spectral_intensity(iteration=iteration, pol=pol) lambd = 2*np.pi*c/info.omega[1:] plt.xlabel(r'Wavelength [$\mu m$]') - plt.title('Laser spectrum [J/m]') - plt.plot( 1.e6*lambd, S[1:]/lambd**2, color='r' ) + plt.title('Laser spectrum [J/\mu m]') + plt.plot( 1.e6*lambd, 1e-6*S[1:]/lambd**2, color='r' ) plt.xlim(0.5,1.2) plt.subplots_adjust(hspace=0.5) From 8872d6381b46baa24a67523a9813e2adde3d6d89 Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Thu, 2 Oct 2025 17:46:03 -0700 Subject: [PATCH 3/4] Use k in new function --- .../staging_injector/templates/analyze_simulation.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/simulation_data/staging_injector/templates/analyze_simulation.py b/simulation_data/staging_injector/templates/analyze_simulation.py index 5b86a20a..3fb7a91e 100755 --- a/simulation_data/staging_injector/templates/analyze_simulation.py +++ b/simulation_data/staging_injector/templates/analyze_simulation.py @@ -190,7 +190,7 @@ def visualize_iteration(iteration): # Plot of the laser spectrum fig.add_subplot(gs[3,1]) S, info = ts.get_laser_spectral_intensity(iteration=iteration, pol=pol) - lambd = 2*np.pi*c/info.omega[1:] + lambd = 2*np.pi/info.k[1:] plt.xlabel(r'Wavelength [$\mu m$]') plt.title('Laser spectrum [J/\mu m]') plt.plot( 1.e6*lambd, 1e-6*S[1:]/lambd**2, color='r' ) From 7b36398287b123247b99c807be4b201833976d65 Mon Sep 17 00:00:00 2001 From: Remi Lehe Date: Thu, 2 Oct 2025 17:49:54 -0700 Subject: [PATCH 4/4] Fix unit formatting --- .../staging_injector/templates/analyze_simulation.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/simulation_data/staging_injector/templates/analyze_simulation.py b/simulation_data/staging_injector/templates/analyze_simulation.py index 3fb7a91e..b5e3fe40 100755 --- a/simulation_data/staging_injector/templates/analyze_simulation.py +++ b/simulation_data/staging_injector/templates/analyze_simulation.py @@ -192,7 +192,7 @@ def visualize_iteration(iteration): S, info = ts.get_laser_spectral_intensity(iteration=iteration, pol=pol) lambd = 2*np.pi/info.k[1:] plt.xlabel(r'Wavelength [$\mu m$]') - plt.title('Laser spectrum [J/\mu m]') + plt.title('Laser spectrum [$J/\mu m$]') plt.plot( 1.e6*lambd, 1e-6*S[1:]/lambd**2, color='r' ) plt.xlim(0.5,1.2)