The code Figure10.py needs the following input files:
driving_ang_pol_u_omega=0.2000_pol1.0000_0.0000_Ltarget=1_d0=0.1000_D0=1.6000_B0=12.0000_Vquad0=10.0000.json
modeprojections_OBC_d0=0.2000_hermTrue_simpFalse_D0=1.6000_Npart=1_B0=12.0000_Vquad0=10.0000_rtrunc=-1.0000.txt

The former contains various information about the steady state of a driven system vs. detuning (e.g. one- and two-photon intensities of the output light, overlap with a static ground state) and can be created by running:
python wfnsteady_manybody_angular_2part_emission_pol.py

The latter contains the f coefficients for intensity rescaling and can be generated by running:
python hex_mb_angular_outputlight_mode2.py

Then, Figure 12 is created by running:
python Figure12.py
