% Optickle input file for LCGT
%
% original file from adVirgo1_sd_hp.m
%

%%% added imperfections 110624
%
% ITM reflectivity difference, BS reflectivity, OMC loss
%
%%%

function opt = lcgt

  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  % PARAMETERS
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  % input power 
%  Pin = 75; % without SBs
  Pin = 88; % temp
  
  % modulation frequencies
  fMod1 =  16.875e6;
  fMod2 =  45e6;
  fMod3 =  39.375e6;
  
  % Lengths
  Lcav  = 3006.69;
  Lprc  =   73.2826;
  Delta =    3.33103; % asymmetry length
  Lsrc  =   73.2826;
  Lprbs =   53.2826; % PRM-BS distance
  
  % Mechanics
  w = 2 * pi * 0.8;      % pendulum resonance frequency
  mI = 30;               % mass of input mirror
  mE = 30;               % mass of end mirror
  mPRM = 10.8;
  mSRM = 10.8;
  mBS  = 30;
  dampRes = [0.001 + 1i, 0.001 - 1i];
  
  % mirror parameters
  lIMX = 41e-6;
  lEMX = 41e-6;
  lIMY = 49e-6;
  lEMY = 49e-6;
  lPRM = 45e-6;
  lSRM = 45e-6;
  lBS  = 45e-6;

  tIMX = 1 - 0.99602 - lIMX; % 0.5% error
  tEMX = 10e-6;
  tIMY = 1 - 0.99598 - lIMY; % 0.5% error
  tEMY = 10e-6;
  tPRM = 0.10;
  tSRM = 1 - 0.92^2 - lSRM;
  tBS  = 0.505; % 1% error
  
  rarBS  = 100e-6;
  rarIMX = 1000e-6 + 600e-6; % AR loss + ITM absorption
  rarIMY = 1000e-6 + 600e-6;
  %rarIMX = 100e-6 + 600e-6;
  %rarIMY = 100e-6 + 600e-6;
  
  % Operating point
  offs_DARM  = -7.5e-13;       % DARM offset in meters
