Number of classes is less than number of spikes and waveforms about combinato HOT 1 CLOSED

#%%  This section works to load files created by combinato and check the spike times
import h5py
import numpy as np
from scipy.io import savemat
import os
from glob import glob

# Directory containing folder output from Combinato
folder = ''

for sfold in os.listdir(folder):
    if sfold.startswith("u"): # Enter micro folders only
        print(folder+sfold)
        sortf = folder+sfold+'/sort_pos_simple/sort_cat.h5' # Find sorting outcomes
        rawf = glob(folder+sfold+'/data_*')[0] # Raw file for waveforms
        if os.path.exists(sortf): # If files don't have any clusters, then none is created
            f = h5py.File(sortf, "r")
            group = f['groups'][:] # overall waveform groups: first index is overall unit index, second is group index (no within group index)
            # Use group to convert matches to neuron group
            classes = f['classes'][:] # each small cluster
            index = f['index'][:] #  seems to index spike times
            matches = f['matches'][:] # orange blue or cyan color in plots
    
            f2 = h5py.File(rawf, "r")
            spike_times = f2['pos']['times'][index] 
            pos_wave = f2['pos']['spikes'][index]
            # spike_timesneg = f2['neg']['times'][index] 
            spike_group = np.zeros(len(classes),dtype='int16')
    
            for i, ic in enumerate(np.unique(classes)): 
                spike_group[classes==ic] = group[i,1]
                
                
            class_times = np.c_[spike_times, classes]
    
            savedict = {'class_times': class_times, 'spike_group':spike_group, 'matches':matches, 'spikes':pos_wave}
    
            savemat(folder+sfold+r'/mat_spikes_'+sfold+'.mat',savedict)