I tried to play an .avi file on allegra. In a normal universe this would be easy, but because its linux I was foiled.

The default video player (Totem) doesn't play .avi or .wmv format. The patches for this work in Suse but not Fedora. Kubuntu but not CentOS, etc.I also tried installing Kplayer, Kaffeine, mplayer, xine, Aktion, Realplay, Helix, etc. They all had compatibility issues with various things but usuallylibdvdread or some gstreamer plugin.So I pressed the BIG update button. This has now started and allegra may never recover. The auto update wouldn't work in default mode becauseof the libdvdread and gstreamer-ugly plugins, so I unchecked those boxes. I think we're going to have this problem as long as we used any kind ofadvanced gstreamer stuff for the GigE cameras (which is unavoidable).

Working well tonight: the handoff of CARM to RF (REFL2I), successful reduction of CARM offset to zero, and transition control of MCL path to the OUT1 from the common mode board.  All that's left in lock acquisition is to try and get the common mode bandwidth up and the boost on.

Tonight I have modified transition steps from als to pdh signals. I have added 1:20 filters to CARM_A and DARM_A filter banks to make them unconditionally stable. These filters made locking more robust -- duty cycle is was ~70% tonight. I have also modified slow/ao crossover to avoid ringing up of lines above 1kHz.

Once AO is engaged with high bandwidth, REFL55 signal looks good and I transition PRCL from 165I to 55I. Optical gain compared to PRMI reduced from 55I/165I = -330 down to 55I/165I = 30 in full lock.

I worked on alignment of ETMs. Looking on the cameras I could improve arm power up to 160 and ifo visibility was 80%. POP22 fluctuated by ~50% and every few minutes we loose lock because POP22 almost touches zero.

Jenne and I interviewd Den this afternoon to make the things clear

- His "duty cycle" is not about the lengths of the lock stretch. He saids, the transition success probability is improved.

- For this improvement, the CARM transition procedure was modified to include turning on 1:20 (Z1P20) filter in CARM_A (i.e. ALS) once CARM_B (i.e. RF) dominates the loop in all frequency.

- I think this transition can be summarized like the attachment. At STEP4, the integration of the ALS is reduced. This actually does not change the stability of the servo as the servo stability is determined by the stability of the CARM_B loop. But this does further allow CARM_B to supress the noise. Or in other word, we can remove the noise coming from the CARM_A loop.

- The POP22 issue: Jenne has the trigger signal that is immune to this issue by adding some amount of POPDC for the trigger.
We can avoid the trigger issue by this technique. But if the issue is due to the true optical gain fluctuation, this may mean that the 11MHz optical gain is changing too much. This might be helped by PRC angular feedforward or RF 22MHz QPD at POP.

I analyzed mode scan data from last week.
Mode matching ratio for Y arm is 86.7 +/- 0.3 %. Assuming we can get rid of TEM01/10 by alignment, this can be improved up to ~ 90%.

Peak search, peak fitting and finnesse calculation:
  I made a python script for doing this. It currently lives in /users/yuta/scripts/modescanresults/analyzemodescan.py.
  What it does is as follows

  1. Read mode scan data(coarse5FSRscan.csv, fine1FSRscan.csv). Each column in the data file should be

[time] [some thing like C1:ALS-BEAT(Y|X)_(COARSE|FINE)_(I|Q)_IN1] [C1:LSC-POY11_I_ERR] [C1:LSC-TRY_OUT]

Each separated by comma. Currently, this script uses only TRY, but it reads all anyway

  2. Find peak in TRY data. For the peak search, it splits data in 1 sec and find local maximum. If the local maximum is higher than given threshold, it recognize it as a peak. If two peaks are very close, it uses higher one. This sometimes fails, because mode scan data we have is not so nice.

  3. Fit each peak with Lorentzian function,

TRY = a*b/(4*(t-c)^2+b^2) + d  (a>0, b>0)

  where a/b is a peak height, b is a linewidth (FWHM), c is a peak position in time, and d is a offset.
  I don't like this, but currently, a/b+c is fixed to the maximum value of TRY data used for fitting. This is because sometimes TRY data is so bad and I couldn't get the peak height correctly. Each points of TRY data doesn't have same error because cavity length is fluctuating and relation between cavity length and TRY is not linear. I think I should use some weighting for the fit, but currently, I just use least squares.

  4. Find TEM00 and calculate FSR in "seconds". I just used "seconds" assuming we did a linear sweep. This script recognize TEM00 from the given threshold.

  5. Calculate finesse using FSR and linewidth of the closest TEM00.

  Below are the result plots from this analysis. Calculated finesse looks quite high (~1000). I think this is from non-linearity in the sweep and error in "measured" line width.



