Merge branch 'master' into 'master'

Adding Hong et al. calculations for Coating Brownian Noise

See merge request gwinc/pygwinc!52

gwinc/ifo/aLIGO/ifo.yaml

@@ -188,7 +188,7 @@ Suspension:
 
 ## Optic Material -------------------------------------------------------
 Materials:
-  MassRadius: 0.17                # m; 
+  MassRadius: 0.17                # m;
   MassThickness: 0.200            # m; Peter F 8/11/2005
 
   ## Dielectric coating material parameters----------------------------------
@@ -203,6 +203,14 @@ Materials:
     Indexhighn: 2.09              # 2020 LMA
     Phihighn: 3.89e-4             # loss angle at 100Hz (Gras 2020)
     Phihighn_slope: 0.1
+    # Following are optional parameters to for providing separate Bulk and
+    # Shear loss angles and to include photoelastic effect in Coatings
+    # Brownian noise calculation.
+    # lossBhighn: 3.89e-4         # Bulk loss angle at 100 Hz
+    # lossBhighn_slope: 0.1       # Bulk loss angle slope
+    # lossShighn: 3.89e-4         # Shear loss angle at 100 Hz
+    # lossShighn_slope: 0.1       # Shear loss angle slope
+    # PEThighn: 0.1095            # Hong et al . PRD 87, 082001 (2013) A.1
 
     ## low index material: silica
     Ylown: 70e9                   # 2020 LMA
@@ -217,6 +225,18 @@ Materials:
                                   # slopes between 0 and 0.3, depending on
                                   # deposition method. Slawek's analysis in
                                   # 10.1103/PhysRevD.98.122001 assumes zero slope.
+    # Following are optional parameters to for providing separate Bulk and
+    # Shear loss angles and to include photoelastic effect in Coatings
+    # Brownian noise calculation.
+    # lossBlown:  2.3e-5         # Bulk loss angle at 100 Hz
+    # lossBlown_slope: 0         # Bulk loss angle slope
+    # lossSlown: 2.3e-5          # Shear loss angle at 100 Hz
+    # lossSlown_slope: 0         # Shear loss angle slope
+    # PETlown: 0.269             # Hong et al . PRD 87, 082001 (2013) A.1
+
+    # Indicator for including effect of amplitude noise due to coatings
+    # brownian noise. Set to 'yes' for including. Otherwise would be ignored.
+    # IncCoatBrAmpNoise: 'yes'     # Hong et al . PRD 87, 082001 (2013) II.E
 
   ## Substrate Material parameters--------------------------------------------
   Substrate:

gwinc/ifo/noises.py

@@ -296,10 +296,21 @@ class CoatingBrownian(nb.Noise):
         w0, wBeam_ITM, wBeam_ETM = arm_cavity(self.ifo)
         dOpt_ITM = coating_thickness(self.ifo, 'ITM')
         dOpt_ETM = coating_thickness(self.ifo, 'ETM')
+        mTi_ITM = None
+        mTi_ETM = None
+        Ic = None
+        if 'IncCoatBrAmpNoise' in materials.Coating:
+            if materials.Coating.IncCoatBrAmpNoise.lower() == 'yes':
+                mTi_ITM = self.ifo.Optics.ITM.Transmittance
+                mTi_ETM = self.ifo.Optics.ETM.Transmittance
+                Ic = ifo_power(self.ifo)[1]   # Circulating Arm Power
+                precomp_mirror(self.freq, self.ifo)
         nITM = noise.coatingthermal.coating_brownian(
-            self.freq, materials, wavelength, wBeam_ITM, dOpt_ITM)
+            self.freq, materials, wavelength, wBeam_ITM,
+            dOpt_ITM, Ic, mTi_ITM)
         nETM = noise.coatingthermal.coating_brownian(
-            self.freq, materials, wavelength, wBeam_ETM, dOpt_ETM)
+            self.freq, materials, wavelength, wBeam_ETM,
+            dOpt_ETM, Ic, mTi_ETM)
         return (nITM + nETM) * 2