function [opt, f0, Q0, m] = optTrapDual(PIR, ratioP, fDetune, T1IR, T1G)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% create the model
% Detune the IR with length and the green with frequency
    
%Deal with args  
if nargin < 1
    PIR = 1;
end

if nargin < 2
    ratioP = 1;
end

if nargin < 3
    fDetune = 0;
end

if nargin < 5
    T1IR = 0.0008; % T1@1064nm
    T1G  = 0.0008; % T1@532nm
end
    

% RF component vector - no sidebands/modulation
fMod   = 20e6;              % Not actually used in a meaningful way
vFrf   = [0; fDetune];
lambda = 1064e-9 * [1 0.5]; % Second one is green
%lambda = 1064e-9 * [1 1]; % Second one is NOT green

% create model
opt = Optickle(vFrf, lambda);
 
% add a source
PG     = ratioP * PIR; % Green power
powerDistribution = [PIR PG];
opt.addSource('Laser', sqrt(powerDistribution));


%Tramsissivity vector
T1Vec = [T1IR T1G];

%Transmissivity values should go [T lambda; T lambda]
T1 = [T1Vec' lambda'];

% Add mirrors
lCav = 0.9;
opt.addMirror('IX', 0, 0, T1);
opt.addMirror('EX', 0, 0.7 / lCav, 0);

opt.addLink('Laser', 'out', 'IX', 'bk', 0);
opt.addLink('IX', 'fr', 'EX', 'fr', lCav);
opt.addLink('EX', 'fr', 'IX', 'fr', lCav);

% Set  mechanical transfer functions
f0 = 172;         % pendulum resonant frequency
w  = 2 * pi * f0; 
Q0 = 3200;        % pendulum Q
m  = 1e-3;        % mass of end mirror

% Get poles
p1    = - w / (2 * Q0) * (1 - sqrt(1 - 4 * Q0^2));
p2    = - w / (2 * Q0) * (1 + sqrt(1 - 4 * Q0^2));
poles = [p1, p2];

% Build zpk model
pendulumModel = zpk([], poles, 1 / m);
opt.setMechTF('EX', pendulumModel);

% Tell Optickle to use this cavity basis
opt.setCavityBasis('IX', 'EX');

% Add REFL optics
opt.addSink('REFL');
opt.addLink('IX', 'bk', 'REFL', 'in', 2);

% Add REFL probes (this call adds probes REFL_DC, I and Q)
phi = 0;
opt.addReadout('REFL', [fMod, phi]);


end