Higher order modes and RF sidebands:
  Assuming the curvature of ITMY/ETMY are flat/57.5 m, Y arm length is 38.6 m(FSR 3.9 MHz), positions of HOMs and RF sidebands(11/55 MHz) in frequency domain should look like the plot below.
  The script for calculating this currently lives in /users/yuta/scripts/modescanresults/HOMRFSB.py, inspired by Yoichi's script for KAGRA


Mode-matching ratio:
  By comparing mode scan data and HOM/RF SB positions in a sophisticated way, you can tell which peak is which.



  From COARSE 5FSR measurement, peak heights are

TEM00 0.884, 0.896, 0.917, 0.905, 0.911
TEM01 0.040, 0.037, 0.051, 0.054, 0.062
TEM02 0.083, 0.078, 0.079, 0.071, 0.078
TEM03 0.018, 0.015, 0.013, 0.015, 0.014

  So the mode-matching ratio is

MMR = 86.2 %, 87.3 %, 86.5 %, 86.6 %, 85.5 %

  From FINE 1FSR measurement, peak heights and mode matching ratio is

TEM00 0.921
TEM01 0.031
TEM02 0.078
TEM03 0.014

MMR = 88.2 %

  Assuming each measurement had same error, mode-matching ratio from these 6 values is

MMR = 86.7 +/- 0.3 %  (error in 1 sigma)

  This can be improved by ~5% by alignment because we still see ~5% of TEM01/10. Study in systematic errors on going.

I changed an ADC channel for GCY_ERR and thus recompiled the c1scy model.

Spectra of BS, PRM, ITMX, ITMY are attached with oplevs ON and OFF (in units of urad). Loops reduce RMS from ~2urad to ~0.3urad but phase margin should be increased. REF traces show loop OFF. <-- really?

Note how PRM pitch and yaw spectra are different in the frequency range 0.5 - 7 Hz; yaw is factor of 50 larger then pitch at 2 Hz.

 All 40m laser safety glasses are cleaned and measured this morning.  Bring your own safety glasses if you have to enter the 40m IFO room.

Glasses were washed in 1% Liquinox water solution and  their transmission measured at  165 mW,  2 mm OD beam of 1064 nm

Roughing down the annuloses required closing V1 for 13 minutes

IFO is 2.2e-5 Torr

Annuloses are not pumped for 30 days, since TP2 failed.  IFO pressure 7e-6 Torr it,  Rga 2.6e-6 Torr

Valve configuration: Vacuum Normal as TP3 is the forepump of Maglev, annuloses are not puped at 1.1 Torr

TP3 50K rpm,  0.15A  24C, foreline pressure 16.1 mTorr

Kiwamu and Koji

Last night we have released PRM from the gluing fixture. All of the six magnets are successfully released from the fixture.

We put SRM on the fixuture and glued a side magnet which we had failed at the last gluing.

We let it cure in the Al house. This should be the last magnet gluing until ETMs are delivered.

[Current status]

ITMX (ITMU03): all of magnets/guiderod/standoffs glued, mirror baked; balance to be confirmed
ITMY (ITMU04): all of magnets/guiderod/standoffs glued, balance confirmed, mirror baked
SRM  (SRMU03): magnets/guiderod/standoff glued; a side magnet gluing in process, balance to be confirmed, last stand off to be glued, mirror to be baked
PRM  (SRMU04): magnets/guiderod/standoff glued; balance to be confirmed, last stand off to be glued, mirror to be baked

TT:            magnets/guiderod/standoff glued; balance to be confirmed, last stand off to be glued, mirror to be baked

1. Pianosa doesn't cache the SVN pwds so you need to re-enter at each SVN up or commit. This is different from the rest of our workstations. We need to determine what behavior we want.
2. Tried to use the netGPIB scripts:

pianosa:gpib 0> ./readSR785.csh rb2
rb2
netgpibdata.py: Command not found.

A new offset-zeroing script has been developed and it is ready to run.

 The motivation is to replace the old zeroing script called offset by a better one because this old script somehow failed to revert the gain settings on a given filter bank.

The new script, named offset2, does the same job, but uses tdsavg instead of using ezcaservo. So it doesn't screw up the gain settings.

Additionally the structure of the script is much simpler than the old offset script, and fewer ezca-functions.

