function p=paramiLCGT()
% parameter set for Advanced LIGO
% see LIGO-T0900043-10
%
% Author: Yuta Michimura

%% UNITS
ppm=1e-6;
MHz=1e6;

%% LASER
p.Pin=125;    %p.Pin=20;

%% RF MODULATION
p.fmod1=9.099471*MHz;
p.fmod2=45.497355*MHz;
p.vMod=[-p.fmod2;-p.fmod1;0;p.fmod1;p.fmod2];
p.g1=0.1i;  % modulation depth (imaginary number for phase modulation)
p.g2=0.1i;

%% FPMI MIRRORS
%BS
p.BSaio = 45;
p.BSChr = 0; % 1/ROC of BS
p.BSThr = 0.5; %BS transmission
p.BSLhr = 0; %BS HR Loss
p.BSRar = 0; %BS AR Reflection (for POB).
p.BSLmd = 0; % Ignore the substrate loss.
p.BSNmd = 1.45; %Index of refraction

TMNmd=1.45;  %reflaction index of test masses
%ITMs
p.ITMXaio = 0;
p.ITMXChr = 1/1934; % 1/ROC of ITMX
p.ITMXThr = 0.014; %ITMX transmission
p.ITMXLhr = 0; %ITMX HR Loss
p.ITMXRar = 0; %ITMX AR Reflection (for POX).
p.ITMXLmd = 0; % Ignore the substrate loss.
p.ITMXNmd = TMNmd; %Index of refraction

p.ITMYaio = 0;
p.ITMYChr = 1/1934; % 1/ROC of ITMY
p.ITMYThr = 0.014; %ITMY transmission
p.ITMYLhr = 0; %ITMY HR Loss
p.ITMYRar = 0; %ITMY AR Reflection (for POY).
p.ITMYLmd = 0; % Ignore the substrate loss.
p.ITMYNmd = TMNmd; %Index of refraction

%ETMs
p.ETMXaio = 0;
p.ETMXChr = 1/2245; % 1/ROC of ETMX
p.ETMXLhr = 0; %ETMX HR Loss
p.ETMXThr = 0; %ETMX transmission
p.ETMXRar = 0; %ETMX AR Reflection (ignore).
p.ETMXLmd = 0; % Ignore the substrate loss.
p.ETMXNmd = TMNmd; %Index of refraction

p.ETMYaio = 0;
p.ETMYChr = 1/2245; % 1/ROC of ETMY
p.ETMYLhr = 0; %ETMY HR Loss
p.ETMYThr = 0; %ETMY transmission
p.ETMYRar = 0; %ETMY AR Reflection (ignore).
p.ETMYLmd = 0; % Ignore the substrate loss.
p.ETMYNmd = TMNmd; %Index of refraction

%% RECYCLING CAVITY MIRRORS
%PRM
p.PRMaio = 0;
p.PRMChr = -1/10.997; % 1/ROC of PRM
p.PRMThr = 0.003; %PRM transmission
p.PRMLhr = 0; %PRM HR Loss
p.PRMRar = 0; %PRM AR Reflection
p.PRMLmd = 0; % Ignore the substrate loss.
p.PRMNmd = 1.45; %Index of refraction

%PR2
p.PR2aio = 0.79;
p.PR2Chr = -1/4.55; % 1/ROC of PR2
p.PR2Thr = 0; %PR2 transmission
p.PR2Lhr = 0; %PR2 HR Loss
p.PR2Rar = 0; %PR2 AR Reflection (ignore).
p.PR2Lmd = 0; % Ignore the substrate loss.
p.PR2Nmd = 1.45; %Index of refraction

%PR3
p.PR3aio = 0.615;
p.PR3Chr = 1/36.00; % 1/ROC of PR3
p.PR3Thr = 0; %PR3 transmission
p.PR3Lhr = 0; %PR3 HR Loss
p.PR3Rar = 0; %PR3 AR Reflection (ignore).
p.PR3Lmd = 0; % Ignore the substrate loss.
p.PR3Nmd = 1.45; %Index of refraction

