TITLE: 01/22 Day Shift: 1530-0030 UTC (0730-1630 PST), all times posted in UTC
STATE of H1: Observing at 154Mpc
OUTGOING OPERATOR: Ryan S
CURRENT ENVIRONMENT:
    SEI_ENV state: CALM
    Wind: 3mph Gusts, 2mph 3min avg
    Primary useism: 0.02 μm/s
    Secondary useism: 0.25 μm/s
QUICK SUMMARY: Locked for 12 hours, useism dropped quickly overnight. Planned commissioning time this morning.

TITLE: 01/22 Eve Shift: 0030-0600 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Observing at 154Mpc
INCOMING OPERATOR: Ryan S
SHIFT SUMMARY:
One lockloss and a successful return to Observing at 3:43 UTC.
Otherwise a quiet night.

LOG:
No Log

TITLE: 01/22 Eve Shift: 0030-0600 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Lock Acquisition
CURRENT ENVIRONMENT:
    SEI_ENV state: USEISM
    Wind: 4mph Gusts, 3mph 3min avg
    Primary useism: 0.03 μm/s
    Secondary useism: 0.64 μm/s
QUICK SUMMARY:

Unknown Lockloss from NLN.

 

The h1cam23 camera process got stuck.  I looked at it with Tony.  We shutdown and restarted the corresponding network switch port to power cycle the camera.  Then restarted the server process via monit.  That didn't work.  I ended up having to go onto h1digivideo2 and manually issue the kill command and send a SIGKILL to the camera process.  After that the camera process started up nicely and gave us a feed.

J. Freed

Before the hoidays as shown in 81658, the Double Mixer had large phase noises at 4096Hz that were thought to be caused by phase offest difference in the Double Mixer. Thanks to Marc Pirello for the suggestion and help to manually adjust the phase difference by extending the wire lengths inside the double mixer. 

DM_PN2.pdf Shows the lower 4096Hz peak in the path adjust(green) as compared to the no path adjust (blue) and both compared to an independant SRS (orange). The harmonics are still large so  DAC output had a ~5kHz low pass applied but with negledgable improvement. Investigation was ongoing when I went on holiday. 

Today I plugged in the double mixer into a network analyzer to see the power spectrum density of the signal 

DM_NA_WIDE.pdf Shows that although the DM has better noise floor, the harmonics of the double mixer are quite a bit worse in the Power spectral desnity than shown in phase noise DM_PN2. The best result (shown) was with NO added length. Unsure why this is. 

DM_NA_NARROW.pdf Besides that the signal has a fairly pure tone, though this is limmited by the resolution.

Investigations ongoing 

