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.
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.
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
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
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