def makeEllipticalFlatTop(self,mx,my,rx,ry) :
print "Intensity:Intensity2D:makeEllipticalFlatTop",mx,my,rx,ry
mx = pl.double(mx)
my = pl.double(my)
rx = pl.double(rx)
ry = pl.double(ry)
self.i = pl.complex64( ((self.xgrid-mx)/rx)**2 + ((self.ygrid-my)/ry)**2 <= 1.0)
def makeEllipticalFlatTop(self, mx, my, rx, ry):
print "Intensity:Intensity2D:makeEllipticalFlatTop", mx, my, rx, ry
mx = pl.double(mx)
my = pl.double(my)
rx = pl.double(rx)
ry = pl.double(ry)
self.i = pl.complex64(((self.xgrid - mx) / rx)**2 +
((self.ygrid - my) / ry)**2 <= 1.0)
def makeGaussian(self, mx, my, sx, sy):
print "Intensity:Intensity2D:makeGaussian", mx, my, sx, sy
self.i = 1 / (sx * sy * 2 * pl.pi) * pl.complex64(
pl.exp(-((self.xgrid - mx)**2 / (2 * sx**2)) -
((self.ygrid - my)**2 / (2 * sy**2))))
def makeGaussian(self,mx,my,sx,sy) :
print "Intensity:Intensity2D:makeGaussian",mx,my,sx,sy
self.i = 1/(sx*sy*2*pl.pi)*pl.complex64(pl.exp(-((self.xgrid-mx)**2/(2*sx**2))-((self.ygrid-my)**2/(2*sy**2))))