Hysteresis model identification by the Flock-of-Starlings Optimization