The code Figure3.py needs the following input files:

'spectrum_OBC_d0=0.0500_hermTrue_simpFalse_D0=2.1000_Npart=3_B0=12.0000_Vquad0=3.0000_rtrunc=-1.0000_mmax9.txt'
'spectrum_OBC_d0=0.3000_hermTrue_simpFalse_D0=1.6000_Npart=2_B0=12.0000_Vquad0=10.0000_rtrunc=-1.0000_mmax7.txt'
'spectrum_OBC_d0=0.3000_hermTrue_simpFalse_D0=0.6000_Npart=2_B0=12.0000_Vquad0=0.0000_rtrunc=-1.0000_mmax5.txt'

The first two files (where the overlap is computed with respect to 2D Laughlin states) can be created by running:
python hex_mb_angular_transform_overlap_manystates2.py   
The code diagonalizes the Hamiltonian and computes the overlap for a given set of parameters. To change the parameters, edit the parameter values inside the code (see also the comments inside)

Similarly, the third file (where the overlap is computed with respect to 1D Laughlin states) can be created by running:
python hex_mb_angular_transform_overlap_manystates3_1d.py

If all the input files are present, then run
python Figure3.py 
to create the figure
