40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
 40m Log, Page 3 of 330 Not logged in
ID Date Author Type Category Subject
16523   Fri Dec 17 22:16:07 2021 YehonathanUpdateBHDSOS assembly

I specifically checked the specification before mounting the mirror. It says clearly "Arrow at the thinnest location pointing towards Side 1". I guess they just ignored it.

As for LO1, I mounted it without noticing the location of the arrow. Later, I checked and the ghost beam was horizontal so I left it as it is. Yeah, I guess I will remount the mirror. Also, what do we do with the pencil markings? It's not vacuum-compatible.

 Quote: We @40m do the convention of the arrow at the thinnest side & pointing the HR side, but nobody says Lambda does the same. We can just remount the mirror without breaking the wires and adjust the pitching if you do it carefully. Does this mean that the LO1 also likely to have the wedge pointing up? Or did you rotate the mirror to have the wedge reflection to be as horizontal as possible?

I had the fear that any mistake in the electronics chain could have been the show stopper.

So I quickly checked the signal assignments for the ADC and DAC chains.

I had initial confusion (see below), but it was confirmed that the electronics chains (at least for LO1) are correct.

Note: One 70ft cable is left around the 1Y0 rack

There are a few points to be fixed:

- It looks like the ADC/DAC card # assignment has been messed up.

CDS DAC0 -> Cable label DAC2 -> AI D2 -> ...
CDS DAC1 -> Cable label DAC0 -> AI D0 -> ...
CDS DAC2 -> Cable label DAC1 -> AI D1 -> ...
(What is going on here... please confirm and correct as they become straight forward)

Once this puzzle was solved I could confirm reasonable connections from the end of the 70 cables to the ADC/DAC.

My system wiring diagram needs to be fixed accordingly too.
This is because the last channel of the first ADC (ADC0) is not available for us and is used for DuoTone.

16521   Fri Dec 17 19:16:45 2021 KojiUpdateBHDSOS assembly

We @40m do the convention of the arrow at the thinnest side & pointing the HR side, but nobody says Lambda does the same.

We can just remount the mirror without breaking the wires and adjust the pitching if you do it carefully.

Does this mean that the LO1 also likely to have the wedge pointing up? Or did you rotate the mirror to have the wedge reflection to be as horizontal as possible?

16520   Fri Dec 17 17:50:17 2021 YehonathanUpdateBHDSOS assembly

I threaded a new wire through a different side block with a magnet and clamped it under a microscope. It was hard, but eventually, I was able to do it by holding the wire on both sides of the side block with weights.

The dented wire was discarded and the side block that was mounted on the AS1 adapter was put aside. I mounted the side block with the new wire on the AS1 adapter.

Anchal and I hanged the AS1 adapter and clamped the wires on the winches of an SOS tower. I balanced the roll and adjusted the height of the magnet with respect to a side OSEM using a camera (attachments 1 & 2).

I shoot the Hene laser on the optic and look at the reflection. I align the laser beam to be as close as possible to the center of the mirror. The OpLev needs to be realigned.

To my surprise, the ghost beam shoots up above the reflected beam! See attachment 3. I check to see that the arrow which marks the thinnest side of the mirror is horizontal (attachment 4). WTF?!

Also, now I realize that the marking on the Lambda optics are pencil markings 😵😵😵.

Anything else? Feel free to edit this entry.

- SUS: AS1 hanging

- SUS: PR3/SR2/LO2 3/4" thick optic hanging

v Electronics chain check (From DAC to the end of the in-air cable / From the end of the in-air cable to the ADC)
[ELOG 16522]

- Long cable installation (4x 70ft)

- In-air cable connection to the flange

- In-vac wiring (connecting LO1 OSEMs)

- LO1 OSEM insertion/alignment

- LO1 Damping test

16518   Thu Dec 16 18:16:36 2021 YehonathanUpdateBHDSOS assembly

Today I glued magnets onto the new 3/4" mirror adapters. I also took the opportunity to make some more side magnets assemblies.

Yesterday I mounted PR3/SR2 3/4" thick mirror onto one of the new adapter. There seem to be no issues for now.

I started the process of suspending AS1 (E2000226-A). The Lambda Optic mirror with the closest specs has Rc = 2 m. I attached side blocks with clamped wires onto adapter number 7 - side block with a magnet on the right.

I then took one of the Lambda Optic mirrors and tried mounting it in the adapter. It was quite difficult to get it right. Unfortunately, I chipped the edge of the substrate (attachment 1) 🤦🏻‍♂️. I put the mirror back in the box and decided to use the spare mirror. I successfully mounted it into the adapter but when I put the clamping screws one of them fell on the mirror 🤦🏻‍♂️🤦🏻‍♂️. There is no visible damage though. I took some pictures (attachment 2-4).

I and Anchal then started suspending the mirror but then we found that one of the wires is dented in the middle 🤦🏻‍♂️🤦🏻‍♂️🤦🏻‍♂️. I'm burned out for today.

Late update: one nice thing that I found yesterday is that the glue is viscous enough to hold the dumbells without a metal sheet from above holding the magnets. This greatly simplifies the gluing process.

16517   Thu Dec 16 17:57:17 2021 AnchalUpdateBHDFinished Coil driver (odd serial number) units tests

S2100619 was fixed by Koji and it passed the test after that.

Quote:
 SERIAL # S2100619 FAIL (CH2 phase)

16516   Thu Dec 16 17:41:12 2021 KojiUpdateBHDCoil driver test failed for S2100619-v1

Good catch. It turned out that the both + and - side of the output stages for CH2 were oscillating at ~600kHz. If I use a capacitance sticks to touch arbitrarily around the components, it stops their oscillation and they stay calm.
It means that the phase margin becomes small while the circuit starts up.

I decided to increase the capacitances C6 and C20 (WIMA 150pF) to 330pF (WIMA FPK2 100V) and the oscillation was tamed. 220pF didn't stop them. After visually checked the signal behavior with an oscilloscope, the unit was passed to Anchal for the TF test.

