LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY - LIGO - CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Document Type LIGO-T070100-00- W 2007/05/17 Validation of the S5 V3 LHO calibrations for GRB 070201 Evan Goetz, Eiichi Hirose, Michael Landry Distribution of this draft: LIGO Scientific Collaboration California Institute of Technology Massachusetts Institute of Technology LIGO Project - MC 18-34 LIGO Project - NW17-161 1200 E California Blvd 175 Albany Street Pasadena CA 91125 Cambridge, MA 01239 Phone (626) 395-2129 Phone (617) 253-4824 Fax (626) 304-9834 Fax (617) 253-7014 E-mail: info@ligo.caltech.edu E-mail: info@ligo.mit.edu WWW: http://www.ligo.caltech.edu/ Processed with L A TEX on 2007/05/17
Abstract We describe validation of the S5 V3 calibration of the LIGO Hanford interferometers during the epoch including GRB070201. As of this writing, the current S5 calibration version is V3. Final V3 epochs for both H1 and H2 were expected to still be representative of the H1 and H2 instruments on the date of GRB070201 (gps time 854378604). Through a series of cross checks of the state of the machines, and comparisons with the quasi-independent calibration procedure autocalibrator, this is shown to be the case. 1 Introduction Both H1 and H2 were locked and in science mode during GRB070201. We have plotted strainequivalent noise curves for H1 and H2, using roughly 55s of DARM ERR data centered on gps time 854378604 (GRB070201). These curves correspond to binary neutron star inspiral ranges of 15.6Mpc (H1) and 6.8Mpc (H2), respectively. The calibration employed here for both H1 and H2 is S5 V3. Strain equivalent noise (Hz 1/2 ) 10 19 10 20 10 21 10 22 GRB070201 Noise curves H2 H1 10 23 Figure 1: Strain equivalent noise curves for H1 and H2 for 55s around the gps event time of GRB070201. 2 H2 Calibration Validation of the H2 calibration for GRB070201 is slightly more complicated than H1, owing to a filter modification prior to the Feb 1 event. We thus deal with H2 first in this document, and then follow with H1. The S5 H2 V3 calibration comprises three distinct calibration epochs, listed in Table 1. The duration of the third epoch is given as nine 9 s, or indefinite. This indicates the calibration is expected then to apply to GRB070201. However, there was an H2 filter modification made on 22 January 2007, necessitating a check on this epoch to ensure its validity at the time of the GRB. page 2 of 7
epoch start (gps) 815155213 824949188 846138794 epoch duration (s) 9793975 21189606 999999999 OLG measurement (gps) 816068773 826542380 849678155 2.1 V3 and V4 model comparisons Table 1: S5 H2 V3 epoch parameters. An official calibration run was performed for H2 on 23 January 2007 during the commissioning break. We made a V4 model from the results of this calibration run and compare it to the V3 model to understand any changes in the DARM loop that were made during the commissioning break period. Figure 2 shows the comparison of V3 and V4 versions of the DARM model. The differences between V3 and V4 versions of the H2 calibration are less than 1 percent in magnitude and less than 1 degree in phase between 10 Hz and 7 khz, a good indication that the filter module change has not significantly impacted the calibration and the V3 calibration is still valid at the time of the Feb 1 GRB. Figure 2: Comparison of the V3 and V4 versions of the DARM model for H2. 2.2 Autocalibrator comparisons As an additional cross-check, we compare autocalibrations (typically made weekly to bi-weekly) against propogated V3 calibration, in order to look for systematic shifts that suggest the V3 calibration is out of date. Autocalibrations rely on the same DC calibrations that are employed by page 3 of 7
the official method, but use a physical transfer function (as opposed to a model), and thus all other components of the calibration are independent. Two autocalibrations bracket the Feb 1 GRB: that of Jan 29, 2007, and Feb 6, 2007 (this for both H2 and H1). Shown in Figures 3 and 5 below are strain curves from the two autocalibration measurements, and the strains produced from the V3 calibration for the same AS Q data use by the autocalibrator. Agreement is more readily judged in the linear ratios Figures 4 and 6. Both ratios show a slight (about 4-5) percent systematic to more sensitive autocalibrations. As this i) is within calibration error, ii) is conservative with respect to the V3 calibration, iii) may demonstrate some problem with the autocalibration sweep templates, we believe the systematic to be small enough such that the offical calibration is acceptable. We have begun testing the autocalibrator procedure in effort to identify the nature of the systematic, and will shortly prove whether or not it resides in the autocalibrator as suspected. 10 17 10 18 H2 Official and Acal strains official V3 cal autocal new f r es srd 10 19 Strain (1/ Hz) 10 20 10 21 10 22 10 23 Figure 3: Comparison of H2 29 Jan 07 autocalibrator spectrum with official calibration method. 2.3 Factors check 2.4 Conlog check We checked the Conlog differences for LSC, ASC, SUS, PSL and IO subsystems between the V3 model time for H2 (gps=849678155) and the time of GRB 070201 (gps=854378604). The state of the instrument shows nothing out of the ordinary at the time of the GRB (i.e. there were no gain changes in the DARM loop or filters not in place). The Conlog differences are mostly due to alignment differences and servo offsets (SUS and ASC). A few example of channels that had changes are H2:SUS RM LLYAW GAIN and H2:LSC-AS3 I OFFSET. The input matrix element (AS Q to DARM ERR) is slightly different between the fiducial calibration time and the time of the GRB, but this is a dynamical parameter in the LSC loop to cancel changes in optical gain over time. page 4 of 7
Ratio of H2 Off. cal/autocal, w0=.764 propogated V3 cal/reformed autocal 1.15 1.1 1.05 1 0.95 0.9 0.85 0.8 Figure 4: Ratio of H2 29 Jan 07 official calibration to autocalibrator output. 3 H1 Calibration There were no changes in filters or gains between the V3 calibration measurement and the time of GRB070201. A new model is not implicitly necessary because the state of the instrument is essentially the same as during the V3 official calibration measurements. Thus the validation is considered simpler and did not necessitate a separate V4 measurement specifically for GRB070201. The four H1 V3 calibration epochs are noted in Table 1 below. epoch start (gps) 815155213 824695694 835044014 843942254 epoch duration (s) 9540481 10348320 8898240 999999999 OLG measurement (gps) 815844976 824791240 835664459 849677446 Table 2: S5 H1 V3 epoch parameters. 3.1 Autocalibrator comparisons Autocalibrator comparisons for Jan 29 and Feb 6 were made for H1, in a similar fashion to that described above for H2. While no significant mean systematic is observed, one as a function of frequency (about ±5%) is indeed present. This may or may not reside within the 4km autocalibrator setup - investigations are also underway. The error is however within the error of the two calibrations, and thus again we claim the V3 calibration sufficient for current analyses of the GRB070201 event. page 5 of 7
10 17 10 18 H2 Official and Acal strains official V3 cal autocal new f r es srd 10 19 Strain (1/ Hz) 10 20 10 21 10 22 10 23 Figure 5: Comparison of H2 6 Feb 07 autocalibrator spectrum with official calibration method. 3.2 Factors check 3.3 Conlog check We checked the Conlog differences between the V3 model time for H1 (gps=849677446) and the time of GRB 070201 (gps=854378604). The state of the instrument shows nothing out of the ordinary at the time of the GRB (i.e. there were no gain changes or filters not in place). As for H2, nearly all the Conlog differences are due to alignment changes and offsets. Examples of channels with different settings are H1:ASC-QPDX 1 OFFSET and H1:SUS-ETMX LRYAW GAIN. Also, as with H2, the intput matrix element for AS Q to DARM ERR is slightly different due to changing optical gain. 4 Response functions For completeness, we plot here the V3 DARM ERR response functions, and the compare them to their V2 counterparts. Figure 11 is a plot of the H1 V3 response function (fourth and final epoch), while Figure 12 compares this function with the V2 response function from its (second and) final epoch. Figure 13 is a plot of the H2 V3 response function (third and final epoch), while Figure 14 compares this function with the V2 response function from its (second and) final epoch. 5 Conclusions The H1 and H2 V3 calibrations (final epochs) are sufficiently representative of the working state of the machines on Feb 1, 2007 to be employed in data analyses of GRB070201. Some sytematics were noted in comparisons with autocalibrator output: we will continue to investigate these effects page 6 of 7
Ratio of H2 Off. cal/autocal, w0=.764 propogated V3 cal/reformed autocal 1.15 1.1 1.05 1 0.95 0.9 0.85 0.8 Figure 6: Ratio of H2 6 Feb 07 official calibration to autocalibrator output. 10 17 10 18 H1 Official and Acal strains official V3 cal reformed autocal, w0=.767 srd 10 19 Strain (1/ Hz) 10 20 10 21 10 22 10 23 Figure 7: Comparison of H1 29 Jan 07 autocalibrator spectrum with official calibration method. The dip at 60Hz is an artifact of the autocalibrator measurement. and they are suspected to reside in the existing LHO autocal setup (such as swept sine template amplitudes which may be saturating the actuation chain). page 7 of 7
Ratio of H1 Off. cal/autocal, w0=.761 propogated V3 cal/reformed autocal 1.15 1.1 1.05 1 0.95 0.9 0.85 0.8 Figure 8: Ratio of H1 29 Jan 07 official calibration to autocalibrator output. 10 17 10 18 H1 Official and Acal strains official V3 cal reformed autocal, w0=.767 srd 10 19 Strain (1/ Hz) 10 20 10 21 10 22 10 23 Figure 9: Comparison of H1 6 Feb 07 autocalibrator spectrum with official calibration method. The dip at 60Hz is an artifact of the autocalibrator measurement. page 8 of 7
Ratio of H1 Off. cal/autocal, w0=.761 propogated V3 cal/reformed autocal 1.15 1.1 1.05 1 0.95 0.9 0.85 0.8 Figure 10: Ratio of H1 6 Feb 07 official calibration to autocalibrator output. Resp. mag. (m/ct) Resp. phase (deg) 10 12 10 14 10 16 200 100 0 100 200 H1 Reference response function Figure 11: V3 response function for H1 DARM ERR. page 9 of 7
Figure 12: Comparison of DARM ERR V3 reference response function to the V2 reference response function for H1. Resp. mag. (m/ct) Resp. phase (deg) 10 12 10 14 10 16 200 100 0 100 200 H2 Reference response function Figure 13: V3 response functions for H2 DARM ERR. page 10 of 7
Figure 14: Comparison of DARM ERR V3 reference response function to the V2 reference response function for H2. page 11 of 7