function [Ns] = ceph(path)

conv = 0.13;

% go into directory ceph and compile alldata
cd(path)
cd images
cd ceph
alldata = compile_data_all('cell',0,0.08,0);
alldata = screen_poles(alldata);

% record data
ct = alldata.midthickness*conv;
cm = mean(alldata.midthickness*conv);
cs = std(alldata.midthickness*conv);
cn = alldata.ncells;

% record 
cd ..
save ceph alldata
    
% do same for noceph dir
cd noceph
alldata = compile_data_all('cell',0,0.08,0);
alldata = screen_poles(alldata);

% record data
nt = alldata.midthickness*conv;
nm = mean(alldata.midthickness*conv);
ns = std(alldata.midthickness*conv);
nn = alldata.ncells;

cd ..
save noceph alldata

% print out data
fprintf('Number of cells in ceph dir: %d\n',cn);
fprintf('Avg mid thickness: %g +- %g\n',cm,cs);

% print out data
fprintf('Number of cells in non-ceph dir: %d\n',nn);
fprintf('Avg mid thickness: %g +- %g\n',nm,ns);

% make histogram normalizing values of the two data sets
if nn>cn
  N = round(nn/cn);
  Nn = nt;
  Nc = repmat(ct,1,N);
  fprintf('The number of ceph cells has been multiplied by %d\n',N);
else
  N = round(cn/nn);
  Nc = ct;
  Nn = repmat(nt,1,N);  
  fprintf('The number of non-ceph cells has been multiplied by %d\n',N);  
end

Ns(1,1:numel(Nn)) = Nn;
Ns(2,1:numel(Nc)) = Nc;

Ns(1,numel(Nn)+1:end) = NaN;
Ns(2,numel(Nc)+1:end) = NaN;

cd ..
close all
figure;
hist([Ns(1,:)',Ns(2,:)'],[0.6:0.01:1.1]);
xlabel('Cell width (microns)');
ylabel('Number of cells');
print('-dpdf','hist-sep.pdf');

cd ..