The modification was also recorded in the DCC S2100619

16515   Thu Dec 16 15:54:08 2021 KojiUpdateElectronicsITMX feedthroughs and in-vac cables installed

Thanks for the installation.

With regard to the connector convention, let's use the attached arrangement so that it will be consistent with the existing flange DSUB configuration. Not a big deal.

16514   Thu Dec 16 15:32:59 2021 AnchalUpdateBHDFinished Coil driver (odd serial number) units tests

We have completed modifications and testing of the HAM Coil driver D1100687 units with serial numbers listed below. The DCC tree reflects these changes and tests (Run/Acq modes transfer functions).

 SERIAL # TEST result S2100609 PASS S2100611 PASS S2100613 PASS S2100615 PASS S2100617 PASS S2100619 FAIL (CH2 phase) S2100621 PASS S2100623 PASS S2100625 PASS S2100627 PASS S2100629 PASS S2100631 PASS S2100633 Waiting for more components S2101649** PASS S2101651** PASS S2101653** PASS S2101655** PASS

** A fix had to be done on the DC power supply for these. The units' regulated power boards were not connected to the raw DC power, so the cabling had to be modified accordingly.

Further, Paco fixed the two even serial number units (S2101648, S211650) that failed the test. The issues were minor soldering mistakes that were easily resolved.

16513   Thu Dec 16 15:04:12 2021 ChubUpdateElectronicsITMX feedthroughs and in-vac cables installed

The ITMX 10" flange with four DSUB-25 feedthroughs has been install with the cables connected at the in-vac side.  See photo; as requested, LO1-1 and LO1-2 are connected to the top row of feedthroughs from left to right respectively and the opposite ends of the cables placed left to right on the laser table.  PR2-1 and PR2-2 are connected to the lower row of feedthroughs from left to right respectively, with the opposite ends placed on the surface below the laser from left to right.  This seemed the easiest way to keep the cable orientation clear.

16512   Thu Dec 16 12:21:16 2021 AnchalUpdateBHDCoil driver test failed for S2100619-v1

Today I found one of the coil driver boards, S2100619 failed the test on CH2. There appears to be an extra phase lag after 10 kHz and some resonant-like feature at 7 kHz. This of course is very high-frequency stuff and maybe we don't care about these deviations. But it could mean something is off with the channel and could potentially lead to failure in the relevant frequency band in the future. I'll need help to debug this. Please see the attachment for details of test failure.

16511   Wed Dec 15 22:09:53 2021 YehonathanUpdateBHDSOS assembly

{Paco, Yehonathan, Anchal}

I cut the wires from the winches and removed the winch adapter plate. I engrave  'LO1' on the SOS tower. Me and Paco wrapped the SOS with foil and transported it to Anchal who put it inside the ITMX chamber.

The transportation seems to be successful. Nothing broke. However, we found that even with the short side OSEM the LO1, as it is now, cannot sit in its designed position since ITMX's side OSEM stands in its way.

If there are magnets on both ITMX sides we can move its side OSEM to the other side. Another option is to resuspend LO1 with a side magnet on its left side.

16510   Wed Dec 15 17:44:18 2021 KojiSummaryBHDPart VIII of BHR upgrade - Placed LO1

If ITMX already has another side magnet, we can migrate the side OSEM of ITMX to the other side. This way, the interference of the OSEMs can be avoided.

16509   Wed Dec 15 16:11:38 2021 AnchalSummaryBHDPart VIII of BHR upgrade - Placed LO1

[Anchal, Yehonathan, Paco]

Today we opened ITMX chamber and removed the following optics and placed them in the Xend flow bench (see attachment 1):

• POPM1
• POPM2

Yehonathan brought his first SOS baby next to ITMX chamber. The suspension was carried by hands throughout. He gave me the suspension over the IMC beam tube from where I placed it on the table. I checked through the OSEMs and the face magnets were still on. I could not verify the side magnet but nothing seemed out of place.

I then moved LO1 near its planned place. I had to bolt it at 1 inch North and 0.5 inch West of its planned position because the side OSEM on ITMX is long and protrudes out of the base footprint. Even if it was small, the current layout would make the OSEM pins of the side OSEMs of ITMX and LO1 very near each other. So we can not place LO1 closer to ITMX from current position. This means the layout needs to be redesigned a bit for the modified position of LO1. I believe it will significantly shift and turn the beam from LO1 to LO2, so we might need to change the beam upstream from TT2 onwards. More discussion is required.

Unfortunately, what I thought was clicking photos was just changing modes between video and image mode, so I have no photos from today but only a video that I recorded in the end.

16508   Wed Dec 15 15:06:08 2021 JordanUpdateVACVacuum Feedthru Install

Jordan, Chub

We installed the 4x DB25 feedthru flange on the North-West port of ITMX chamber this afternoon. It is ready to go.

16507   Wed Dec 15 13:57:59 2021 PacoUpdateComputersupgraded ubuntu on zita

[Paco]

Upgraded zita's ubuntu and restarted the striptool script.

16506   Tue Dec 14 19:29:42 2021 PacoUpdateBHD1Y0 rack work for LO1

[Paco]

Two coil drivers have been installed on 1Y0 (slots 6, 7, for LO1 SOS). All connections have been made from the DAC, AI board, DAC adapter, Coil driver, Sat Amp box. Then no SOS load installed, all return connections have been made from Sat Amp box, ADC adapter, AA board, and to ADC. We will continue this work tomorrow, and try to test everything before closing the loop for LO1 suspension.

16505   Tue Dec 14 14:02:33 2021 YehonathanUpdateBHDSOS assembly

{Yehonathan, Paco}

