Suspension thermal noise sub-budgets

Adds suspension thermal noise sub-budgets. This addresses 1/3 of #51 (closed). Suspension thermal noise can now be broken up into horizontal and vertical noise of each of the suspension stages. With one exception, which does not change the total noise, the calculation is the same as the current calculation but the contributions of each stage and direction are tracked separately.

In the current calculation the horizontal and vertical spring constants of each stage are broken up into an upper and lower joint and are stored separately. The total suspension thermal noise is then calculated by adding the contributions of each of these joints using the appropriate temperature of each joint. Instead of adding the contributions of each joint all at once, this MR makes it possible to track each stage with separate Noise classes and add them all to a sub-budget.

The noise for each stage is the sum of the noise from that stage's lower joint and the noise from the upper joint of the stage below. In the case of the top mass, the contribution from the upper joint is also included. In the current calculation, for the upper stages, the imaginary part of the horizontal spring constant is equally distributed between the upper and lower joints while it is only attributed to the upper joint for the test mass. The one change this MR makes to the calculation itself is to distribute this equally between the top and bottom joints as is done for the upper stages.

The imaginary part of the vertical spring constant is only attributed to the upper joint of each stage. The blade springs are the only contribution to the vertical spring constant for the upper stages while the lower stage also includes the continuum fibers or ribbons. It's unclear how to break this up into an upper and lower stage as is done in figure 14 of https://arxiv.org/abs/2001.11173, so this MR does not have vertical noise of the test mass. (The noise is included, it's just attributed entirely to the PUM.)

The sub-budgets are only implemented for the Cosmic Explorer interferometers since these are the only ones that have any other sub-budgets. The CE1 suspension thermal noise budget is shown below as an example.

gwinc/ifo/noises.py

@@ -124,57 +124,8 @@ def ifo_power(ifo, PRfixed=True):
     power.Tpr = Tpr
     return power
 
-
-############################################################
-# precomp functions
-############################################################
-
-def precomp_suspension(f, ifo):
-    pc = Struct()
-    pc.VHCoupling = Struct()
-    if 'VHCoupling' in ifo.Suspension:
-        pc.VHCoupling.theta = ifo.Suspension.VHCoupling.theta
-    else:
-        pc.VHCoupling.theta = ifo.Infrastructure.Length / const.R_earth
-    hForce, vForce, hTable, vTable = suspension.suspQuad(f, ifo.Suspension)
-    pc.hForce = hForce
-    pc.vForce = vForce
-    pc.hTable = hTable
-    pc.vTable = vTable
-    return pc
-
-
 ##################################################
 
-############################################################
-# calibration
-############################################################
-
-def dhdl(f, armlen):
-    """Strain to length conversion for noise power spetra
-
-    This takes into account the GW wavelength and is only important
-    when this is comparable to the detector arm length.
-
-    From R. Schilling, CQG 14 (1997) 1513-1519, equation 5,
-    with n = 1, nu = 0.05, ignoring overall phase and cos(nu)^2.
-    A small value of nu is used instead of zero to avoid infinities.
-
-    Returns the square of the dh/dL function, and the same divided by
-    the arm length squared.
-
-    """
-    c = const.c
-    nu_small = 15*pi/180
-    omega_arm = pi * f * armlen / c
-    omega_arm_f = (1 - sin(nu_small)) * pi * f * armlen / c
-    omega_arm_b = (1 + sin(nu_small)) * pi * f * armlen / c
-    sinc_sqr = 4 / abs(sin(omega_arm_f) * exp(-1j * omega_arm) / omega_arm_f
-                       + sin(omega_arm_b) * exp(1j * omega_arm) / omega_arm_b)**2
-    dhdl_sqr = sinc_sqr / armlen**2
-    return dhdl_sqr, sinc_sqr
-
-
 
 def precomp_suspension(f, ifo):
     pc = Struct()
@@ -200,7+151,7 @@
     pc.ASvac = noise_dict['ASvac']
     pc.SEC = noise_dict['SEC']
     pc.Arm = noise_dict['arm']
     pc.Injection = noise_dict['injection']
     pc.PD = noise_dict['pd']
 
     # FC0 are the noises from the filter cavity losses and FC0_unsqzd_back