Problem:

Diagnostic test tools was starting with errors.

Cause:

After the reboot of the frame builder machine yesterday by Alex, the diagconfd daemon was not getting started by xinetd.  There was a sequence error in the startup where xinetd was being called before mounting drives from linux1.

Important Note:

If you do not see the "nds" line you would not have diagnostic tests enabled in the DTT:

[controls@rosalba apps]$ diag -i | grep nds
nds * * 192.168.113.202 8088 * 192.168.113.202

Solution:

Alex changed /etc/xinetd.d/diagconfd file to point to /opt/apps/gds/bin/diagconfd instead of /opt/apps/bin/diagconf.  He also ensured xinetd started after mounting from linux1.

Alex's Suggestion:
My feeling is we should get rid of this feature and have an NDS address
entry box in the "Online" tab in the DTT with the default "nds". I
mentioned this to Jim Batch and he greed with me, so maybe he is going to
implement this. So maybe you guys want to request the same thing too, send
the request to Rolf and Jim, so we can have the last demon exercised.

The diagnostic script I've written is available in the 'caltech/users/kiwamu/work/20091026_OMC-LSC-diag/src'.

To run the script, you can just execute 'run_OMC_LSC.sh' or just call the m-file ' OMC_LSC_timinig.m'  from matlab.

NOTES:

The script destructs the lock of DARM and OMC, be careful if you are working with IFO.

The input matrix of ITMX has been diagonalized.

The evaluation of this diagonalisation  will be done tonight by freely swinging ITMX again.

