def METHOD_NAME(cls):
    cls.system.part.add(
        pos=np.random.random((2 * cls.N0, 3)) * cls.system.box_l,
        type=cls.N0 * [cls.types["A-"], cls.types["H+"]])
    cls.RE.add_reaction(
        gamma=cls.Ka,
        reactant_types=[cls.types["HA"]],
        product_types=[cls.types["A-"], cls.types["H+"]],
        default_charges=cls.charges_dict)