%SRM
%No SRM installed for iLCGT

%SR2
p.SR2aio = 0.87;
p.SR2Chr = -1/6.428; % 1/ROC of SR2
p.SR2Thr = 0; %SR2 transmission
p.SR2Lhr = 0; %SR2 HR Loss
p.SR2Rar = 0; %SR2 AR Reflection (ignore).
p.SR2Lmd = 0; % Ignore the substrate loss.
p.SR2Nmd = 1.45; %Index of refraction

%SR3
p.SR3aio = 0.785;
p.SR3Chr = 1/36.00; % 1/ROC of SR3
p.SR3Thr = 0; %SR3 transmission
p.SR3Lhr = 0; %SR3 HR Loss
p.SR3Rar = 0; %SR3 AR Reflection (ignore).
p.SR3Lmd = 0; % Ignore the substrate loss.
p.SR3Nmd = 1.45; %Index of refraction

%% OTHER MIRRORS (using LCGT parameters)
% partially reflecting BS for REFL
p.INaio=45;
p.INChr=0;
p.INThr=1-100*ppm;
p.INLhr=0;
p.INRar=0;
p.INLmd=0;
p.INNmd=1.45;

% pick off mirror for AS
p.ASSPLITaio=45;
p.ASSPLITChr=0;
p.ASSPLITThr=0.05;
p.ASSPLITLhr=0;
p.ASSPLITRar=0;
p.ASSPLITLmd=0;
p.ASSPLITNmd=1.45;

% half mirrors at REFL and AS (for splitting A and B)
p.HALFaio=45;
p.HALFChr=0;
p.HALFThr=0.5;
p.HALFLhr=0;
p.HALFRar=0;
p.HALFLmd=0;
p.HALFNmd=1.45;

%% LENGTHS 
%Michelson part
p.Las=50e-3;  %Schnupp asymmetry
p.LMIavg=(5.3828+5.3328)/2;  %Average length of the Michelson arms 
p.LBS_ITMX=p.LMIavg+p.Las/2; % Michelson X arm
p.LBS_ITMY=p.LMIavg-p.Las/2; % Michelson Y arm

%Arm Cavity Length
p.Larm=3994.5;

%PRC lengths
p.LPRM_PR2 = 16.6037; %Distance between PRM and PR2
p.LPR2_PR3 = 16.1558; %Distance between PR2 and PR3
p.LPR3_BS = 19.5384; %Distance between PR3 and BS

%SRC lengths
p.LSR3_BS = 19.368; %Distance between SR3 and BS
p.LSR2_SR3 = 15.4607; %Distance between SR2 and SR3
p.LSRM_SR2 = 15.726; %Distance between ASSPLIT and SR2

%% Mechanical TFs (using LCGT parameters)
mTM = 30;
rTM = 0.25/2;
dTM = 0.15;
QTM = 1e5; %Mechanical Q of the TM pendulum
wTM = sqrt(9.8/0.4); %Pendulum freq.

InTM = (rTM^2/4 + dTM^2/12)*mTM; %Moment of innertia
wTMPit = 2*pi*1.0; %Pit resonant freq.
QTMPit = 100; %1e5;

p.tfTM = zpk([], -wTM*[1/(2*QTM) + sqrt((1/(2*QTM))^2 - 1), ...
                      1/(2*QTM) - sqrt((1/(2*QTM))^2 - 1)], 1/mTM);
                  
p.tfTMPit = zpk([], -wTMPit*[1/(2*QTMPit) + sqrt((1/(2*QTMPit))^2 - 1), ...
                      1/(2*QTMPit) - sqrt((1/(2*QTMPit))^2 - 1)], 1/InTM);

%% tickle01
p.ftickle01=10; %tickle01-ing frequency