[Hirshfeld partitioning](https://doi.org/10.1007/BF00549096) could be easily implemented in IntegratorXX. Here, the weight function is defined by the ratio of electron densities $$ w_A({\bf r}) = \frac {n_A({\bf r})} {\sum_B n_B({\bf r})} $$