I fixed the overall height of the adapter (attachment 1). I put an OSEM next to the side magnet. I positioned a camera in front of the SOS and connected it to my laptop for live streaming. I painted a line indicating the height of the OSEM plates and a line in between with the mean height. I discarded the wire clamp on the suspension block I released the wires from the winches and pulled on them until the magnet was roughly in the right position. I clamped the wires back on the winches and adjusted them until the magnet was on the middle line. I also verified that the roll of the adapter is aligned as before by making sure that the horizontal features on the adapter are parallel to the horizontal features on the SOS tower.

The wires were clamped to the suspension block using a new wire clamp.

I Found that locking the counterweight setscrew changes the alignment. Today we verified this effect. We released the setscrew and pre-compensated by adjusting the counterweight such that when the setscrew was locked the mirror was aligned.

We measured mechanical resonances (attachment 2). This time the yaw motion was very quiet so we got a smaller peak for the yaw. The peaks are the same as before. Y readout has peaks around the pitch and pos resonances that don't appear in the X readout. I'm not sure what they are. Maybe coming from the QPD electronics.

We locked the adapter using the EQ stops. We made sure the alignment stays close to ideal.

We installed OSEMs on the SOS. We chose suboptimal OSEMs because LO1 will only be used for steering. I made a spreadsheet copying the OSEM catalog into it. There we mark which OSEM goes where.

I cleaned the optic using the ion gun with a pressure of 30 PSI.

The next steps are:

1. Engrave the SOS tower.

2. Cut the wire at the winches and remove the winch adapter plate.

3. Wrap the SOS with foil.

4. Install the SOS in the vacuum chamber.

16504   Tue Dec 14 11:33:29 2021 TegaUpdatePEMgit repo for temp sensor and sus medm

[Temperature sensor]

Added new temp EPICs channels to database file (/cvs/cds/caltech/target/c1pem1/tempsensor/C1PEMaux.db)

Added new temp EPICs channels to slow channels ini file (/opt/rtcds/caltech/c1/chans/daq/C0EDCU.ini)

[SUS medm screen]

Moved new SUS screen to location : /opt/rtcds/userapps/trunk/sus/c1/medm/templates/NEW_SUS_SCREENS

Place button on the new screen to link to the old screen and replace old screens link on sitemap.

Fixed Load Coefficient button location issue

Fixed LOCKIN flow diagram issue

Fixed watchdog labelling issue

Linked STATE WORD block to FrontEnd STATUS screen

Replaced the 2x1 pit/yaw filter screens for LOCK and DAMP fliters with 3x1 LPY filter screen

*Need some more time to figure out the OPTLEV red indicator

16503   Mon Dec 13 15:05:47 2021 TegaUpdatePEMgit repo for temp sensor and sus medm

[temperature sensor]

git repo: https://git.ligo.org/40m/tempsensor.git

todo

Update the temp sensor channels to fit with cds format, ie. "C1:PEM-TEMP_EX", "C1:PEM-TEMP_EY", "C1:PEM-TEMP_BS"

- Use FLOAT32_LE data format for the database file (/cvs/cds/caltech/target/c1pem1/tempsensor/C1PEMaux.db) to create the new channels.

- Keep the old datadase code and channels so we can compare with new temp channels afterwards. Also we need a 1-month overlap b4 deleting the old channels.

[sus medm screen]

git repo: https://git.ligo.org/40m/susmedmscreen.git

todo (from talk with Koji)

- Damp filter and Lock filter buttons should open a 3x1 filter screen so that the 6 filters are opened by 2 buttons compared to the old screen that has 3 buttons connected to 2X1 filter screen

- Make the LOCKIN signla modulation flow diagramlook more like the old 40m screen since that is a better layout

- Move load coefficient button to top of sus medm screen (beside stateword)

- The rectangular red outline around the oplev display is confusing and needs to be modified for clarity

- COMM tag block should not be 3D as this suggests it is a button. Make it flat and change tag name to indicate individual watchdog control as this better reflect its functionality. Rename current watchdog switch to watchdog master is it does what the 5 COMM switches do at once.

- Macro pass need to be better documented so that when we call the sus screens from locations other than sitemap, we should know what macro variables to pass in, like DCU_ID etc.

- Edit sitemap.adl to point only to the new screens. Then create a button on the new screen that points to the old screen. This way, we can still access the old screen without clogging sitemap.

- Move the new screen location to a subfolder of where the current sus screens reside, /opt/rtcds/userapps/trunk/sus/c1/medm/templates

- Rename the overview screen (SUS_CUST_HSSS_OVERVIEW.adl) to use the SUS_SINGLE nomenclature, i.e. SUS_SINGLE_OVERVIEW.adl

- Keep an eye of the cpu usage of c1pem as we add BLRMS block for other optics.

https://dcc.ligo.org/LIGO-S2101689

https://dcc.ligo.org/LIGO-S2101690

https://dcc.ligo.org/LIGO-S2101691

https://dcc.ligo.org/LIGO-S2101692

The units are completely passive right now and has option to extend to have a dewhitening board added inside.
So the power switch does nothing.

Some of the components for the dewhitening enhancement are attached inside the units.

16501   Fri Dec 10 19:22:01 2021 KojiUpdateVACPumping down the RGA section

