function [] = onsetParameterSweepPaper()
currentFile = mfilename( 'fullpath' );
[pathstr,~,~] = fileparts( currentFile );
filename = [pathstr,filesep,'..',filesep,'Data',filesep,'parameterSweep3d.dat'];
data = readtable(filename);
dimInput = data{:,1};
nInput = data{:,2};
phiEqInput = data{:,3};
seedInput = data{:,4};
zeroFInput = data{:,5};
iterationFactorInput = data{:,6};
alphaInput = data{:,7};
phiJInput = data{:,8};

dim = 3;
nP = 1024;
phiEq = 0.4;
numTry = 32;
alphaP = [0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.95 0.99]';
betaP = [3 6 9 12]';
phiM = NaN(size(alphaP)); phiS = NaN(size(alphaP)); success = NaN(size(alphaP));
for i = 1:numel(alphaP)
    thisAlpha = alphaP(i);
    phiJ = phiJInput(alphaInput==thisAlpha & iterationFactorInput==6);
    phiM(i) = mean(phiJ);
    phiS(i) = std(phiJ)/sqrt(numel(phiJ));
    success(i) = numel(phiJ);
end

success = success/numTry;
    
phiM = NaN(size(betaP)); phiS = NaN(size(betaP)); success = NaN(size(betaP));
for i = 1:numel(betaP)
    thisBeta = betaP(i);
    phiJ = phiJInput(alphaInput==0.9 & iterationFactorInput==thisBeta);
    phiM(i) = mean(phiJ);
    phiS(i) = std(phiJ)/sqrt(numel(phiJ));
    success(i) = numel(phiJ);
end

success = success/numTry;



figure(4)
clf; hold all
errorbar(betaP/3,success,sqrt(success.*(1-success)/numTry),'.r','MarkerSize',10,'Linewidth',1.5)
xlabel('\beta')
ylabel('Fraction Jammed')
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
axis([0 5 0.5 1])

% Gardner Plots
filename = [pathstr,filesep,'..',filesep,'Data',filesep,'parameterSweepGardner.dat'];
data = readtable(filename);
dimInput = data{:,1};
nInput = data{:,2};
phiEqInput = data{:,3};
seedInput = data{:,4};
zeroFInput = data{:,5};
iterationFactorInput = data{:,6};
alphaInput = data{:,7};
phiJInput = data{:,9};

phiM = NaN(size(alphaP)); phiS = NaN(size(alphaP));
for i = 1:numel(alphaP)
    thisAlpha = alphaP(i);
    phiJ = phiJInput(alphaInput==thisAlpha & iterationFactorInput==6);
    phiM(i) = mean(phiJ);
    phiS(i) = std(phiJ)/sqrt(numel(phiJ));
end

figure(5)
clf; hold all
errorbar(alphaP,phiM,phiS,'.r','MarkerSize',10,'Linewidth',1.5)
xlabel('\alpha')
ylabel('\phi_G')
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
axis([0 1 0.625 0.64])

phiM = NaN(size(betaP)); phiS = NaN(size(betaP));
for i = 1:numel(betaP)
    thisBeta = betaP(i);
    phiJ = phiJInput(alphaInput==0.9 & iterationFactorInput==thisBeta);
    phiM(i) = mean(phiJ);
    phiS(i) = std(phiJ)/sqrt(numel(phiJ));
end

figure(6)
clf; hold all
errorbar(betaP/3,phiM,phiS,'.r','MarkerSize',10,'Linewidth',1.5)
xlabel('\beta')
ylabel('\phi_G')
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
axis([0 5 0.632 0.642])

% Now, do phiJError from filestructure
basedir = [pathstr,filesep,'..',filesep,'Data',filesep,'on3params',filesep,'3',filesep,'1024',filesep,'0.4'];