I will modify some scripts which use the old offset script so that all the offset-zeroing is done by offset2.

P.S.

Useful scripts are listed on the 40m wiki

http://blue.ligo-wa.caltech.edu:8000/40m/Computers_and_Scripts/All_Scripts

I wanted to restart the c1oaf model. As usual, the first time the model was restarted, it came back online with a 0x2bad error. This isn't even listed in the diagnostics manual as one of the recognized error states (unless there is a typo and they mean 0x2bad when they say 0xbad). The fix that has worked for me is to stop and start the model again, but of course, there is some chance of taking all the vertex FEs down in the process. No permutation of mxstream and daqd process restarts have cleared this error. We need some CDS/RCG support to look into this issue and fix it, it is not reasonable to go through reboots of all the vertex FEs every time we want to make a model change.

Another hysteresis test has begun at 1:50 PT, Dec/19.

It will finish after 3 or 4 hours. During the measurement the PSL mechanical shutter will be kept closed.

The heat drives the ants to the lab!  Make sure the light doors are tight on the chambers.

Do not leave organic trash or food boxes in the 40m to attrack ants !

We observed one or two ants climbing over PMC optics without booties and safety glasses.

The floor was mopped with strong Bayer Home Pest Control solution in the Vertex area.

Do not work inside the 40m lab if you are sensitive to chemicals!

I spotted around 4 within 30 minutes working at the PSL table even after the deathly spray. They seem to be running down from the cables on the oscilloscope rack to the table and the optics.

 The floor is mopped again with strong PEST CONTROL SOLUTION in water in the Vertex area.

Do not plane to work in the IFO-room till noon if you are sensitive to chemicals!

 Ants are back on the PSL table. We'll be mopping the floor with pest control solution tomorrow morning.

So, near 2 of the trashcans in the control room and underneath a desk there are hundrends of ants. Is this normal?

apache has been restarted.
How to: search "apache" on the 40m wiki

CDS changed the suffix for all aquired channel names from _DAQ to _DQ.  When we rebuilt the sus models, described in the previous log, the channel names were changed and consequently the channel files were completely rewritten.

To fix the issue, the latest archived channel file was copied back into the chans directory, and the suffixes were changed, as so:

cd /opt/rtcds/caltech/c1/chans
cp archive/C1SUS_110602_155403.ini  C1SUS.ini
sed -i 's/DAQ/DQ/g' C1SUS.ini

We then restarted the models and the framebuilder.

I ran Zach's arbcav on our SRC with curved TTs and the situation looks much worse than the PRC.

I used the following parameters

SRM: RoC = 142 m, T = 10%
ITM: RoC = 83.1e3 m, T = 1.4%
SRC length: 5.37 m

In this case, with TT RoC of -600, the combined cavity g-factor = 0.9986, and astigmatism from SR3 makes the cavity patently not stable.  You have to go up to an RoC of -710 before the cavity is just over the edge.

This is just a simple rerun of arbcav from #7995 but with the PR2/3 RoCs set to 600, instead of -600.  Overall g-factor = 0.922, and the modes are well separated:

This doesn't take into account the effect of traveling through the substrates (still working on it).  It assumes the PR2/3 have been moved such that the cavity fold lengths remain the same.

This is something that we need to keep in mind: we will need to adjust the position of the PR2/3 to keep the fold lengths the same.

Arbcav with half PRC (flat temporary mirror in front of BS), PR2 RoC = 600m, PR3 RoC = -600m:

NOTE: this does NOT include the affect of the PR2 substrate in the cavity.  Arbcav does not handle that.  It would have to be modified to accept arbitrary ABCD matrices.

NOTE: I added to the mode plot the frequency separation of the first HOMs from the carrier (\omega_{10/01}), in units of carrier FSR.

We need expected finesse and g-factor to compare with mode-scan measurement. Can you give us the g-factor of the half-PRC and what losses did you assumed to calculate the finesse?

Also, flipped PR2 should have RoC of - R_HR * n_sub (minus measured RoC of HR surface multiplied by the substrate refractive index) because of the flipping.
According to Jenne dictionary, HR curvature measured from HR side is;

PRM: -122.1 m
PR2: -706 m
PR3: - 700 m
TM in front of BS: -581 m

Please use these values to calculate expected g-factor so that we don't get confused.

I noticed that Koji used a high reflector for the ITMs for his full PRC arbcav calculation. I just redo it here with the correct ITM transmission and RoC for completeness.

In this case the finesse is 95, instead of 121.