%  offs_DARM  = 15.0e-13;       % DARM offset in meters
  tuning_SRC = 2*(pi/2-1.5088);        % SRC tuning in radians
  %tuning_SRC = 0.14;        % SRC tuning in radians
  
  % Demodulation phases, tuned later
  SYM_SB1_PHASE = 0;
  SYM_SB2_PHASE = 0;
  POBS_SB1_PHASE = 0;
  POBS_SB2_PHASE = 0;

  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  % create model
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  vFrf = [-fMod2, -fMod3, -fMod1, 0, fMod1, fMod3, fMod2];
  opt = Optickle(vFrf);
  
  % laser source (no IMC)
  opt = addSource(opt, 'Laser', sqrt(Pin) * (vFrf == 0));
  
  % Modulators
  gamma = 0.1;
  opt = addRFmodulator(opt, 'Mod1', fMod1, 1i * gamma);
  opt = addRFmodulator(opt, 'Mod2', fMod2, 1i * gamma);
  opt = addRFmodulator(opt, 'Mod3', fMod3, 1i * gamma * 0);
  
  opt = addLink(opt, 'Laser', 'out', 'Mod1', 'in', 0.1);
  opt = addLink(opt, 'Mod1',  'out', 'Mod2', 'in', 0.1);
  opt = addLink(opt, 'Mod2',  'out', 'Mod3', 'in', 0.1);
  
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%% X cavity
  
  % input mirror with pick-off port
  opt = addBeamSplitter(opt, 'IMX_AR', 1, 0, 1-rarIMX, 0, 0);
  %                          angle, 1/Rcurv, THR,  LHR
  opt = addMirror(opt, 'IMX', 0,    1/500000, tIMX, lIMX);
  %                          angle, 1/Rcurv, THR,  LHR
  opt = addMirror(opt, 'EMX', 0,    1/7000, tEMX, lEMX);

  opt = addLink(opt, 'IMX_AR', 'bkA', 'IMX', 'bk', 0);
  opt = addLink(opt, 'IMX', 'bk', 'IMX_AR', 'bkB', 0);
  
  opt = addLink(opt, 'IMX', 'fr', 'EMX', 'fr', Lcav);
  opt = addLink(opt, 'EMX', 'fr', 'IMX', 'fr', Lcav);
  opt = setCavityBasis(opt, 'IMX', 'EMX');
  
  
  opt = setMechTF(opt, 'IMX', zpk([], -w * dampRes, 2 / mI)); % changed by Kentaro 110621
  opt = setMechTF(opt, 'EMX', zpk([], -w * dampRes, 2 / mE));

  %%%%%%%%%%%%%%%%%%%%%%%%%%% Y cavity
  % input mirror with pick-off port
  %                          angle, 1/Rcurv, THR,  LHR
  opt = addBeamSplitter(opt, 'IMY_AR', 1, 0, 1-rarIMY, 0, 0);
  opt = addMirror(opt, 'IMY', 0,    1/500000, tIMY, lIMY);
  opt = addMirror(opt, 'EMY', 0,    1/7000, tEMY, lEMY);

  opt = addLink(opt, 'IMY_AR', 'bkA', 'IMY', 'bk', 0);
  opt = addLink(opt, 'IMY', 'bk', 'IMY_AR', 'bkB', 0);
  
  opt = addLink(opt, 'IMY', 'fr', 'EMY', 'fr', Lcav);
  opt = addLink(opt, 'EMY', 'fr', 'IMY', 'fr', Lcav);
  opt = setCavityBasis(opt, 'IMY', 'EMY');
  
  opt = setMechTF(opt, 'IMY', zpk([], -w * dampRes, 2 / mI)); % changed by Kentaro 110621
  opt = setMechTF(opt, 'EMY', zpk([], -w * dampRes, 2 / mE));
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%% Recycling cavities
  % power recycling
  opt = addMirror(opt, 'PRM', 0, 0, tPRM, lPRM);
  % signal recycling
  opt = addMirror(opt, 'SRM', 0, 0, tSRM, lSRM);
  
  % beam splitter
  opt = addBeamSplitter(opt, 'BS', 45, 0, tBS, lBS, rarBS);
  
  opt = setMechTF(opt, 'PRM', zpk([], -w * dampRes, 2 / mPRM)); % changed by Kentaro 110621
  opt = setMechTF(opt, 'SRM', zpk([], -w * dampRes, 2 / mSRM));
  opt = setMechTF(opt, 'BS',  zpk([], -w * dampRes, 2 / mBS));
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%% Links

  opt = addLink(opt, 'Mod3',  'out', 'PRM', 'bk', 1);
  
  opt = addLink(opt, 'PRM',  'fr',  'BS',  'frA', Lprbs);
  opt = addLink(opt, 'BS',   'frB', 'PRM', 'fr', Lprbs);
  
  opt = addLink(opt, 'SRM',  'fr',  'BS', 'bkA', Lsrc - (Lprc - Lprbs));
  opt = addLink(opt, 'BS',   'bkB', 'SRM', 'fr', Lsrc - (Lprc - Lprbs));
  
  opt = addLink(opt, 'BS',  'frA', 'IMY_AR', 'frA', (Lprc - Lprbs) - Delta/2);
  opt = addLink(opt, 'IMY_AR', 'frB',  'BS', 'frB', (Lprc - Lprbs) - Delta/2);
  
  opt = addLink(opt, 'BS',  'bkA', 'IMX_AR', 'frA', (Lprc - Lprbs) + Delta/2);
  opt = addLink(opt, 'IMX_AR', 'frB',  'BS', 'bkB', (Lprc - Lprbs) + Delta/2);
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%% Read-out optics (to reduce power)
  
  opt = addMirror(opt, 'SYM_BS', 0, 0, 0.122);
  opt = addSink(opt, 'SYM');
  opt = addLink(opt, 'PRM', 'bk', 'SYM_BS', 'fr', 2);
  opt = addLink(opt, 'SYM_BS', 'bk', 'SYM', 'in', 2);
  
  opt = addMirror(opt, 'POBS_BS', 0, 0, 1);
  opt = addSink(opt, 'POBS');
  opt = addLink(opt, 'BS', 'piB', 'POBS_BS', 'fr', 2);
  opt = addLink(opt, 'POBS_BS', 'bk', 'POBS', 'in', 2);
    
  %%%%%%%%%%%%%%%%%%%%%%%%%%% ASY optics and OMC
  
  %%% 1 percent loss for ASYp
  opt = addSink(opt, 'ASY');
  opt = addSink(opt, 'ASYs');
  opt = addSink(opt, 'ASYp');
  opt = addMirror(opt, 'ASY_BS', 0, 0, 0.99);