(Somehow I couldn't get any data for ETMX from the DAQ channels. I will try it again tonight.)

(details)

For solving the matrix, I used Yuta's python code called inmartixoptimizer.py.

I took the transfer functions of UL->UR, UL->LL and UL->LR as described in this entry.

In the measurement, the frequency bin was set to 0.001 Hz and the data were 50 times averaged on dtt.

Here is the new input matrix.

[[ 0.87059649  1.14491977  1.07992057  0.90456317]
 [ 0.64313916  0.55555661 -1.44997325 -1.35133098]
 [ 1.13979571 -1.19186285 -0.89606597  0.77227546]]

This matrix should give a better performance than before.

I went inside to align the beam on WFS and noticed that oscillations in yaw are ~10 times stronger then in pitch. I've plot rms of pitch and yaw measured by LID sensors and saw that MC3 yaw rms motion is a few times larger then pitch.

Also MC1 input diag matrix does not diagonalize signals for pitch and yaw. In the spectrums of these signals all 4 resonance are equally seen during the free swinging. I think we should rediagonalize MC1.

Another thing is that if MC1 and MC3 are on the same stack,  pitch and yaw spectrums of these mirrors should be comparable. But MC1 signal is ~2-3 times larger then of MC3. I think we should correct calibration.

What are "LID" sensors?  Do you mean the OSEM shadow sensors?  I'm pretty sure that's what you meant, but I'm curious what "LID" means.

I've written MC123 input matrixes to the front-end.

MC1 diagonalization is poor, better then before, but still pitch is seen in pos and yaw. Either smth is malfunctioning or flags touch sensors and do not move freely. On the plot mc1_new black lines - before, red - after rediagonalization.

 I've actuated on MC1 with UL, UR, LR, LL coils in turn and measured sensor readings. All coils separately work fine from the first look.

On the plot: black - free mirror, blue - UL coil actuation, green - UR, grey - LR, red - LL.

 I've manually corrected MC1 input matrix by looking at UL, UR, LL, LR transfer functions between each other. This improved pos significantly and slightly yaw.

Motivation:
  MC is aligned from the A2L measurement, but to do the beam centering more precisely, we need coils to be balanced.
  There are several ways to balance the coils, like using oplev or WFS QPD RF channels.
  But oplev takes time to setup, especially for MC3. Also, c1ioo WFS channels were newly setup and haven't been checked yet.
  So, I decided to use OSEM sensors.
  An OSEM sensor itself is sensitive to every DOF of an optic motion, but we can diagonalize them using 4 OSEM sensors and proper input matrix.

Method:
  1. Measure transfer functions between
     ULSEN and URSEN (H_UR(f))
     ULSEN and LRSEN (H_LR(f))
     ULSEN and LLSEN (H_LL(f))

  2. Make a matrix A.

    A =  [[ 1           1           1          ]
          [ H_UR(f_pos) H_UR(f_pit) H_UR(f_yaw)]
          [ H_LR(f_pos) H_LR(f_pit) H_LR(f_yaw)]
          [ H_LL(f_pos) H_LL(f_pit) H_LL(f_yaw)]]

    where f_dof are resonant frequencies.

  3. A is

    s = Ad

   where vectors s^T=[ULSEN URSEN LRSEN LLSEN] and d^T=[POS PIT YAW].
   So,

    d = Bs = (A^TA)^(-1)A^Ts

   where A^T is transpose of A.

   B is the input matrix that diagonalizes 3 DOFs.

What I did:
  1. Measured the TFs using diaggui and exported as ASCII.

  2. Made a script that reads that TF file, calculates and sets a new input matrix B.
    /cvs/cds/caltech/users/yuta/scripts/inputmatrixoptimizer.py
   You need to set resonant frequencies to use the script.

   New input matrices for MCs are;

C1:SUS-MC1_INMATRIX
[[ 1.17649712  0.94315611  0.85065054  1.02969624]
 [ 0.55939288  1.28066594 -0.85235358 -1.3075876 ]
 [ 1.23467139 -0.74521928 -1.29394051  0.72616882]]

C1:SUS-MC2_INMATRIX
[[ 1.12630748  1.01451545  0.9013457   0.95783137]
 [ 1.03043025  0.67826036 -1.37270598 -0.91860341]
 [ 0.83546271 -1.26311029 -0.6456881   1.2557389 ]]

C1:SUS-MC3_INMATRIX
[[ 1.18212117  1.26419447  0.77744155  0.77624281]
 [ 0.79344415  0.84959646 -1.10946339 -1.247496  ]
 [ 1.00225331 -0.84807863 -1.21772132  0.93194674]]

  I ignored SIDE this time.

Result:
  Spectra of each SUSDOF_IN1_DAQ before diagonalization (INMATRIX elements all 1 or -1) were


  After diagonalization, spectra are


  As you can see, each SUSDOF has only single peak (and SIDE peak) after the diagonalization.
  SUSSIDE still has 4 peaks because SIDE is not included this time.

  For MC2, POS to SUSPIT and POS to SUSYAW got worse. I have to look into them.

Effect of resonant frequency drift:
  As you can compare and see from the spectra above, resonant frequencies of MC1 are somehow drifted(~0.5%) from Nov 9 to Nov 13.
  If resonant frequency you expected was wrong, calculated input matrix will be also wrong.
  The effect of 0.5% drift and wrong input matrix can be seen from this spectra. DOFs are not clearly separated.
 

Plan:
 - learn how to use diaggui from command line and fully automate this process
 - balance the coils using these diagonalized SUSPOS, SUSPIT, SUSYAW

There has been some input matrix diagonalization in the past by Yuta and Kiwamu, but I find the automation to be not totally satisfactory.

It would be better if we could automatically fit the data to find the Suspended optic eigenfrequencies and then use that to get the matrix. So I wrote a peak fitter to get the matrix.

It gets the data from mafalda with NDS2, then it makes the PSDs, and then starts with some initial guesses (based on looking at the plots) and them uses fminsearch to get the peak frequencies and Q's.

Using the output of this, we can use Yuta's method and take the passive transfer functions with the free swing data (from April 30, so we got do do it quick) to get the input matrix.

Doing the SUS input matrix is nice for having good damping (as long as we remember to include SIDE), but my motivation is to produce a good null stream from the 4 face sensors so that we can estimate the sensor noises at all times.

 I've moved all of my SOS peak fitting stuff into the scripts area so that Leo can make it better:

/cvs/cds/rtcds/caltech/c1/scripts/SUS/peakFit

findPeaks.m gets the data and makes the fitted spectra that I put in the previous entry.

findMatrix.m is the barely started script that ought to take the TF data and output the matrix to the MEDM screen.

I have updated the scripts/SUS/peakFit/ directory so that it now finds the SUS input matrix coefficients in addition to just finding the free-swinging peaks.

Procedure:

1. Get OSEM sensors data via NDS2 from a time when the optics have been kicked and then left free swinging.
2. Downsample the data to 64 Hz and save.
3. Make power spectra with a 1 mHz resolution (i.e. we need a few hours of data) and  ~10 averages.
4. use the fminsearch lorentzian peak fitter -> save the peak frequencies
5. Make Transfer Function estimate matrix at the peak frequencies between all OSEMs (this makes a 5x4 complex matrix)
6. The matrix should be real, so make sure its mostly real and then take the real part only
7. Normalize (height of biggest peak for each f_DOF should be 1)
8. Add a Butterfly mode vector. This makes the sensing matrix go from 5x4 to 5x5. (Butterfly a.k.a. Pringle)
9. Invert
10. Normalize so that the biggest element in each Sensor2DOF column is 1.
11. Load values into MEDM screen and then verify by another free swinging data run.

The attached PDF shows how much rejection of the unwanted DOFs we get between the existing diagonal input matrix and this new empirical matrix. Previously, the decoupling was only a factor of a few for some of the modes. Now the decoupling is more like orders of magnitude (at least according to this calculation). It will be worse when we load it and then try another free swinging run. However, the fact that the suppression can be this good means that the variation in the coefficients at the ~hours time scale is at least this small (~< 0.1%)

That's the basic procedure, but there are a lot of important but mainly technical details:

1. Free swinging data must be taken with the angle bias ON. Otherwise, we are not measuring the correct sensing gain (i.e. the magnets are not in their nominal place within the OSEM-LED beam)
2. Data must be checked so that the shadow sensor outputs are in their linear regime: if they are exploring the cubic part, then the fundamental is being suppressed.
3. Instead of just using the peak frequency, I average a few points around the peak to get better SNR before inversion. I think this will make the results more stable.
4. All previous input matrix diagonalization efforts (Buckley, Sakata & Kawamura, Black, Barton, Gonzalez, Adhikari & Lawrence, Saulson,...) for the past ~15 years have been using the spectra's peak height data. Today's technique uses the TF and so is more precise. The coherent transfer function is always better than just using the magnitude data.
5. This method is now fairly automatic - there's no human intervention in fudging values, choosing peak heights, frequencies, etc.
6. We'll have to rerun this, of course, after the mirrors are aligned and after the OSEM whitening fiasco is cleaned up somewhat.

I'll set the optics to be aligned and then swing tonight.

I've finally completed the SUS/peakFit/ scripts which find the new input matrix for the SUS. MC1, MC2, MC3, and ITMX have been matrix'd.

I tried to do the BS, but it came out with very funny matrix elements. Also the BS is missing its DAQ channels again (JAMIE !) so we can't diagnose it with the free swinging method.

To continue, we have to get some good data and try this again. Right now there are some weird issues with a lot of the optics. I've also set the damping gains for the optics with the new matrices.

 Ex.

new_matrix = findMatrix('ITMX')

writeSUSinmat('ITMX', new_matrix)

writeSUSinmat.m

this script writes the values to the MEDM input SUS matrix. To do the writing, I used the low level 'caput' command instead of ezcawrite since the ezca libraries are getting deprecated.

caput doesn't really have good diagnostics, so I use matlab to check the return status and then display to the terminal. You can just rerun it if it gives you an error.

A coupled of normalization notes:

1) The POS/PIT/YAW rows are scaled so that the mean of abs(FACE elements) = 1. Previously, I had the max element = 1.

