Umap embedding when working with Harmony on Multiome datasets

[sylestiel]

@rcorces

When working with multiple multiome datasets, if I choose to run Harmony after the step
proj <- addCombinedDims(proj, reducedDims = c("LSI_ATAC", "LSI_RNA"), name = "LSI_Combined")

1. How can I generate/render Harmony based embedding for LSI_ATAC and LSI_RNA

2. How do I ensure that the getMarkerFeatures is based upon the Harmony correction?