%centroid idx is the first
clust_sample_centroid_idx = clust_dist2mu_sorted_idx(1);

centroid_FM_idx = not_labeled_idx(aux_idx(aux_idx2(clust_sample_centroid_idx)));

centroid_QRS_idx = (centroid_FM_idx);

aux_hb_seq =    max(1, QRS_locations(centroid_QRS_idx) - round(2*sampling_rate)): ...
                min(ECG_header.nsamp, QRS_locations(centroid_QRS_idx) + round(2*sampling_rate));

aux_hb_seq2 =   max(1, centroid_FM_idx - 5):min(cant_QRS_locations, centroid_FM_idx + 2);

aux_hb_seq3 =    max(1, QRS_locations(centroid_QRS_idx) - round(0.2*sampling_rate)): ...
                min(ECG_header.nsamp, QRS_locations(centroid_QRS_idx) + round(0.5*sampling_rate));


cCentroid = {   (QRS_locations(centroid_QRS_idx) - aux_hb_seq(1) + 1) ... %relative position of the heartbeat
                ECG_total(aux_hb_seq,:) ... % ecg excerpt
                centroid_FM_idx - aux_hb_seq2(1) + 1 ... % relative position of the heartbeat 
                RR_intervals(aux_hb_seq2)/sampling_rate ... % RR interval sequence
                (QRS_locations(centroid_QRS_idx) - aux_hb_seq3(1) + 1) ... %relative position of the heartbeat                
                ECG_total(aux_hb_seq3,:) ... % morph excerpt
                clust_data_m ... % cluster size
                };
            
%closer samples
quartile_size = round(clust_data_m * 0.25);
if( quartile_size > CantSamplesClose )
    clust_samples_close_idx = clust_dist2mu_sorted_idx( 1 + randsample(quartile_size, CantSamplesClose ) );
else
    if( clust_data_m < 2 )
        % the centroid itself
        clust_samples_close_idx = clust_dist2mu_sorted_idx;
    else
        clust_samples_close_idx = clust_dist2mu_sorted_idx( 2:(2+quartile_size-1) );
    end
end

clust_samples_close_idx = sort(clust_samples_close_idx);
clust_samples_close_FM_idx = not_labeled_idx(aux_idx(aux_idx2(clust_samples_close_idx)));
clust_samples_close_QRS_idx = (clust_samples_close_FM_idx);


aux_ECGslices = [];
aux_RRslices = [];
aux_MorphSlices = [];
aux_rel_posQRS = [];
aux_rel_posRR = [];
aux_rel_MorphposQRS = [];

   
    aux_seq = colvec(1:length(clust_samples_close_QRS_idx));
    
    aux_hb_seq = arrayfun(@(a)( max(1, QRS_locations(clust_samples_close_QRS_idx((a))) - round(2*sampling_rate)): ...
                                min(ECG_header.nsamp, QRS_locations(clust_samples_close_QRS_idx((a))) + round(2*sampling_rate)) ) , ...
                           aux_seq, 'UniformOutput', false);

    aux_hb_seq2 = arrayfun(@(a)(    max(1, clust_samples_close_FM_idx((a)) - 5): ...
                                    min(cant_QRS_locations, clust_samples_close_FM_idx((a)) + 2) ) , ...
                           aux_seq, 'UniformOutput', false);

    aux_hb_seq3 = arrayfun(@(a)( max(1, QRS_locations(clust_samples_close_QRS_idx((a))) - round(0.2*sampling_rate)): ...
                                min(ECG_header.nsamp, QRS_locations(clust_samples_close_QRS_idx((a))) + round(0.5*sampling_rate)) ) , ...
                           aux_seq, 'UniformOutput', false);

    %ECG slices
    aux_ECGslices = [aux_ECGslices; cellfun(@(a)(ECG_total(a,:)), aux_hb_seq, 'UniformOutput', false)];
    %RR slices
    aux_RRslices = [aux_RRslices; cellfun(@(a)(RR_intervals(a)/sampling_rate), aux_hb_seq2, 'UniformOutput', false)];
    %Morph slices
    aux_MorphSlices = [aux_MorphSlices; cellfun(@(a)(ECG_total(a,:)), aux_hb_seq3, 'UniformOutput', false)];
    
    %relative positions
    aux_rel_posQRS = [aux_rel_posQRS; cellfun(@(a,b)(QRS_locations(clust_samples_close_QRS_idx((a))) - b(1) + 1), num2cell(aux_seq), aux_hb_seq)];
    aux_rel_posRR =  [aux_rel_posRR; cellfun(@(a,b)(clust_samples_close_FM_idx((a)) - b(1) + 1), num2cell(aux_seq), aux_hb_seq2)];
    aux_rel_MorphposQRS = [aux_rel_MorphposQRS; cellfun(@(a,b)(QRS_locations(clust_samples_close_QRS_idx((a))) - b(1) + 1), num2cell(aux_seq), aux_hb_seq3)];
                       

