Capturing beam images with camera objects

The finesse.detectors.camera module includes a number of Camera objects for probing single pixels, slices or full images of beam(s). These cameras are split into categories according to these dimensions as well as another two categories corresponding to whether it is detecting the field amplitude and phase (i.e. a sub-class of ComplexCamera) or the pixel intensity (ie. a sub-class of CCDCamera).

Equations describing the quantities computed by each of these two latter categories can be found in Camera Equations.

Complex (field) cameras

Listed below are the camera objects which are used for probing the per-pixel field amplitude and phase for a beam at a specific frequency.

finesse.detectors.camera.FieldCamera(name, ...)

Camera for detecting the full image of the beam in terms of amplitude and phase.

finesse.detectors.camera.FieldScanLine(name, ...)

Camera for detecting a slice of the beam in terms of amplitude and phase.

finesse.detectors.camera.FieldPixel(name, node)

Camera for detecting a single pixel of the beam in terms of the amplitude and phase.

These are listed in order of decreasing dimensions.

CCD cameras

Listed below are the camera objects which are used for probing the per-pixel beam intensity of all (non-audio) fields present at the detection node.

finesse.detectors.camera.CCD(name, node, ...)

Camera for detecting the full image of the beam in terms of the intensity.

finesse.detectors.camera.CCDScanLine(name, ...)

Camera for detecting a slice of the beam in terms of the intensity.

finesse.detectors.camera.CCDPixel(name, node)

Camera for detecting a single pixel of the beam in terms of the intensity.

These are listed in order of decreasing dimensions.

Note

All camera classes take the optional argument w0_scaled in their constructors. This flag is True by default and indicates that all \(x\) and \(y\) coordinates should be scaled by the waist size of the beam. Any outputs (including plots) will then be in units of the beam waist. By setting this flag to False, you can prevent this behaviour and instead have the \(x\) and \(y\) axes in units of metres.

Detecting a single pixel

A single pixel of a beam profile can be detected with Pixel type cameras. The return type of this detector will be a single number (complex for a field, real for a CCD). A simple example below shows how one can use a FieldPixel to detect, for example, the central pixel (at \(x = y = 0\)) of a beam distribution on reflection from a mirror whilst tilting the mirror in the x-direction.

import finesse
finesse.configure(plotting=True)

IFO = finesse.Model()
IFO.parse("""
l L0 P=1
s s0 L0.p1 ITM.p1 L=1
m ITM R=0.99 T=0.01
gauss gL0 L0.p1.o q=(-2+3j)

# a FieldPixel detector - defaults to x = y = 0 and f = 0
fpx refl_px ITM.p1.o

modes(maxtem=2)
xaxis(ITM.xbeta, lin, 0, 100u, 100)
"""
)

out = IFO.run()
out.plot();
../../_images/imaging_beams_0_0.svg

Note that, as implied in the code above, the constructor for FieldPixel uses default values of zero for x and y (so that the central pixel is the default), and it probes at a frequency offset of 0 Hz (carrier) by default too.

Probing a slice of a beam

Slices of a beam profile are detected with ScanLine type cameras. These detectors return a 1D numpy.ndarray where the scanning axis is specified at construction time via the direction argument. An example is shown below, where the x-axis of a beam image, in terms of the intensity, is computed and plotted using a CCDScanLine detector.

IFO = finesse.Model()
IFO.parse("""
l L0 P=1
s s0 L0.p1 ITM.p1 L=1
m ITM R=0.99 T=0.01
gauss gL0 L0.p1.o q=(-2+3j)

# A CCDScanLine detecting x in [-3, 3] with y = 0 by default
ccdline refl_x ITM.p1.o xlim=3 npts=100

modes(maxtem=2)
"""
)

out = IFO.run()
out.plot();
../../_images/imaging_beams_1_0.svg

A single value can be specified for x, y when either of these is the scanning axis - this will then create the scanning axis array ranging from, e.g., \(-|x|\) to \(+|x|\) if \(x\) is the axis to be scanned. The value for the non-scanned axis will take on a default of zero if not specified.

Capturing a full beam image

Full two-dimensional beam profiles can be captured with Image type cameras. Detectors of this type return a 2D numpy.ndarray. An example is shown below, using the same model as above but now taking the full beam distribution image in terms of the intensity with a CCD detector.

# Remove the CCDScanLine detector from above
IFO.remove(IFO.refl_x)

# Add a CCD with x in [-1, 1] and y in [-1, 1] with 300 pixels per axis
IFO.parse("""
ccd refl_im ITM.p1.o xlim=1 ylim=1 npts=80
""")

out = IFO.run()
out.plot();
../../_images/imaging_beams_2_0.svg

Producing animated beam images

Animated plots of beam distributions can be created by adding a full image detector (i.e. a CCD or FieldCamera) to a model, then scan over a model parameter (using the xaxis command for example) to produce 3D data which will be animated using the parameter scan axis when plotting via ArraySolution.plot().

An example of this is shown below, where a distorted beam is injected into a Fabry-Perot cavity. A cavity scan is then performed whilst a CCD captures the beam image, at each mirror detuning, to produce a mode scan animation.

IFO = finesse.Model()
IFO.parse(
"""
l L0 P=1

s s0 L0.p1 ITM.p1

m ITM R=0.99 T=0.01 Rc=-10
s CAV ITM.p2 ETM.p1 L=1
m ETM R=0.99 T=0.01 Rc=10

cav FP ITM.p2.o

tem(L0, n=0, m=1, factor=0.9)
tem(L0, m=1, n=0, factor=0.9)
tem(L0, n=0, m=2, factor=0.8)
tem(L0, n=1, m=1, factor=0.4)
tem(L0, n=2, m=0, factor=0.5)

ccd trns ETM.p2.o xlim=2.5 ylim=2.5 npts=80

modes(maxtem=2)
xaxis(ITM.phi, lin, 0, 180, 60)
"""
)

out = IFO.run()
figures, animations = out.plot(cmap="hot")