The scan result was ~x10 higher than the previously reported scan on 2020/9/15 (https://nodus.ligo.caltech.edu:8081/40m/15570), which was sort of high from the reference taken on 2018/7/18.

This just could mean that the vacuum level at the RGA was x10 high.
We'll just go ahead with the vacuum repair and come back to the RGA once we return to "vacuum normal".

Meanwhile, I asked Jordan to turn off the RGA to make it cool down. I shut off RGA section and turned TP2 off.

16500   Fri Dec 10 18:55:58 2021 TegaUpdateCDSNew SUS medm screen update

Turns out the BLRMS monitoring channels for MC1, MC2, MC3, ITMY and SRM already exist in c1pem. So I modified the new SUS screen to display the BLRMS info for the aforementioned optics. Next step is to add the BLRMS monitor for PRM, ITMX, ETMX and ETMY. This would require extending the number of inputs for the "SUS" block in c1pem to accomodate the additional inputs from the remaining optics.

16499   Fri Dec 10 15:59:23 2021 PacoUpdateBHDFinished Coil driver (even serial number) units tests

[Paco, Anchal]

We have completed modifications and testing of the HAM Coil driver D1100687 units with serial numbers listed below. The DCC tree reflects these changes and tests (Run/Acq modes transfer functions).

 SERIAL # TEST result S2100608 PASS S2100610 PASS S2100612 PASS S2100614 PASS S2100616 PASS S2100618 PASS S2100620 PASS S2100622 PASS S2100624 PASS S2100626 PASS S2100628 PASS S2100630 PASS S2100632 PASS S2101648** FAIL (Ch1, Ch3 run mode) S2101650** FAIL (Ch3 run mode) S2101652** PASS S2101654** PASS

** A fix had to be done on the DC power supply for these. The units' regulated power boards were not connected to the raw DC power, so the cabling had to be modified accordingly (see Attachment #1)

16498   Fri Dec 10 13:02:47 2021 Ian MacMillanSummaryComputersQuantization Noise Calculation Summary

I am trying to replicate the simulation done by Matt Evans in his presentation  (see Attachment 1 for the slide in particular).

He defines his input as $x_{\mathrm{in}}=sin(2\pi t)+10^{-9} sin(2\pi t f_s/4)$ so he has two inputs one of amplitude 1 at 1 Hz and one of amplitude 10^-9 at 1/4th the sampling frequency  in this case: 4096 Hz

For his filter, he uses a fourth-order notch filter. To achieve this filter I cascaded two second-order notch filters (signal.iirnotch) both with locations at 1 Hz and quality factors of 1 and 1e6. as specified in slide 13 of his presentation

I used the same procedure outlined here. My results are posted below in attachment 2.

Analysis of results:

As we can see from the results posted below the results don't match. there are a few problems that I noticed that may give us some idea of what went wrong.

First, there is a peak in the noise around 35 Hz. this peak is not shown at all in Matt's results and may indicate that something is inconsistent.

the second thing is that there is no peak at 4096 Hz. This is clearly shown in Matt's slides and it is shown in the input spectrum so it is strange that it does not appear in the output.

My first thought was that the 4kHz signal was being entered at about 35Hz but even when you remove the 4kHz signal from the input it is still there. The spectrum of the input shown in Attachment 3 shows no features at ~35Hz.

The Input filter, Shown in attachment 4 shows the input filter, which also has no features at ~35Hz. Seeing how the input has no features at ~35Hz and the filter has no features at ~35Hz there must be either some sort of quantization noise feature there or more likely there is some sort of sampling effect or some effect of the calculation.

To figure out what is causing this I will continue to change things in the model until I find what is controlling it.

I have included a Zip file that includes all the necessary files to recreate these plots and results.

16497   Thu Dec 9 21:57:35 2021 YehonathanUpdateBHDSOS assembly

{Yehonathan, Tega}

We took the optic out of the SOS tower and removed the side blocks. We mounted new side blocks with wires already clamped in them in the reverse order.

The Adapter was placed back into the SOS and the wires were threaded through the wire clamp and suspended on the winches. The roll of the optic was balanced using a camera (attachment 1).

The pitch was balanced. this time I used 2 counterweights instead of 1 in order to not have to take so much of the weight out.

The mechanical resonances were measured by taking a 100 sec time series of QPD readout and doing PSD estimation (attachment 2). The mirror motion was damped as much as possible before taking the measurement.

3 peaks below 1.5Hz can be seen with frequencies of 755mHz (Yaw), 942mHz, 1040mHz (Pitch + Pos). The pitch/pos peaks are a bit close to each other, I bet if we go back to 1 counterweight the situation will be better.

While inserting the side OSEM I realized I didn't check the overall height of the adapter. The magnet was too high. I will fix it tomorrow and repeat the roll balancing.

16496   Thu Dec 9 18:22:36 2021 TegaUpdateCDSNew SUS medm screen update

Work on the medm screen for SUS RMS monitor is ongoing. The next step would be to incorporate this into the SUS medm screen, add the BLRMS model to the SUS controller model, recompile, check that the channels are being correctly addressed, then load the appropriate bandpass and lowpass filters.

16495   Thu Dec 9 00:32:56 2021 TegaUpdateCDSNew SUS medm screen update

The new SUS screen can be reached via sitemap -> IFO SUS button -> NEW ETMX dropdown menu link. Please use and provide feedback. Not sure exactly if we need/want the display screens after the IOP model on the right of the medm screen. I have not been able to locate the corresponding channels but did not want to remove them until I was sure that we don't plan to add these features to our screens. When all bugs have been ironed out, we can use appropriate macro substitution for the other optics.

The next feature to add is the BLRMS to the coil and PD channels. I plan to combine the PEM BLRMS medm implementation with the sus_single_BLRMS model block (located in  /opt/rtcds/userapps/release/cds/c1/models). This way we use the latest BLRMS block in "/opt/rtcds/userapps/release/cds/common/models/BLRMS.mdl" whilst also leveraging the previous work done on the sus_single_BLRMS model, which neatly fits into our current SUS model.

16494   Wed Dec 8 10:14:43 2021 JordanUpdateVACPumping down the RGA section

After an overnight pumpdown/RGA warm up, I took a 100 amu scan of the RGA volume and subsequent pumping line. Attached is a screenshot along with the .txt file. Given the high argon peak (40) and the N2/O2 ratio, it looks like there is a decent sized air leak somehwere in the volume.

Are we interested in the hydrocarbon leak rates of this volume? That will require another scan with one of the calibrated leaks opened.

Edit: Added a Torr v AMU plot to see the partial pressures

 Quote: So that Jordan can run the RGA scan this afternoon, I ran TP3 and started pumping down the RGA section. Procedure: - Same 1~4 - Same 5 - 6 Opened only the backing path for TP3 - 7 Turned on TP3 only - TP3 reached the nominal full speed @75kRPM - 11 Opened V5 to pump the pump spool -> Immediately reached P3<1e-4 - 12 Opened VM3 to pump the RGA section -> Immediately reached P4<1e-4 The pumps are kept running. I'll come back later to shut down the pumps. => Jordan wants to heat the filament (?) and to run the scan tomorrow. So we decided to keep TP3 running overnight. I switched TP3 to the stand-by mode (= lower rotation speed @50kRPM)

16493   Tue Dec 7 13:12:50 2021 KojiUpdateVACPumping down the RGA section

So that Jordan can run the RGA scan this afternoon, I ran TP3 and started pumping down the RGA section.

Procedure:
- Same 1~4
- Same 5
- 6 Opened only the backing path for TP3
- 7 Turned on TP3 only

- TP3 reached the nominal full speed @75kRPM

- 11 Opened V5 to pump the pump spool -> Immediately reached P3<1e-4
- 12 Opened VM3 to pump the RGA section -> Immediately reached P4<1e-4

The pumps are kept running. I'll come back later to shut down the pumps.
=> Jordan wants to heat the filament (?) and to run the scan tomorrow.
So we decided to keep TP3 running overnight. I switched TP3 to the stand-by mode (= lower rotation speed @50kRPM)

16492   Tue Dec 7 10:55:25 2021 Ian MacMillanSummaryComputersQuantization Noise Calculation Summary

[Ian, Tega]

Tega and I have gone through the IIR Filter code and optimized it to make sure there aren't any areas that force high precision to be down-converted to low precision.

For the new biquad filter we have run into the issue where the gain of the filter is much higher than it should be. Looking at attachments 1 and 2, which are time series comparisons of the inputs and outputs from the different filters, we see that the scale for the output of the Direct form II filter shown in attachment 1 on the right is on the order of 10^-5 where the magnitude of the response of the biquad filter is on the order of 10^2. other than this gain the responses look to be the same.

I am not entirely sure how this gain came into the system because we copied the c code that actually runs on the CDS system into python. There is a gain that affects the input of the biquad filter as shown on this slide of Matt Evans Slides. This gain, shown below as g, could decrease the input signal and thus fix the gain. However, I have not found any way to calculate this g.

With this gain problem we are left with the quantization noise shown in Attachment 4.

State Space:

I have controlled the state space filter to act with a given precision level. However, my code is not optimized. It works by putting the input state through the first state-space equation then integrating the result, which finally gets fed through the second state-space equation.

This is not optimized and gives us the resulting quantization noise shown in attachment 5.

However, the state-space filter also has a gain problem where it is about 85 times the amplitude of the DF2 filter. Also since the state space is not operating in the most efficient way possible I decided to port the code chris made to run the state-space model to python. This code has a problem where it seems to be unstable. I will see if I can fix it

16491   Mon Dec 6 18:23:04 2021 YehonathanUpdateBHDSOS assembly

I installed OSEMs on the LO1 SOS. To my surprise, the side magnet is not in the center of the side OSEM. It completely misses the LED as can be seen in the attachment.

Looking at the CAD model, it turns out the position of the OSEM on the left side plate is different from the position of the OSEM on the right side plate in the SOS tower.

We need to take the optic down, swap the right and left side blocks, and resuspend it.

There is not enough 🤦🏻‍♂️in the world

16490   Mon Dec 6 14:26:52 2021 KojiUpdateVACPumping down the RGA section

Jordan reported that the RGA section needs to be pumped down to allow the analyzer to run at sufficiently low pressure (P<1e-4 torr).
The RGA section was pumped down with the TP2/TP3. The procedure is as listed below.
If the pressure go up to P>1e-4 torr, we need to keep the pump running until the scan is ready.

----
### Monitor / Control screen setup ###
1. On c1vac: cd /cvs/cds/caltech/target/c1vac/medm
3. RGA section (P4) 3.6e-1 torr / P3/P2 still atm.

### TP2/TP3 backing ###
5. Turn on AUX RP with the circuit breaker hanging on the AC.
6. Open manual valve for TP2/3/ backing / PTP2/3 ~ 8torr

### TP2/TP3 starting ###
7. Turned on TP2/TP3 with the Standby OFF

### Pump down the pump spool ###
8. Connect manual RP line (Quick Connect)
9. Turned on RP1/RP3 -> quickly reached 0.4 torr
10. Open V6 for pump spool pumping -> Immediately go down to sufficiently low pressure for TP2/TP3.
(10.5 I had to close V6 at this point)
11. Open V5 to start pumping pump spool with TP3 (TP2 still stand by) -> P3 immediately goes down below 1e-4 torr. This automatically closed V6 because of the low pressure of P3 (interlocking)

### Pump down the RGA section ###
12. Open VM3 to pump down RGA section -> P4 goes down to <1e-4 torr
13. P2 is still 2e-3. So decided to open V4 to use TP2 (now it's ready) too. -> Saturated at 1.7e-3

### Shutting down ###
14. Close VM3
15. Close V4/V5 to isolate TP2/TP3
16. Stop TP2/TP3 -> Slowing down
17. Stop RP1/RP3
18. Close the manual valves for TP2/3/ backing
19. Stop AUX RP with the circuit breaker hanging on the AC.

16489   Wed Dec 1 12:57:08 2021 YehonathanUpdateBHDSOS assembly

{Yehonathan, Tega}

We glued some magnets onto modified side blocks. We followed pretty much the same procedure form last time. The music wires were clamped on the side blocks away from the optic adapter. The clamps were screwed down using the low profile screws the Jordan obtained from UC components to allow maximum clearance for the coils (Attachment 2).

The modified side blocks with wires already clamped in them were screwed onto the adapter. We put the adapter on the EQ stops and do rough adjustments, making sure the side magnet is roughly on the center of a coil we inserted to the side. The wires were threaded through the wire clamp on the suspension block and clamped on the winches. We realigned the Oplev beam such that it is parallel to the table using the quad photodiode.

We balanced the height of the adapter. This time we used a camera facing the adapter. The tilt of the camera was set by drawing a straight line (red line in attachment 1) such that the bottom clamps on the SOS are parallel to it.

Then, we adjusted the winches such that the screws on the side blocks are on the same green line on attachment 1.

Once the height was adjusted, we clamped the wire on the suspension block and cut it at the winches.

We balanced the optic. We had to take out the front counterweight to get the balance (attachment 3).

We checked whether the wire is touching anything. We confirmed that it doesn't. The wire goes nicely through the new hole on the side blocks (attachment 4, 5).

We measure the resonance frequency for both yaw (attachment 6) and pitch (attachment 7). They both seem to be sub-Hz. The pitch trace also shows that the oscillations are centered around 0 meaning the optic is balanced.

There's still the issue of what to do with the lower back EQ stop not touching the adapter.

16488   Tue Nov 30 17:11:06 2021 PacoUpdateGeneralMoved white rack to 1X3.5

[Paco, Ian, Tega]

We moved the white rack (formerly unused along the YARM) to a position between 1X3, and 1X4. For this task we temporarily removed the hepas near the enclosures, but have since restored them.

16487   Tue Nov 30 11:03:44 2021 YehonathanMetaphysicsGeneralToilet tank broken

a plumber came in yesterday and fixed the issue.

 Quote: The toilet tank in the big bathroom stopped refilling. I contacted PPService@caltech.edu and put up an "Out of Order sign".

16486   Mon Nov 29 15:24:53 2021 HangHowToGeneralFisher matrix vs length of each FFT segment

We have been discussing how does the parameter estimation depends on the length per FFT segment. In other words, after we collected a series of data, would it be better for us to divide it into many segments so that we have many averages, or should we use long FFT segments so that we have more frequency bins?

My conclusions are that:

1). We need to make sure that the segment length is long enough with T_seg > min[ Q_i / f_i ], where f_i is the resonant frequency of the i'th resonant peak and the Q_i its quality factor.

2). Once 1) is satisfied, the result depends weakly on the FFT length. There might be a weak hint preferring a longer segment length (i.e., want more freq bins than more averages) though.

