Model-independent description of<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>B</mml:mi><mml:mo>→</mml:mo><mml:mi>π</mml:mi><mml:mi>l</mml:mi><mml:mi>ν</mml:mi></mml:math>decays and a determination of<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo></mml:math>

We propose a new parametrization of the $B\ensuremath{\rightarrow}\ensuremath{\pi}$ vector form factor ${f}_{+}({q}^{2})$ as an expansion in powers of a conformal mapping variable, which satisfies unitarity, analyticity and perturbative QCD scaling. The unitarity constraint is used also for defining the systematic error of the expansion. We fit with the new parametrization the available experimental and theoretical information on exclusive $B\ensuremath{\rightarrow}\ensuremath{\pi}l\ensuremath{\nu}$ decays, making a conservative estimate of the effects of correlations in the systematic and statistical errors of the lattice results. With four parameters to describe ${f}_{+}({q}^{2})$, the systematic error is negligible in the whole semileptonic region. We also obtain $|{V}_{ub}|=(3.50\ifmmode\pm\else\textpm\fi{}0.24)\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}3}$ where, in our approach, the uncertainty is predominantly statistical.