function p=parambLCGT(rse,gfactor,PitYaw)
% Parameter set for bLCGT
% rse: 'BRSE' or 'DRSE'
% gfactor: 'positive' for positive g-factor, 'negative' for negative g-factor
% 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

%% MODEL NAME
p.modelName = ['bLCGT-',rse,'-',gfactor,'-',PitYaw];
if strfind(rse,'BRSE')>0
    p.DRSE=0;
elseif strfind(rse,'DRSE')>0
    p.DRSE=1;
end
if strfind(gfactor,'positive')>0
    ITMChr=1/7125;
    ETMChr=1/7125;
    PRMChr=1/-42.67;
    SRMChr=1/-46.61;
elseif strfind(gfactor,'negative')>0
    ITMChr=1/1900;
    ETMChr=1/1900;
    PRMChr=1/458.128519465;
    SRMChr=1/458.128519465;
end

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

%% LASER
p.Pin=77.5;

%% RF MODULATION
p.fmod1=16.880962*MHz;
p.fmod2=45.015898*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
Nsilica = 1.44967;
TMNmd = 1.754;  %reflaction index of test masses

%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 = Nsilica; %Index of refraction

%ITMs
p.ITMXaio = 0;
p.ITMXChr = ITMChr; % 1/ROC of ITMX
p.ITMXThr = 0.004; %ITMX transmission
p.ITMXLhr = 45*ppm; %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 = ITMChr; % 1/ROC of ITMY
p.ITMYThr = 0.004; %ITMY transmission
p.ITMYLhr = 45*ppm; %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 = ETMChr; % 1/ROC of ETMX
p.ETMXThr = 10*ppm; %ETMX transmission  (for TRX)
p.ETMXLhr = 45*ppm; %ETMX HR Loss
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 = ETMChr; % 1/ROC of ETMY
p.ETMYThr = 10*ppm; %ETMY transmission  (for TRY)
p.ETMYLhr = 45*ppm; %ETMY HR Loss
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 = PRMChr; % 1/ROC of PRM
p.PRMThr = 0.1; %PRM transmission (1 for iLCGT)
p.PRMLhr = 0; %PRM HR Loss
p.PRMRar = 0; %PRM AR Reflection
p.PRMLmd = 0; % Ignore the substrate loss.
p.PRMNmd = Nsilica; %Index of refraction

%PR2
p.PR2aio = 0.6860;
p.PR2Chr = 1.0/-3.0764084715; % 1/ROC of PR2
p.PR2Thr = 500*ppm; %PR2 transmission  (for POP)
p.PR2Lhr = 0; %PR2 HR Loss
p.PR2Rar = 0; %PR2 AR Reflection (ignore).
p.PR2Lmd = 0; % Ignore the substrate loss.
p.PR2Nmd = Nsilica; %Index of refraction

%PR3
p.PR3aio = 0.6860;
p.PR3Chr = 1.0/24.9164838708; % 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 = Nsilica; %Index of refraction

%SRM
p.SRMaio = 0;
p.SRMChr = SRMChr; % 1/ROC of SRM
p.SRMThr = 0.1536; %SRM transmission (1 for iLCGT)
p.SRMLhr = 0; %SRM HR Loss
p.SRMRar = 0; %SRM AR Reflection
p.SRMLmd = 0; % Ignore the substrate loss.
p.SRMNmd = Nsilica; %Index of refraction

%SR2
p.SR2aio = 0.6860;
p.SR2Chr = 1.0/-2.98718007727; % 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 = Nsilica; %Index of refraction

%SR3
p.SR3aio = 0.6860;
p.SR3Chr = 1.0/24.9164838708; % 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 = Nsilica; %Index of refraction

%% OTHER MIRRORS
% input steering mirrors
p.IN1aio=45;
p.IN1Chr=0;
p.IN1Thr=0;
p.IN1Lhr=0;
p.IN1Rar=0;
p.IN1Lmd=0;
p.IN1Nmd=Nsilica;

p.IN2aio=45;
p.IN2Chr=0;
p.IN2Thr=0;
p.IN2Lhr=0;
p.IN2Rar=0;
p.IN2Lmd=0;
p.IN2Nmd=Nsilica;

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

% half mirrors 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=Nsilica;

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