2) The SIDE row is scaled so that the SIDE element = +1.

3) I then normalized the ROWS according to the geometrical factors that Jamie has calculated and almost put into the elog.

All these scripts have been added to the SVN. I've removed the large binary data files from the directory though. You can just rsync them in to your laptop if you want to run this stuff remotely.

Besides the purpose of correctly tuning the suspensions, my hidden goal in the input matrix diagonalization has been to figure out what the 'true' sensing noise of the OSEMs is so that we can accurately predict the noise impact on the OAF.

The attached plot shows the DOFs of ITMX calibrated into microns or microrad as per Jamie's ethereal input matrix calculations.

The main result is in the ratio of POS to BUTTER. It tells us that even at nighttime (when this data was taken) we should be able to get some reduction in the arms at 1 Hz.

Whether we can get anything down to 0.1 Hz depends on how the arm control signal compares to the POS signal here. I leave it to Jenne to overlay those traces using a recent Arm lock.

Dichroic mirror quotes are in the wiki.

ATF is pricy.

We got a good price from Laseroptik, but the wedges are 5 arcminutes. The fused silica grade is 0F, meaning the  homogeneity is 5 ppm instead of 1ppm.  I requested an other large wedge quote on the substrates.We may have to get substrates from somebody else and ship it to Germany