I have calibrated ETMX and ETMY actuators and added a template armSpectra.xml into /users/Templates directory.

Template shows control and error signals of both arms. Procedure is standard: calibrate control to meters and match error based on UGF measurement. XARM UGF: 200 Hz, YARM UGF 210 Hz.

Noise level at high frequencies (>100 Hz) for YARM is 3*10^-15 and is factor of 3 better then for XARM. Servo gains are in the same ratio. I think there is less light on POX than on POY RF PD because I checked phase rotation and analog gain. I assume transimpedances are the same.

ITMX Pitch: 142   microrad/counts
ITMX Yaw:   145   microrad/counts
ITMY Pitch: 257   microrad/counts
ITMY Yaw:   206   microrad/counts

ETMX Pitch: 318   microrad/counts
ETMX Yaw:   291   microrad/counts
ETMY Pitch: 309   microrad/counts
ETMY Yaw:   299   microrad/counts

BS Pitch:    70.9 microrad/counts
BS Yaw:      96.3 microrad/counts

PRM Pitch:   78.5 microrad/counts
PRM Yaw:     79.9 microrad/counts

SRM Pitch:  191   microrad/counts
SRM Yaw:    146   microrad/counts

 Since we don't have an arm length precision measurement (i.e. better than centimeters), why not just do as Koji suggests and use the ALS to get the frequency spacing between a few red FSR and then we have the measurement solid ?

Arm lengths measured using ALS. Both the arms were estimated to have the same length (to the order of a centimeter) 37.51 m

How
I used ALS error signal to lock the arms and scanned the arm to find 4 consecutive IR resonances. From the beat note frequencies measured using the spectrum analyser during IR resonance, the FSR, and hence the length of the arms were calculated.

The estimated lengths are not very precise down to a mm given the resolution of the spectrum analyser. We have brought out the rubidium clock to use as reference for the spectrum analyser to challenge the measurements.

Arm lengths were measured using ALS

X arm length = 37.79 +/- 0.05 m

Y arm length = 37.81 +/- 0.01 m

Whats and whys

We want to measure the arm length with an accuracy of say a mm.

This would mean a measurement precision of 1e-3/40=25ppm. (1mm in 40m)

So the required measurement resolution on the spectrum analyser is 25ppm*4MHz=100Hz (assuming the cavity FSR is roughly 4MHz). 

Although the spectrum analyser does not limit the measurement precision, we are limited by the noise in ALS at 1000Hz rms. So we can use ALS only to measure arm length precise to the order of a few mm.

RXA: Not that we really need to right now, but even with an ALS noise of 1000 Hz, we can can do better just by averaging at each resonance point. And fitting a line as you have already done gets even better:

http://en.wikipedia.org/wiki/Propagation_of_uncertainty

The Spectrum analyser was reference locked to the rubidium clock @10MHz for these measurements.

The FSRs of the arms

X arm = 3.9671e+06 +/- 4.8535e+03 Hz

Y arm = 3.9648e+06 +/- 1.1064e+03 Hz

Attachments:

1&2. Plots representing the arm scans showing the beat frequency for which IR resonates in the arm vs the ALS offset (position of the ETM).

3. Data and code (zip file)

P.S. We had trouble scanning the arms using ALS. This was because the slow servo was not enabled. Hence ALS was losing its PDH lock everytime we scanned past a couple of FSRs.

I'm continuing the arm loss measurements Yuki was making. I'm first familiarizing myself with the procedures for the measurement Johannes describes.

I'm not very familiar with the medm screens, so I'm just kind of poking around and checking with Gautam. I do the following:

1. Turned Xarm ASS dither on, then off.
2. Turned X and Y ALS on, then off shortly after
   1. Realizing I needed some guidance, I found this page on lock acquisition on the wiki
   2. Gautam showed me how to align/lock the IFO so I could take some notes, and we locked the Y arm, misaligned X.
3. I put the PD back in the AS beam path to get the ASDC signal, and approximately centered the beam. This PD is on channel 1 of the scope, which is at 192.168.113.24.
4. I centered the beam onto the MC2 PD that Yuki had installed. This PD is on channel 2 of the scope.
   1. Both scope channels are set to 1V scale (I also had tried 500mV, and it didn't seem to make a difference) and 10s time axis spacing (maximum integration time, since we're looking for a DC effect. Is this what we want?)
   2. The impedance for both channels is 1MOhm.
