The code Figure7.py needs the following input files:
modeprojections_OBC_d0=0.1000_hermTrue_simpFalse_D0=1.6000_Npart=1_B0=12.0000_Vquad0=3.0000_rtrunc=-1.0000.txt
output_OBC_d0=[d0]_D0=1.6000_Npart=2_B0=12.0000_Vquad0=3.0000_Lgs2.json
where [d0] stands for 0.1000, 0.2000, 0.3000, 0.4000


The file modeprojections_OBC_d0=0.1000_hermTrue_simpFalse_D0=1.6000_Npart=1_B0=12.0000_Vquad0=3.0000_rtrunc=-1.0000.txt, containing the f coefficients for rescaling of intensities, can be generated by running:
python hex_mb_angular_outputlight_mode2.py

The files output_OBC_d0=[d0]_D0=1.6000_Npart=2_B0=12.0000_Vquad0=3.0000_Lgs2.json, containing various information on eigenstates and emitted light, can be generated using:
python hex_mb_angular_transform_orbbasis.py 
(the value of d0 can be changed by editing the code, see comments inside)

Then. Figure 7 is created by running:
python Figure7.py
