function p=paramiLCGT(PitYaw)
% Parameter set for iLCGT
% uses parambLCGT.m for common parameters
% PitYaw: 'pitch' for pitch, 'yaw' for yaw
% See http://gw.icrr.u-tokyo.ac.jp/JGWwiki/LCGT/subgroup/ifo/MIF/OptParam  (最終更新日時 2011-04-19 17:43:58 更新者 YoichiAso)
% Check demodulation phases and Gouy phases!
% Important vectors and matrices for ASC are defined below!
% Author: Yuta Michimura

%% COMMON PARAMETERS
p=parambLCGT('BRSE','negative',PitYaw);

%% UNITS
ppm=1e-6;
MHz=1e6;
lambda=1.0640e-06;
pm=1e-12;

%% MODEL NAME
p.modelName = ['iLCGT-',PitYaw];

%% LASER
p.Pin=2;

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

%% FPMI MIRRORS
TMNmd = 1.44967;

ITMChr=1/14000;
ETMChr=1/7500;
p.ITMXChr = ITMChr;
p.ITMYChr = ITMChr;
p.ETMXChr = ETMChr;
p.ETMYChr = ETMChr;

%ITMs
p.ITMXNmd = TMNmd; %Index of refraction
p.ITMYNmd = TMNmd; %Index of refraction

%ETMs
p.ETMXNmd = TMNmd; %Index of refraction
p.ETMYNmd = TMNmd; %Index of refraction


%% RECYCLING CAVITY MIRRORS
%PRM
p.PRMThr = 1; %PRM transmission
p.PRMNmd = 1; %Index of refraction

%PR3
p.PR3Thr = 0; %PR3 transmission  (no POP for iLCGT)

%SRM
p.SRMThr = 1; %SRM transmission
p.SRMNmd = 1; %Index of refraction

%% OTHER MIRRORS

%% LENGTHS
%Michelson part
p.Las=3.3299;  %Schnupp asymmetry
p.LMIavg=25+26.718;  %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=3000-26.718-34.932;


%% ATTENUATORS
% Tuned (about 50mW at each probe (except AS))
p.AttPOP=0;  % no POP for iLCGT
p.AttREFL=0.9945;
p.AttAS=0;
p.AttTRX=0;
p.AttTRY=0;

%% MECHANICAL TRANSFER FUNCTIONS
% fitted from iLCGT suspension simulation data by T. Sekiguchi (see fitTFbypeaks.m)
if strfind(PitYaw,'pitch')>0
    ZQ=[0.3837 1e3];
    PQ=[0.3516 1e3
        1.3243 1e3];
    p.tfTMPit = filterdesign('zpk',ZQ,PQ,19.6161);
elseif strfind(PitYaw,'yaw')>0
    ZQ=[0.0210 1e3
        0.0695 1e3];
    PQ=[0.0175 1e3
        0.0216 1e3
        0.6281 1e3];
    p.tfTMPit = filterdesign('zpk',ZQ,PQ,19.6243);
end

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

%% DEMOD AND GOUY PHASES
% Tuned using LCGT_ASC.m(and DemodGouy.m)
% Check again when lengths or modulation frequencies are changed.
% Some constants for calculating demodulation phases
c=299792458;
omegamod1=2*pi*p.fmod1;
omegamod2=2*pi*p.fmod2;

% Demod phases are in degrees
p.demodphasePOP1=0;      % demod phase for POP f1 demodulation I phase   (no POP for iLCGT)
p.demodphasePOP2=0;     % demod phase for POP f2 demodulation I phase
L_REFL=2*(p.LPRM_PR2+p.LPR2_PR3+p.LPR3_BS+p.LMIavg);
p.demodphaseREFL1=L_REFL*omegamod1/c/pi*180;  % fine-tuned
p.demodphaseREFL2=L_REFL*omegamod2/c/pi*180+0.3;  % fine-tuned
L_AS=p.LPRM_PR2+p.LPR2_PR3+p.LPR3_BS+2*p.LMIavg+p.LSR3_BS+p.LSR2_SR3+p.LSRM_SR2;
p.demodphaseAS1=L_AS*omegamod1/c/pi*180+0.2;      % fine-tuned

% Gouy phases are in radians
p.GouyPOPA=0/180*pi;
p.GouyPOPB=0/180*pi;
p.GouyREFLA=69.7/180*pi;
p.GouyREFLB=-4.3/180*pi;
p.GouyASA=70.5/180*pi;
p.GouyASB=p.GouyASA+0.5*pi;
p.GouyTRA=64.5/180*pi;
p.GouyTRB=p.GouyTRA+0.5*pi;

%% SET VECTORS and MATRICES for ASC
% Mirrors to drive (put BSM probes for each of them!)
p.driveNames={'ETMX','ETMY','ITMX','ITMY','BS','PR3','PR2'};

% BSM probes (pay attention to the order! same as driveNames) (A and B for BeamSplitters)
p.BSMprobeNames={'ETMX_BSM','ETMY_BSM','ITMX_BSM','ITMY_BSM','BS_BSMA','BS_BSMB','PR3_BSMA','PR3_BSMB','PR2_BSMA','PR2_BSMB'};  

% DOFs to control
p.cDrvNames={'CS','CH','DS','DH','BS','PR3','PR2'};

% Matrix for transforming MIRROR basis to DOF basis (see LIGO-T0900511 section 3.2, LIGO-T080186 eqs. (3.1)-(3.2))
p.gITM=1-p.Larm*p.ITMXChr;
p.gETM=1-p.Larm*p.ETMXChr;
r=2/((p.gITM-p.gETM)+sqrt((p.gITM-p.gETM)^2+4));
p.mDrv = eye(length(p.cDrvNames),length(p.driveNames));
                %  EX  EY  IX  IY 
p.mDrv(1:4,1:4) = [ 1,  1,  r,  r;  %CSOFT
                    r,  r, -1, -1;  %CHARD
                    1, -1,  r, -r;  %DSOFT
                    r, -r, -1,  1];  %DHARD

% selected probes for getting error signals
p.selectedprobeNames = {'REFL_A1I','TRX_ADC','AS_A1Q','TRY_ADC','REFL_B1Q','REFL_B1I'};

% DC probes for monitoring DC power at selected probes
p.selectedprobeDC = p.selectedprobeNames;
for kk=1:length(p.selectedprobeNames)
    p.selectedprobeDC{kk}=[p.selectedprobeNames{kk}(1:end-2),'DC'];
end

%% SUSPENSION DATA
% TM actuator TF data
if strfind(PitYaw,'pitch')>0
    p.TM_Hact_TFdata='../iLCGT_ASC/seism_angle_110513/TypeA/typeA_tf_wRM_pitch_110513.txt';
elseif strfind(PitYaw,'yaw')>0
    p.TM_Hact_TFdata='../iLCGT_ASC/seism_angle_110513/TypeA/typeA_tf_wRM_yaw_110513.txt';
end