%this file is for setting parameters of MZI 
%% LCGT model parameters

%(BRSE), MZI

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

p.Name = 'MZI_yamamoto';

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

% Set this to 1 for DRSE  
% this is uploaded by main code (see figure "SRCL Sweep")
p.DRSE = 0;

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

%Whether to use MZI or not %have to add correspondent sentence to other cases??
p.useMZI = 1;

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

%whether to use RF AM modulator, for calculating RAM;
p.RAM = 0;

%whether to split beam to two for cancel out RF AM modulator
p.ParaRAM = 0;

% Check the difference in result b/w MZIM1 and MZIM2 for Mid-Fringe
p.MZIM2 = 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;

p.doTickle = 1;

%% additional parameters only for MZI setting

%% DC Readout Setting

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

%% 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= round(16880962.0755/3);

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

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

% some parameters, for carrier "p.lambda" (see COMMON_PARAMS.m)
p.lambda1 = c/p.fmod1;
p.lambda2 = c/p.fmod2;
p.lambda3 = c/p.fmod3;


% set the phase difference[degree] b/w EOM typeA and B
%p.Phi1 = 96.352;
p.Phi1 = 0;
p.Phi3 = 0;

% set offset on path so as to realize phase difference mentioned above
p.f1Off = p.lambda1*p.Phi1/360;
p.f3Off = p.lambda3*p.Phi3/360;

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

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

% RC Lengths   
p.Lprc = 66.5913279884;
p.Lsrc = 66.5913279884;
p.Lsrc = p.Lsrc - 7e5/p.fmod1; %Adjustment to compensate for the DARM offset.
p.Las = 3.32985084757;

% Supposed DC offset from Mid-Fringe
%p.MidOffset = 1.659e-9; % for 1%
%p.MidOffset = 1.706e-10; % for 0.1%
p.MidOffset = 1.6934e-10; % for 0.1% (29Jan2018, ISC)


%% SRC Detuning  %actually I want to set detune phase by the output of MZI.
% Detuning of SRC
% 86.5 deg in Buonanno & Chen convention (i.e. 90deg for non detune)
% this value is uploaded by main code.(see figure "SRCL Sweep")


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}
if ~abs(p.g3)
    % whether to use f3 RF modulator for MZM
    p.f3EOM = 0;
else
    p.f3EOM = 1;
end


if p.f3EOM
    p.signalPorts = {'AS_DC', 'REFL_1I', 'REFL_1DmQ', 'REFL_2DmI', 'REFL_1DmI'};
else
    if p.DRSE
        p.signalPorts = {'AS_DC', 'REFL_2I', 'REFL_1Q', 'POP_2I', 'POP_1I'};
    else
        p.signalPorts = {'AS_DC', 'REFL_1I', 'REFL_1Q', 'POP_2I', 'POP_1I'};
        p.Phi3 = 0;
    end
end
    

% REFL port
%                    |
%                   ATTN
%                    |
%                    v
%                    |
%   PD4 <--/----/----/
%          |    |    |
%          v    v    v
%         PD3  PD2  PD1
%
% PDs are assigned to the following frequencies.
% PD1 -> REFL_1, PD2 -> REFL_2, PD3 -> REFL_1D, PD4 -> REFL_2D

if p.f3EOM
    %use Dm Ports(Lock acquisition phase)
    p.REFLBS1_R = 2.5/3;
    p.REFLBS2_R = 0;
    p.REFLBS3_R = 1.5/2;
else
    if p.DRSE
        p.REFLBS1_R = 0.2;
        p.REFLBS2_R = 1;
        p.REFLBS3_R = 0;
    else
        p.REFLBS1_R = 0;
        p.REFLBS2_R = 0;
        p.REFLBS3_R = 0;
    end
end


%% Scattered Light Parameters

%Whether to include scattered light or not
p.IncludeSCL = 0;

%{
%MICH PO transmissivity
p.MICHX_POThr =1-50*ppm;

%REFL PO transmissivity
p.REFL_POThr =1-50*ppm;
%}
%% Load common parameters 2
COMMON_PARAMS2;