Causal inference with Approximate Bayesian Computation for Multivariate Hawkes Processes