We present here a new method to model the shape of the 3-{mu}m absorption band in the reflectance spectra of meteorites and small bodies. The band is decomposed into several OH/H2O components using Exponentially Modified Gaussian (EMG) profiles, as well as possible organic components using Gaussian profiles when present. We compare this model to polynomial and multiple Gaussian profile fits and show that the EMGs model returns the best rendering of the shape of the band, with significantly lower residuals. We also propose as an example an algorithm to estimate the error on the band parameters using a bootstrap method. We then present an application of the model to two spectral analyses of smectites subjected to different H2O vapor pressures, and present the variations of the components with decreasing humidity. This example emphasizes the ability of this model to coherently retrieve weak bands that are hidden within much stronger ones.