%% LCGT model parameters


%% Units & Constants for convenience
ppm=1e-6;
nm=1e-9;
pm=1e-12;
MHz = 1e6;
percent = 1/100;

c=299792458;
lambda=1.0640e-06;
p.lambda = lambda;
 
%% ==== Edit the following parameters and flags ===========

%% Name
%This string is appended to the result folder name

p.Name = 'SDMforOfs';
%% Flags
% Set this to 1 to include imperfections (like arm asymmetry)
p.notperfect = 1;

% Set this to 1 to have AR reflections
p.ARreflective = 1;

% Set this to 1 to put arm offset for DC readout
p.DCReadout = 1;

% Set this to 1 to fix all optics (to remove optical spring effects)
p.fixedOptics = 1;

% Set this to 1 for DRSE
p.DRSE = 1;

% Whether to use MZ or not
p.useMZ = 0;

if ~p.useMZ
    p.REFL_ATTN = 0.5;
else
    p.REFL_ATTN = 0;
end

% Whether to generate OSA like plots
p.osaplots = 0;

% Whether to generate mirror sweep plots
p.sweepplots = 1;

% Whether to generate vector sensing matrix
p.vectorSensM = 1;

%% Tickle setting
p.ftickle = 1e-9;

%% DC Readout Setting
% Arm offset for DC readout
p.armOffset = p.DCReadout * 1.5*pm;

%% SRC Detuning
% Detuning of SRC
% 86.5 deg in Buonanno & Chen convention (i.e. 90deg for non detune)
p.detunePhase = 86.5;
p.dTune = lambda * (90-p.detunePhase)/360;
if p.DRSE
    p.posOffsetSRM = p.dTune;
else
    p.posOffsetSRM = 0;
end

%% MICH offset
p.MICHoffset = 0*pm;

%% Arm Asymmetries

% Asymmetry in the reflectivities
p.armAsym = -1*percent;

% Asymmetry of the round trip loss
p.ArmLossAsym = -p.notperfect*10*ppm;

p.BSasym = p.notperfect*p.armAsym; %BS reflectivity asymmetry
p.ITMasym = p.notperfect*p.armAsym; %ITM reflectivity asymmetry
p.ETMasym = p.notperfect*p.armAsym; %ETM reflectivity asymmetry

%% Signal extraction ports
%
% Signal ports used for {DARM, CARM, MICH, PRCL, SRCL}

p.signalPorts = {'AS_DC', 'REFL_2I', 'POP_1Q', 'POP_2I', 'POP_1I'};
%% Control UGFs

p.DARMUGF = 200;
p.CARMUGF = 10e3;
p.MICHUGF = 10;
p.PRCLUGF = 50;
p.SRCLUGF = 50;
%% Feed forward

% Error in the feed forward transfer function
p.FeedForwardError = 1*percent;
%% Displacement Noise Requeirements

% Contributions of mirror displacement noises (seismic and thermal) should
% be smaller than the DARM quantum noise by this factor.

p.DispNoiseSafetyFactor = 10;

%% Laser Parameters
p.Pin=82;  %Laser Power

%MC FSR
fMC=11.25*MHz/2;

% RF Sideband frequencies
p.fmod1=fMC*2;  %First modulation Frequency
p.fmod2=fMC*8; %Second modulaiton
p.fmod3=fMC*7;

%Modulation Depth (imaginary number for phase modulation)
p.g1=i*0.1; %f1 
p.g2=i*0.1; %f2
p.g3=0;   %f3

% Number of higher oreder sidebands to be computed 
p.sbOrder=3;

p.sbOrder2 = p.sbOrder -1;
if p.sbOrder2 < 0
    p.sbOrder2 = 0;
end
p.sbOrder3 = p.sbOrder -2;
if p.sbOrder3 < 0
    p.sbOrder3 = 0;
end

%Modulation vector
vMod1=generateRFFrequencyVector(p.fmod1, p.sbOrder, ...
                            p.fmod2, p.sbOrder, p.fmod3, p.sbOrder);
vMod2=generateRFFrequencyVector(...
      p.fmod2-p.fmod1, p.sbOrder2,p.fmod2+p.fmod1, p.sbOrder2);
vMod3=generateRFFrequencyVector(...
      p.fmod3-p.fmod1, p.sbOrder2,p.fmod3+p.fmod1, p.sbOrder2);
vMod4=generateRFFrequencyVector(...
      (p.fmod2-p.fmod1)-p.fmod3, p.sbOrder3,(p.fmod2+p.fmod1)-p.fmod3, p.sbOrder3);
vMod5=generateRFFrequencyVector(...
      (p.fmod2-p.fmod1)+p.fmod3, p.sbOrder3,(p.fmod2+p.fmod1)+p.fmod3, p.sbOrder3); 
vMod6=generateRFFrequencyVector(...
      p.fmod2-p.fmod1, p.sbOrder2,p.fmod2+p.fmod1, p.sbOrder2);
p.vMod = unique([vMod1;vMod2;vMod3;vMod4;vMod5;vMod6]);
%% ===== You usually don't need to edit below =====
% Load common parameters
COMMON_PARAMS;