%PRC lengths
p.LPRM_PR2 = 14.7614883609; %Distance between PRM and PR2
p.LPR2_PR3 = 11.0660636313; %Distance between PR2 and PR3
p.LPR3_BS = 15.7637759961; %Distance between PR3 and BS

%SRC lengths
p.LSRM_SR2 = 14.7412236675; %Distance between SRM and SR2
p.LSR2_SR3 = 11.1115048847; %Distance between SR2 and SR3
p.LSR3_BS = 15.7385994361; %Distance between SR3 and BS

% other lengths are all 0

%% ATTENUATORS
% Tuned (about 50mW at each probe (except AS))
p.AttPOP=0.73;
p.AttREFL=0.76;
p.AttAS=0;
p.AttTRX=0.974;
p.AttTRY=0.974;

%% OPERATING POINTS
% 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;
p.posOffsetPRM=0;
p.posOffsetSRM=lambda/4+p.DRSE*p.dTune;    % depending on RSE configuration
p.armOffset=1.3*pm;        % DC readout

%% MECHANICAL TRANSFER FUNCTIONS (torque on TM to angle of TM)
% fitted from bLCGT suspension simulation data by T. Sekiguchi (see fitTFbypeaks.m)
if strfind(PitYaw,'pitch')>0
    ZQ=[0.8188 4e2
        0.8900 4e2
        1.0225 4e2
        1.2983 4e2
        1.8017 4e2
        5.1185 4e2
       16.5220 4e2];
    PQ=[0.7618 4e2
        0.8560 4e2
        0.8999 4e2
        1.1683 4e2
        1.3349 4e2
        1.8320 4e2
        8.7234 4e2
       26.4898 4e2];
    p.tfTMPit = filterdesign('zpk',ZQ,PQ,9.0914);
    ZQ=[0.4420 4e2
        0.6134 4e2
        0.6449 4e2
        0.6893 4e2];
    PQ=[0.4372 4e2
        0.4673 4e2
        0.6272 4e2
        0.6742 4e2
        1.1814 4e2];
    p.tfBSPit = filterdesign('zpk',ZQ,PQ,5.3191);
    ZQ=[0.4420 4e2
        0.6134 4e2
        0.6449 4e2
        0.6893 4e2];
    PQ=[0.4372 4e2
        0.4673 4e2
        0.6272 4e2
        0.6742 4e2
        1.1814 4e2];
    p.tfPRMPit = filterdesign('zpk',ZQ,PQ,19.6079);
    p.tfSRMPit = p.tfPRMPit;
    p.tfPR2Pit = p.tfPRMPit;
    p.tfSR2Pit = p.tfPRMPit;
    p.tfPR3Pit = p.tfPRMPit;
    p.tfSR3Pit = p.tfPRMPit;
elseif strfind(PitYaw,'yaw')>0
    ZQ=[0.0420 0.5e1
        0.0764 4e2
        0.1377 4e2
        0.6134 4e2
        1.7044 4e2];
    PQ=[0.0272 0.5e1
        0.0554 0.5e1
        0.0799 4e2
        0.3481 4e2
        1.5000 4e2
        2.0472 4e2];
    p.tfTMPit = filterdesign('zpk',ZQ,PQ,9.0909);
    ZQ=[0.4182 4e2
       1.2628 4e2];
    PQ=[0.0410 0.5e1
        1.0225 4e2
        1.4752 4e2];
    p.tfBSPit = filterdesign('zpk',ZQ,PQ,5.3191);
    ZQ=[0.4182 4e2
       1.2628 4e2];
    PQ=[0.0410 0.5e1
        1.0225 4e2
        1.4752 4e2];
    p.tfPRMPit = filterdesign('zpk',ZQ,PQ,19.6078);
    p.tfSRMPit = p.tfPRMPit;
    p.tfPR2Pit = p.tfPRMPit;
    p.tfSR2Pit = p.tfPRMPit;
    p.tfPR3Pit = p.tfPRMPit;
    p.tfSR3Pit = p.tfPRMPit;
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=60.7;      % demod phase for POP f1 demodulation I phase
p.demodphasePOP2=12.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+0.17;  % fine-tuned
p.demodphaseREFL2=L_REFL*omegamod2/c/pi*180+0.62;  % 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
if strfind(gfactor,'positive')>0
    p.GouyPOPA=-45.22/180*pi;
    p.GouyPOPB=38.76/180*pi;
    p.GouyREFLA=89.10/180*pi;
    p.GouyREFLB=88.41/180*pi;
    p.GouyASA=2.42/180*pi;
    p.GouyASB=-87.59/180*pi;
    p.GouyTRA=62.65/180*pi;
    p.GouyTRB=p.GouyTRA+0.5*pi;