5. I ran the script to start the loss measurement on the Y arm.
   1. python2 armloss_dcrefl_asdcpd_scope.py 192.168.113.24 1 2 5 YARM
   2. I'm reading ~15 (au?) for the MC channel and ~5% of that out the AS, which seems to make sense to me and looked to be about what Yuki the ratios when I checked the log files. However, I'm a bit confused by the normalization, because the maximum output of the MC PD is 10V, and indeed the scope's display is reading under 10V.

I've left the script running.

After running this script Friday night, i noticed Saturday that the data hadn't saved. Scrolling up inthe terminal, I couldn't see where I'd run the script, so I thought I'd forgotten to run it as I was making last minute changes to the scope settings Friday before leaving.

Monday it turns out I hadn't forgotten to run the script, but the script itself was getting hung up as it waited for ASS to settle, due to the offset on the ETM PIT or YAW setpoints. The script was waiting until both pitch and yaw settled to below 0.7, but yaw was reading ~15; I think this is normal, and it looks like Yuki had solved this problem by waiting for the DEMOD-OFFSET to become small, rather than just the DEMOD signal to be small. Since this is a solved problem, I think I might be using an old script, but I'm pretty sure I'm running the one in Johannes' folder that Yuki is referencing for example here. The scripts in /yutaro_scripts/ have this DEMOD-OFFSET functionality commented out, and anyway those scripts seem to do the 2D loss maps rather than 1D loss measurements.

In the meantime I blocked the beams and ran the script in DARK mode. The script is saving data in /armloss/data/run_20181105/, and runs with no exceptions thrown.

However, when I try to dither align the YARM, I get an error that "this is not a degree of freedom that has an ASS". I'm alsogetting some exceptions from MEDM about unavailable channels. It must have been something about donatella not initializing, because it's working on pianosa. I turned on YARM ADS from pianosa. Monitoring from dataviewer, I see that LSC-TRY_OUT has some spikes to 0.5, but it's mostly staying near 0. I tried returning to the previous frozen outputs, and also stepping around ETMY-[PIT/YAW] from the IFO_ALIGN screen, but didn't see much change in the behavior of LSC-TRY. I missed the other controls Gautam was using to lock before, and I've also made myself unclear on whether ASS is acting only on angular dof, or also on length.

I unblocked the beams after the DARK run was done.

I'm checking out the data this morning, running armloss_AS_calc.py using the parameters Yuki used here.

I made the following changes to scripts (measurement script and calculator script)

• Included the 'hour' of the run in the armloss_dcrefl_* script. This way, we can run more than once a day without overwriting data.
• Changed the calculator script to loop over all iterations of locked/misaligned states, and calculate the loss for adjacent measurements.
  • That is, the measurement script will make a measurement with the arm locked, then with it misaligned, and repeat that N times
  • The calculator now finds the loss for the nth iteration using *_n_locked and *_n_misaligned, and finds N separate loss measurements
  • The dark signal is also computed N times, though all of the dark measurements are made before running the arm scripts, so they could be all integrated together.
  • All of these are saved in the same directory that the data was grabbed from.

I repeated the 'dark' measurements, because I need 20 files to run the script and the measurements before had the window on the scope set larger than the integration time in the script, so it was padded with bad values that were influencing the calculation.

On running the script again, I'm getting negative values for the loss. I removed the beamstops from the PDs, and re-centered the beams on the PDs to repeat the YARM measurements.

Please check your data file and compare with those Johannes made last year. I think the power in your data file may have only three-disits and flactuate about 2%, which brings huge error. (see elog: 40m/14254)

I did the measurement of the arm loss on both X and Y arm by running the armLoss script.

The results will be posted later.

Here are the results of the arm loss measurements, which I have done before the vent.

I ran the existing matlab script, called 'armLoss.m', to estimate the loss. The script resides in /scripts/LSC.

(Y arm)

 Round trip loss =  154.668624 +/- 11.343204 ppm

The figure above is a time series of the measurement.

In the lower plot the power in the ASDC_PD are plotted. The green dotted-curve is the power when the Y arm is unlocked.

The blue dotted-curve is the one when the Y arm is locked.

In the upper plot the estimated loss from each combination of locked/unlocked power are plotted.

(X arm)

Round trip loss = ????? 50 ppm ?????

The obtained time series looked wired because difference in the ASDC power when the arm was locked/unlocked were small.

This small difference results in such a small loss.

To see what was going on I will look at the trend data.

I've plotted TRX, TRY, PD12I and PD11Q.  Arm powers after locking increase for a few tens of minutes, peak out, and then decrease before lock is lost.