I am using mathematica to generate plots which is supposed to be the sum of a lot of dirac delta functions centered at different positions. The function smoothkerneldistribution can replace these delta functions with very narrow Gaussians and the graph can be obtained very quickly. But now I want to put a "weight" on each of these Gaussian kernels so that each gaussian has a multiplicative factor associated to it and its contribution in the graph can be tuned. This does not seem to be an option for smoothkerneldistribution and calculating the graph directly as a sum of gaussians takes much longer than using the function for the case without weights (they should be approximately the same because it's adding just a prefactor). Does anyone know any convenient ways to do it?