Yes, you can extend this computation to include other curvature tensors. You just need to add more replacement rules changing the abstract tensor (say RiemannCD[-a, -b, -c, d]) by the corresponding explicit tensor (say CTensor[array of components, bases][-a, -b, -c, d]). The Ricci and RicciScalar rules are already present in the notebook I sent, but not the Riemann and Einstein rules.
Jose.