We systematically explore the parametrization dependence of the Parton Distribution Functions (PDFs) to better quantify the true uncertainty from global QCD analyses. To achieve this, we employ a novel technique that automates the generation of polynomial parametrizations for PDFs using Bézier curves. This technique is implemented in a C++ module, named Fantômas, which is now included within the xFitter program. As an example, we examine the charged pion PDF. Our analysis reveals that the sea and gluon distributions in the pion are strongly correlated, and PDF solutions featuring a vanishing gluon and a large quark sea are still experimentally allowed.