phiJE = NaN(numel(alphaP),1);
phiJEE = NaN(numel(alphaP),1);
ngood = zeros(numel(alphaP),1);
phiM = NaN(numel(alphaP),1); phiS = NaN(numel(alphaP),1);
for i = 1:numel(alphaP)
    thisAlpha = alphaP(i);
    thisPhiJE = NaN(32,1);
    thisPhiJ = NaN(32,1);
    for j = 0:31
        try
            phiList = load([basedir,filesep,num2str(j),filesep,'2',filesep,'6',filesep,num2str(thisAlpha),filesep,'phiList.dat']);
            thisPhiJ(j) = phiList(end);
            thisPhiJE(j+1) = load([basedir,filesep,num2str(j),filesep,'2',filesep,'6',filesep,num2str(thisAlpha),filesep,'phiJError.dat']);
            ngood(i) = ngood(i) + 1;
        catch
        end
    end
    phiJE(i) = geomean(thisPhiJE(~isnan(thisPhiJE)));
    phiJEE(i) = geostd(thisPhiJE(~isnan(thisPhiJE)));
    phiM(i) = nanmean(thisPhiJ);
    phiS(i) = nanstd(thisPhiJ)/sqrt(32);
end

alphaP = alphaP(ngood>0);
phiJE = phiJE(ngood>0);
phiJEE = phiJEE(ngood>0)./sqrt(ngood(ngood>0));
ngood = ngood(ngood>0);

figure(1)
clf; hold all
errorbar(alphaP,phiM,phiS,'.r','MarkerSize',10,'Linewidth',1.5)
xlabel('\alpha')
ylabel('\phi_J')
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
axis([0 1 0.638 0.645])

figure(7)
clf; hold all
errorbar(alphaP,phiJE,(phiJE - phiJE./phiJEE),(phiJE.*phiJEE - phiJE),'.r','MarkerSize',10,'Linewidth',1.5)
%errorbar(alphaP,phiJE,(phiJE - phiJE./phiJEE)./sqrt(ngood),(phiJE.*phiJEE - phiJE)./sqrt(ngood),'.r','MarkerSize',10,'Linewidth',1.5)
xlabel('\alpha')
ylabel('\Delta\phi_{RS}/\phi_J')
axis([0 1 0.000001 0.02])
set(gca,'yscale','log')
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
set(gca,'YMinorTick','Off')
yticks([1e-6 1e-5 1e-4 1e-3 1e-2])
% axis([0 1 0.625 0.64])

phiJE = NaN(numel(betaP),1);
phiJEE = NaN(numel(betaP),1);
ngood = zeros(numel(betaP),1);
phiM = NaN(numel(betaP),1); phiS = NaN(numel(betaP),1);
for i = 1:numel(betaP)
    thisBeta = betaP(i);
    thisPhiJE = NaN(32,1);
    for j = 0:31
        try
            phiList = load([basedir,filesep,num2str(j),filesep,'2',filesep,num2str(thisBeta),filesep,'0.9',filesep,'phiList.dat']);
            thisPhiJ(j) = phiList(end);
            thisPhiJE(j+1) = load([basedir,filesep,num2str(j),filesep,'2',filesep,num2str(thisBeta),filesep,'0.9',filesep,'phiJError.dat']);
            ngood(i) = ngood(i) + 1;
        catch
        end
    end
    phiJE(i) = geomean(thisPhiJE(~isnan(thisPhiJE)));
    phiJEE(i) = geostd(thisPhiJE(~isnan(thisPhiJE)));
    phiM(i) = nanmean(thisPhiJ);
    phiS(i) = nanstd(thisPhiJ)/sqrt(32);
end

betaP = betaP(ngood>0);
phiJE = phiJE(ngood>0);
phiJEE = phiJEE(ngood>0)./sqrt(ngood(ngood>0));
ngood = ngood(ngood>0);

figure(6)
% clf; hold all
errorbar(betaP/3,phiM,phiS,'.b','MarkerSize',10,'Linewidth',1.5)
xlabel('\tau')
ylabel('\phi_J')
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
axis([0 5 0.638 0.645])
% yticks([0.641 0.642 0.643 0.644])

figure(8)
clf; hold all
errorbar(betaP/3,phiJE,(phiJE - phiJE./phiJEE),(phiJE.*phiJEE - phiJE),'.r','MarkerSize',10,'Linewidth',1.5)
xlabel('\tau')
ylabel('\Delta\phi_{RS}/\phi_J')
axis([0 5 0.000001 0.02])
plotSpecs()
set(gca,'Linewidth',2)
set(gca,'FontSize',20)
set(gca,'yscale','log');
set(gca,'YMinorTick','Off')
yticks([1e-6 1e-5 1e-4 1e-3 1e-2])
% axis([0 1 0.625 0.64])


end