calculate FFT based on an oscillogram
If it helps, you have a sine wave multiplied by a square wave added to a phase shifted sine wave multiplied by another phase shifted square wave.

Addition translates to addition in fourier space, multiplication maps to convolution.

So what you're going to get is the FT of a square wave shifted along the w axis and with an imaginary component, I think.