def METHOD_NAME(y_true_reference: pd.Series, y_pred_reference: pd.Series) -> Tuple:
    """
    Calculate sampling error components for Root Mean Squared Error (RMSE) using reference data.
    Parameters
    ----------
    y_true_reference: pd.Series
        Target values for the reference dataset.
    y_pred_reference: pd.Series
        Predictions for the reference dataset.
    Returns
    -------
    (std,): Tuple[np.ndarray]
    """
    squared_error = (y_true_reference - y_pred_reference) ** 2
    squared_error_std = np.std(squared_error)
    squared_error_mean = np.mean(squared_error)
    return squared_error_std, squared_error_mean