Printing to std::out suggests the following particle sets exist somewhere: ionset0 ionset01 ionset12 ionset03. Code should only construct ionset0 and ionset1.