%  opt = addMirror(opt, 'ASY_BS', 0, 0, 0.9999);
  opt = addLink(opt, 'SRM', 'bk', 'ASY_BS', 'fr', 2);
  opt = addLink(opt, 'ASY_BS', 'fr', 'ASYp', 'in', 2);
  
  %%% short triangle OMC cavity, preliminary
  opt = addMirror(opt, 'OMCin',  0, 0, 0.014, 10e-6); % preliminary
  opt = addMirror(opt, 'OMCout', 0, 0, 0.014, 10e-6);
%  opt = addMirror(opt, 'OMCin',  0, 0, 1); % temp
%  opt = addMirror(opt, 'OMCout', 0, 0, 1);
  opt = addMirror(opt, 'OMCend', 0, 0, 0, 10e-6);
  
  opt = addLink(opt, 'ASY_BS', 'bk', 'OMCin', 'bk', 1);
  opt = addLink(opt, 'OMCin', 'fr', 'OMCout', 'fr', 1e-3);
  opt = addLink(opt, 'OMCout', 'fr', 'OMCend', 'fr', 0.1);
  opt = addLink(opt, 'OMCend', 'fr', 'OMCin', 'fr', 0.1);
  
  opt = addMirror(opt, 'ASY_PDloss', 0, 0, 0.90);
  opt = addLink(opt, 'OMCout', 'bk', 'ASY_PDloss', 'fr', 1);
  opt = addLink(opt, 'ASY_PDloss', 'bk', 'ASY', 'in', 1);
  
  nOMC = getDriveIndex(opt, 'OMCend');
  opt = setDriveOffset(opt, nOMC, 1.064e-6 / 4);
  
  %%%%%%%%%%%%%%%%%%%%%%%%%%% Detectors
  
  opt = addProbeAt(opt, 'ASY_DC',    'ASY',  'in',  0, 0);

  opt = addProbeAt(opt, 'SYM_SB1_P', 'SYM', 'in',  fMod1, SYM_SB1_PHASE); 
  opt = addProbeAt(opt, 'SYM_SB2_P', 'SYM', 'in',  fMod2, SYM_SB2_PHASE); 
  opt = addProbeAt(opt, 'SYM_SB1_Q', 'SYM', 'in',  fMod1, SYM_SB1_PHASE + 90); 
  opt = addProbeAt(opt, 'SYM_SB2_Q', 'SYM', 'in',  fMod2, SYM_SB2_PHASE + 90); 
  
  opt = addProbeAt(opt, 'POBS_SB1_P', 'POBS', 'in',  fMod1, POBS_SB1_PHASE); 
  opt = addProbeAt(opt, 'POBS_SB2_P', 'POBS', 'in',  fMod2, POBS_SB2_PHASE); 
  opt = addProbeAt(opt, 'POBS_SB1_Q', 'POBS', 'in',  fMod1, POBS_SB1_PHASE + 90); 
  opt = addProbeAt(opt, 'POBS_SB2_Q', 'POBS', 'in',  fMod1, POBS_SB2_PHASE + 90); 
      
  % Set operating points
  nSRM = getDriveIndex(opt, 'SRM');
  nEMX = getDriveIndex(opt, 'EMX');
  nEMY = getDriveIndex(opt, 'EMY');
  opt = setDriveOffset(opt, nSRM, 1.064e-6 * (1 / 4 + tuning_SRC/(4*pi)));
  opt = setDriveOffset(opt, nEMX,  offs_DARM / 2);
  opt = setDriveOffset(opt, nEMY, -offs_DARM / 2);

save lcgt.mat