Following up on log  82252,  I made movies of HAM3 from the PR2 camera viewport last Thursday as I moved CPy over its pitch range and then, at nearly maximum pitch up, I ran it over its yaw range. In the movie (https://youtu.be/0f0F7cWqpwQ ), the beam spot appears to stay on the ballast baffle. At this highest pitch there appears to be a second spot on the support structure of the baffle. This second spot would be a little lower at the nominal setting. A steeply angled baffle that stuck down a little below the edge of the table might mitigate both of these stray beams.

TITLE: 01/22 Eve Shift: 0030-0600 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Lock Acquisition
OUTGOING OPERATOR: TJ
CURRENT ENVIRONMENT:
    SEI_ENV state: USEISM
    Wind: 5mph Gusts, 3mph 3min avg
    Primary useism: 0.03 μm/s
    Secondary useism: 0.62 μm/s
QUICK SUMMARY:
H1 is Locked and Observing.
Microseism is falling back down.

Noting random Issues :
MC2 PR2 Cameras still unplugged likely due to robert's work in the the LVEA. 
ITMY camera not working for some reason.

H1:PEM-CS_DUST_LAB2 Still out sick

TITLE: 01/22 Day Shift: 1530-0030 UTC (0730-1630 PST), all times posted in UTC
STATE of H1: Observing at 155Mpc
INCOMING OPERATOR: Tony
SHIFT SUMMARY: Maintenance day today, recovery was straight forward in terms of alignment, but we had a lock loss after DRMI ASC and then one at Transition_From_ETMX, which made recovery take a bit longer.
LOG:

• 16:08 Lock loss
• 1838 All CS SEI tripped from ISC-C4 power supply work. Those in the control room at the time didn't realize that there was a connection between this rack and all of the CS ISIs. This power supply supplies the 71MHz RF that goes to a fanout for all of the these ISI's CPSs.
• 2349 Observing
• 0000 Dropped observing for 8 seconds due to some SDF diff in calcs.

[Keita, Mayank]

We optimized the ALSY PDH demod phase of the CM board.

We measured the OLTF of the PDH loop using the setup detailed in https://alog.ligo-wa.caltech.edu/aLOG/index.php?callRep=82166 at 3142Hz. (Note that the raw number measured is 40dB smaller than the actual gain due to electronics gain in the generic DAQ filter in the denominator.)

We changed the delay line between the local oscillator and the demod board to maximize the OLTF amplitude. 
Before we changed anything, we measured the ALSY OLTF and it was about 9 dB smaller than ALSX at 3142Hz. As we changed the delay line we eventually had to decrease the IN1 LOCKED gain of the CM board (H1:ALS-Y_REFL_LOCK_LOCKEDGAIN) by 3dB so that ALSY OLTF becomes about the same as ALSX. ALSY locked and ran fine with that setting.

Initial values: Delay line Phase: 223 Degree (421 steps), CM board LOCK IN1 gain +10dB, TF magnitude measured -37 dB.
Tuned value: Delay line Phase: 165 Degree (311 steps), CM board LOCK IN1 gain +7dB, TF magnitude measured -27 dB. 

Attached is the screenshot of before (blue) and after (red) the tuning as well as ALSX (black dotted)

Here're the bare-bones instructions to run the DARM loop critique -- an automated set of plots that were inspired by the presentations given over the years that compare the transfer functions of the DARM loop between the two sites, e.g. G1501372, G1700316, G2000149.

(0) Conda activate the latest pydarm environment.
    ~$ conda activate /ligo/groups/cal/conda/pydarm
This is a symbolic link to the latest tagged version of pydarm.
    () ~$ which pydarm
        /ligo/groups/cal/conda/pydarm/bin/pydarm
    () ~$ pydarm --version
        20240821.0
    () ~$ ipython
    In [1]: import pydarm
    In [2]: pydarm.__version__
    Out [2]: '20240821.0' 
Note the version of pydarm is the same on the command line and within the ipython session *because* we activated the correct pydarm conda environment. The command line pydarm internally activates this environment automatically, but any old ipython session may not have the right version of pydarm if you just import "from scratch" without making sure you're in the right environment first. So -- any time you're running script or python session pyDARM these days -- make sure that you've activated the conda environment.

Is this the right version? check https://git.ligo.org/Calibration/pydarm/-/tags see if that version is the latest, i.e. the top of the list.

It is. And this will *always* be true, by process of deploying the latest tag (see deploy script).
To compare previous versions of tags, go to https://git.ligo.org/Calibration/pydarm/-/compare?from=master&to=master and change the drop-down menu to compare one tag to the other.
If you want to activate an environment with an older version of pydarm, then head to the folder /ligo/groups/cal/conda/ and activate the environment you like.


(1) Grab the DARM loop model parameter sets that you want to plot and/or compare. If at LHO, wanting H1 parameter files, then (back out on the command line, not in yet in a python session)
Find the list of reports and tags:
   () ~$ pydarm ls -r
From that list, find the model parameters installed in the front-end:
   () ~$ pydarm ls -r | exported

Once you find what you like, go into that report directory, e.g.
   () ~$ cd /ligo/groups/cal/H1/reports/20240927T211612Z
Look for pydarm_H1.ini, and copy it somewhere you like, or at least store where it lives.

If you're at LHO, want L1 parameter files, then the easiest thing is to download it from the web. Head to the https://ldas-jobs.ligo-la.caltech.edu/~cal/ and scroll down to the list of all reports, and find the entry that's been exported most recently e.g. 20240408T200652Z, and download the pydarm_L1.ini and save it to a known location.

(2) Start critiquing:
~$ ipython
In [1]: from pydarm.darm import DARMModel

In [2]: d = DARMModel('/ligo/groups/cal/H1/reports/20240927T211612Z/pydarm_H1.ini')

# Plot a single parameter file's model:
In [3]: d.plot(filename='/ligo/home/jeffrey.kissel/2025-01-21/20240927T211612Z_pydarm_H1_critique.pdf',label=['H1'])

# Compare against another, say, previous version of the model
d_prenewdarm = DARMModel('/ligo/groups/cal/H1/reports/20231027T203619Z/pydarm_H1.ini')

from pydarm.plot import critique

critique(d,d_prenewdarm,filename='/ligo/home/jeffrey.kissel/2025-01-21/20240927T211612Z_vs_20231027T203619Z_pydarm_H1_critique.pdf',label=['H1 20240927T211612Z', 'H1 20231027T203619Z'])

# Compare two IFOs:
In [9]: d_H1 = DARMModel('/ligo/groups/cal/H1/reports/20240927T211612Z/pydarm_H1.ini')

In [10]: d_L1 = DARMModel('/ligo/home/jeffrey.kissel/2025-01-21/20240408T200652Z_pydarm_L1.ini')

critique(d_H1,d_L1,filename='/ligo/home/jeffrey.kissel/2025-01-21/20240927T211612Z_H1_vs_20240408T200652Z_L1_critique.pdf',label=['H1 20240927T211612Z','L1 20240408T200652Z'])

# Compare one IFO's parameters, while updating one of the parameters on the fly
In [9]: d_H1 = DARMModel('/ligo/groups/cal/H1/reports/20240927T211612Z/pydarm_H1.ini')
In [10]: d_H1_moreAA = DARMModel('/ligo/groups/cal/H1/reports/20240927T211612Z/pydarm_H1.ini') #Intentionally the same
# Update the parameter you want,
In [6]: d_H1_moreAA.sensing.omc_filter_noncompensating_modules
Out[6]: [[9, 10], [9, 10]]

In [7]: d_H1_moreAA.sensing.omc_filter_noncompensating_modules = [[9,10,10],[9,10,10]]

In [8]: critique(d_H1,d_H1_moreAA,filename='/ligo/home/jeffrey.kissel/2025-01-21/20240927T211612Z_H1_vs_20240927T211612Z_H1_2x16kDigitalAA_critique.pdf',label=['Nominal','1x More 16k AA'])

LVEA Zone 2A temperature rate alarm limit was relaxed from 0.001 Celsius/minut to 0.0012.  This will reduce daily alarms when the sky is clear and the sun is heating that part of the LVEA.

FAMIS 31069

Due to the chiller swap (alog82379) and PMC/RefCav alignment adjustments (alog82377), several trends show changes from this morning. Most notably is the large drop in RefCav transmission, which likely shows a need to tune up the FSS path on-table.

I've added some logic to ISC_LOCK in DRMI and PRMI to check if we've already tried PRMI, and or CHECK_MICH. Two counters are created and = 0 in each acquisition. On our second pass through PRMI or DRMI in an aquisition the POPAIR timers will be doubled and the threshholds will be halved. The changes have been svn commited.

J. Oberling, R. Short

While recovering the PSL after the chiller swap this morning (alog82379), we noticed the alignment of the PMC and RefCav looked poor according to the cameras, so we took the time to touch these up while still in the maintenance period. Starting with the PMC, I was able to slightly improve the alignment using the picomotors:

Moving on to the RefCav, once the FSS was locked I was able to improve the alignment here as well; the TPD signal increased from 620mV to 737mV. Not as much gain here as we hoped, and the RefCav Refl camera spot still looks like two vertical lobes, so the FSS path may need to be adjusted on-table.

Finally, since the output of the PMC had increased, I changed the ISS RefSignal to bring the diffracted power back down to our desired setpoint of around 4%.

All maintenance activities have completed and we are almost done with initial alignment.

[Ryan, Matthew]

Sweep done at 12:51 PSC, only things to note were two view ports (HAM2 and HAM3) had foil wrap instead of coverings. Unplugged a few scissor lifts, left one by the high bay plugged in.

J. Kissel

Pushing forward and making filter file changes while out of observing, I've changed the arrangement and amount of filtering in the OMC DCPD A1, A2, B1, and B2 banks again. 

Last time, I configured the banks to be analog channel agnostic, so I used up the spare filter banks to include copies of both the A and B analog electronics compensation in each of the four paths (LHO:82313). Talking with Louis and Evan this morning, in light of the discovery that more digital anti-aliasing would be prudent, I instead use up the spare banks for copies of the existing anti-aliasing filters (much like Louis and Evan had done over the weekend).

Check out screenshot of the configuration setup for one test we plan to do: using the 1st attachment in LHO:82375, comparing the performance of the RED, Dec65k*Dec16k filtered nominal trace to GREEN Dec16k*Dec16k and MAGENTA Dec65k*Dec16k*Dec16k. These options, which have lesser phase impact that the straight 2x copy of the existing filters, will also tell us whether it's the frequency content in the signal between 7 and 32 kHz or that of 32 kHz and above that matters in the GW band.

E. Goetz, J. Kissel, L. Dartez

Summary:
There is evidence that some of the lines found in the gravitational wave band are actually aliased lines from higher frequencies. So far it is unclear exactly how many of the lines in the run-averaged list are due to this problem and if the lines are stationary or if violin ring-ups may induce more lines in band due to aliasing artifacts. Further investigation is needed, but this investigation suggests that the current level of anti-aliasing is insufficient to suppress the artifacts at high frequency aliasing into the gravitational wave band.

Details:
We used the live, test-point acquisition of the 524 kHz sampled DCPD data in DTT, channel H1:OMC-DCPD_B1_OUT (equivalent to the nominal H1:OMC-DCPD_B0_OUT channel used in loop). This channel had the same 524k-65k and 65k-16k decimation digital AA filtering applied. The time series was exported from DTT and processed by stitching together the time segments into a single time series. Then one can process the time series using the scipy.signal.welch() function of 1) the full 524 kHz sampled data, 2) 524 kHz sampled data decimated (no additional AA filtering) by a factor of 32 to get the 16 kHz sampled frequency data, 3) 524 kHz sampled data decimated using additional AA filtering by using the scipy.signal.decimate() function which has a built-in anti-aliasing filter.

