Lilli S, Jeff K
We have processed all valid measurements for the observing period from 0611 (powered up to 37W) to 0828 (fixed SRC deturning), including high freq data. New GPR hdf5 files are generated, and the uncertainty budget is reestimated.
The actuation and sensing GPR fitting results are attached in the pdf files:
A_GPR_allmeas0327-0904_model0416.pdf
C_GPR_allmeas0328-0828_model0416.pdf
The scripts used to generate the GPR fittings are:
^trunk/Runs/O3/H1/Scripts/Uncertainty/process_allmeas_writeGPRHDF5_model20190416-A.py
^trunk/Runs/O3/H1/Scripts/Uncertainty/process_allmeas_writeGPRHDF5_meas20190328-20190828_withHF_model20190416-C.py
The resulting hdf5 files are:
^trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_GPR_20190416_meas0327-0904.hdf5
^trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_GPR_20190416_meas0328-0828_withHFafter0611_nofsQcorr_fmin30Hz_lengthscaleZp5.hdf5
The uncertainty plot was generated using this command line:
python3 RRNom.py --IFO=LHO --HDF5_A_MCMCresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_MCMC_20190404.hdf5 --HDF5_A_GPRresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_GPR_20190416_meas0327-0904.hdf5 --HDF5_C_MCMCresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_MCMC_20190404.hdf5 --HDF5_C_GPRresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_GPR_20190416_meas0328-0828_withHFafter0611_nofsQcorr_fmin30Hz_lengthscaleZp5.hdf5 --modelPath=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/params/ --modelFilename=modelparams_H1_20190416 --IFOmodel=modelPars --outDir=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/ --plot1SigmaUncs --sampleNumber=1000 --seed=1234 --version=ref --gpsTime=1247687618
NOTE: Please see the 3rd pdf file showing the GPR fitting with all available HF meas (three additional sets 0515, 0531, and 0610 compared to the results above). The 0610 (lasting from 0610 to 0618) data set (green dots with large systematics in the upper right corner) is obviously problematic because it includes the power up on 0611. Hence we decided to throw it away. We did not include 0515 and 0531 data sets either since they were taken before the power up and did not match the IFO configuration during 0611-0828 (in fact including them or not does not impact the final results).
For chunk 2 (0611-0828), I've also added the 0925 A meas to A GPR since actuation does not change over the full period. The fitting results are in the pdf file.
The script used to generate the GPR fitting is updated:
^trunk/Runs/O3/H1/Scripts/Uncertainty/process_allmeas_writeGPRHDF5_model20190416-A.py
The resulting hdf5 file is:
^trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_GPR_20190416_meas0327-0925.hdf5
The uncertainty plot was generated using this command line:
python3 RRNom.py --IFO=LHO --HDF5_A_MCMCresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_MCMC_20190404.hdf5 --HDF5_A_GPRresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_GPR_20190416_meas0327-0925.hdf5 --HDF5_C_MCMCresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_MCMC_20190404.hdf5 --HDF5_C_GPRresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_GPR_20190416_meas0328-0828_withHFafter0611_nofsQcorr_fmin30Hz_lengthscaleZp5.hdf5 --modelPath=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/params/ --modelFilename=modelparams_H1_20190416 --IFOmodel=modelPars --outDir=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/ --plot1SigmaUncs --sampleNumber=1000 --seed=1234 --version=ref --gpsTime=1247687618