https://github.com/BlackrockMicrosystems/Python-Utilities/blob/2b93f0d10ff8fd9dcb627e3eb579ca88e5f39b87/brpylib.py#L501
With the original codes one may lose 2 temporal samples for each spike waveform. It seems that the right codes should be wfs = np.ndarray((nPackets,int((self.basic_header['BytesInDataPackets']-(tsBytes+6))/2)), '<h', rawdata, tsBytes+6, (self.basic_header['BytesInDataPackets'], 2))