rowlesmr
I'm starting to look into some Fourier transform things for some peak shape things, and have a couple of (initial) questions about how it works in Topas.
0
Is FT_K=0 taken to be at the peak position?
Looking at the Gaussian and Lorentzian examples in test_examples\ft:
1
X-space -> G=(2 Sqrt(Ln(2)/Pi)/w) Exp(-4 Ln(2) (X/w)^2) goes to
FT_K-space -> Gt = Exp(-(Pi w k)^2/(4 Ln(2)))
If I do the same transformation in Maple, I get a slightly different expression
k-space -> Gt = Exp(-(w k)^2/(16 Ln(2)))
which can be taken care of if (Topas k) = 2 Pi (Maple k). Is this correct? Tech ref is unclear on this.
2
X-space -> L=(2/(Pi w)) / (1+4(X/w)^2) goes to
FT_K-space -> Lt = Exp(-Pi w k)
If I do the same transformation in Maple, I get a slightly different expression
k-space -> Gt = Exp(-(1/2) w k) Heaviside(k) + Exp((1/2) w k) Heaviside(-k)
where Heaviside(k) is a piecewise function =0, k < 0, =1 k>0, undef k=0.
Plotting the Topas FT_K-space function over k = -5..5 gives a single exponential, whereas the Maple version gives back-to-back exponentials.
How are convolutions defined in Topas in this instance? If I wanted to write down the fourier transform function, which is correct?
3
WPPM_ft_conv is said to be done on data interpolated to s-space.
How is s-space defined? is it Sin(Th)/Lam? or some other definition?
.
I think that's all for now.
Matthew