We also plotted in DTT the individual channels against the 16k H1:OMC-DCPD_B_OUT_DQ channel, showing that some of the lines are visible in the in-loop DCPD 16 kHz channel, but not visible in the test point 524 kHz channels.

Figure 1: ASD of raw 524 kHz data (blue), decimated data without any extra anti-aliasing filter applied (orange), and decimated data with additional anti-aliasing filtering (green). Orange peaks are visible above the blue and green traces above ~750 Hz.
Figure 2: ASD ratio of the decimated data without anti-aliasing filtering to the raw data showing the noise artifacts
Figure 3: Zoom of figure 2 near 758 Hz
Figure 4: ASD computed from DTT showing DCPD B Ch 1, 5, 9, 13 and and the H1:OMC-DCPD_B_OUT_DQ channel at the same time as the 9 and 13 channels were acquired (limitations of the front end handling of 524 kHz test points to DTT)
Figure 5: Zoom of figure 4 near 758 Hz

We were also interested in the time-variability of these artifacts and watched the behaviour of H1:OMC-DCPD_B_OUT_DQ, and saw amplitude variations on the order of factors of a few and frequency shifts on the order of 0.1 Hz, at least for the artifacts observed near 758 Hz.

Figures 4 and 5 indicate that there are more artifacts not necessarily directly caused by aliasing; perhaps these are non-linearity artifacts? This needs further study.

A count of the number of 0.125 Hz frequency bins from 0 to 8192 Hz of the ratio between downsampling without additional anti-aliasing filtering and the raw 524 kHz ASD indicates that ~4900 bins are above a threshold of 1.1 (though most of those bins are above 2 kHz, as indicated by figure 2).