=================================================================

To reach the conclusion, I performed the following numerical experiment.

I considered a simple pendulum with resonant frequency f_1 = 0.993 Hz and Q_1 = 6.23. The value of f_1 is chosen such that it is not too special to fall into a single freq bin. Additionally, I set an overall gain of k=20. I generated T_tot = 512 s of data in the time domain and then did the standard frequency domain TF estimation. I.e., I computed the CSD between excitation and response (with noise) over the PSD of the excitation. The spectra of excitation and noise in the readout channel are shown in the first plot.

In the second plot, I showed the 1-sigma errors from the Fisher matrix calculation of the three parameters in this problem, as well as the determinant of the error matrix \Sigma = inv(Fisher matrix). All quantities are plotted as functions of the duration per FFT segment T_seg. The red dotted line is [Q_1/f_1], i.e., the time required to resolve the resonant peak. As one would expect, if T_seg <~ (Q_1/f_1), we cannot resolve the dynamics of the system and therefore we get nonsense PE results. However, once T_seg > (Q_1/f_1), the PE results seem to be just fluctuating (as f_1 does not fall exactly into a single bin). Maybe there is a small hint that longer T_seg is better. Potentially, this might be due to that we lose less information due to windowing? To be investigated further...

I also showed the Fisher estimation vs. MCMC results in the last two plots. Here each dot is an MCMC posterior. The red crosses are the true values, and the purple contours are the results of the Fisher calculations (3-sigma contours). The MCMC results showed similar trends as the Fisher predictions and the results for T_seg = (32, 64, 128) s all have similar amounts of scattering << the scattering of the T_seg=8 s results. Though somehow it showed a biased result. In the third plot, I manually corrected the mean so that we could just compare the scattering. The fourth plot showed the original posterior distribution.

