function [detx, dety] = read_dhits_f(SST, n, File_name_base)

%this function reads in detector hit positions from TAPSIM results_data
%SST is the run number
%n is the number of files to stitch

% Author: Andrew Breen, MPIE, 09/07/2018

%File_name_base = 'TAPSIMSimID';

%SST = 80; %this is the experiment number

file_names = {};
file_names{1} = [File_name_base num2str(SST) '.00000001'];


%number of files to stitch
%n = 30;

for i = 1:n
    file_names{i+1} = [File_name_base num2str((SST+i*10^-4), '%.08f')];
end


% %compile list of file names
% file_names = {'results_data.00000001';
%     'results_data.00010000';
%     'results_data.00020000';
%     'results_data.00030000';
%     'results_data.00040000';
%     'results_data.00050000';
%     'results_data.00060000';
%     'results_data.00070000';
%     'results_data.00080000';
%     'results_data.00090000';
%     'results_data.00100000';
%     'results_data.00110000';
%     'results_data.00120000';
%     'results_data.00130000';
%     'results_data.00140000';
%     'results_data.00150000';
%     'results_data.00160000';
%     'results_data.00170000';
%     'results_data.00180000';
%     'results_data.00190000';
%     'results_data.00200000';
%     'results_data.00210000';
%     'results_data.00220000';
%     'results_data.00230000';
%     'results_data.00240000';
%     'results_data.00250000';
%     'results_data.00260000';
%     'results_data.00270000';
%     'results_data.00280000';
%     'results_data.00290000';
%     'results_data.00300000';
%     'results_data.00310000';
%     'results_data.00320000';
%     'results_data.00330000';
%     'results_data.00340000';
%     'results_data.00350000';
%     'results_data.00360000';
%     'results_data.00370000';
%     'results_data.00380000';
%     'results_data.00390000';
%     'results_data.00400000';
%     'results_data.00410000';
%     'results_data.00420000';
%     'results_data.00430000';
%     'results_data.00440000';
%     'results_data.00450000'};

detx = [];
dety = [];

for i = 1:size(file_names,2);
[data,result]=readtext(file_names{i},'[,\t]');

if strcmp(data(40), 'ASCII')
    detx= [detx cell2mat(data(41:end,8))'];
    dety= [dety cell2mat(data(41:end,9))'];
else
    detx= [detx cell2mat(data(30:end,8))'];
    dety= [dety cell2mat(data(30:end,9))'];
end
end

detx = detx';
dety = dety';

% detx_clip = detx(sqrt(detx.^2.*dety.^2)<0.1);
% dety_clip = dety((sqrt(detx.^2.*dety.^2)<0.1);

    