MLT quote is outrageously high

REO is not interested in this low volume job.

 The Laseroptik quote is here.The 2 degrees wedge cost is $40 on each optics!  See wiki

 Instead of doing RCG stuff, I went to Millikan to work on data analysis as I couldn't stand the fumes from the construction.  (this morning, 8am) 

 Diesel fumes are pumped away from control room AC intakes  with the help of newly installed  reflector boxes on the CES wall  fans.........see # 2272

Noise below 10 Hz became larger again compared with the data before (see here #4352)

Note that the Y-axis is in MHz.

Today we tried the Schmitt trigger DFD, and while it works it does not improve the noise performance.  At least part of our problem is coming from the discrete nature of our DFD algorithm, so I would propose that an industrious day job person codes up a new DFD which avoids switching.  We can probably do this by mixing the input signal (after high-passing) with a time-delayed copy of itself... as we do now, but without the comparator.  This has the disadvantage of giving an amplitude dependent output, but since we are working in the digital land we can DIVIDE.  If we mix the signal with itself (without delay) to get a rectified version, and low-pass it a little, we can use this for normalization.  The net result should be something like:

output = LP2[ s(t) * s(t - dt) / LP1[ s(t) * s(t) ]],

where s(t) is the high-passed input and LP is a low-pass filter.  Remember not to divide by zero.

I modified the c1gfd.mdl simulink model.  I made a backup as c1gfd_20110325.mdl. 

The first change was to use a top_names block to put everything in.  The block is labeled ALS.  So all the channels will now be C1:ALS-GFD_SOMETHING.  This means medm channel names will need to be updated.  Also, the filter modules need to be updated in foton because of this.

I then proceeded to add the suggested changes made by Matt.  To avoid a divide by zero case, I added a saturation part which saturates at 1e-9 (note this is positive) and 1e9.

I looked once again to the local OSEM sensors and MC length signals. Then I replaced 1e-20 to 1e-50 in the if-statement of the iir_filter function. Here I report about the difference of the signals in question.

First we look at the MC2 OSEM local sensor. In the figure below the psd of the signal is presented in three cases - with a free MC2 mirror without feedback, with a feedback signal and with a feedback signal with corrected if-statement. We can see that without FB the wire resonances are high and dumped when OSEMs are on. However we can see that below 1 Hz the psd of the sensor signal with 1e-20 in the if-statement is higher then psd of the sensor signal from free mirror. FB with 1-50 in the if-statement fixes this problem. 

If we take a look on the plot of the coherence between GUR1_X and SENSOR signals we can see that coherence is corrupted when 1e-20 is used in the is-statement and is good when 1e-50 is used.

 Next we look at the psd of the MC length. We can see how strongly these curves diverge below 1 Hz. The MC_F signal was also corrupted at higher frequencies.

 The coherence between MC_F and GUR1_X is also improved.

I made a LISO fit of the measured TF of the daughter board, so that I can digitally invert the daughter board whitening. Results attached. (Inverse) Filters have been uploaded to the ALS X Foton filter banks.

Channels C1:X05-MADC0_TP_XXX with XXX 2-9, 14-19, 21-27, 29-31 showed digital zeros.

Some of these channels are used in c1scy.mdl, e.g. for OSEM stuff. I guess this is not optimal...

 Last night I found that there were many dust particles on the second steering mirror in the REFL path on the AS table.

Looking at it through an IR viewer, I saw the REFL beam hitting one of the biggest dust particles on that mirror.

This dust particle maybe causing the glitches or maybe not.

Anyway because it's always better to have clean mirrors, I will wipe the steering mirror in this evening and check the presence of the glitches again.

 {Rana and Kiwamu}

Yesterday we disabled the sendmail daemon and the isdn daemon on allegra because we don't need these daemons always running.

-   How to disable/enable daemons:

sudo ntsysv

While investigating the BIO situation with the LSC machine and the iscaux2 processor last night, we wondered if maybe the Anti-Aliasing filters were mistakenly disabled. But why do we need these anyway?

Our ADCs digitize at 64 kHz and there is a digital lowpass in the IOP at 5 kHz before we downsample to 16 kHz. So mainly we're trying to prevent some aliasing at the 64 kHz IOP rate. But our analog AA filter is a 8th order ELP at 7570 Hz, so its overkill.

So, I propose that we bypas the AA via hardwiring the board and implement a 10 kHz pole in the whitening board (D990694) before the whitening by turning R127, etc. into a 0.1 uF cap. Along with the 100 Ohm series resistor, this will make a pole at ~15 kHz. Probably ought to check that the input resistor is metal film. Also, if we replace C158/C159, etc. with a 0.47 nF cap, we'll get 2 poles at 35 kHz to limit the higher frequencies from saturating.

Bought dish soap and scrub sponges today and placed them under the sink with the other dish supplies.

Backed up all the wikis. Theyr'e in wiki_backups/*.tar.xz (because xz -9e gives better compression than gzip or bzip2)

Moved old user directories in the /users/OLD/

I disabled the feedback to the PZT1 PITCH in the Y arm dithering scripts so that it won't push the beam away from the good point.

Currently one has to do a manual alignment only for the PZT PITCH but the rest of DOFs are still able to be automatically aligned with the script.

We made a model for the dither angular stabilization system c1ass.mdl. The attached file shows the diagram.

The idea is to dither a combination of 6 optics (ETMs, ITMs, PZTs) at different frequencies and demodulate three PDs (TRX, TRY, REFL11I). Then form the DOFs from demodulted signals, filter, and send each DOF to a combination of optics.

This is enough to get started with arm cavities alignment (we may need to add the BS for the Y arm). More optics and PD can be added as they become available and/or needed.

The DAC for the fast PZT  are not connected and have to be commissioned.

I modified our existing c1ass model to include alignment of input steering TT1 and TT2 for YARM and BS for XARM. Corresponding medm screens are also created.

Dithering:

ETM_PIT: frequency = 6 Hz, amplitude = 100 cnts
ETM_YAW: 8 Hz, 400 cnts
ITM_PIT: 11 Hz, 800 cnts
ITM_YAW: 14 Hz, 1200 cnts

These values were chosen by looking at cavity transmission and length signals - excitation peaks should be high enough but do not shake the optics too much.

Demodulation:

LO for each degree of freedom is mixed with cavity length and transmission signals that are first bandpassed at LO frequency. After mixing low-pass filter is applied. Phase rotation is chosen to minimize Q component

Sensing matrix:

8 * 8 matrix was measured by providing excitation at 0.03 Hz to optics and measuring the response in the demodulated signals. Excitation amplitude was different for each optics to create cavity transmission fluctuations of 25%

Though coherence was > 0.95 during the measurement for each element (except for TT -> Length signals), after inverting and putting it to control servo, loops started to fight each other. So I decided to try a simple diagonal matrix:

TT1_PIT -> ETM_PIT_TRANS, TT1_YAW -> ETM_YAW_TRANS, TT2_PIT -> ITM_PIT_TRANS, TT2_YAW -> ITM_YAW_TRANS,

ITM_PIT -> ETM_PIT_LENGTH, ITM_YAW -> ETM_YAW_LENGTH, ETM_PIT -> ITM_PIT_LENGTH, ETM_YAW -> ITM_YAW_LENGTH

And this matrix worked much better.

Control loops:

8 loops are running at the same time. UGF for input steering loops is 20 mHz, for cavity axis loops - 80 mHz. Slower loop is stronger at low frequencies so that cavity axis servo follows input steering alignment.

Results:

When I started experiment the cavity was misaligned, transmission was ~0.4. Servo was able to align the cavity in ~30 seconds. This time depends on mirrors misalignment as well as input optics and cavity axis misalignment relative to each other.

When servo converged I disturbed ETMY, ITMY, TT1 and TT2. Servo was able to compensate for this.

Excitation lines seen by transmission and length of the cavity are suppressed as shown on the attached as pdf figures.

Note:

Though the servo is able to align the cavity during my tests, this does not mean it will work perfectly any time. So please, if you lock, try to use the servo for alignment. If something goes wrong we'll fix it. This is better then to align IFO by hands every time.

I've put 4 scripts into ASS directory for YARM alignment. They should be called from !Scripts YARM button on c1ass main medm screen.

Scripts configure the servo to align the cavity and then save computed offsets. If everything goes right, no tuning of the servo is needed.

Call TRANS MON script to monitor YARM transmission, then "ON" script for aligning the cavity, then "SAVE OFFSETS" and "OFF" for turning the servo off.

ON script:

• sets demodulation gains that I used during OL measuments
• sets LO oscillator frequency and amplitude for each optic
• sets demodulation phase rotation
• sets sensing matrix
• sets servo gains for each degree of freedom
• sets up limits for servo outputs
• gently increases the common gain from 0 to 1

SAVE OFFSETS script:

• holds servo outputs
• sets servo common gain to 0 and clears outputs
• reads old optics DC offsets
• computes new DC offsets
• writes new offsets to C1:SUS-OPTIC_ANGLE_OFFSET channel
• holds off servo outputs

OFF script:

• sets LO amplitudes to 0
• blocks servo outputs

Notes:

SAVE OFFSET script writes DC offsets to C1:OPTIC_ANGLE_OFFSET channel, not to _COMM channel!

LIMITS are set to 500 for cavity axis degrees of freedom and to 0.5 for input steering. Usually servo outputs is ~30% if these numbers. But if something goes wrong, check this for saturation.

DC offsets of all 8 degrees of freedom are written one by one but the whole offset of put at the same time. This works fine so far, but we might change it to ezcastep in future.

I wanted to center beams on the XARM cavity mirrors using c1ass model. I've run XARM setup script and then turned dithering on. Cavity went out of lock because calculated offsets were incorrect.

I was using TRX only and calculated rotation phases for ITM and ETM pitch and yaw. For this I've added a low pass filter into Q-quadrature bank and made DC value at the output to be zero by adjusting the phase. I've put gains (+1 or -1)  in the I quadrature such that output was positive.

Then I've set the sensing matrix to identity as I decided to deal with separate loops. Of coarse, they are mixed by the cavity, but at least in the control system they are distinguished. Old matrix summed error signals in one degree of freedom from both mirrors. This makes more sense but still not precise because coils are not ideally diagonalized.

Then I've adjusted gains for control loop for every degree of freedom. I've ended up with (0.1; 0.1; 0.1; -0.1). I did not use large gains as I wanted slow convergence because of the demodulation low-pass filter time response constant of 20 sec. Coupling (I quadrature) was reduced from (0.9, 0.3, 2.4, 1.2) to zeros (0-0.1) in ~5 minutes, TRX increased from 0.73 to 0.90.

There is one thing that I do not understand yet. I think controllers should minimize angle -> length coupling that is proportional to I-quadrature if phase is correct. But phase depends on alignment and when the feedback loops are on, phase drifts. I could see it during my measurement. But I did not find any script that smoothly tunes phase such that coupling is all in I-quadrature. I guess this is not hard to set a gradient descent algorithm that minimizes DC value of Q-quadrature. Or how this is usually done?

The first column is distance from the front face of the laser in cm.
The second column is beam diameter in the horizontal direction in microns.
The third column is the beam diameter in the vertical direction in microns.

The SP table was found open this morning. Please, do not make optics dirty!

I cleaned up the tops of the SP table.

Stop storing your junk, boxes, laptops, etc. on the optical tables. This includes the big SP table. Please move all of that junk into racks or shelves, etc.

ETMX optical table was left open. Burned toast award goes to ......?

 me

Type: How not to

Please do not leave optical tables open! You will be held responsible for creating dirty optics.

Specially heavy items: old analoge scope or hardware loaded boxes......etc

 The table cover section holding crossbars are not evenly spaced.

You have to center each cover section on the cross bar so it is supported on both sides !

I will clean up on this table tomorrow

Yesterday evening I took nearly all of the masks, gloves, gowns, alcohol wipes, hats, and shoe covers. These were the ones in the cleanroom cabinets at the east end of the Y-arm, as well as the many boxes under the yarm near those cabinets.

This photo album shows the stuff, plus some other random photos I took around the same time (6-7 PM) of the state of parts of the lab.

Control room to outside door was realigned.

                                                   It is self closing now.

Control room to IFO door lock optimized to soft closing.

All other doors lubricated by Alex of the key shop.

[steve,ericq,gautam]

We re-checked IMC locking, arm alignments (we were able to lock and dither align both arms today, and also made the michelson spot look reasonable on the camera) and made sure that the AS and REFL spots were in the camera ballpark. We then proceeded to remove the heavy doors off ITMY and BS/PRM chambers. We also quickly made sure that it is possible to remove the side door of the OMC chamber with the current crane configuration, but have left it on for now.

The hunt for clipping now begins.

Here for the downsampling process we use a low-pass Bessel digital filter of order 6, normalized cut-off frequency = 0.1. In the plot presented below we compare the results with downsampling ratio = 1, 2, 4.

We can see that increasing the downsampling ratio, we increase the error of the filter. Moreover, the error at some particular frequency f seems to depend on the ratio f/Fs, where Fs - sampling frequency (2048 Hz) devided by downsampling ratio. Error is the same for all curves below 1 Hz but then begins to increase as we increase the sownsampling ratio. In order to figure out what the problem is - mistake in the filter code, inaccurate upsample algorithm or this is NLMS particularity, I've changed sampling frequency in the chans/daq/C1PEM.ini and C1IOO.ini files from 2048 Hz to 512 Hz for corresponding channels. Now, we compare the error from the filter working with 2048 Hz frequency, downsampling ratio = 4, low-pass filter = Bessel of order 6, normalized cut-off = 0.1 and filter working with 512 Hz sampling frequency, without downsampling and with corresponding Bessel low-pass filter with normalized frequency 0.4.

MC_F measurement at 2048 Hz was done during the day, for that reason red curved is slightly higher then green in the resonance frequencies. But still we can see that these two cases are very much alike. For this reason, it seems that NLMS filter works better with higher sampling frequencies.

All PEM and IOO DQ channels disappeared. These channels were commented in C1???.ini files though I've uncommented them a few weeks ago. It happened after these models were rebuild, C1???.ini files also changed. Why?

I added the channels back. mx_stream died on c1sus after I pressed DAQ reload on medm screen. For IOO model it is even worse. After pressing DAQ Reload for C1IOO model DACQ process dies on the FB and IOO machine suspends.

I rebooted IOO, restarted models and fb. Models work now, but there might be an easier way to add channels without rebooting machines and demons.

manufacturer info

Jamie and Shuresh moved in Jenne's 11 drawers cabinet and relocated old note book boxes on the inside of the vac tube.

 Barring other chores for next Wednesday, we're going to spend Wednesday afternoon populating the new cabinet with all of the optics hardware: posts, forks, dogs, everything!  It's going to be so organized and awesome!!

 Early morning drift in pitch. This plot is meaningless because there is no real  light on IP-Ang

The beam is clipping on the pick off mirror at ETMY chamber. The beam is half beam size too high. Yaw is perfect

http://www.ligo.caltech.edu/~misi/summary/day/20150504/

I left the arms locked last night. Looks like the drift in the Y arm power is related to the Y arm control signal being much bigger than X.

Why would it be that Y > X  ?

Why is the Y arm drifting so much?

The " PSL FSS Slow Actuator Adjust " was brought back to range from 1.5 to 0.3 yesterday as ususual. Nothing else was touched.

I'm not sure if the timing scale is working correctly on theses summery plots. What is the definition of today?

The y-arm became much better as I noticed it at 5pm