def METHOD_NAME(self):
    for domain in ["negative", "positive"]:
        self.submodels[
            f"{domain} open-circuit potential"
        ] = pybamm.open_circuit_potential.SingleOpenCircuitPotential(
            self.param, domain, "lead-acid main", self.options, "primary"
        )
        self.submodels[
            f"{domain} oxygen open-circuit potential"
        ] = pybamm.open_circuit_potential.SingleOpenCircuitPotential(
            self.param, domain, "lead-acid oxygen", self.options, "primary"
        )