40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log, Page 333 of 336  Not logged in ELOG logo
ID Date Author Type Categorydown Subject
  8418   Fri Apr 5 01:28:56 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : populating the table

I started populating the end table; the TRY path to start with. I found that I need to redo the cables/electronics layout around the table as we have only one cable feedthrough hole with the new box right now. I need another hand with this and will have Annalisa help me tomorrow.

P.S. I misaligned PRM and restored ETMY to get TRY flashes. I tweaked ETMY to see strong TEM00 flashes.

Old slider positions on medm screen in case we need to restore them:

        TT1    TT2     ITMY    ETMY

P   -1.3586  0.8443   0.9114   -3.7693

Y   0.3049   1.1507  -0.2823  -0.2761

 

  8419   Sat Apr 6 09:21:36 2013 ranaUpdate40m UpgradingNot a fan of the new plastic box yet

 

 1) We still need to drill and install the thumbscrew latches which secure the lids to the table. We cannot use the tables as an acoustic enclosure with loose lids.

2) For the camera issue, the idea is to put the longpass filters on the front of the cameras: then they are only sensitive to light with wavelength > 800 nm.

3) Whenever any interferometer work is happening the light switches must be in the positions which have been marked on them (and which most everyone ignores foolishly). We have never been insensitive to the room lights; the black table enclosures just give us a false sense of security. Room lights impact the interferometer noise.

  8423   Mon Apr 8 16:37:26 2013 SteveUpdate40m UpgradingETMY enclosure wall transmission

 

 Enclosure cover #1 transmission measured  in 1064 nm,  156 mW,  P polarization and beam size ~ 1 mm

As condition:  fully assembled, protective layer removed, tinted- adhesive activated on yellow acrylic on top of each other.

T = 1.2 % in 20 minutes exposure test. This agrees with the test measurement of 6-18-2012

There is a reflected  2-3 cm circular glare that is barely visible on sensor card. It is well below 1 mW level

As we are installing the NPRO with ~350 mW of power we have to address what additional shield should be installed.

The June 2012 test with 1W power burned through of the 3 layer  IR coated films in 3-4 hours.

We 'll use Aluminum shields in the high power path till we come up with better solution.

  8425   Tue Apr 9 00:15:18 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : populating the table

[Den, Annalisa, Manasa]

The Alberto laser was moved from the PSL table. The height of the heat sink rendered a beam height of only 3 inches. I did not want to deal with changing beam height at the table. So, we went ahead and used the old heat sink. I used the beam scan to make measurements of the beam width to match my mode-matching calculations and found some mismatch with the measurements done earlier. So I will measure the beam width again before alignment.

I will also have to change the layout because of the supporting posts that have come up with the new box. Annalisa is doing a COMSOL model to check what the thickness of these supporting posts should be so that the box stays stiff.

  8427   Tue Apr 9 00:32:57 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : TRY temporarily in place

The TRY path on the end table is temporarily in place to help IFO locking.
The Y arm transmission was steered to get TRY back on the PD and the camera. I found that TRY is a couple of inches off in yaw at the end table (comparing to the CAD layout and the earlier layout) and I believe it is because of the changes in input pointing.
I've used a Y1 mirror to steer the Y transmission to an R98% BS. The reflected beam falls on PDA520 and the transmitted beam is steered to the camera. The earlier normalization of TRY is no more valid as the power distribution at the PD has changed.

  8430   Tue Apr 9 08:37:54 2013 SteveUpdate40m UpgradingEndtable upgrade for auxiliary green laser : TRY temporarily in place

Quote:

The TRY path on the end table is temporarily in place to help IFO locking.
The Y arm transmission was steered to get TRY back on the PD and the camera. I found that TRY is a couple of inches off in yaw at the end table (comparing to the CAD layout and the earlier layout) and I believe it is because of the changes in input pointing.
I've used a Y1 mirror to steer the Y transmission to an R98% BS. The reflected beam falls on PDA520 and the transmitted beam is steered to the camera. The earlier normalization of TRY is no more valid as the power distribution at the PD has changed.

