I find the answer of this question and writing it here in case that anyone needs it in the future:
When we think of low pass or band pass in the frequency domain, a gaussian low pass filter must just pass the low frequency infomation and what about its counterpasrt in the bandpass? that will be a shift in frequency domain around the f0. Shift in frequency domain is a multiplication in time domain(e^(-jw0t)). I did this buy multiplying my low pass gaussian filter in cos(w0t).