I have updated the correction filter for the dtt DARM spectra that are projected in the wall and screens in the control room. I hope Evan will like it.
[The DARM spectra]
Here is a comparison of the newly corrected and previous DARM spectra.
As shown in the attached screenshot, the high frequency excess wing that has been shown above 3 kHz is now gone. The noise level at around 7 kHz now seems a bit too low compared with the f-shaped GWINC curve (shown in green). This maybe a real calibration error because I have applied the best and latest knowledge about the sensing chain for correcting the spectrum.
The correction filter is made by a matlab script:
/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/PreER8/H1/Scripts/ControlRoomCalib/H1DARM_FOM_correction.m
This script creates a text file which contains the correction filter in [freq, dBmag, phasedeg] format from 1 Hz to 7444 Hz. The text file can be found at
/ligo/svncommon/CalSVN/aligocalibration/trunk/Runs/PreER8/H1/Scripts/ControlRoomCalib/DARM_FOM_calibration.dat
Both script and text files are checked in to the svn. Also I updated the latest FOM dtt template, H1_DARM_FOM_20150722.xml
so that it now uses the new correction filter.
[Corrections newly included]
Here is a list of what I newly included.
Note that the previous correction was made only from the IOP down sampling filter and the numerical whitening zpk([1,1,1,1,1], [100,100,100,100,100]) which are still included in the new correction.
For those of us who like to look at the DCPD streams, I modified the script to also produce the appropriate correction TF (using the uncompensated preamp poles, the AA filter, and the downsampling filters only).
There is still a little bit of sloping at high frequencies, even in the null stream (which we expect to be flat above 100 Hz). Maybe the AA model needs to be tuned?
I happened to come across this code again to reproduce the calibration correction filter for O1. I then found some bugs in the code which are now fixed.