The code Figure6.py needs the following input files:
output_OBC_d0=[d0]_D0=0.6000_Npart=2_B0=12.0000_Vquad0=0.0000_Lgs2.json
spectrum_OBC_d0=[d0]_hermTrue_simpFalse_D0=0.6000_Npart=2_B0=12.0000_Vquad0=0.0000_rtrunc=-1.0000_mmax5.txt
where [d0] stands for d0=0.1000, 0.2000, 0.3000, 0.4000, 0.5000

The files output_OBC_d0=[d0]_D0=0.6000_Npart=2_B0=12.0000_Vquad0=0.0000_Lgs2.json, containing various information on eigenstates and output light, are created by running:
python hex_mb_angular_transform_orbbasis.py

while the files spectrum_OBC_d0=[d0]_hermTrue_simpFalse_D0=0.6000_Npart=2_B0=12.0000_Vquad0=0.0000_rtrunc=-1.0000_mmax5.txt, containing the eigenspectra and overlap with model states, are created by running:
python hex_mb_angular_transform_overlap_manystates3_1d.py

The value of d0, as well as other parameters, can be changed in the code itself (see comments inside).

When all the input files are present, the figure can be create by running:
python Figure6.py
