def measure_raw(state_vect, shots=1_000_000, random_seed=1):
    """
    Make measurements along the z basis, return raw data

    :param state_vect: qiskit.quantum_info.Statevector object
    :param shots: int, representing # of measurements
    :param random_seed: int, for setting the 'randomness'
    :return: a list of np.arrays containing the measurement results shape = (shots, num_qbits)
    """

    state_vect.seed(random_seed)
    data_lst = state_vect.sample_memory(shots)
    results = []
    for measurement in data_lst:
        results.append([float(j) for j in list(measurement)])

    return results