The code Figure15.py needs following input files:
spectrum_OBC_d0=[d0]_hermTrue_simpFalse_D0=0.6000_Npart=[Npart]_B0=12.0000_Vquad0=0.0000_rtrunc=-1.0000_mmax5.txt
where [d0] stands for 0.1000, 0.2000, 0.3000, 0.4000, 0.5000, 1.0000, 1.7000 and [Npart] stands for 1,2,3
These files contain energy spectra with decay rates and overlaps with model states

They can be generated by running:
python hex_mb_angular_transform_overlap_manystates3_1d.py

By editing the code one can change the d0 and Npart, as well as turn on or off the calculation of decay rate and overlap with model states

Note Figure15.py code assumes that overlap is NOT calculated in Npart=1 case (when overlap calculation makes no sense), and that the plotted quantity (decay rate for Npart=1, overlap for Npart=2,3) is always in the second column of the input file. 


With the input files present, the figure is generated by running
python Figure15.py