if( iscell(aux_ECGslices) )
    cCloserExamples = { aux_rel_posQRS ... %relative position of the heartbeat
                        aux_ECGslices ... % ecg excerpt
                        aux_rel_posRR ... % relative position of the heartbeat 
                        aux_RRslices ... % RR interval sequence
                        aux_rel_MorphposQRS ... % Morph relative position of the heartbeat 
                        aux_MorphSlices ... % Morph interval sequence
                        };
else
    cCloserExamples = { aux_rel_posQRS ... %relative position of the heartbeat
                        {aux_ECGslices} ... % ecg excerpt
                        aux_rel_posRR ... % relative position of the heartbeat 
                        {aux_RRslices} ... % RR interval sequence
                        aux_rel_MorphposQRS ... % Morph relative position of the heartbeat 
                        {aux_MorphSlices} ... % Morph interval sequence
                        };
    
end

%distant samples
if( quartile_size > CantSamplesFar )
    clust_samples_far_idx = clust_dist2mu_sorted_idx( clust_data_m - randsample(quartile_size, CantSamplesFar ) + 1  );
else
    if( clust_data_m < 3 )
        % all examples
        clust_samples_far_idx = clust_dist2mu_sorted_idx;
    else
        clust_samples_far_idx = clust_dist2mu_sorted_idx(3*quartile_size:clust_data_m);
    end
end


clust_samples_far_idx = sort(clust_samples_far_idx);
clust_samples_far_FM_idx = not_labeled_idx(aux_idx(aux_idx2(clust_samples_far_idx)));
clust_samples_far_QRS_idx = (clust_samples_far_FM_idx);

aux_ECGslices = [];
aux_RRslices = [];
aux_MorphSlices = [];
aux_rel_posQRS = [];
aux_rel_posRR = [];
aux_rel_MorphposQRS = [];

    
    aux_seq = colvec(1:length(clust_samples_far_QRS_idx));
    
    aux_hb_seq = arrayfun(@(a)( max(1, QRS_locations(clust_samples_far_QRS_idx((a))) - round(2*sampling_rate)): ...
                                min(ECG_header.nsamp, QRS_locations(clust_samples_far_QRS_idx((a))) + round(2*sampling_rate)) ) , ...
                           aux_seq, 'UniformOutput', false);

    aux_hb_seq2 = arrayfun(@(a)(    max(1, clust_samples_far_FM_idx((a)) - 5): ...
                                    min(cant_QRS_locations, clust_samples_far_FM_idx((a)) + 2) ) , ...
                           aux_seq, 'UniformOutput', false);

    aux_hb_seq3 = arrayfun(@(a)( max(1, QRS_locations(clust_samples_far_QRS_idx((a))) - round(0.2*sampling_rate)): ...
                                min(ECG_header.nsamp, QRS_locations(clust_samples_far_QRS_idx((a))) + round(0.5*sampling_rate)) ) , ...
                           aux_seq, 'UniformOutput', false);
                       
    %ECG slices
    aux_ECGslices = [aux_ECGslices; cellfun(@(a)(ECG_total(a,:)), aux_hb_seq, 'UniformOutput', false)];
    %RR slices
    aux_RRslices = [aux_RRslices; cellfun(@(a)(RR_intervals(a)/sampling_rate), aux_hb_seq2, 'UniformOutput', false)];
    %Morph slices
    aux_MorphSlices = [aux_MorphSlices; cellfun(@(a)(ECG_total(a,:)), aux_hb_seq3, 'UniformOutput', false)];
    
    %relative positions
    aux_rel_posQRS = [aux_rel_posQRS; cellfun(@(a,b)(QRS_locations(clust_samples_far_QRS_idx((a))) - b(1) + 1), num2cell(aux_seq), aux_hb_seq)];
    aux_rel_posRR =  [aux_rel_posRR; cellfun(@(a,b)(clust_samples_far_FM_idx((a)) - b(1) + 1), num2cell(aux_seq), aux_hb_seq2)];
    aux_rel_MorphposQRS = [aux_rel_MorphposQRS; cellfun(@(a,b)(QRS_locations(clust_samples_far_QRS_idx((a))) - b(1) + 1), num2cell(aux_seq), aux_hb_seq3)];


if( iscell(aux_ECGslices) )
    cDistantExamples ={ aux_rel_posQRS ... %relative position of the heartbeat
                        aux_ECGslices ... % ecg excerpt
                        aux_rel_posRR ... % relative position of the heartbeat 
                        aux_RRslices ... % RR interval sequence
                        aux_rel_MorphposQRS ... % Morph relative position of the heartbeat 
                        aux_MorphSlices ... % Morph interval sequence
                        };
else
    %force cell contents.
    cDistantExamples ={ aux_rel_posQRS ... %relative position of the heartbeat
                        {aux_ECGslices} ... % ecg excerpt
                        aux_rel_posRR ... % relative position of the heartbeat 
                        {aux_RRslices} ... % RR interval sequence
                        aux_rel_MorphposQRS ... % Morph relative position of the heartbeat 
                        {aux_MorphSlices} ... % Morph interval sequence
                        };
end

clear aux