The code Figure5.py needs the following input files:
modeprojections_OBC_d0=0.3000_hermTrue_simpFalse_D0=0.6000_Npart=1_B0=12.0000_Vquad0=0.0000_rtrunc=-1.0000.txt
output_OBC_d0=0.3000_D0=0.6000_Npart=2_B0=12.0000_Vquad0=0.0000_Lgs2.json

The former, containing the f coefficients for intensity rescaling, are created with:
python hex_mb_angular_outputlight_mode2.py  

The latter, containing various information about the eigenspectrum and the output light (some of which are used when plotting the figure), can be created by running:
python hex_mb_angular_transform_orbbasis.py

Then, the figure is created by running: 
python Figure5.py
