Add CO2 and Ar gas noise functions to noises.py

gwinc/ifo/noises.py

@@ -833,6 +833,42 @@ class ExcessGasScatteringH2O(nb.Noise):
         return n * 2 / sinc_sqr
 
 
+class ExcessGasScatteringCO2(nb.Noise):
+    """Excess gas scattering for CO2
+
+    """
+    style = dict(
+        label='CO2 scattering',
+        color='xkcd:slate blue'
+    )
+
+    def calc(self):
+        cavity = arm_cavity(self.ifo)
+        species = self.ifo.Infrastructure.ResidualGas.CO2
+        n = noise.residualgas.residual_gas_scattering(
+            self.freq, self.ifo, cavity, species)
+        dhdl_sqr, sinc_sqr = dhdl(self.freq, self.ifo.Infrastructure.Length)
+        return n * 2 / sinc_sqr
+
+
+class ExcessGasScatteringAr(nb.Noise):
+    """Excess gas scattering for Ar
+
+    """
+    style = dict(
+        label='Ar scattering',
+        color='xkcd:violet'
+    )
+
+    def calc(self):
+        cavity = arm_cavity(self.ifo)
+        species = self.ifo.Infrastructure.ResidualGas.Ar
+        n = noise.residualgas.residual_gas_scattering(
+            self.freq, self.ifo, cavity, species)
+        dhdl_sqr, sinc_sqr = dhdl(self.freq, self.ifo.Infrastructure.Length)
+        return n * 2 / sinc_sqr
+
+
 class ExcessGasDampingH2(nb.Noise):
     """Residual gas damping for H2
 
@@ -865,9 +901,6 @@ class ExcessGasDampingN2(nb.Noise):
     @nb.precomp(sustf=precomp_suspension)
     def calc(self, sustf):
         species = self.ifo.Infrastructure.ResidualGas.N2
-        # n = noise.residualgas.residual_gas_damping(
-        #     self.freq, self.ifo, species, sustf)
-        # return 4 * n
         return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
 
 class ExcessGasDampingH2O(nb.Noise):
@@ -883,12 +916,40 @@ class ExcessGasDampingH2O(nb.Noise):
     @nb.precomp(sustf=precomp_suspension)
     def calc(self, sustf):
         species = self.ifo.Infrastructure.ResidualGas.H2O
-        # n = noise.residualgas.residual_gas_damping(
-        #     self.freq, self.ifo, species, sustf)
-        # return 4 * n
         return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
 
 
+class ExcessGasDampingCO2(nb.Noise):
+    """Excess gas damping for CO2
+
+    """
+    style = dict(
+        label='CO2 damping',
+        color='xkcd:slate blue',
+        linestyle='--',
+    )
+
+    @nb.precomp(sustf=precomp_suspension)
+    def calc(self, sustf):
+        species = self.ifo.Infrastructure.ResidualGas.CO2
+        return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
+
+
+class ExcessGasDampingAr(nb.Noise):
+    """Excess gas damping for Ar
+
+    """
+    style = dict(
+        label='Ar damping',
+        color='xkcd:violet',
+        linestyle='--',
+    )
+
+    @nb.precomp(sustf=precomp_suspension)
+    def calc(self, sustf):
+        species = self.ifo.Infrastructure.ResidualGas.Ar
+        return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
+
 def _calc_residual_gas_damping(f, ifo, species, sustf):
     """Redundant calculation of res gas damping for each species