model.loaded_spectra += [input_spectrum] plt.plot(model.loaded_spectra[0].x, model.loaded_spectra[0].lineshape) filename2 = os.path.join(directory, r'data\scatterers2.json') model.loadScatterers(filename2) model.unscattered_spectrum = model.loaded_spectra[-1] #%% model.scattering_medium.d_through_gas = 10000 model.scattering_medium.pressure = 40 model.setCurrentScatterer('Shirley') model.changeAlgorithm(1) model.calculation.n_iter = 50 model.algorithm_option = 'bulk' model.scattering_medium.scatterer.norm_factor = 0.6 shirley_params = {'edge': 0, 'exponent': 0.05, 'fermi_width': 0.01} model.scattering_medium.scatterer.loss_function.editComponent( 0, shirley_params) model.scatterSpectrum() fig, ax = plt.subplots(1, 2) ax[0].plot(model.simulated_spectrum.x, model.simulated_spectrum.lineshape) ax[1].plot(model.scattering_medium.scatterer.loss_function.x, model.scattering_medium.scatterer.loss_function.lineshape) ax[1].set_xlim(0, 200)