In this talk, I will discuss a novel strategy to fit experimental data using an amplitude ansatz satisfying the constraints of Analyticity, Crossing, Unitarity, and UV completeness. The fit strategy requires both the use of S-matrix Bootstrap methods and non-convex Particle Swarm Optimization (PSO) techniques. As a proof of principle, I will focus on $\pi\pi$ scattering. Using this procedure, I will show how to construct numerically a full-fledged scattering amplitude that fits the available experimental and lattice data, and that features all the known QCD spectrum with quantum numbers $I^G=0^+,1^+$ below 1.4 GeV, plus an additional surprise.