elseif strfind(gfactor,'negative')>0
% default
%     p.GouyPOPA=-14.0/180*pi;
%     p.GouyPOPB=-73.2/180*pi;
%     p.GouyREFLA=88.7/180*pi;
%     p.GouyREFLB=-86.8/180*pi;
%     p.GouyASA=5.8/180*pi;
%     p.GouyASB=-85.4/180*pi;
%     p.GouyTRA=-64.1/180*pi;
%     p.GouyTRB=p.GouyTRA+0.5*pi;

% PRC615_SRC220
%     p.GouyPOPA=42.0/180*pi;
%     p.GouyPOPB=-35.6/180*pi;
%     p.GouyREFLA=86.7/180*pi;
%     p.GouyREFLB=87.1/180*pi;
%     p.GouyASA=5.3/180*pi;
%     p.GouyASB=-83.0/180*pi;
%     p.GouyTRA=-64.1/180*pi;
%     p.GouyTRB=p.GouyTRA+0.5*pi;

% PRC200_SRC770
% PRM: REFL_B1I for this Gouy phase pair
%     p.GouyPOPA=16.0/180*pi;
%     p.GouyPOPB=-73.2/180*pi;
%     p.GouyREFLA=88.7/180*pi;
%     p.GouyREFLB=-1.3/180*pi;
%     p.GouyASA=0.1/180*pi;
%     p.GouyASB=-89.9/180*pi;
%     p.GouyTRA=-64.1/180*pi;
%     p.GouyTRB=p.GouyTRA+0.5*pi;

% PRC735_SRC760
%     p.GouyPOPA=61.0/180*pi;
%     p.GouyPOPB=-26.0/180*pi;
%     p.GouyREFLA=84.6/180*pi;
%     p.GouyREFLB=84.9/180*pi;
%     p.GouyASA=0.1/180*pi;
%     p.GouyASB=-89.0/180*pi;
%     p.GouyTRA=-64.2/180*pi;
%     p.GouyTRB=p.GouyTRA+0.5*pi;

% PRC165_SRC175
%     p.GouyPOPA=-13.0/180*pi;
%     p.GouyPOPB=-76.3/180*pi;
%     p.GouyREFLA=88.7/180*pi;
%     p.GouyREFLB=-85.6/180*pi;
%     p.GouyASA=6.8/180*pi;
%     p.GouyASB=-84.1/180*pi;
%     p.GouyTRA=-64.1/180*pi;
%     p.GouyTRB=p.GouyTRA+0.5*pi;
    
% PRC350_SRC250    
%     p.GouyPOPA=11.0/180*pi;
%     p.GouyPOPB=-59.1/180*pi;
%     p.GouyREFLA=88.3/180*pi;
%     p.GouyREFLB=-89.2/180*pi;
%     p.GouyASA=4.4/180*pi;
%     p.GouyASB=-85.6/180*pi;
%     p.GouyTRA=-64.1/180*pi;
%     p.GouyTRB=p.GouyTRA+0.5*pi;
% 20120612
%    p.GouyPOPA=-17.0/180*pi;
%    p.GouyPOPB=-76.4/180*pi;
%    p.GouyREFLA=88.4/180*pi;
%    p.GouyREFLB=-85.7/180*pi;
%    p.GouyASA=6.7/180*pi;
%    p.GouyASB=-83.6/180*pi;
%    p.GouyTRA=-61.4/180*pi;
%    p.GouyTRB=p.GouyTRA+0.5*pi;   

% final design, armoffset 1.3pm
    if p.DRSE;
        p.GouyPOPA=-38.0/180*pi;
        p.GouyPOPB=-76.5/180*pi;
        p.GouyREFLA=88.3/180*pi;
        p.GouyREFLB=85.0/180*pi;
        p.GouyASA=-22.9/180*pi;
        p.GouyASB=-83.4/180*pi;
        p.GouyTRA=-62.7/180*pi;
        p.GouyTRB=p.GouyTRA+0.5*pi;
    else
        p.GouyPOPA=-8.0/180*pi;
        p.GouyPOPB=-76.4/180*pi;
        p.GouyREFLA=88.3/180*pi;
        p.GouyREFLB=-88.4/180*pi;
        p.GouyASA=6.7/180*pi;
        p.GouyASB=-83.7/180*pi;
        p.GouyTRA=-61.4/180*pi; %+0.33*pi;
        p.GouyTRB=p.GouyTRA+0.5*pi;
    end
