Several updates on matlab scripts modifications I recently made during data processing

1- generate_QUAD_Model_Production.m. Produces a matlab model for damped and undamped quad suspension  (/ligo/svncommon/SusSVN/sus/trunk/QUAD/Common/MatlabTools/QuadModel_Production)

The filter file used in the script to compute the modeled closed loop was out of date from the filters installed on our quads, so I updated it with the most recent damping filters : "dampingfilters_QUAD_2013-06-14.mat" living in /ligo/svncommon/SusSVN/sus/trunk/QUAD/Common/FilterDesign/MatFiles

design_damping_QUAD_20130614.m (that creates dampingfilters_QUAD_2013-06-14.mat) won't work if susTag is set to "L1:SUS-ETMY" since there's no 'etmy' variable in the ISI2SUS_projection_file.mat

2-plotquad_matlabtfs.m. Processes the measurements and plots them against our model (/ligo/svncommon/SusSVN/sus/trunk/QUAD/Common/MatlabTools)

The option of plotting the measurements with the damping active against the closed loop model has been added to the script, with a new flag "isDamped=true/false" at the beginning of the script

3-Matlab_TFs.m. Launches automated transfer functions (/ligo/svncommon/SusSVN/sus/trunk/Common/MatlabTools/SchroederPhaseTools)

The script has been modified in order to take transfer functions on several suspensions at the same time. The drive signal txt file used by awgstream to send an excitation is now suspension specific, as well as the diary file. Those two files were creating conflicts and error messages when two similar tf scripts were beeing run at the same time.

Also, in order to make the script easier to read, and for future debugging purposes, I moved the sensor and actuator channel list out of the code, and packed it in a subfunction "generate_channel_list.m" located in /ligo/svncommon/SusSVN/sus/trunk/Common/MatlabTools/SchroederPhaseTools/

4-collect_schroeder_TF_v9DQ_sma.m. Subfunction of Matlab_TFs.m. It creates a schroeder phased drive signal

From now, when Matlab_TFs.m is ran, the one time excitation signal, as well as the full drive signal will be displayed and plotted, as a function of time before the signal is sent to the DAC.

Yesterday this IOC lost communication to the dust monitor in the diode room. I restarted it today and it has reconnected. I suspect the Comtrol serial to Ethernet converter (10.105.0.100) is on the network switch that was remounted yesterday.

The SATABOY QFS disk system used by h1fw1 and h1nds1 suffered a failure at 23:31 last night local time. I have turned off h1fw1 while we investigate the problem.

Please use h1nds0 for archived data until this problem has been resolved. h1nds1 will only support real time data.

Attached are plots of dust counts requested from 4 PM August 12 to 4 PM August 13.

Attached are plots of dust counts requested from 4 PM August 11 to 4 PM August 12.

Out of the 4 DoF (arm optical axis pos and angle, green input beam pos and angle, 2x for pitch and yaw not counted), 3 DoF will likely need a frequent initial alignment - the input beam position at the ETM is the one DoF that remains fixed.

Below is a short description of the suggested alignment loops. Attached is a draft RTCDS model block for e.g. X-arm pitch (we need 4 of them for X/Y P/Y).

Signals to be demodulated:
L: cavity PDH length error signal (green)
R: cavity reflected power (green)
T: cavity transmitted power (green)

Actuation basis (in rad):
(ITM, ETM, TMS)

g-factors:
g1=1-L/R1;
R1= RoC of ITM
g2=1-L/R2;
R2= RoC of ETM

Dither drive frequencies:
f1: ITM dither frequency
fT: TMS dither frequency

Control loops:
1) Fix ITM actuation node
L --->x ---> CTRL_FILT -----> (g1,1,1)
      ^
      |
      f1

Remark:
- Actuation at 1.65-2Hz or above 4Hz seems possible in pitch
- Pitch centering offset at 1.7Hz: ~2.3mm (slightly freq. dependent)
- Pitch centering offset above 4Hz: ~4.7mm
- Could provide a slight additional length drive to compensate this offset
- no such issues for yaw
- Actuation leaves the following quantities unchanged:
  - arm axis position on ETM
  - angle between arm axis and input beam


2) Input beam to arm axis angle overlap at ETM focal point
R or T ---> x ---> CTRL_FILT -----> (1,g2,g2)
            ^
            |
            f1
Remark:
- Actuation leaves the following quantities unchanged:
  - arm axis position on ITM
  - arm axis to input beam separation at ETM focal point
  - input beam position at ETM (trivial)

3) Input beam to arm axis angle overlap at ETM
R or T ----> x ---> CTRL_FILT ----->  (0,0,1)
             ^
             |
             fT
Remark:
- Actuation leaves the following quantities unchanged:
  - arm axis position and angle
  - input beam position at ETM (trivial)

Alternatives:
- GigE camera centering on ITM (instead of 1)
- ETM dither possible, but due to the steep drop-off of the Quad actuation, I expect reasonable drive amplitudes (producing a decent dip in power buildup) to only be possible at very low frequencies. As a result, since the dither frequency separation dictates the achievable bandwidth, it might be difficult to run all dither loops at the same time.

General Remarks:
- The dither drive is done in the optics basis, because any other basis would require significant attention to the relative phase of different actuation functions.
- To be able to run the loops in sequence, an orthogonal actuation function is desirable.
- repeat procedure of X- and Y -arm, pitch and yaw
- Once both arms are aligned, turn this feed-back off.

Work for today:
- Kiwamu and Karen, PSL enclosure/room cleaning, and PSL switched from Science to commisioning.
- Hugh, cable strain relieve at HAM2 and HAM3.
- Richard, end station cleanroom work, electrical.
- Jodi and Bubba, end station work.
- Hugh, to HAM6 work.
- David, DAQ reboot.
- Cyrus, network switch re-rack.