16485   Wed Nov 24 17:13:31 2021 YehonathanMetaphysicsGeneralToilet tank broken

The toilet tank in the big bathroom stopped refilling. I contacted PPService@caltech.edu and put up an "Out of Order sign".

16484   Wed Nov 24 14:34:15 2021 YehonathanUpdateBHDSaving on SUSAUX slow channels

Koji found out that the stock for BIO Acromag modules is very low and that the lead time for ordering new ones is ~ 1-year X-o.

We figure we might need to minimize the number of modules but still keep the Acromag chassis functional.

Looking at the new C1AUXEY feed-throughs spreadsheet one can see that we actually normally need only 1 BIO (not 2) module since there are 16 suspensions related bios + 1 green shutter which is unrelated to SUSAUX so there is no room to cut back here.

There are 16 analog input channels, 5 for PDMONs and 5 VMONs, and 6 spares which require 2 ADCs. Removing the spares and 2 monitoring channels will be enough to get us to 1 ADC.

16483   Wed Nov 24 14:15:15 2021 YehonathanUpdateBHDSOS assembly

Late update. We got 2 modified side blocks from Jordan a few days ago. Yesterday, I glued a side magnet to one of the modified side blocks.

I took the opportunity to reglue some magnets that were knocked off from the adapters. I did this for 2 adapters only since w need 4 shallow adapters and we already had 2 complete ones.

Today, Jordan gave us the rest of the modified side blocks clean and baked. We are ready to suspend a mirror today.

16482   Wed Nov 24 13:44:19 2021 ranaSummaryComputersQuantization Noise Calculation Summary

