Lilli S, Jeff K,

Lately we have tracked down some bugs and misunderstandings -- see more details in 48272, 48276. Now we have made a new model file /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/params/modelparams_H1_20190404.py, using the measurements taken on 2019-04-03. This new model has been validated and installed in CAL-CS.

--------------------------------------------------

The measurements used to create the model live in /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Measurements
The processing scripts live in /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Scripts
See attached pdfs for the fitting results. See the detailed fitting output below.

To validate the model, we processed the broadband + sweeps measurements taken last night -- /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Measurements (2019-04-10)

The results are attached as BB_deltaL/pcal and SS_deltaL/pcal files. With the response correction factors and a 40-microsecond delay applied, the results (red dots) look good. (Scripts are/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Scripts/CALCS_FE/process_broadband_pcal_20190410.py and /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Scripts/FullIFOSensingTFs/process_sweep_pcal2darm_20190410.py)

--------------------------------------------------

We also did an estimate of the uncertainty for the reference model (i.e. all kappas set to 1).

The GPR fitting using multiple measurements are shown in multi-meas-GPR.pdf (4 meas for sensing and 2 meas for actuation); Script -- /ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Scripts/Uncertainty/process_allmeas_writeGPRHDF5_20190404.py

The uncertainty plot for the reference model is attached.

Fitting details are as below:

---------------------------------------------------------------------
Sensing
---------------------------------------------------------------------
Inverse Sensing FOTON values: [NB: SRCD2N zpk gain based on sensing sign in parameters file]
SRCD2N: zpk([410.5556;0.0428+i*4.4689;0.0428-i*4.4689],[0.1;0.1;7000],1,"n")gain(1997.42)
Gain: gain(3.077e-07)

Inverse Sensing without cavity pole FOTON values for CFTD path: [NB: SRCD2N zpk gain based on sensing sign in parameters file]
SRCD2N: zpk([0.0428+i*4.4689;0.0428-i*4.4689],[0.1;0.1],1,"n")gain(1997.42)
Gain: gain(3.077e-07)

Parameter                                | Quantiles (0.15, 0.50, 0.84)
---------------------------------------------------------------------
Optical gain, H_c (ct/m)                 | 3.249e+06, 3.25e+06, 3.252e+06
Cavity pole, f_cc (Hz)                   | 409.8, 410.6, 411.3
Detuned SRC spring frequency, f_s (Hz)   | 4.425, 4.469, 4.513
Detuned SRC spring quality factor, 1/Q_s | 64.37, 52.19, 43.88
Residual time delay, tau_c (usec)        | 0.5522, 0.9963, 1.438
------------------------------- OR ----------------------------------
Optical gain, H_c (ct/m)                 | 3.25e+06 (+1430,-1408) or (+0.04399%,-0.04333%)
Optical gain, H_c (mA/pm)                | 4.34 (+0.001909,-0.00188) or (+0.04399%,-0.04333%)
Cavity pole, f_cc (Hz)                   | 410.6 (+0.7268,-0.7248) or (+0.177%,-0.1765%)
Detuned SRC spring frequency, f_s (Hz)   | 4.469 (+0.04354,-0.04401) or (+0.9743%,-0.9847%)
Detuned SRC spring quality factor, Q_s   | 52.19 (+275.5,-275.9) or (+18.94%,-18.92%)
Residual time delay, tau_c (usec)        | 0.9963 (+0.4416,-0.4441) or (+44.32%,-44.58%)

---------------------------------------------------------------------
UIM
---------------------------------------------------------------------
The Force Coefficient "answer" in terms of ...
    Force per count of longitudinal DAC request (just down stream of DRIVEALIGN bank):
    Gain = 7.67e-08 (N/ct)
    Force per actuator signal (current or voltage):
    Gain = 1.634 (N/A)
 
All MCMC Results (w/ Uncertainty) for UIM Parameters:
    Parameter                              | Quantiles (0.15, 0.50, 0.84)
    ---------------------------------------------------------------------
    Actuator gain, H_c (N/ct)              | 7.666e-08, 7.67e-08, 7.673e-08
    Residual time delay, tau_A (usec)      | 58.27, 61.05, 63.85
    ------------------------------- OR ----------------------------------
    Actuator gain, H_c (N/ct)              | 7.67e-08 (+3.684e-11,-3.709e-11) or (+0.04803%,-0.04835%)
    Residual time delay, tau_A (usec)      | 61.05 (+2.794,-2.786) or (+4.576%,-4.563%)

---------------------------------------------------------------------
PUM
---------------------------------------------------------------------
The Force Coefficient "answer" in terms of ...
    Force per count of longitudinal DAC request (just down stream of DRIVEALIGN bank):
    Gain = 6.036e-10 (N/ct)
    Force per actuator signal (current or voltage):
    Gain = 0.02947 (N/A)
 
All MCMC Results (w/ Uncertainty) for PUM Parameters:
    Parameter                              | Quantiles (0.15, 0.50, 0.84)
    ---------------------------------------------------------------------
    Actuator gain, H_c (N/ct)              | 6.033e-10, 6.036e-10, 6.039e-10
    Residual time delay, tau_A (usec)      | 38.6, 40.88, 43.17
    ------------------------------- OR ----------------------------------
    Actuator gain, H_c (N/ct)              | 6.036e-10 (+3.143e-13,-3.131e-13) or (+0.05207%,-0.05187%)
    Residual time delay, tau_A (usec)      | 40.88 (+2.296,-2.276) or (+5.617%,-5.569%)

---------------------------------------------------------------------
TST
---------------------------------------------------------------------
The Force Coefficient "answer" in terms of ...
    Force per count of longitudinal DAC request (just down stream of DRIVEALIGN bank):
    Gain = 4.727e-12 (N/ct)
    Force per actuator signal (current or voltage):
    Gain = 4.427e-11 (N/V**2)
 
All MCMC Results (w/ Uncertainty) for TST Parameters:
    Parameter                              | Quantiles (0.15, 0.50, 0.84)
    ---------------------------------------------------------------------
    Actuator gain, H_c (N/ct)              | 4.724e-12, 4.727e-12, 4.729e-12
    Residual time delay, tau_A (usec)      | 5.352, 5.947, 6.538
    ------------------------------- OR ----------------------------------
    Actuator gain, H_c (N/ct)              | 4.727e-12 (+2.411e-15,-2.393e-15) or (+0.051%,-0.05062%)
    Residual time delay, tau_A (usec)      | 5.947 (+0.5917,-0.5942) or (+9.951%,-9.993%)