Unfortunately, there is no software with such function built-in, nevertheless, consider that in terms of velocity the function is just a power of the velocity multiplied by a gaussian function. If you multiply your data by an opportune power of the velocity (which depends on the type of the detector you are using), you end up with something that could be deconvoluted with simple gaussian functions.

More info are available in http://www.sciencedirect.com/science/article/pii/S0169433212015759

I don't know too much about this Coulomb Boltzmann function, but if you have access you should try Mathematica. It has amazing features including convolutions, deconvolutions, fourier transforms, etc.

Using Origin we have to define a user function which is not a simple one. But there is another possibility which is described in a paper entitled Gaussian energy distribution of fast ions emitted by laser-produced plasmas, Applied Surface Science 272 (2013) 46– 49. The trick consists in multiplication of the ion collector signal with a term v-3 (velocity), i.e. time-of-flight t^3. Then a Gaussian function can be applied for deconvolution.

I am very much thankful to Domenico Delle Side, Dragos Palade and Josef Krasa for helping me with valuable information .