This looks great. I think what we want to see mainly is just the noise in the 32 bit IIR filtering subtracted from the 64 bit one.

It would be good if Tega can look through your code to make sure there's NO sneaky places where python is doing some funny casting of the numbers. I didn't see anything obvious, but as Chris points out, these things can be really sneaky so you have to be next level paranoid to really be sure. Fox Mulder level paranoia.

And, we want to see a comparison between what you get and what Denis Martynov put in an appendix of his thesis when comparing the Direct Form II, with the low-noise form (also some slides from Matt Evans on thsi from a ~decade agoo). You should be able to reproduce his results. He used matlab + C, so I am curious to see if it can be done all in python, or if we really need to do it in C.

And then...we can make this a part of the IFOtest suite, so that we point it at any filter module anywhere in LIGO, and it downloads the data and gives us an estimate of the digital noise being generated.

16481   Wed Nov 24 11:02:23 2021 Ian MacMillanSummaryComputersQuantization Noise Calculation Summary

I added mpmath to the quantization noise code. mpmath allows me to specify the precision that I am using in calculations. I added this to both the IIR filters and the State-space models although I am only looking at the IIR filters here. I hope to look at the state-space model soon.

Notebook Summary:

I also added a new notebook which you can find HERE. This notebook creates a signal by summing two sine waves and windowing them.

Then that signal is passed through our filter that has been limited to a specific precision. In our case, we pass the same signal through a number of filters at different precisions.

Next, we take the output from the filter with the highest precision, because this one should have the lowest quantization noise by a significant margin, and we subtract the outputs of the lower precision filters from it. In summary, we are subtracting a clean signal from a noisy signal; because the underlying signal is the same, when we subtract them the only thing that should be left is noise. and since this system is purely digital and theoretical the limiting noise should be quantization noise.

Now we have a time series of the noise for each precision level (except for our highest precision level but that is because we are defining it as noiseless). From here we take a power spectrum of the result and plot it.

After this, we can calculate a frequency-dependent SNR and plot it. I also calculated values for the SNR at the frequencies of our two inputs.

This is the procedure taken in the notebook and the results are shown below.

Analysis of Results:

The first thing we can see is that the precision levels 256 and 128 bits are not shown on our graph. the 256-bit signal was our clean signal so it was defined to have no noise so it cant be plotted. The 128-bit signal should have some quantization noise but I checked the output array and it contained all zeros. after further investigation, I found that the quantization noise was so small that when the result was being handed over from mpmath to the general python code it was rounding those numbers to zero. To overcome this issue I would have to keep the array as a mpmath object the entire time. I don't think this is useful because matplotlib probably couldn't handle it and it would be easier to just rewrite the code in C.

The next thing to notice is sort of a sanity check thing. In general, low precision filters yield higher noise than high precision. This is a good quick sanity check. However, this does not hold true at the low end. we can see that 16-bit actually has the highest noise for most of the range. Chris pointed out that at low precisions that quantization noise can become so large that it is no longer a linearly coupled noise source. He also noted that this is prone to happen for low precision coefficients with features far below the Nyquist frequency like I have here. This is one explanation that seems to explain the data especially because this ambiguity is happening at 16-bit and lower as he points out.

Another thing that I must mention, even if it is just a note to future readers, is that quantization noise is input dependent. by changing the input signal I see different degrees of quantization noise.

Analysis of SNR:

One of the things we hoped to accomplish in the original plan was to play around with the input and see how the results changed. I mainly looked at how the amplitude of the input signal scaled the SNR of the output. Below I include a table of the results. These results were taken from the SNR calculated at the first peak (see the last code block in the notebook) with the amplitude of the given sine wave given at the top of each column. this amplitude was given to both of the two sine waves even though only the first one was reported. To see an example, currently, the notebook is set up for measurement of input amplitude 10.

 0.1 Amplitude of input 1 Amplitude 100 Amplitude 1000 Amplitude 4-bit SNR 5.06e5 5.07e5 5.07e5 5.07e5 8-bit SNR 5.08e5 5.08e5 5.08e5 5.08e5 16-bit SNR 2.57e6 8.39e6 3.94e6 1.27e6 32-bit SNR 7.20e17 6.31e17 1.311e18 1.86e18 64-bit SNR 6.0e32 1.28e32 1.06e32 2.42e32 128-bit SNR unknown unknown unknown unknown

As we can see from the table above the SNR does not seem to relate to the amplitude of the input. in multiple instances, the SNR dips or peaks in the middle of our amplitude range.

16480   Tue Nov 23 18:02:05 2021 AnchalUpdateIMCMC autolocker shifted to python3 script running in docker

I finished copying over the current autolocker bash script functionality into a python script which runs using a simple configuration yaml file. To run this script, one needs to ssh into optimus and :

controls@optimus|~> cd /opt/rtcds/caltech/c1/Git/40m/scripts/MC
controls@optimus|MC> sudo docker-compose up -d
Creating mc_AL_MC_1 ... done

That's it. To check out running docker processes, one can:

controls@optimus|MC> sudo docker ps

And to shut down this particular script, in the same directory, one can

controls@optimus|MC> sudo docker-compose down
Removing mc_AL_MC_1 ... done

If the docker image requires to be rebuild in future, go to the directory where Dockerfile is present and run:

controls@optimus|MC> sudo docker build -t pyep .

I had to add PyYAML package in the pyepics docker image already present on docker hub, thanks to Andrew.

For now, I have disabled the MCautolocker service on Megatron. To start it back again, one would need to ssh into megatron and do following:

~> sudo systemctl enable MCautolocker
~> sudo systemctl start MCautolocker

Let's see for a day how this new script does. I've left PSL shutter open and autolocker engaged.

To do: Fix the C1:IFO-STATE epics channel definition so that it takes its bits from separate lock status channels instead of scripts writign the whole word arbitrarily.

