%% iLCGT model parameters

%BRSE, SDM, use REFL_1I for CARM, with BAE

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

p.Name = 'iLCGT';

%% 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 = 0;

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

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

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

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

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

% Whether to generate vector sensing matrix
p.vectorSensM = 1;
 
%% Load common parameters
% First load the common parameters
% Some of the parameters defined in the file will be overwritten below

COMMON_PARAMS;

%% Tickle setting
p.ftickle = 1e2;

%% DC Readout Setting

%Homodyne phase [deg]
%Arm offset is automatically adjusted to achieve this HD phase.
p.HDphase = 58.2; % = -(121.8-180)

%% 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.posOffsetSRM + p.dTune;
end

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

p.signalPorts = {'ASPO_1I', 'REFL_1I', 'REFL_1Q', 'POP_2I', 'POP_1I'};

p.REFLBS1_R = 0;

%% Control UGFs

p.DARMUGF = 200;
p.CARMUGF = 10e3;
p.MICHUGF = 50;
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

% Power at BS
% The input laser power will be adjusted to achieve this power.
p.Pbs = 780;

%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*10;

%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=4;

% Generate RF vector
p.vMod = generateRFFrequencyVector([p.fmod1,p.fmod2],p.sbOrder,1);

%% FPMI

% PRM is transpalent
p.PRMThr = 1-45*ppm;

% SRM is transpalent
p.SRMThr = 1-45*ppm;

p.ASPOMThr = 0.01; %Pick up all the AS light