Robust Recursive Partitioning for Heterogeneous Treatment Effects with Uncertainty Quantification