TITLE: 08/15 Eve Shift: 2330-0500 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Lock Acquisition
INCOMING OPERATOR: Oli
SHIFT SUMMARY: Quiet shift with H1 observing for almost the whole duration except for a lockloss right at the end. Initial alignment just finished, now starting main locking sequence.
LOG:
Lockloss @ 04:45 UTC after 18+ hours locked - link to lockloss tool
Perhaps caused by the smallest of ETMX glitches starting in L3 right before the lockloss? Doesn't exactly look to me like a usual one of these, though. Otherwise, no obvious cause.
Jumping right into an alignment given the long lock stretch.
TITLE: 08/14 Day Shift: 1430-2330 UTC (0730-1630 PST), all times posted in UTC
STATE of H1: Observing at 151Mpc
INCOMING OPERATOR: Ryan S
SHIFT SUMMARY: Locked for 13 hours. We had calibration and comissioning this morning, where we rode through a few earthquakes with help from some new configuration (alog86364). Ryan S noticed ITMX mode 2 has been growing slowly for this entire lock. It was never damped because it was almost non existant at the start of the lock. He will address it when practical.
LOG:
| Start Time | System | Name | Location | Lazer_Haz | Task | Time End |
|---|---|---|---|---|---|---|
| 18:49 | LASER | LVEA is LASER HAZARD | LVEA | YES | LVEA IS LASER HAZARD | 09:49 |
| 19:22 | SPI | Jeff | Opt Lab | yes | Inventory | 19:25 |
TITLE: 08/14 Eve Shift: 2330-0500 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Observing at 153Mpc
OUTGOING OPERATOR: TJ
CURRENT ENVIRONMENT:
SEI_ENV state: CALM
Wind: 25mph Gusts, 16mph 3min avg
Primary useism: 0.02 μm/s
Secondary useism: 0.12 μm/s
QUICK SUMMARY: H1 has been locked for 12.5 hours and observing since commissioning wrapped up this morning.
I noticed the ~500Hz violin mode peak is high in DARM and saw that indeed ITMX mode 2 is elevated and has been growing this entire lock stretch. It looks like the gain was never turned on for this mode, which happens if the monitor for a mode is low enough when the full power damping is turned on. Since the violin damping Guardian sees the mode is growing, it's constantly zeroing the gain in the damping filters, so we've been unable to turn on the nominal gain of -15. To fix this, we would have to take H1 out of observing, take the VIOLIN_DAMPING node to 'DAMPING_ON_SIMPLE', then back to 'DAMP_VILINS_FULL_POWER' so that the node resets its "node growing" checks and sees the mode needs damping. I plan to do this if either L1 drops out of observing or if the 500Hz peak on DARM gets up to 10-16 m/Hz-1/2
I ran a bruco yesterday that showed lots of low-ish broadband coherence of SRCL, PRCL and MICH.
I have a previous good measurement of the PRCL coupling, so I used that to iteratively fit an improved feedforward. Based on the last test I did of the PRCL feedforward, I am not certain that just adjusting the gain is good enough right now to reduce the coupling significantly. screenshot
I also measured and tried to fit a better SRCL feedforward. The challenge continues to be fitting the low and high frequency behaviors of the SRCL coupling simultaneously. The misfit of the high frequency behavior continues to create a broad bump of SRCL noise around 100-300 Hz.
I made some small improvements of the SRCL coupling, but there is still room for more improvement. screenshot
I still need to remeasure and refit the MICH feedforward.
The new SRCL FF is in FM6 and the PRCL FF is in FM8. I updated the guardian and the SDF.
Follow up bruco taken after these changes shows large reduction in PRCL coherence, some reduction of SRCL coherence. Minimal change in MICH coherence.
Edgard, Ivey
This Tuesday Oli ran into issues installing the new blend filters and new SR3 Yaw transfer function fits for the OSEM estimator (see LHO: 86319). An error from the cplxpair function, which pairs complex conjugates, appeared only when converting the fit from zpk into state space. Somehow, this would introduce a small but significant difference between any two complex pairs, which would cause the computer to reject it as a complex pair. To resolve this issue, we ran a function to average each complex pair. We tested it with our mock PR3 model in the lab, and this seemed to fix the error cplxpair was encountering.
We were not able to reproduce the error that Oli noticed with the blend filters in our lab, but we will look more into this, and hopefully the issue will be fixed by next Tuesday when Oli tests the estimator again with the new fits.
The new cleaned SR3 yaw fits were added to SusSVN and can be found here '~/SusSVN/sus/trunk/HLTS/Common/FilterDesign/Estimator/fits_H1SR3_Y_2025-08-05.mat'. More on the fits LHO: 86233.
The new blend filter with updated fits can be found here '~/SusSVN/sus/trunk/HLTS/Common/FilterDesign/Estimator/Estimator_blend_skinnynotch_SR3yaw_20250723.mat'. More on the blend filter LHO: 86265.
DQ shifters: Riley McNeil and Emil Lofquist-Fabris
Daily observing duty cycle: 57.05%, 67.36%. 99.68%, 59.91%, 91:23%, 69.16%, 100.00%. Week average: 77.77%.
The BNS range this week sat around 150 Mpc, however there were a few sudden drops in the range by about 5-10 Mpc. These were seen on Tuesday, Thursday, and Friday.
There were several issues with the squeezer this week that resulted in sudden increases in noise and glitching, as well as a couple drops from observing.
There was a recurring spike in the y-manifold beam tube motion [X] at 17:00 every day. This was seen in the previous week as well, this is caused by the daily dewar fill noise at the Y-manifold cryopump. (alog)
There was also recurring noise from ~11:00-14:00 in several corner station accelerometers also seen in previous shifts.
Just above 20 Hz (~22-25 Hz) there has been a line of glitching that we believe to be light scattering. This has been seen since at least the start of July 2025.
There has also been a lot of glitching/noise from ~20-40 Hz, consistent with previous weeks, likely from an AC unit housed inside the VPW. (alog)
There was a 34 hour lock segment this week, from Tuesday August 5th 22:19 to August 7th 08:31.
The detector also spent all of Sunday observing as well.
See full report here: https://wiki.ligo.org/DetChar/DataQuality/DQShiftLHO20250804
This morning a 6.2 magnitude earthquake near Vanuatu happened during the commissioning window, so we tested the asc gain reversion scripts Elenna came up with in this alog. It took a little longer than I would have liked to get the scripts working, so we didn't actually transition until basically the peak ground velocities, but the transition went pretty smoothly. I think this eq is tied for the biggest we stayed locked through for all of O4.
Attached image shows ground velocities on the top row (both peakmon and the ITMY Z STS 30-100mhz blrms, second row is the ISC_LOCK state, third and fourth show channels related to the asc transition. Fourth row is the first gain changed by the script CHARD Y, the third row is the SRCL FF gain, which is the last thing changed in the script. The dashed white vertical lines on the CHARD trace and the top row ground traces show the time the transition took to complete.
Even if the transition was a little late I think that this test likely still saved the lock, I have often seen on eqs like this that the IFO stays locked through the peak ground motion, but loses lock some time shortly after, while the eq is still ringing down.
There is still a lot of work to be done as the transition is currently handled by a couple of scripts on the ISI_CONFIG overview medm and the transition will knock us out of Observe, still unclear what the best way to automate this would be, but still a pretty good test.
It's been remarkable and wonderful to see how well we've been surviving earthquakes lately!
This is such a big win, that I suggest we not worry too much about if this pops us out of Observing. If we're able to automate saving the lock, and then going back to Observe when the ground motion is compatible with our usual loops, that's already a huge leap forward on improving our duty cycle.
I have updated ISC LOCK guardian and LSCPARAMS to have the gains for the lownoise ASC and length feedforwards as parameters in LSCPARAMS. Then, TJ helped me figure out how to import LSCPARAMS into these two scripts so it draws those gains from the parameters instead of having them hardcoded. This way, if we decide to update a feedforward or ASC gain, it will also be updated for the earthquake script.
This script also changes ASC filters, but an update to have those filters as parameters in LSCPARAMS is a much more involved change to ISC_LOCK, so I have to think about that one a bit. That means if we change ASC filters, we have to remember to update this script, for now.
Thu Aug 14 10:08:28 2025 INFO: Fill completed in 8min 24secs
/ligo/home/camilla.compton/Documents/sqz/templates/dtt/20250814higher_order_modes.xml screenshot attached.| Type | Time (UTC) | Angle | DTT Ref | Notes |
| SQZ | 15:54:00 - 15:59:00 | (-)133 | ref 0 | While taking CAL measurements |
| FDS Mid - SQZ | 16:01:30 - 16:04:00 | (-)105 | ref 1 | |
| FDS Mid SQZ, -50urad pit +3deg | 16:07:30 - 16:09:30 | (-)109 | ref 2 | |
| FDS Mid - SQZ | 16:10:30- 16:12:30 | (-)111 | ref 3 | Redid at 4dB ASQZ |
| FDS Mid SQZ, -50urad yaw | 16:14:30 - 16:16:30 | (-)114 | ref 4 | |
| FDS Mid SQZ, +50urad yaw | 16:17:45 - 16:19:45 | (-)109 | ref 5 | |
| FDS Mid SQZ, +100urad yaw | 16:21:30 - 16:23:30 | (-)106 | ref 6 | |
| FDS Mean SQZ, +100urad yaw | 16:26:00 - 16:28:00 | N/A | ref 7 | |
| FDS Mean SQZ, | 16:29:00 - 16:31:00 | N/A | ref 8 | |
|
No SQZ
|
16:37:00 - 16:57:00
|
N/A
|
ref 9
|
took 1500 avg (~10mins)
|
| OPO Setpoint | Amplified Max | Amplified Min | UnAmp | Dark | NLG | Note |
| 80 | 0.1121155 | 0.001849 | 0.0068192 | -1.16e-5 | 16.4 | Without Optimizing Temp |
| 90 | 0.257827 | 0.001813 | 37.7 |
| Type | Time (UTC) | Angle | DTT Ref | SQZ |
|
No SQZ
|
16:37:00 - 16:57:00
|
N/A
|
ref 0
|
|
| ASQZ NLG 37 | 17:43:30-17:46:30 | (-)81 | ref 1, 10 | 19dB |
| SQZ NLG 37 | 17:51:00-17:54:00 | (-)141 | ref 2, 11 | -4.5dB |
Took more 5kHz/10kHz data with the NLG at 37 (effects of HOM easier to see). DTT saved as /ligo/home/camilla.compton/Documents/sqz/templates/dtt/20250814higher_order_modes.xml
| Type | Time (UTC) | Angle | DTT Ref | Notes |
|
No SQZ
|
16:37:00 - 16:57:00
|
N/A
|
ref 0
|
|
| ASQZ NLG 37 | 17:43:30-17:46:30 | (-)81 | ref 1, 10 | |
| SQZ NLG 37 | 17:51:00-17:54:00 | (-)141 | ref 2, 11 | |
| Mid SQZ | 18:00:30 - 18:02:30 | (-)118 | ref 12 | |
| Mid SQZ +4cts DHARD PIT | 18:06:30 - 18:08:30 | (-)118 | ref 13 | Difference only at 10k? Repeat |
| Mid SQZ | 18:09:30 - 18:11:30 | (-)118 | ref 14 | |
| Mid SQZ +4cts DHARD PIT | 18:12:30 - 18:14:30 | (-)118 | ref 15 | |
| Mid SQZ -4cts DHARD PIT | 18:15:45- 18:17:45 | (-)118 | ref 16 | |
| Mid SQZ +4cts DHARD YAW | 18:19:30 - 18:21:30 | (-)118 | ref 17 | No difference seen |
| Mid SQZ SRM YAW 1urad (offset 0.3) | 18:28:00 - 18:30:00 | (-)120 | ref 18 | Big difference at 5kHz, none at 10kHz |
| Mid SQZ SRM YAW -1urad (offset -0.3) | 18:33:30 - 18:35:50 | (-)114 | ref 19 | Got better |
| Mid SQZ SRM YAW -2urad (offset -0.6) | 18:41:30 - 18:43:30 | (-)108 | ref 20 | Overshot, got worse, mode didn't flip |
| Mid SQZ SRM YAW -1.5urad (offset -0.45) | 18:46:30 - 18:48:30 | (-)108 | ref 21 | |
| Mid SQZ | 18:52:00 - 18:52:00 | (-)118 | ref 22 | Back to nominal |
| OPO Setpoint | Amplified Max | Amplified Min | UnAmp | Dark | NLG | Note |
| 80 | 0.153078 | 0.002067 | 0.0068192 | -1.16e-5 | 22.4 | With Optimizing Temp |
Running another calibration measurement today following the usual broadband and simulines.
Simulines Start:
PDT: 2025-08-14 08:36:33.095985 PDT
UTC: 2025-08-14 15:36:33.095985 UTC
GPS: 1439221011.095985
Simulines End:
PDT: 2025-08-14 08:59:53.303431 PDT
UTC: 2025-08-14 15:59:53.303431 UTC
GPS: 1439222411.303431
Files:
2025-08-14 15:59:53,136 | INFO | File written out to: /ligo/groups/cal/H1/measurements/DARMOLG_SS/DAR
MOLG_SS_20250814T153633Z.hdf5
2025-08-14 15:59:53,145 | INFO | File written out to: /ligo/groups/cal/H1/measurements/PCALY2DARM_SS/
PCALY2DARM_SS_20250814T153633Z.hdf5
2025-08-14 15:59:53,150 | INFO | File written out to: /ligo/groups/cal/H1/measurements/SUSETMX_L1_SS/
SUSETMX_L1_SS_20250814T153633Z.hdf5
2025-08-14 15:59:53,156 | INFO | File written out to: /ligo/groups/cal/H1/measurements/SUSETMX_L2_SS/
SUSETMX_L2_SS_20250814T153633Z.hdf5
2025-08-14 15:59:53,161 | INFO | File written out to: /ligo/groups/cal/H1/measurements/SUSETMX_L3_SS/
SUSETMX_L3_SS_20250814T153633Z.hdf5
TITLE: 08/14 Day Shift: 1430-2330 UTC (0730-1630 PST), all times posted in UTC
STATE of H1: Observing at 153Mpc
OUTGOING OPERATOR: Tony
CURRENT ENVIRONMENT:
SEI_ENV state: CALM
Wind: 2mph Gusts, 0mph 3min avg
Primary useism: 0.01 μm/s
Secondary useism: 0.10 μm/s
QUICK SUMMARY: Observing for 4 hours, automated relock. Environment is calm. The only active alarm is the known vacuum PT242B pressure. Planned calibration and commissioning time today 1530-1930 UTC.
TITLE: 08/14 Eve Shift: 2330-0500 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Observing at 153Mpc
INCOMING OPERATOR: Tony
SHIFT SUMMARY: Very quiet shift with H1 observing throughout, despite a couple hours of high winds blowing through with gusts briefly over 40mph. There's a M5.6 quake inbound from the South Pacific, but the EQ response tool just places it in the lower bound of the "earthquake" region. H1 has now been locked for 14.5 hours.
J. Kissel Recall the SPI pathfinder is pathfinding more than just the future fully assembled sensor array's roll in improving seismic isolation, but also -- at a low-level -- pathfinding some new individual types optomechanical components for LIGO. This aLOG covers the latest in the path for the Schaefter + Kirchhoff (SuK) fiber collimators. The first of SPI's Schaefter + Kirchhoff (SuK) fiber collimators (DCC:S0272502, ICS:S0272502) is now Class B clean after having gone through one round of massaging with isopropyl alcohol and air-baking (see CNB:2187), at the planned time and levels for the future Class-A vacuum-bake (peak temp 85 [deg C]; ramp up 6 hours, hold for 48 hours, ramp down for 6 hours). The first open question was "will the lens remain intact and undamaged through the bake." the definitive answer: YES -- the lens remained whole and, at least visually, without defect. While I didn't use a optical fiber microscope, the lens appears very much intact with no obvious cracks or glinting. See the attached Visual Inspection. Next up was to check it's optical performance as (hopefully) a more quantitative measure of performance change (hopefully not peformance degradation). Prior to the bake, I'd set the lens position on the collimator to achieve some level of collimation, and characterized the beam profile -- see LHO:84825. After the bake, I used the same characterization setup -- augmented only to keep the now Class-B fiber collimator clean (see LHO:86299) -- to remeasure the beam profile to see if the collimator still projected the same beam quantitatively. The beam remains as I had collimated it, with the waist within 100 [cm] of the prebake position. I did *not* adjust the tuned "before" pre-bake lens position at all prior to taking the "after" data. Check out the 2025-08-07_spifc_S0272502_beamprofile_fit.pdf attachment. - First page compares the two data sets, X axis (parallel to the optical table surface plane) on the top panel, and Y-axis (perpendicular to the optical table plane) on the bottom panel. You can see that the change in waist position / size is Using Delta = (2025-08-06)' - (2025-06-03), and % Diff = [(2025-08-06)' - (2025-06-03)]/(2025-06-03) z0x' - z0x = +0.0544 [m] (+4% difference) z0y' - z0y = -0.1103 [m] (-8% difference) w0x' - w0x = -12.14 [um] (1.3% difference) w0y' - w0y = -10.17 [um] (1.1% difference) Excellent. Because the X waist position moved in +Z, and the Y waist position moved in -Z, the beam has become a bit more astigmatic. This is evident in the "Far Field" projection of the model on pages 3 and 5. But, given that the measurement setup is limited with the furthest data point being z(meas_max) - z0 ~ 5.41 [m] - 1.4 [m] ~ 4.0 [m] away from the waist, I wouldn't claim that the measurements perfectly constrain the far-field behavior. Recall we *want* the waist to be at z=0, at the fiber collimator and this collimator's lens position was *not* tuned to that simply because of user error. I'm fully confident I *can* set the lens position such that the waist *is* at the fiber collimator, and after doing so the 5.41 [m] NanoScan position will get a bit more "in to" the far field, which will hopefully better constrain the model, and thus get a better handle on how astigmatic the beam gets after baking (or even *if* it gets consistently astigmatic). If we define the dimensionless astigmatism parameter, A, as the (zRx - zRy) difference in Rayleigh range, divided by the (zRx + zRy) sum (with the Rayleigh range defined by the fit waist, zR = pi * w0^2 / lambda), then the change (% Difference) in astigmatism is only 2025-06-03 A = +2.3983 2025-08-06 A' = +2.3409 (A' - A)/A = 0.02494 = 2.5% which seems totally tolerable from pre- to post-bake. Regardless, for the remaining to collimators that we've yet to bake, we'll be setting the collimation (lens position) post-bake anyways, so how it changes across a bake is moot. Whether the absolute value of A ~ +2.35 is tolerable an open question, that we'll work to answer in the mean time. Pages 2 and 4 show the model zoomed in to the data between z = 0 and 6 [m]. Here you can also see that the fit doesn't *perfectly* match the data their either, so there's another systematic grain of salt to take with the assessment of change. Minor note: I was not consistent with the orientation of the collimator w.r.t. to the optical table surface; I oriented the collimator 90-deg from the 2025-06-03 vs 2025-08-06 data (because I found out / rediscovered that the 2025-08-06 position is how you align the p-pol transmission with the optical table; see Figure 5 of T2400413). I've flip-flopped the X and Y axes data for the waist size in the 2025-06-03 data to account for this (which is why the careful reader would notice a difference between this entry's version of the 2025-06-03 results from that in LHO:86342). But -- all in all -- this looks good enough! I've taken these results to indicate that we can move forward with the full Class-A clean-and-bake. All three collimators, (S0272502, S0272503, S0272504) in the queue now -- see CNB:2243.
In consultation with Sheila, I've added a new flag in lscparams.py called "ignore_sqz" (currently set to False) which when True, removes SQZ_MANAGER from ISC_LOCK's managed nodes list and sets up ISC_LOCK such that no requests are made of SQZ_MANAGER. I've additionally removed all instances of ISC_LOCK requesting SQZ_MANAGER to 'NO_SQUEEZING.'
The motivation for this stems from an issue encountered earlier this week (alluded to in Monday's shift summary) where after H1 dropped observing from the SQZ PMC having to relock, SQZ_MANAGER eventually stalled (which is not unexpected) and ISC_LOCK then requested it to 'NO_SQUEEZING' for a then-not-well-understood reason. After looking into the frequently used "unstall nodes" decorator in ISC_LOCK, I learned that the revive function it calls simply requests the stalled subordinate node to whatever its last requested state was. The catch here is that the last requested state refers only to what the manager's last request was to the subordinate, not whatever request a user, other node, or standalone script may have made, as the manager node has no way of knowing about requests outside of its own. A discrepancy between a manager node's last request to a subordinate and a different request to that subordinate can be seen with a notification on the manager saying the subordinate's request changed.
My understanding of the sequence of events that led to the confusion on Monday is that commissioners had been working with the squeezer and wanted its Guardian manager node to remain in 'NO_SQUEEZING' while H1 was relocking following their work. In its design at the time, in a few different states such as 'INJECT_SQUEEZING', ISC_LOCK would check the status of SQZ_MANAGER, and if it was in 'NO_SQUEEZING', move on with the state and re-request SQZ_MANAGER to 'NO_SQUEEZING'. This makes SQZ_MANAGER's manager's last request 'NO_SQUEEZING', so when SQZ_MANAGER later stalled, it was requested back to 'NO_SQUEEZING' even though a user had set SQZ_MANAGER to its nominal state sometime later. This is why I've removed requests to 'NO_SQUEEZING' in ISC_LOCK; I believe it's a sound assumption that if SQZ_MANAGER is in 'NO_SQUEEZING', someone wants it there and Guardian should just move on. Further, this makes it so that the only request ISC_LOCK ever makes to SQZ_MANAGER is its nominal state (except in 'DOWN', of course, which would then later get updated to the nominal), meaning there should be no confusion as to what its last request was.
Changes have been saved and committed to svn, but ISC_LOCK has not yet been loaded. This should be done at the next drop from observing.
We loaded this and tested going to and from ISC_LOCK's Inject_Squeezing state with the SQZ manager in No_Squeezing and FREQ_DEP_SQZ. All worked well.
If we ever need to go to observing without squeezing, we should keep this in mind as I'm not 100% confident we won't get some manager confusion depending on when we do the transition. We'll cross that bridge when we get to it though.
Today, we measured the calibration at three different ESD biases. First, we measured at the current bias of 269 V, and then our O4 standard bias of 136 V. Then, I stepped up to a higher bias of 409 V.
| ESDAMON value | Bias Offset | L3 Drivealign gain | Calibration report | Notes |
| 269 V | 6.0542453 | 88.28457 |
alog: 86337 report: 20250813T153848Z |
only 1 hour thermalized at measurement time current operating bias |
| 136 V | 3.25 | 198.6643 |
alog: 86339 report: 20250813T162026Z |
nominal O4 bias, calibration model fit at this bias |
| 409 V | 8.89 | 57.587 |
alog: 86341 report: 20250813T174921Z |
ESD saturation warnings while at this bias Took 5 minutes of quiet time, cal lines on, at this new bias start: 18:12:54 UTC end: 18:18:00 UTC |
The attached plot compares the three broadband measurements at each ESD bias. It seems like the overall systematic error decreases as we increase the ESD bias.
To step up the ESD bias, I used guardian code that Sheila attached to this alog. Another relevant alog comparing simulines results at different biases is here.
Attaching figures comparing the sensing function, the actuation function, and the open loop gain (olg). All the figures are formatted in the same way, where the left side shows the bode plot from each report and the right shows the bode plot from of each measurement ratio to a reference. I used the latest exported calibration measurement "20250719T225835Z" as a reference. From 10 Hz to 1 kHz the sensing and the actuation function residuals are within 5%. The OLG is within 10% with one outlier at 410 Hz.
We are trying to understand how the systematic error is changing at each bias voltage, even though we think we are correcting the drivealign gain to account for the actuation change.
Francisco and I made some plots of how the modeled error changes. First, we pulled the model from the 20250719T225835Z report, since that is our current calibration model. Then, we pulled the kappas at the time of the lowest bias voltage measurement, since that is the bias voltage that our model is based on. We applied the kappas from that measurement time, and then calculated a new response function assuming an additional TST actuation change, ranging from no change (0%) to 1.5% change. Then, we compared each of these response functions to the kappa-corrected model.
To be clear, we are calculating the new response function as:
R = 1/C_model + (error_factor*TST_model + PUM_model + UIM_model) * D_model
The "model" in this case also has the kappa-corrected values applied, which are:
{'c': 0.98335475,
'f_c': 447.65558,
'uim': 1.0052187,
'pum': 1.0012773,
'tst': 1.0183398}
Looking at our results side-by-side with the broadband pcal measurement, we see some similarites. However, it's not exactly the same, since the frequency dependence appears slightly different in the measurement than the model.
There are some other comparisons to be made, but we can start with these. The script I used to make this plot is saved in /ligo/groups/cal/H1/ifo/scripts/compare_models_tst_err.py
Leo, Jennie, Camilla, WP 12694
Jennie followed instructions for set up in 85775. We removed the SQZ beam iris at the bottom of the LPM (added for alignment capture during OFI vent work). Then took beam profiles in this SQZ path of the SEED beam with various PSAMS settings, adjusted PSAMS settings as in 85775 and used the servo for nominal settings. All data is attached. Jennie then reverted settings back to nominal.
Leo, Jennie W., Camilla The attached pdf contains all the beam q parameters fitted to the collected beam width data. Only the 13.5% data was fitted, as the D4S data was too inconsistent to obtain confident q values. Fitting was performed with the a la mode beamPath.fitBeamWidth function. The attached q parameters were individually plotted using a la mode and verified for their data-fitting accuracy. As mentioned in the document, all q parameters are located immediately after the interaction with ZM5 (through the view of BM4 -> ZM4 -> ZM5 beam travel).
Leo, Jennie W., Camilla Attached is a plot of the q manifold from the q parameter data, which allows for characterizing the beam smoothly with respect to ZM4/5 strain gauge voltage values. The image is taken from the presentation uploaded to T2500228. The real plot will likely have slightly different labels to axes. Link to git code for plotting: https://git.ligo.org/leendert.schrader/alm-beam-simulation-for-sqz/-/tree/main