Code covered by the BSD License  

Highlights from
OSCAR

image thumbnail

OSCAR

by

 

08 Jul 2008 (Updated )

An optical FFT code to simulate Fabry Perot cavities with arbitrary mirror profiles

Example_cavity_eigenmodes.m
clear all; close all; clear classes
addpath(genpath('Classes'));

disp('---------------------------------------------------------------------------')
disp('                  OSCAR V3.14')
disp('  ')


% Define the grid for the simulation: 64 X 64,10 cm X 10 cm
G1 = Grid(64,0.35);

% Define the incoming beam, it is a dummy beam since it is not used to
% calculate the cavity eigen modes. However it is necessary to define the
% cavity
E_input = E_Field(G1,'w',0.02);

% Define the 2 mirrors, RofC = 2500m, 9 cm in diameter, transmission 1%,
% no loss

% Advanced Virgo
IM = Interface(G1,'RoC',1420,'CA',0.33,'T',0.01);
EM = Interface(G1,'RoC',1672,'CA',0.33,'T',0.01);
C1 = Cavity1(IM,EM,3000,E_input);

% Advanced LIGO
% IM = Interface(G1,'RoC',1934,'CA',0.33,'T',0.01);
% EM = Interface(G1,'RoC',2245,'CA',0.33,'T',0.01);
% C1 = Cavity1(IM,EM,4000,E_input);

% Calculate the kernel of the cavity
C1 = Calculate_RT_mat(C1);

% Display the cavity eigen modes
Display_cavity_modes(C1,'N',20);
  
% Display the Airy peak of the first 30 (lowest clipping loss) eigenmodes:
% Display_cavity_modes(C1,'N',30,'Airy',1)

Contact us