def METHOD_NAME(self, system: System):
    """
    Initialize the system according to the instructions given in _setup_momenta.
    Args:
        system (object): System class containing all molecules and their replicas.
    """
    if self.temperature.shape[0] != 1:
        if self.temperature.shape[0] != system.n_molecules:
            raise InitializerError(
                "Initializer requires either a single temperature or one per molecule."
            )
    if self.remove_com:
        system.remove_center_of_mass()
    if self.wrap_positions:
        system.wrap_positions()
    self._setup_momenta(system)
    if self.remove_translation:
        system.remove_translation()
    if self.remove_rotation:
        system.remove_com_rotation()