get EICs for adducts, isotopes and neutral loss masses
multiEICplus(adducts = c(0, 1, 2, 3), mz, rt, ...)
adducts
numeric() of mass shifts
...
all other arguments passed on to multiEIC