16479   Mon Nov 22 17:42:19 2021 AnchalUpdateGeneralConnected Megatron to battery backed ports of another UPS

[Anchal, Paco]

I used the UPS that was providing battery backup for chiara earlier (a APS Back-UPS Pro 1000), to provide battery backup to Megatron. This completes UPS backup to all important computers in the lab. Note that this UPS nominally consumes 36% of UPS capacity in power delivery but at start-up, Megatron was many fans that use up to 90% of the capacity. So we should not use this UPS for any other computer or equipment.

While doing so, we found that PS3 on Megatron was malfunctioning. It's green LED was not lighting up on connecting to power, so we replaced it from the PS3 of old FB computer from the same rack. This solved this issue.

Another thing we found was that Megatron on restart does not get configured to correct nameserver resolution settings and loses the ability to resolve names chiara and fb1. This results in the nfs mounts to fail which in turn results in the script services to fail. We fixed this by identifying that the NetworkManager of ubuntu was not disabled and would mess up the nameserver settings which we want to be run by systemd-resolved instead. We corrected the symbolic link: /etc/resolv.conf -> /run/systemd/resolve/resolv.conf. the we stopped and diabled the NetworkManager service to keep this persistent on reboot. Following are the steps that did this:

> sudo rm /etc/resolv.conf
> ln -s /etc/resolv.conf /run/systemd/resolve/resolv.conf
> sudo systemctl stop NetworkManager.service
> sudo systemctl disable NetworkManager.service

16478   Mon Nov 22 16:38:26 2021 TegaSummaryComputer Scripts / ProgramsSUS Plant Plan for New Optics

[Tega, Ian]

TODO

1. Investigate cross-coupling btw the various degrees of freedom (dof) - turn on noise for each dof in the plant model and measure the transfer function of the other dofs.

2. Get a closed-loop transfer function using noise injection and give a detailed outline of the procedure in elog - IN1/IN2 for each TM_RESP filter while the others are turned off.

3. Derive analytic model of the closed-loop transfer functions for comparison.

4. Adapt control filters to fit optimized analytical solutions.

16477   Thu Nov 18 20:00:43 2021 Ian MacMillanSummaryComputer Scripts / ProgramsSUS Plant Plan for New Optics

[Ian, Raj, Tega]

Here is the comparison between the results of Raj's python model and the transfer function measurement done on the plant model by Tega and me.

As You can see in the graphs there are a few small spots of disagreement but it doesn't look too serious. Next we will measure the signals flowing through the entire plant and controller.

For a nicer (and printable) version of these plots look in the zipped folder under Plots/Plant_TF_Individuals.pdf

16476   Thu Nov 18 15:16:10 2021 AnchalUpdateGeneralMoved Chiara to 1X7 above nodus powered with same UPS

[Anchal, Paco]

We moved chiara to 1X7 above nodus and powered with same UPS from a battery backed port. The UPS is at 40% load capacity. The nameserver and nfs came back online automatically on boot up.

16475   Thu Nov 18 14:29:01 2021 KojiSummaryBHDBHD invac optics / opto-mechanics

I went through the optics list (in the BHD procurement google spreadsheet) and summarized how to build them.

The red ones are what we need to purchase. Because of the strange height of the LMR mounts, the post needs to have none half-integer inch heights.

They need to be designed as the usual SS posts are not designed to be vac compatible (not because of the material but the design like screw hole venting).

We also need to check how many clean forks we have.
-> The components were ordered except for the custom posts.

ssome partssss
Name Optic Mount Mount OH Post Post OH Fork / Base Base OH Total Height Notes
POP_SM5  Previous POYM1 / 2" Y1-2037-0 LMR2V Thorlabs 1.36 Custom Post 4.14 SS Fork 0 5.5
POP_SM4 New CM254-750-E03 Thorlabs LMR1V Thorlabs 0.87 Newport 9953+PLS-T238 3.88 BA1V / BA2V 0.75 5.5
BSOL1 New 2" VIS BB2-E02 LMR2V Thorlabs 1.36 Custom Post 4.14 SS Fork 0 5.5
ITMYOL1 New 2" VIS BB2-E02 LMR2V Thorlabs 1.36 Custom Post 4.14 SS Fork 0 5.5
ITMYOL2 New 2" VIS BB2-E02 LMR2V Thorlabs 1.36 Custom Post 4.14 SS Fork 0 5.5
SRMOL1 New 2" VIS BB2-E02 LMR2V Thorlabs 1.36 Custom Post 4.14 SS Fork 0 5.5
ASL LA1779-C Thorlabs or KPX217AR.33 Newport LMR2V Thorlabs 1.36 Custom Post 4.14 SS Fork 0 5.5
GRY_SM1 Y2-2037-0 (in hand) DLC   DLC Post   DLC Fork   5.5
BHDBS CVI (In hand) DLC 2 DLC Post   DLC Fork   5.5 (3" post for BHD)
LO3 Lambda (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)
LO4 Lambda (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)
AS3 Lambda (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)
OMC1R3 Y1-1025-45P (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)
OMC1R4 Y1-1025-45P (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)
OMC2R3 Y1-1025-45P (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)
OMC2R4 Y1-1025-45P (in hand) POLARIS-K1-2AH Thorlabs 1 Custom Post 4.5 SS Fork 0 5.5 (3" post for BHD)

OMC1R1 Y1-1025-45P (in hand) LMR1V Thorlabs 0.87 Custom Post 4.63 SS Fork 0 5.5 (3.13" post for BHD)
OMC2R1 NB1-K14 Thorlabs LMR1V Thorlabs 0.87 Custom Post 4.63 SS Fork 0 5.5 (3.13" post for BHD)

16474   Wed Nov 17 17:37:53 2021 AnchalUpdateGeneralPlaced Nodus and fb1 on UPS power

Today I placed nodus and fb1 on UPS battery backed supply. Now power glitches should not hurt our cds system.

ELOG V3.1.3-