Deliveries/vendors:
- Praxair, LN2 delivery to CP1, level now at 93.2.
- Unifirst, change out mats.

Set the HAM6 Optical Table to -220.1mm Gz per E1000403 & D1100398.  Used Mon 802.

The table started about 1 to 1.5mm low.  Finished at average -220.3 +-0.2mm; still need to do horizontal alignment so this is close enough on the vertical for now.

Went through all the Feed Thrus on these chambers and added strain relief.  This consists of the rod (bolt, all-thread) & zip tie.  I also fussed with cable routing some to get cables off the HEPI Crossbeams (Floating.)  I had to power down the HAM2 CPS racks for a time.
Couple photo attached if you care to see what I'm talking about.

The workmen from Craftsman arrived this morning to install the balance of the flooring at Y End. However, there was not enough flooring on-site to complete the task. (Yes, we had checked before hand but the guy who really knew the answer was on vacation...) The Craftsman work crew was called off and additional flooring ordered. The ETA for the flooring to Seattle is 23 Aug but Craftsman will expedite if possible. 

In the meantime, we will work on populating the chamber cleanroom, assembling the work platform and modifying the CPB cleanroom.

Created H0EDCU_TIDAL.ini, added to DAQ

Restarted all DAQ systems newly built against RCG2.7.1. Added all DCU systems back into DAQ exceeding initial 85 DCU limit. Verified that dataValid flag in frames are correct (unity no longer written). Patrick verified no corruption of DUST channels.

New broadcaster channel list installed (list from JZ).

Create first version of h1pemex model (copy of h1pemey). Built model against RCG2.7.1 (all other models still at 2.7). Started model and added to DAQ.

H1EDCU_DAQ.ini expanded to include new DAQ PV channels. Created DC0 Front end channel list from rtsystab list.

I enabled port monitoring for port 0/42 (sush34) on sw-msr-h1fe, mirroring traffic to port 0/6 for troubleshooting purposes.  It is not known what side effects can occur on the Netgear switches in this configuration.

Built new daqd and nds from tags/advLigoRTS-2.7.1, installed on x1dc0, x1nds1, x1fw1, then restarted the daq system.

I re-racked sw-lvea-aux0 in the EE room today, as specified under WP 4086.  It should be in it's now (hopefully) final location, in RU34 below the patch panel feeding the LVEA.

Here is the current status of the aLIGO SEI work regarding the HAM chambers. Everything in red is a change from last week's status. Everything in green is available.

The PSL was switched to the commissioning mode this morning for approximately one hour for the weekly cleaning. Karen and I entered the PSL enclosure. After the cleaning the PSL was then switched back to the science mode again.

Note: when we followed Robert's instruction to switch it to the commissioning mode the HEPAs in the laser room didn't come up for some reason. We could notice it because the laser room was too quiet. This might be some sort of bug. Toggling the HEPAs off and on on the touch panel eventually made it running.

Beginning measurements of IMC Crossovers. 

Not sure why, but I found the IMC locked, but with the top (M1) stage turned OFF, gain sert to zero, with FMs 1 and 2 on. In order to engage M1, I
- turned off FM1 (the boost), turned on inputs and outputs
- used a TRAMP of 20, and ramped to a gain of 1.0 (probably didn't need to go that slow, but I didn't want to blow the lock)
- Turned on the boost filter.

Also found FM9 (CLP300) of M3 LOCK OFF, so I turned that ON.

The LOCK Gains are now [M1,M2,M3] = [+1,+0.2,-300].

Finally, all of the notch filters (FMs 6-10 in M3_ISCINF Bank) were OFF. Turning them on (both the R3 and V3 notches [FMs 6 & 7], the 1st, 2nd, and 3rd harmonic violin mode notches [FMs 8-10] independently) makes the loop unstable, as indicated by the MC2 output at all stages increasing wildly hitting the limiters. We'll have to explore / model that later.

Save the notches, and the factor of 2 difference in M2 and M3 gain, the configuration should be the same as what's at LLO.

Results of ETMX main and reaction chain top mass transfer functions for Phase 2c, after final alignement, are attached below. Most of the measurements are showing good agreement with the model, except for pitch, as on last thursday. The first two modes of the main chain @~0.55Hz and ~1.30Hz appears to be shifted down in frequency (page 5 of first pdf). The second pitch mode of the reaction chain @ ~ 1.33Hz is shifted up, and we can see vertical to pitch as well as transverse to pitch cross couplings from the extra 2Hz and 2.22Hz resonnances (page 5 of the second pdf). After plotting ETMX pitch transfer function against other fiber quads (3rd pdf), it seems like a common feature between quads. The results will be discussed during the sus meeting, but are approved for now since livingston ITMY looks identical.

Plots attached are described below :

1- H1 M0-M0 ETMX phase 2c main chain top mass transfer functions (in red) against "fiber" model (in blue), for the 6 degrees of freedom

2- H1 R0-R0 ETMX phase 2c reaction chain top mass transfer function (in red) against "erm (end reaction mass)" model (in blue), for the 6 degrees of freedom

3-a) Pitch to pitch main chain comparison between model and other fiber quads : L1 ITMX (June 2013, orange), L1 ITMY (Mai 2013, black), H2 ETMY (Mai 2013, pink), H1 ETMX (last week during alignement, cyan) and H1 ETMX (end of last week after alignement, green)

3-b) Pitch to pitch reaction chain comparison between erm model, H1 ETMX (most recent) and H1 ETMY (Mai 2013)