CC_Scan_ALL_Stack.m

%CCCP
%Scan All Stacked Cross Correlations
%Scans stacks of cross correlations from all components, all stations
%Read in written/generated template data from file
%all input parameters are modifiable in template_data_input.m
%all general parameters are modifiable in general_settings.m


%Load data about templates
template_data_input
%Load control Settings
general_settings





%Adjust time data to MATLAB readable format
start_date = sprintf('%04d-%02d-%02d %02d:%02d.%d',start_year,start_month,start_day,start_minute,start_second);
end_date = sprintf('%04d-%02d-%02d %02d:%02d.%d',end_year,end_month,end_day,end_minute,end_second);

%Check to ensure the proper flow of time
start_date = datenum(start_date);
end_date = datenum(end_date);
delta_time = end_date - start_date;
is_real = delta_time > 0;
if is_real == 0;
    fprintf('The Campaign for Real Time requests that you please ensure that your ending date follows your starting date. \n');
end

%Check for directory structure
directory_check = sprintf('./%s/%s',base_folder,all_stack_folder);
if exist(directory_check,'dir') ~= 7;
    mkdir(directory_check);
    fprintf('Cross Correlation directory created \n');
end
directory_check = sprintf('./%s/%s',base_folder,waveform_folder);
if exist(directory_check,'dir') ~= 7;
    mkdir(directory_check);
    fprintf('Waveform directory created \n');
end

directory_check = sprintf('./%s/%s',base_folder,result_folder);
if exist(directory_check,'dir') ~= 7;
    mkdir(directory_check);
    fprintf('Result directory created \n');
end



%Failure Control
failure = 0;




%Loop over data for requested time and stack cross correlation functions
All_event_listing_savename = sprintf('./%s/%s_to_%s_events_All_stacks',base_folder,datestr(start_date),datestr(end_date));
events = {'Time Index', 'CC Value','Time (UTC)','MATLAB Time','Threshold Value','Template'};
row = 2;
total_detections = 0;


for time = start_date:CC_increment:(end_date-CC_increment)
    
    start_time = time;
    end_time = time + CC_increment;
    if end_time >= end_date
        end_time = end_date;
    end
    
    for template_count = 1:length(template_list(:,1))
        single_template = template_list(template_count,:);
        template = single_template(template_count).template;
        ALL_CC_Stacked_savename = sprintf('%s/%s/ALL_CC_Stacked_%s_%s.mat',base_folder,all_stack_folder,template,datestr(start_time,30));
        
        repeat = 1;
        while repeat == 1;
            if exist(ALL_CC_Stacked_savename,'file') == 0
                while exist(ALL_CC_Stacked_savename,'file') == 0
                    fprintf('%s %s\n',datestr(time),template);
                    fprintf('Previous ALL Stacked Cross Correlation NOT found\n');
                    fprintf('Running Stacking Program\n');
                    CC_All_Stack
                end
            elseif exist(ALL_CC_Stacked_savename,'file') == 1
                load(ALL_CC_Stacked_savename);
                if isempty(ALL_Stacked_CC) == 1
                    while isempty(ALL_Stacked_CC) == 1
                        fprintf('%s %s\n',datestr(time),template);
                        fprintf('Previous ALL Stacked Cross Correlation NOT found\n');
                        fprintf('Running Stacking Program\n');
                        CC_All_Stack
                    end
                end
            else
                fprintf('%s %s\n',datestr(time),template);
                fprintf('Previous ALL Stacked Cross Correlation found\n');
                
                
                data = get(ALL_Stacked_CC,'data'); %Review CC for prospects
                threshold = 9*mad(data);
                freq = get(ALL_Stacked_CC,'freq');
                [PeakCorr,PeakIndex] = getpeaks(data,'NPEAKS',candidates);
                TimeIndex = PeakIndex/freq;
                DisplayTime = time;
                detections = 0;
                for z = 1:candidates
                    if PeakCorr(z) >= threshold
                        DisplayTime = time + (TimeIndex(z)/(60*60*24));
                        events{row,1} = TimeIndex(z);
                        events{row,2} = PeakCorr(z);
                        events{row,3} = datestr(DisplayTime,'dd mmmm yyyy HH:MM:SS.FFF');
                        events{row,4} = DisplayTime;
                        events{row,5} = threshold;
                        events{row,6} = template;
                        row = row + 1;
                        detections = detections + 1;
                        total_detections = total_detections + 1;
                    end
                end
                fprintf('%d events detected\n',detections);
                fprintf('%d total events detected\n',total_detections);
                repeat = 0;
                
                
            end
        end
    end
end

save(All_event_listing_savename,'events');
%Use dlmcell to spit out a text file
text_savename = sprintf('%s/%s/%s_%s_to_%s_multiple_station.csv',base_folder,result_folder,template,start_date,end_date);
dlmcell(text_savename,events,',');
fprintf('%s saved\n',All_event_listing_savename);