def METHOD_NAME(self):
    """
    Get the cumulative work array for the nonequilibrium switching.
    The arrays have the indexing [lambda_state, protocol_index, step_index].
    Returns
    -------
    cumulative_works : np.array [2, n_protocols, n_steps_per_protocol]
        The cumulative work for the protocols that were run.
    """
    return self._cumulative_work