Temporary  acrylic wind guard added between enclosure and ETMY transmission window to help IFO locking

  8432   Tue Apr 9 21:27:48 2013 JenneUpdate40m Upgrading TRY temporarily in place

Quote:

I've used a Y1 mirror to steer the Y transmission to an R98% BS. The reflected beam falls on PDA520 and the transmitted beam is steered to the camera. The earlier normalization of TRY is no more valid as the power distribution at the PD has changed.

 To take this into account, last night, I reduced the TRY gain by a factor of 2.  This is not exactly correct - when the layout is finalized we need to figure out what the pickoff situation used to be (we think, based on the Xend, that it could have been 0.5*0.9), and do the correct normalization.

  8436   Wed Apr 10 14:25:04 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : LWE Alberto laser beam profile

I measured the beam profile of the Alberto laser (that will be the ETMY end laser) as I found the data inconsistent with the measurements done earlier.

Method

The laser was set to nominal current (ADJ =0) and the output measured 330mW. I used a 99% BS and measured the beam profile of the transmitted light at several points along propagation using Beamscan. I am attaching the data and matlab script for the fit for future reference.

LWE_Alberto_BP.png

Attachment 2: beam_waist.zip
  8458   Wed Apr 17 02:20:13 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

 

Assembly progress:

1. ETMY oplev setup has been put together. Because of the shift in the TRY path, I had to modify the oplev path on the table as well.

TRY_OPL.jpg

2. Green laser setup coming together:
    (i) Used a HWP after the NPRO to convert s-polarization to p-polarization. (Verified by introducing a PBS after the HWP and then removed later).
    (ii) Lens focuses the beam into the Faraday.
    (iii) Used steering mirrors to align the beam to the faraday. With 320mW before the Faraday, I was able to get 240mW after the output aperture. The spec sheet for the faraday specifies a 93% transmission; but what I measure is only 75%.

GRY.jpg

 

  8472   Mon Apr 22 17:43:09 2013 SteveUpdate40m Upgrading ETMY optical table & enclosure is ready for optics

Quote:

Quote:

 Enclosure is at the east end. It has it's bottom o-ring in place. It will be ready for optics tomorrow around 5pm

I have to shim out the enclosure, finish leveling the table and cut surgical tubing O-ring for the top.

 

 Glued surgical latex tubing with super glue into O-ring shape. The existing in place tubing K-100, OD 0.125" (actual size 0.140"), wall 0.031", ID 0.062".

I have just found out that tolerances on tubing OD are + - 0.026" by the manufacturer. I'm getting larger tubing for better fit.

The table is ready for optics.

Things left to do:

1, finalize o-ring size  2, finish cable feedthrough  3, finalize window connection 4, IR-Thermashield strips for bridge sides 5, replace bridge support post with solid one

 

 I'm working on to improve the quality of the enclosure.

The short comings are: more cable feedthroughs needed, latches to anchor top covers air tight and posts to support bending bridges.

Red triangles are compression latches at 10 places to hold the top air tight on surgical tubing

Green lines represent 4  posts of Al 1" OD to support the covers and maximize their eigenfrequencies.

Black crosses are 4 spring loaded push-bottom quick release pins to anchor the top covers to the bridges. This connection will  not be air tight.

(quarter-turn wing head fastener have the same problem) I'm thinking of some solution to minimize the leak.    

