def get_quad_bias_act(ff):
zeros = []
poles = [1.8, 1.8, 1590]
gain = 40 * (2 * esd_coeff * esd_bias)
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('quadTstTstLL.txt', ff)
return ff, abs(driver_tf * sus_resp) * ligo_dac_noise(ff)
def get_quad_tst_act(ff):
zeros = [50, 50, 3250]
poles = [2.2, 2.2, 152]
gain = 2 * (2 * esd_coeff * esd_bias)
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('quadTstTstLL.txt', ff)
return ff, abs(driver_tf * sus_resp) * ligo_dac_noise(ff)
def get_BS_Mar_act(ff):
zeros = []
poles = [300]
# 28mN/A per pair, 30 ohms (resistor+coil), noise adds incoherently
gain = sqrt(2) * 28e-3 / 30
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('BS_Mar_resp.txt', ff)
return ff, abs(driver_tf * sus_resp) * virgo_dac_noise / 200
def get_TM_Mir_act(ff, res=10):
zeros = []
poles = [300]
# 3.4mN/A per pair, 20 ohms (resistor+coil), noise adds incoherently
gain = sqrt(2) * 3.4e-3 / (res + 20)
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('TM_Mir_resp.txt', ff)
return ff, abs(driver_tf * sus_resp) * virgo_dac_noise
def get_bs_pum_act(ff):
# T1100479
zeros = [10, 20, 10]
poles = [1, 200, 100]
gain = 0.32e-3 * 0.963
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('bsfmM2M3LL.txt', ff)
return ff, abs(driver_tf * sus_resp) * ligo_dac_noise(ff) / 280
def get_quad_pum_act(ff):
# T1100378
zeros = [6, 20, 1.35]
poles = [0.5, 250, 110]
gain = 0.268e-3 * 0.0309
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('quadPumTstLL.txt', ff)
return ff, 2 * abs(driver_tf * sus_resp) * ligo_dac_noise(ff)
def get_BS_Mir_act(ff):
zeros = []
poles = [300]
# 1.7mN/A per pair, 30 ohms (resistor+coil), noise adds incoherently
#gain = sqrt(2)*1.7e-3 / 30
# factor of 6 increase for mirror from combination of magnets, coils, DAC actuation changes
gain = sqrt(2) * 1.7e-3 * 6 / 30
_, driver_tf = nbu.freqresp((zeros, poles, gain), ff)
sus_resp = nbu.get_complex_interp('BS_Mir_resp.txt', ff)
return ff, abs(driver_tf * sus_resp) * virgo_dac_noise / 200