function []=figure3b()
close all
figure()
hold on
%% Data parameters
sample_index=62;
CP_Indices=[12,13,14,15,16,17,18];
%%
NCP=size(CP_Indices,2);
for iCP_dummy=1:NCP
    CP_Index=CP_Indices(iCP_dummy);
    %% Loading data
    data_name=sprintf('q_raw_binned_%d_%d.dat',sample_index,CP_Index);
    fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
    Ndata=round(size(numbers,1)/2);
    q_bin_mids_bare=numbers(2*(1:Ndata)-1);
    q_binned=numbers(2*(1:Ndata));
    
    %% Spline fitting
    qs=(0:0.001:1).';
    q_binned = spline(q_bin_mids_bare,q_binned,qs);
    %% Plotting
    plot(qs,q_binned)
    xlim([0,1])
    ylim([0,3.3])
    xlabel('$Q$','Interpreter','latex')
    ylabel('$P(Q)$','Interpreter','latex')
end
print(gcf, '-dpdf', 'figure3b.pdf');
end