Reports until 14:39, Wednesday 25 September 2019
H1 CAL
ling.sun@LIGO.ORG - posted 14:39, Wednesday 25 September 2019 - last comment - 14:36, Thursday 03 October 2019(52132)
Process all valid LF+HF meas for 0909 model and update uncertainty estimate

Lilli S, Jeff K,

All valid low-fre and high-freq meas for the new 0909 model are processed together (including 0821, 0828, 0904x2, 0909x2, 0916 LF meas, and 0903 HF meas).

The sensing GPR fitting results are attached in the pdf files.

The script used to generate the GPR fitting is:
^trunk/Runs/O3/H1/Scripts/Uncertainty/process_allmeas_includingPCALXHFdata_writeGPRHDF5_model20190909-C.py

The resulting hdf5 file is:
^trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_GPR_20190909_nofsQcorr_fmin20Hz_lengthscaleZp5_withPCALXHFdata.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_20190909_multi.hdf5 --HDF5_C_MCMCresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_MCMC_20190904_nomconfig_fmin20Hz_finaltrial.hdf5 --HDF5_C_GPRresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_GPR_20190909_nofsQcorr_fmin20Hz_lengthscaleZp5_withPCALXHFdata.hdf5 --modelPath=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/params/ --modelFilename=modelparams_H1_20190909 --IFOmodel=modelPars --outDir=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/ --plot1SigmaUncs --sampleNumber=1000 --seed=1234 --version=ref --gpsTime=1251687618

The online uncertainty code will start to use this new hdf5 for 0909 model.

Images attached to this report
Non-image files attached to this report
2019-09-16_H1_allO3BMeas_nofsQcorr_fmin20Hz_lengthscaleZp5_withPCALXHFdata_sensingFunction_referenceModel_vs_allMeasurements.pdf
H1_allO3BMeas_nofsQcorr_fmin20Hz_lengthscaleZp5_withPCALXHFdata_model2019-09-09_meas2019-09-16_sensingFunction_GPR.pdf
Comments related to this report
ling.sun@LIGO.ORG - 10:03, Thursday 26 September 2019 (52146)
Link

The last LF measurement before the break taken on 0925 was added to the GPR fitting. The fitting results are shown in the attached pdfs. The uncertainty budget slightly improves compared to the results posted yesterday. Here only one set of HF data is included (only the 0903 set is valid for the new model 0909). The HF part can be further improved when more HF data sets are added.

Images attached to this comment
Non-image files attached to this comment
jeffrey.kissel@LIGO.ORG - 10:28, Friday 27 September 2019 (52160)
Link 

I've re-run the RRNom.py command that Lilli indicates above, but have added the "--saveSummaries flag such that text files representing the curves above are saved:
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_20190909_multi.hdf5 --HDF5_C_MCMCresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_MCMC_20190904_nomconfig_fmin20Hz_finaltrial.hdf5 --HDF5_C_GPRresults=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_C_GPR_20190909_nofsQcorr_fmin20Hz_lengthscaleZp5_withPCALXHFdata.hdf5 --modelPath=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/params/ --modelFilename=modelparams_H1_20190909 --IFOmodel=modelPars --outDir=/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/ --plot1SigmaUncs --sampleNumber=1000 --seed=1234 --version=ref --gpsTime=1251687618 --saveSummaries

The text files are now committed to the repo here:
/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Results/Uncertainty/
    2019-09-26_O3_LHO_GPSTime_1251687618_ref_RelativeResponseUncertainty_FinalResults.txt
    2019-09-26_O3_LHO_GPSTime_1251687618_ref_RelativeResponseUncertainty_MinMax.txt
ling.sun@LIGO.ORG - 14:59, Friday 27 September 2019 (52164)
Link

Now I've added the last set of A meas (0925) to the hdf5 file: ^trunk/Runs/O3/H1/Results/Uncertainty/O3_H1_A_GPR_20190909_multi.hdf5. The resulting uncertainty budget is attached.

The GPR script is: ^trunk/Runs/O3/H1/Scripts/Uncertainty/process_allmeas_writeGPRHDF5_model20190909-A.py

The RRNom command is the same as above.

The txt summary files are ^trunk/Runs/O3/H1/Results/Uncertainty/2019-09-27_O3_LHO_GPSTime_1251687618_ref_RelativeResponseUncertainty_FinalResults.txt and 2019-09-27_O3_LHO_GPSTime_1251687618_ref_RelativeResponseUncertainty_MinMax.txt

Images attached to this comment
jeffrey.kissel@LIGO.ORG - 14:36, Thursday 03 October 2019 (52286)
Link 

To be consistent with the file name convention established in Chunks 1 and 2, I've changed the name of the script that produces the GPR estimate and committed the file name change:
/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/O3/H1/Scripts/Uncertainty/
    process_allmeas_writeGPRHDF5_meas20190828-20191001_withHF_model20190909-C.py