end

%% 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','PRM','SR3','SR2','SRM'};

% BSM probes (pay attention to the order! same as driveNames)
p.BSMprobeNames={'ETMX_BSM','ETMY_BSM','ITMX_BSM','ITMY_BSM','BS_BSMA','BS_BSMB','PR3_BSMA','PR3_BSMB','PR2_BSMA','PR2_BSMB','PRM_BSM','SR3_BSMA','SR3_BSMB','SR2_BSMA','SR2_BSMB','SRM_BSM'};  

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

% 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
if strfind(gfactor,'positive')>0
    p.selectedprobeNames = {'REFL_A2I','TRX_ADC','AS_A1Q','TRY_ADC','POP_B1Q','POP_A2Q','POP_ADC','REFL_BDC','POP_A1I','AS_BDC'};
%    p.selectedprobeNames = {'REFL_A2I','TRX_ADC','AS_A1Q','TRY_ADC','POP_B1Q','REFL_B1I','POP_ADC','REFL_ADC','POP_A1I','AS_BDC'};
elseif strfind(gfactor,'negative')>0
                            % CS,      CH,        DS,       DH,      BS,      PR3,      PR2,      PRM,       SR3,      SRM
    p.selectedprobeNames = {'TRX_ADC','REFL_A2I','TRY_ADC','AS_A1Q','POP_A1Q','POP_A2Q','POP_BDC','REFL_BDC','POP_B1I','AS_BDC'};
%    p.selectedprobeNames = {'TRX_ADC','REFL_A2I','TRY_ADC','AS_A1Q','POP_A1Q','POP_A2Q','POP_BDC','REFL_BDC','POP_B1I','AS_ADC','AS_BDC'};
end

% 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_Hrad_TFdata='./suspensionTF/oplev_120314/120314_TypeA_TMp.txt';  % torque on TM to angle of TM
    p.TM_Hact_TFdata='./suspensionTF/oplev_120314/120314_TypeA_RMp.txt';  % torque from RM to angle of TM
    p.TM_seism_TFdata='./suspensionTF/oplev_120314/120313_typeA_p.txt'; % seismic noise of TM
    p.BS_Hrad_TFdata='./suspensionTF/oplev_120314/120314_TypeB_TMp.txt';
    p.BS_Hact_TFdata='./suspensionTF/oplev_120314/120314_TypeB_RMp.txt';
    p.BS_seism_TFdata='./suspensionTF/oplev_120314/120311_TypeB_p.txt';
    p.RC_Hrad_TFdata='./suspensionTF/oplev_120314/120314_TypeB_TMp.txt';
    p.RC_Hact_TFdata='./suspensionTF/oplev_120314/120314_TypeB_RMp.txt';
    p.RC_seism_TFdata='./suspensionTF/oplev_120314/120311_TypeB_p.txt';
elseif strfind(PitYaw,'yaw')>0
    p.TM_Hrad_TFdata='./suspensionTF/oplev_120314/120314_TypeA_TMy.txt';  % torque on TM to angle of TM
    p.TM_Hact_TFdata='./suspensionTF/oplev_120314/120314_TypeA_RMy.txt';  % torque from RM to angle of TM
    p.TM_seism_TFdata='./suspensionTF/oplev_120314/120313_typeA_y.txt'; % seismic noise of TM
    p.BS_Hrad_TFdata='./suspensionTF/oplev_120314/120314_TypeB_TMy.txt';
    p.BS_Hact_TFdata='./suspensionTF/oplev_120314/120314_TypeB_RMy.txt';
    p.BS_seism_TFdata='./suspensionTF/oplev_120314/120311_TypeB_y.txt';
    p.RC_Hrad_TFdata='./suspensionTF/oplev_120314/120314_TypeB_TMy.txt';
    p.RC_Hact_TFdata='./suspensionTF/oplev_120314/120314_TypeB_RMy.txt';
    p.RC_seism_TFdata='./suspensionTF/oplev_120314/120311_TypeB_y.txt';
end