function L1 = L1IFOparams; % Parameter file for Interferometer LSC loop model % % version 1.0.1 Updated AA & AI filter to match board measurements % version 1.2.3 Cavity pole -> 82.6 Andri elog ~August FOTON_dir = ''; %%%%%%%% frequencies of interest %%%%%%%%%%%%%%%% darm.fl = 9; % lower frequency of band, Hz darm.fu = 8000; % upper frequency of band, Hz darm.npt = 301; % number of points in band darm.fs = 16384; % sampling frequency, Hz darm.ugf = 171.5; % unity gain frequency of loop, Hz %%%%%%%%% parameters of the plant %%%%%%%%%%%%%%% darm.cavpole = 82.6; % cavity pole, Hz darm.pendf0 = 0.76; % pendulum eigenfrequency, Hz darm.pendQ = 10; % pendulum Q darm.ETMXcal = 0.38e-9; % DC calibration of % ETMX, m/count, from % elog 03/11/12 darm.ETMYcal = 0.49e-9; % DC calibration of % ETMY, m/count, from % elog 03/11/12 darm.tdelay = 159e-6; % time delay in loop, sec darm.armlength = 3995.15; % arm length in meters darm.hflowpass = ... zpk([],-2*pi*[33e3 33e3 33e3],... (2*pi*33e3)^3); % RC lowpass in L1 after mixer %%%%%%%%%%% digital filters %%%%%%%%%%%%%%%%%%%%%%%% %% specify vector of DARM digital filters engaged %%% %% Note: FM1 = 0, FM2 = 1, etc %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% darm.filterfile = [FOTON_dir 'L1LSC.txt']; darm.digitalfilters = [0,1,2,3]; darm.susfilterfile1 = [FOTON_dir 'L1SUS_ETMX.txt']; darm.susdigitalfilters1 = [1,3,4]; darm.susfilterfile2 = [FOTON_dir 'L1SUS_ETMY.txt']; darm.susdigitalfilters2 = [1,3,4]; %%%%%%%%%%% digital gains %%%%%%%%%%%%%%%%%%%%%%%%%% darm.darm2etmx = 2.5; % Output matrix: DARM to ETMX darm.darm2etmy = -2.5; % Output matrix: DARM to ETMY darm.itmtrx = 0.00129; % Input matrix: AS_Q to DARM darm.lsc_gain = -1.3; % DARM filter module gain darm.DCcal = darm.ETMYcal * darm.darm2etmy -... % DARM_CTRL cal darm.ETMXcal * darm.darm2etmx; %%%%%%%%%%% snubber component values %%%%%%%%%%%%%%% %% snubber is a series RC (R_snub & C_snub), in %%%% %% parallel with the coil %%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% darm.R_ser = 3000; % Coil Driver series resistor darm.R_snub = 680; % snubber series resistor darm.C_snub = 0.022e-6; % snubber series cap darm.C_cabl = 800e-12; % ribbon cable capacitance darm.R_coil = 22; % OSEM coil resistance darm.H_coil = 3.3e-3; % OSEM coil inductance % Misc Info (mostly unused) -------------------------------------------------------- darm.opgain = 1.2e9; % Calculated watts/meter darm.efficiency = 0.9; % PO fraction * 90/10 split darm.eta = 0.8; % Amps / Volt darm.Zrf = 140; % Tank circuit impedance darm.Grf = 10; % Pre-amp gain darm.Cable = 0.67; % 3.5 dB of loss in 100' of RG-405 darm.PS = 0.707; % 3 dB loss in power splitter darm.Mixer = 0.5; % 6 dB loss in mixer darm.WG = 10^(24/20); % 24 dB of whitening gain darm.AA = 2; % Gain of 2 in single-diff conv darm.ADC = 32768/10; % 16-bit Analog to Digital conversion darm.AS1_Q_GAIN = 0.005; % Compensates some whitening gain darm.electronics_gain = darm.eta *... % Counts / Watt darm.Zrf *... darm.Grf *... darm.Cable *... darm.PS *... darm.Mixer *... darm.WG *... darm.AA *... darm.ADC *... darm.AS1_Q_GAIN; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %----------------------------------------------------------------------------------- % Parameter file for MICH loop model % %%%%%%%% frequencies of interest %%%%%%%%%%%%%%%% mich.fl = 0.9; % lower frequency of band, Hz mich.fu = 1000; % upper frequency of band, Hz mich.npt = 301; % number of points in band mich.fs = 16384; % sampling frequency, Hz mich.ugf = 10; % unity gain frequency of loop, Hz %%%%%%%%% parameters of the plant %%%%%%%%%%%%%%% mich.cavpole = 1e6; % Something must happen here mich.pendf0 = 0.75; % pendulum eigenfrequency, Hz mich.pendQ = 10; % pendulum Q mich.RMcal = 0.38e-9; % DC calibration of RM, m/count mich.BScal = 0.8e-9; mich.tdelay = 150e-6; % time delay in loop, sec mich.schnupp = 0.31; % (ly-lx) mich.hflowpass = ... zpk([],-2*pi*[100e3],... (2*pi*100e3)); % RC lowpass in L1 after mixer %%%%%%%%%%% digital filters %%%%%%%%%%%%%%%%%%%%%%%% %% specify vector of MICH digital filters engaged %%% %% Note: FM1 = 0, FM2 = 1, etc %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% mich.filterfile = [FOTON_dir 'L1LSC.txt']; mich.digitalfilters = [0,1,2,3,5]; %%%%%%%%%%% digital gains %%%%%%%%%%%%%%%%%%%%%%%%%% mich.mich2rm = -10.5; % Output matrix: MICH to RM mich.mich2bs = 7.4; % Output matrix: MICH to BS mich.itmtrx = -0.666; % Input matrix: POB_Q to MICH mich.lsc_gain = -0.08; % MICH filter module gain mich.DCcal = (sqrt(2) * mich.BScal * mich.mich2bs +... mich.RMcal * mich.mich2rm) -... mich.RMcal * mich.mich2rm; % MICH_CTRL cal %%%%%%%%%%% snubber component values %%%%%%%%%%%%%%% %% snubber is a series RC (R_snub & C_snub), in %%%% %% parallel with the coil %%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% mich.R_ser = 3000; % series resistor mich.R_snub = 680; % snubber series resistor mich.C_snub = 0.022e-6; % snubber series cap mich.C_cabl = 800e-12; % cable capacitance mich.R_coil = 22; % coil resistance mich.H_coil = 3.3e-3; % coil inductance % Misc Info (mostly unused) -------------------------------------------------------- mich.opgain = 2; % Calculated watts/meter mich.efficiency = 0.9 * 80e-6; % EO Shutter, clipping, etc. mich.eta = 0.8; % Amps / Volt mich.Zrf = 380; % Tank circuit impedance mich.Grf = 10; % Pre-amp gain mich.Cable = 0.67; % 3.5 dB of loss in 100' of RG-405 mich.PS = 0.707; % 3 dB loss in power splitter mich.Mixer = 0.5; % 6 dB loss in mixer mich.WG = 10^(36/20); % 24 dB of whitening gain mich.AA = 2; % Gain of 2 in single-diff conv mich.ADC = 32768/10; % 16-bit Analog to Digital conversion mich.POB_Q_GAIN = 0.125; % Compensates some whitening gain mich.electronics_gain = mich.eta *... % Counts / Watt mich.Zrf *... mich.Grf *... mich.Cable *... mich.PS *... mich.Mixer *... mich.WG *... mich.AA *... mich.ADC *... mich.POB_Q_GAIN; %------------------------------------------------------------------------------- %******************************************************************************* % Parameter file for PRC loop model % %%%%%%%% frequencies of interest %%%%%%%%%%%%%%%% prc.fl = 9; % lower frequency of band, Hz prc.fu = 1000; % upper frequency of band, Hz prc.npt = 301; % number of points in band prc.fs = 16384; % sampling frequency, Hz prc.ugf = 30; % unity gain frequency of loop, Hz %%%%%%%%% parameters of the plant %%%%%%%%%%%%%%% prc.cavpole = 100e3; % PRC pole?? prc.pendf0 = 0.75; % pendulum eigenfrequency, Hz prc.pendQ = 10; % pendulum Q prc.RMcal = 0.38e-9; % DC calibration of RM, m/count prc.tdelay = 170e-6; % time delay in loop, sec prc.rclength = 9.204; % (l1+l2)/2 prc.hflowpass = ... zpk([],-2*pi*[100e3],... (2*pi*100e3)); % RC lowpass in L1 after % mixer %%%%%%%%%%% digital filters %%%%%%%%%%%%%%%%%%%%%%%% %% specify vector of PRC digital filters engaged %%% %% Note: FM1 = 0, FM2 = 1, etc %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% prc.filterfile = [FOTON_dir 'L1LSC.txt']; prc.digitalfilters = [0,1,2,3,4,5]; %%%%%%%%%%% digital gains %%%%%%%%%%%%%%%%%%%%%%%%%% prc.prc2rm = 7.4; % Output matrix: PRC to RM prc.itmtrx = 0.17; % Input matrix: POB_I to PRC prc.lsc_gain = -0.125; % PRC filter module gain prc.DCcal = -prc.RMcal * prc.prc2rm; % PRC_CTRL cal %%%%%%%%%%% snubber component values %%%%%%%%%%%%%%% %% snubber is a series RC (R_snub & C_snub), in %%%% %% parallel with the coil %%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% prc.R_ser = 3000; % series resistor prc.R_snub = 680; % snubber series resistor prc.C_snub = 0.022e-6; % snubber series cap prc.C_cabl = 800e-12; % cable capacitance prc.R_coil = 22; % coil resistance prc.H_coil = 3.3e-3; % coil inductance % Misc Info (mostly unused) -------------------------------------------------------- prc.opgain = 1.2e9; % Calculated watts/meter prc.efficiency = mich.efficiency; % PO fraction * 90/10 split prc.eta = mich.eta; % Amps / Volt prc.Zrf = mich.Zrf; % Tank circuit impedance prc.Grf = mich.Grf; % Pre-amp gain prc.Cable = mich.Cable; % 3.5 dB of loss in 100' of RG-405 prc.PS = mich.PS; % 3 dB loss in power splitter prc.Mixer = mich.Mixer; % 6 dB loss in mixer prc.WG = 10^(36/20); % 36 dB of whitening gain prc.AA = mich.AA; % Gain of 2 in single-diff conv prc.ADC = mich.ADC; % 16-bit Analog to Digital conversion prc.POB_I_GAIN = 0.125; % Compensates some whitening gain prc.electronics_gain = prc.eta *... % Counts / Watt prc.Zrf *... prc.Grf *... prc.Cable *... prc.PS *... prc.Mixer *... prc.WG *... prc.AA *... prc.ADC *... prc.POB_I_GAIN; %********************************************************************************** % Parameter file for CARM loop model % %%%%%%%% frequencies of interest %%%%%%%%%%%%%%%% carm.fl = 9; % lower frequency of band, Hz carm.fu = 8000; % upper frequency of band, Hz carm.npt = 301; % number of points in band carm.fs = 16384; % sampling frequency, Hz carm.ugf = 150; % unity gain frequency of loop, Hz %%%%%%%%% parameters of the plant %%%%%%%%%%%%%%% carm.cavpole = 1; % cavity pole, Hz carm.pendf0 = 0.75; % pendulum eigenfrequency, Hz carm.pendQ = 10; % pendulum Q carm.ETMXcal = 0.38e-9; % DC calibration of ETMX, m/count carm.ETMYcal = 0.38e-9; % DC calibration of ETMY, m/count carm.tdelay = 100e-6; % time delay in loop, sec carm.armlength = 3995.15; % arm length in meters carm.hflowpass = ... zpk([],-2*pi*[33e3 33e3 33e3],... (2*pi*33e3)^3); % RC lowpass in L1 after mixer %%%%%%%%%%% digital filters %%%%%%%%%%%%%%%%%%%%%%%% %% specify vector of CARM digital filters engaged %%% %% Note: FM1 = 0, FM2 = 1, etc %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% carm.filterfile = [FOTON_dir 'L1LSC.txt']; carm.digitalfilters = [0,1,2]; %%%%%%%%%%% digital gains %%%%%%%%%%%%%%%%%%%%%%%%%% carm.carm2etmx = -2.5; % Output matrix: CARM to ETMX carm.carm2etmy = -2.5; % Output matrix: CARM to ETMY carm.itmtrx = 0.05; % Input matrix: REFL_I to CARM carm.lsc_gain = -1.5; % CARM filter module gain carm.DCcal = carm.ETMXcal * carm.carm2etmx +... carm.ETMYcal * carm.carm2etmy; %%%%%%%%%%% snubber component values %%%%%%%%%%%%%%% %% snubber is a series RC (R_snub & C_snub), in %%%% %% parallel with the coil %%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% carm.R_ser = 3000; % series resistor carm.R_snub = 680; % snubber series resistor carm.C_snub = 0.022e-6; % snubber series cap carm.C_cabl = 800e-12; % cable capacitance carm.R_coil = 22; % coil resistance carm.H_coil = 3.3e-3; % coil inductance % Misc Info (mostly unused) -------------------------------------------------------- carm.opgain = 1.2e9; % Calculated watts/meter carm.efficiency = 0.125; % PO fraction * 90/10 split carm.eta = 0.8; % Amps / Volt carm.Zrf = 320; % Tank circuit impedance carm.Grf = 10; % Pre-amp gain carm.Cable = 0.67; % 3.5 dB of loss in 100' of RG-405 carm.PS = 0.707; % 3 dB loss in power splitter carm.Mixer = 0.5; % 6 dB loss in mixer carm.WG = 10^(30/20); % 36 dB of whitening gain carm.AA = 2; % Gain of 2 in single-diff conv carm.ADC = 32768/10; % 16-bit Analog to Digital conversion carm.REFL_I_GAIN = -1.0; % Compensates some whitening gain carm.electronics_gain = carm.eta *... % Counts / Watt carm.Zrf *... carm.Grf *... carm.Cable *... carm.PS *... carm.Mixer *... carm.WG *... carm.AA *... carm.ADC *... carm.REFL_I_GAIN; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ % - - - Miscellaneous Parameters - - - %_____________________________________________________________________________ los.m = 10.5; % LOS mass in kg los.phi = 1e-3; % Loss angle of steel wire los.wirelength = 0.442; % LOS wire length los.fp = 0.75; % LOS pend freq %%%%%%%%%%% AA & AI filtering %%%%%%%%%%%%%%% [z,p,k] = ellip(4,4,60,2*pi*7570,'s'); misc.ai = zpk(z,p,k*10^(4/20)) * zpk([],-2*pi*13e3,2*pi*13e3); % Fudged Anti-Imaging Filter [z,p,k] = ellip(8,0.001,80,2*pi*7570,'s'); misc.aa = zpk(z,p,k*10^(0.001/20)) *... zpk([],-2*pi*32768,2*pi*32768); % Schematic AA Filter %[z,p,k] = ellip(8,.035,80,2*pi*7570,'s'); %misc.aa = zpk(z,p,k*10^(0.035/20)); misc.c = 299792458; misc.ec = 1.6022e-19; misc.lambda = 1064e-9; misc.nu = misc.c / misc.lambda; misc.alpha = 1/137.0360; %____________________________________________________________________________ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L1.darm = darm; L1.mich = mich; L1.prc = prc; L1.carm = carm; L1.los = los; L1.misc = misc; return