Violet _ steel plate (1" wide, 15" long, 0.125" thick)  between the two posts will anchor the quick release pins and make bridge rigid.

Blue rectangle is an other cable feedthrough exiting on the chamber side.

Planning to substitute window with soft - air tight ( Aluminized thin wall hose )  connection to vacuum view port where white circle is representing the Al adaptor ring.

Updated after Wednesday meeting 4-24-2013

 

Attachment 1: more2.jpg
more2.jpg
  8473   Mon Apr 22 19:48:56 2013 JenneUpdate40m Upgrading4 pins enough?

Are 4 of these spring loaded pins enough?  I'm not sure how one pin can hold 2 lids at each point.  It seems like we need 8 pins.

  8475   Tue Apr 23 15:00:20 2013 JenneUpdate40m Upgrading4 pins enough?

Quote:

Are 4 of these spring loaded pins enough?  I'm not sure how one pin can hold 2 lids at each point.  It seems like we need 8 pins.

 Steve has explained to me that the pins will go in between the 2 lids, with a big washer, so that one pin holds both lids at the same time.  4 is the right number.

  8486   Wed Apr 24 15:27:59 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress: New layout

Layout that will be improved upon over the next few days.

Things that need to be updated:
1. Waist size at all optics
2. Beam size at detectors and choice of lenses
3. IPANG & green PD proposed positions

 ETMY_endtable_New.png

  8522   Thu May 2 01:19:41 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

I started to put together optics at the endtable. I am attaching the layout with the green blocks showing the optics that are assembled and will not be moved henceforth unless somebody contradicts.

1. Power after HWP = 314mW
    Power before faraday = 310mW
    Power after faraday = 300mW (the power loss while aligning the faraday earlier was due to the AR coating on the focusing lens before the faraday - it was AR coated for visible and that accounted to the power lost)

2. Since we do not know the length of TGG inside faraday, I measured the beam profile after the faraday so that I can trace the beam without any errors to calculate exact mode matching solutions.

3. The NPRO beam seems to be obviously elliptical as seen on the IR card and also from beam profile measurement. So we cannot skip including cylindrical lenses in the layout.

ETMY_0502.png

 

  8543   Tue May 7 19:06:45 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

I have updated the waists (W)  and beam diameters (D) at 1064nm optics on the endtable.

I am not able to locate the characteristics of QPD-Y and oplev PD and hence took the beam diameter to be half of the detector surface area to determine their positions.

Beam diameter on PDA520 used for TRY was calculated using the transimpedance and responsivity of the PD from an old elog in 2004.

ETMY_endtable_New.png 

  8556   Thu May 9 01:36:32 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

Progress with end table:

Parts in green show assembled optics that will not require any changes. Parts in yellow are in place but will need either change of lenses in their optical path or change in position.

0509.png 

  8560   Thu May 9 22:29:23 2013 ManasaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

[Annalisa, Manasa]

More optics have been put on the table. Direction of the rejected beam from the 532nm faraday estimated to be ~1.7 deg along -y axis.

Transmon QPD, TRY and camera have beams on them for locking Y arm. Oplev configuration is waiting for it's lens to arrive.

 

  8561   Fri May 10 20:05:21 2013 AnnalisaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

I rotated some mounts along the green beam path, and I started aligning the beam again.

The beam is aligned up to the waveplate just before the doubler crystal, even if I couldn't reach more than 88% transmission for the Faraday. Next week I will finish the alignment and I'll put the lenses that Manasa already ordered.

 

  8565   Mon May 13 21:35:55 2013 AnnalisaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

Yend table - Current status

OPLEV

Today the 2m focal length lens along the oplev path (just after the laser) has been added. In Manasa's layout it allows to have a beam waist of 3.8mm on the OPLEV QPD, even if it seems to be smaller.

The laser is closer to the box wall than the layout shows (it's on the line n.1 instead of line n.9), so maybe it has to be moved in the position shown in the layout, as Steve suggests, to leave empty space just before the window.

Rana suggests a 2mm diameter beam on the QPD, so a new calculation has to be done to add a second lens.

GREEN

The beam has been aligned until the doubler, but after the crystal it it has a small tilt, so a better alignment has to be done.

Moreover, the beam waist has to be measured after the Faraday for the green, in way to choose the focal length of the lenses necessary for the mode matching.

Then the three steering mirrors to send the beam into the arm have to be put.

TRANSMON PATH

A lens which has to be put on the Transmon path (already ordered) has to be added, and the beam alignment on the QPD-y and on the PDA520 has to be done.

  8568   Tue May 14 01:13:35 2013 JenneUpdate40m UpgradingTRY realigned

Koji noticed that earlier this afternoon the Yarm ASS was working, but then after dinner it was no longer working.  I saw that the ETMY trans camera beam was clipped.  These things precipitated a visit to the Yend station. 

I saw that the beam on the optic that steers the camera beam to the camera was very, very low, almost falling off the optic.  The only mirror which steers to this optic is the harmonic separator which reflects the IR, and transmits the green.  I turned the pitch knobs on the harmonic separator until the beam was roughly centered on all 3 optics between the separator and the camera (BS to QPD, BS to TRYDC and Y1 for camera).  The yaw was fine, so I didn't touch it.

I then adjusted the steering mirror to the camera, and the BS pointing to the DC PD.  I have not touched the BS pointing to the QPD.  Once the beam was on the TRY PD, Koji ran the ASS script, and I recentered the beam on the DC PD.  During this time, Koji had the Yarm triggering using -1 in the POYDC element of the matrix.

The harmonic separator is not mounted in a nice way (I'm assuming that Annalisa is in the middle of things, and she'll get back to it after the green work), so the TRY PD and camera will need to be aligned again, so I didn't do any ASS-recentering-ASS iteration tonight.

The Yarm ASS works nicely again, getting TRY to ~0.89 . 

  8571   Tue May 14 14:24:56 2013 SteveUpdate40m Upgradingenclosure removal

 

 I'm planning to remove the ETMY optical table enclosure and move it over to CES Shop 8am Thursday morning.

We'll install spring loaded lathes, hooks and quick release pins.

The bridge will be reinforced with steel plate to support release pins on posts.

There will be an other cut out for cable feedtrough as it is shown on elog #8472

Let me know if this timing does not fit your work.

  8576   Tue May 14 21:03:15 2013 AnnalisaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

 

GREEN

The new lenses arrived, and I put the right 250mm before the doubler. I'm still not so confident with the alignment, because I cannot get more than 11-12 uW out from the "green" Faraday, with more than 200uW going in.

TRANSMON

I replaced the Y1 mirror with an HR1064-HT532. The alignment has to be done. Today the 50cm focal length lens arrived, and I'm going to put in tomorrow.

 

  8584   Wed May 15 21:27:39 2013 AnnalisaUpdate40m UpgradingEndtable upgrade for auxiliary green laser : progress

 

GREEN

I still have problems in maximizing the power out from the doubler. I realized that the real green power I obtain is about 30 uW, and it is the power which really enters the Faraday.

Before I was measuring it just after the Harmonic separator, and there was some residual IR beam which increased the power on the power meter, that's why I obtained about 200 uW.

I also tried to slightly vary the position of the mode matching lens, but I was not able to get more than 30 uW on the power meter.

TRANSMON PATH

The 50 cm focal length lens has been added in the position shown on Manasa's layout, and the beam has been focused on the PD.

 

 

  31   Tue Oct 30 16:55:40 2007 tobinRoutine Drag-wiping perfected
Steve, Tobin

Steve procured an assortment of syringes from the bio storeroom and we practiced drag-wiping the SOS in the flow bench. Using a 50 microliter Hamilton syringe to deliver 16 microliters of methanol seems perfect for drag-wiping the small optics. Drag-wiping in the downward direction seems to work very well, since we can squirt the optic directly in the center, and the (half) piece of kodak lens tissue fits easily between the bottom two earthquake stops.
  108   Thu Nov 15 18:36:48 2007 JohnSummary PSL table work
I've rotated the lambda/2 plate to 340deg (from 6 deg) and blocked one arm of the Mach-Zender. Undo both if you need to.
  130   Wed Nov 28 12:43:53 2007 AndreyBureaucracy Here was the PDF-file of my presentation

I was making a report with powerpoint presentation during that Wednesday's 40-m meeting.

Here was the pdf-file, but LATER IN THE EVENING I CREATED A WIKI-40M-page describing the algorithm, and now the pdf-file is ON THAT WIKI-40M PAGE.


NOTE ADDED AFTER THE PRESENTATION: I double checked, I am indeed taking the root-mean-square of a difference, as we discussed during my talk.

My slide #17 "Calculation of differential length" was wrong, but now I corrected it.
  215   Thu Dec 27 12:12:02 2007 pkpUpdate Update on GigE Camera
So I finally got the linux software to compile on mafalda. I got the software to dump all the information regarding the camera onto a file. I tried to take a tiff snap and came up empty. So I looked at the configuration file and realized that the camera thinks that the frame-rate is a nan. Am reading up the manual to fix the frame-rate manually and then will attempt to take another snap.

All the files are in a folder called Prosilica in /home/controls/ on mafalda. All the executables are in /home/controls/Prosilica/bin-pc/x86/* . On another note, I am looking for a name for the camera. Any suggestions are welcome.
  254   Wed Jan 23 09:27:55 2008 steveUpdate etmy sus damping restored
Seismis event trips etmy
Attachment 1: etmysus.jpg
etmysus.jpg
  548   Fri Jun 20 02:20:33 2008 YoichiUpdate PMC transmittance degradation
The PMC transmitted light power has been degrading constantly for last two weeks (see the attachment 1).
I went down to 2.55V.
The output of the MOPA is constant during this period. More strangely, the reflected power from the PMC is also constant.
One possible explanation is the contamination of the PMC mirrors. But I don't know why it started two weeks ago.

I tweaked the alignment of PMC and was able to recover the transmitted power to above 2.7V (attachment 2).
I will keep eye on this issue.
Attachment 1: pmc_trans_long_trend.pdf
pmc_trans_long_trend.pdf
Attachment 2: pmc_trans_improvement.pdf
pmc_trans_improvement.pdf
  563   Wed Jun 25 09:46:45 2008 SharonUpdate Adaptive Filters
I have been learning about different methods for applying adaptive filters to improve the Mode Cleaner lock in specific, and other LIGO systems in general.
Finding the exact number of coeffs we would like to apply for our FIR adaptive filter is very important to the efficiency of the filter. Getting this number higher might improve the accuracy of the filter, but costs time we do not have. Another important number to find is the step size. The step size is the variable that controls how far back we want to look into our data for finding the new coeffs. To understand more about the step size it is necessary to learn about the standard deviation of our input and output signals. By getting the step size too big, we are considering long term behavior, but might be missing out on a short term one.
In the near future I will be learning about the meanings of these variables and their contribution to the over all accuracy of our filters.
Results will be posted.
  577   Thu Jun 26 18:29:44 2008 JenneUpdate MC Back online
Jenne, Rob, Yoichi

I was playing with the Mode Cleaner earlier today, working on measuring the effect of the new filter (see next post for loop gain measurements etc.), and bumped something which made it so the Mode Cleaner would not lock.

After much poking around by Rob and Yoichi, we found that the TNC-to-2 pin LEMO from Out1 of the MC Servo Board to Ch1 of the Pentek Generic board to the right of the MC Servo Board was bad. If we touched the cable, the MC would lock, but as soon as we let go, the MC would fall out of lock. The LEMO end of the cable was not heat shrink wrapped, so the 2 wires could have been touching. I replaced the cable, and the MC locked immediately after plugging it in, so I think that has fixed the problem.
  680   Wed Jul 16 11:26:47 2008 Max JonesUpdate This Week
Baffles.

I got a battery for the magnetometer today which is slightly too large (~2 mm) in one dimension. Not sure what I'm going to do.

I'm attempting to calibrate the magnetometer but I'm having a hard time calibrating the axis that I cannot simply put through a coil parallel to the coils length. I have attempted to use the end fields of the solenoid but the measurements from the magnetometer are significantly different from the theoretical calculations.

I would appreciate suggestions. - Max.
  704   Mon Jul 21 09:52:05 2008 SharonUpdate Adaptive code changes
Thanks to Alex, we now save the coefficients of the adaptive filter every cycle. When we choose ENABLE: OFF on the MEDM screen, suppressing the signal to the MC1, we stop and save the last coefficients. When enabling it again, it starts from the last coefficients saved. I will take some measurements today to check how this contributes to the adaptation rate. If you change the number of taps or the number of AUX channels, the coefficients are again set to zero.
  709   Mon Jul 21 19:48:57 2008 SharonUpdate how tp restart C1ASS
How to restart C1ASS:

1. reboot
2. as root: caltech/target/c1ass:> ./startass
3. no need for root: burtgooey

that's it...
  720   Wed Jul 23 10:47:05 2008 SharonUpdate Weekly update
This week I spent some time with Alex who updated the adaptive code to save the filter's coeffs all the time, stop when I open the loop, and reload the latest coeffs. when I start it again.
The point was to minimize the adaptation rate. Unfortunately, seems it is making some filters go wild, so it is not in use yet.
After taking some more measurements with the adaptive filter running, I have noticed a new peak in the signal around 22-23 Hz. My first assumption was that this is caused due to internal resonance of MC1 (which is driven when the adaptive code is running, and not when it's not). Therefore, I drove MC1 without the adaptive filter looking for the same peak... which wasn't there.
This sent me back to the adaptive code... Seems there is a matrix in the simulink file of the adaptive filter which doesn't have an MEDM screen. I am now working on making this screen. Once I am done with that, and make sure there is correlation between simulink and MEDM, I'll keep on chasing the peak in the code.
  723   Wed Jul 23 13:52:26 2008 SharonUpdate MEDM changes
There is a new MEDM screen now when you go from c1ass>top>pem.
Instead of having 12 "mini filters" screens go to 8 outputs (with the wrong correlation impression from the table), we have a 24X8 matrix.
This matrix is there so you could choose which noise signals you want to send to the adaptive code. When you indicate the number of noise channels you are going to use
on the nAUX option on the screen top, you are choosing the channels 1 to nAUX. Channels 15-22 are the seismic and accelerometers we are now using. (you can see the order in Jenne's post 673).
Hope this will make things clearer.
Attachment 1: matrix
  739   Fri Jul 25 13:30:53 2008 SharonUpdate Changes in ASS computer
I editted the simulink diagram of the ASS computer so it now has 2 more channels reading 2 sets of the FIR coefficients to match Alex's changes in the C code.
The new simulink has already been compiled and can be found in /cvs/cds/caltech/users/alex/cds/advLigo/src/epics/simLink/ass.mdl
I backed up the old file and it's also in that folder under ass_BAK_24_jul.mdl

There is also a backup of the old ASS.ini file in caltech/chans/daq/C1ASS_BAK_24_jul.ini

Will update once it's all set and running
  740   Fri Jul 25 17:32:46 2008 SharonUpdate ASS computer
So, it seems a bit too complicated getting the coefficients the way I wanted it to happen (simulink-.ini...).
I returned everything to the way it was and it's all working. The new plan is to choose the specific channel I want to find its instantanous coefficients, let the adaptive code run for a while, setting mu and tau to zero (freezing the coefficients), and exciting the noise signal channel taking the transfer function. This way I can find the filter I want to simulate with an IIR filter.
Once I have the mode cleaner to myself, I'll start posting results.
  747   Mon Jul 28 12:02:32 2008 SharonUpdate accelerometers settings
Jenne, Sharon


We looked again at the channels of the accelerometers and there are some updates. Last time when we reported, we crossed the ADAP channels and the accelerometer. Now that there is a new MEDM screen, with which you can control which channels goes to which adaptive channels, this has no meaning...
Therefore, the channels that go with the noise source channels are:

PEM 15 MC1_X
PEM 16 MC1_Y
PEM 17 MC1_Z
PEM 18 MC2_X
PEM 19 MC2_Y
PEM 20 MC2_Z
PEM 21 SEIS

disregard the last post regarding these channels by Jenne, since I am changing the ADAP channels all the time...
  750   Mon Jul 28 17:58:05 2008 SharonUpdate TOP screen changes
I wanted to test the adaptive code with a downsampling rate of 32 instead of 16. To do this I entered a 32 Hz ((2048/32)/2 - match Nyquist Freq.) low pass filter on the ERROR EMPH, MC1 and the relevant PEM channels.
  760   Tue Jul 29 21:04:55 2008 SharonUpdate OSEM's Power Spectrum
From 16:30 this afternoon
Attachment 1: ITMY2.JPEG
ITMY2.JPEG
Attachment 2: ITMY.JPEG
ITMY.JPEG
Attachment 3: ITMX2.JPEG
ITMX2.JPEG
Attachment 4: ITMX.JPEG
ITMX.JPEG
Attachment 5: ETMY2.JPEG
ETMY2.JPEG
Attachment 6: ETMY.JPEG
ETMY.JPEG
Attachment 7: ETMX2.JPEG
ETMX2.JPEG
Attachment 8: ETMX.JPEG
ETMX.JPEG
  765   Wed Jul 30 12:36:19 2008 SharonUpdate Weekly update
This week included many computer's issues. I tested Alex's new C code (the one that saves the FIR coefficients and restores them when you start running the code again). Seems there is an improvement in the adaptation time, but not a significant one (more details on the coming report). I had to recompile simulink and the FB whenever I wanted to find a solution for taking the record of those coefficients. This is so I could simulate the adaptive filter with a regular IIR filter and compare the two.
After Rob tried to help and it seems to be an impossible to a huge hassle mission, we thought of a different method to do this. I re-compiled the old simulink file and restored the .ini file and all should be back in place. Instead of finding the FIR coefficients, I am going to use one noise source in the adaptive filter, stop the adaptation (by setting mu and tau to 0), and put excitation instead of the noise signal. The transfer function I will get between the exc. and MC1_IN1 is the filter I am looking for.

Also seems that whenever I get the MC unstable, and the adaptive code stops itself, it doesn't come back. Setting the reset flag to a different number (anything other than 0) and pressing the reset button will get it working again, but the CPU will always flip and the ASS computer needed a restart. Still haven't found a problem in the C code, but that's the plan. Moreover, I want to change Alex's code, so that instead of starting from zero like in Matt's code, or starting from the old coefficients like in Alex's, it is going to calculate a Wiener Filter as the first set of coefficients. This will hopefully reduce the adaptation time.

I have also been working on my progress report, and stood in line for the MC... Still standing...
  814   Fri Aug 8 11:04:34 2008 SharonUpdate MCL Wiener filter
I took some old data from Rana and converted the units of the Weiner filter to m/m so to make the effect of the seismometer and accelerometers more obvious.

The data is in counts, and so to convert to m this is what I did:

%%% MC_L calibration
v_per_counts = 5/32768;
v_per_v = 3;
amp_per_N = 1/0.064;

%%% Accelerometers calibration
v_per_counts_acc = 61.045e-6;
g_per_v = 9.8/100;

%%% Seismometer calibration
v_per_counts_seis = 61.045e-6;
m_per_s_per_s_per_volt = 9.8/100;
m_per_v_per_s = 1/345;



for jj=1:6
hfmatm(:,jj)=hfmat(:,jj).*(v_per_counts.*v_per_v.*amp_per_N.*f)./(v_per_counts_acc*g_per_v); %%% accelerometers' data
end
hfmatm(:,7)=hfmat(:,7).*(v_per_counts.*v_per_v.*amp_per_N)./(v_per_counts_seis*m_per_v_per_s); %%% Seismometer data
Attachment 1: m_per_m.png
m_per_m.png
  841   Fri Aug 15 16:45:43 2008 SharonUpdate Converting from FIR to IIR
I have been looking into different techniques to convert from FIR to IIR. This is so we can see how effective the adaptive FIR filter is in comparison to an IIR one with fewer taps.
For now I tried 2 matlab functions: Prony and stmcb (which works on LMS method).
I used the FIR wiener code, MC1_X, (c1wino) and applied the FIR to IIR algorithm.
Seems like the stmcb one works a bit better, and that there isn't much effect for having 1000 and not 400 taps.

Will keep updating on more results I have, and hopefully have the MC in time to actually check this live.
Attachment 1: fir2iir50.png
fir2iir50.png
Attachment 2: fir2iir400.png
fir2iir400.png
Attachment 3: fir2iir1000.png
fir2iir1000.png
  850   Tue Aug 19 10:36:34 2008 SharonUpdate Calibrating accelerometers
I took apart the accelerometers near MC1 and MC2.
The 2 sets of 3 accelerometers are now covered by a box on the floor. Please try not to move them... I will place it all back once I am done calibrating.
  854   Tue Aug 19 17:00:19 2008 SharonUpdate Wiener TF calibration - update
This is an update for post 814

I added the calibration gains I got for the accelerometers (I realize I am just calibrating the accelerometers to themselves and this is not m/m exactly since we don't really know which accelerometer is doing exactly what we want it to do. However, since we are talking on relative small numbers, this shouldn't really change much).


I also added another missing gain for the seismometer. Rana has previously installed a 4300 ohm resistor in the seismometer, which changed the gain to 4300/(4300+5000) = 0.46 (this is from the manual). Moreover, there is a gain of 100 on the SR560. This comes up to an extra gain of 46, meaning multiplying the seismometer's counts by 1/46.
Attachment 1: m_per_m.png
m_per_m.png
  855   Tue Aug 19 17:15:34 2008 SharonUpdate MEDM
I plugged in the gains I got for the accelerometers in the accelerometers' filters in the PEM screen of the adaptive filter
  863   Wed Aug 20 17:02:01 2008 SharonUpdate More FIR to IIR
I tested another method for converting from FIR to IIR other than the 2 mentioned in post 841.
I got this one from Yoichi, called poles fitting, you can read about it more if you want at: http://www.rssd.esa.int/SP/LISAPATHFINDER/docs/Data_Analysis/DA_Six/Heinzel.pdf.

Seems it's not doing much good for us though.

I am attaching a PDF file with the plots, which have N=50,100,600,1000, respectivaly.
Attachment 1: polefit1.pdf
polefit1.pdf polefit1.pdf polefit1.pdf polefit1.pdf
  870   Fri Aug 22 13:58:39 2008 SharonSummary Trend of the Wiener TF
In order to understand if we really need an adaptive filter, I used old data of MC_L and the accelerometers and seismometer to see if the Wiener (ideal) TF between MC_L and the others really changes all the time.
Two tests I made:
  • Compare the TF after different segments of time, starting from the same point. Meaning, measuring the TF after 5,10,15,20... minutes, looking when and if the TF stablizes (stops changing).
  • Compare the TF between same-length segments, from different times. Meaning, comparing for example 2 segments of 10 minutes taken from different times.

Results:
  • As you can see in the attached PDF, the changes start being minor after 200,000 data points, which correspont to 200,000/256 s, which is approximately 13 minutes.
    If you look at the PDF file, it is arranged from shorter times to longer in the order of: 3, 6, 13, 26 and 39 minutes.

  • As expected, the TF between different segmants of the same length is not completely the same. Again, you can look at the attached PDF.
    Sorry the titles are the same. Each 2 consecutive pages represent the same length of segment in different times. The order of segment's lengths is: 3, 13, 26 and 39 minutes

How do I explain what's going on?

Since the Wiener filter finds the correlation matrix between the data and the noise signals, it will maintain some kind of familiar shape when we don't add a significant amount of unusual data. I am assuming that if I had looked at longer time periods, we could see a more significant change in the TF in time. When looking at different times, the average noise is likely to be different which can explain the change in the correlation matrix and the TF.

To sum up

I think we should give adaptive filtering a go.
Attachment 1: Same_start_time1.pdf
Same_start_time1.pdf Same_start_time1.pdf Same_start_time1.pdf Same_start_time1.pdf Same_start_time1.pdf
Attachment 2: same_length.pdf
same_length.pdf same_length.pdf same_length.pdf same_length.pdf same_length.pdf same_length.pdf same_length.pdf same_length.pdf
  916   Wed Sep 3 18:45:01 2008 AlbertoConfiguration PD3 gain
Alberto, Yoichi,

We found that the PD3 servo was unstable with a gain of 1, so we switched it to 0.5
ELOG V3.1.3-