| ID |
Date |
Author |
Type |
Category |
Subject |
|
17314
|
Sun Nov 27 15:30:22 2022 |
Chris | Update | Optimal Control | IMC alignment controller testing |
Five more mode cleaner alignment controllers were tested this morning (remotely). These were designed to run in tandem with the standard controller, instead of supplanting it. Before the test, c1ioo was burt restored back to the settings of the previous test on Oct 28, and in MC TRANS PIT/YAW filter banks the 80 dB gain filters were disengaged and outputs were enabled. Subsequently, all settings were returned to the original values. Each test consisted of five minutes with pitch alignment uncontrolled, five minutes with the standard controller only, and twenty minutes with both controllers enabled. GPS times for each phase of testing are the following:
- musgo
- OL start 1353602764
- CL start 1353603074
- policy start 1353603410
- musgo_ghost
- OL start 1353604697
- CL start 1353605007
- policy start 1353605355
- musgo_stumble
- OL start 1353606574
- CL start 1353606884
- policy start 1353607229
- musgo_goldfish
- OL start 1353608446
- CL start 1353608756
- policy start 1353609099
- musgo_late
- OL start 1353610321
- CL start 1353610631
- policy start 1353610971
|
|
17333
|
Sun Dec 4 10:03:02 2022 |
Chris | Update | Optimal Control | IMC alignment controller testing |
Another five mode cleaner alignment controllers were tested last night (remotely), running in tandem with the standard controller. As before, c1ioo was burt restored to Oct 28 and the MC TRANS PIT/YAW 80dB gain filters were disabled before the test. Each test consisted of five minutes with pitch alignment uncontrolled, five minutes with the standard controller only, and twenty minutes with both controllers enabled.
The first four tests went smoothly, but the last controller (goldfish_short) repeatedly broke the lock. Eventually I got it running with an output gain of 0.5, strong enough to see the misbehavior without unlocking the mode cleaner.
GPS times for each phase of testing were the following:
- ichabod
- Open loop 1354165556
- Closed loop 1354165866
- Policy 1354166241
- ichabod_2
- Open loop 1354167462
- Closed loop 1354167772
- Policy 1354168113
- ichabod_3
- Open loop 1354169357
- Closed loop 1354169667
- Policy 1354170006
- goldfish_long
- Open loop 1354171297
- Closed loop 1354171607
- Policy 1354172022
- goldfish_short
- Open loop 1354173255
- Closed loop 1354173565
- Policy 1354173924 (output gain 1.0, immediately unlocked the cavity)
- Policy 1354175008 (output gain 0.1, 2 min)
- Policy 1354175189 (output gain 0.3, 2 min)
- Policy 1354175376 (output gain 0.5, 20 min)
|
|
17371
|
Wed Dec 21 12:38:32 2022 |
Chris | Update | Optimal Control | IMC alignment controller testing |
Three additional mode cleaner alignment controllers were tested Sunday night (remotely). They were run in tandem with the (recently improved) standard controller. Each test consisted of five minutes with pitch alignment uncontrolled, five minutes with the standard controller only, and twenty minutes with both controllers enabled.
GPS times for each phase of testing were the following:
- slug_functional
- Open loop 1355466269
- Closed loop 1355466579
- Policy 1355466987
- slug_hippocamp
- Open loop 1355468210
- Closed loop 1355468520
- Policy 1355468849
- slug_hippocamp_slow
- Open loop 1355470093
- Closed loop 1355470403
- Policy 1355470855
|
|
17
|
Fri Oct 26 09:10:17 2007 |
steve | Routine | PEM | PEM &PSL trend |
The fires are out, lab particle counts are up.
Psl HEPAs are at 100% and mobel HEPAs are just turned on
20 days plot and 5 hrs plot below |
| Attachment 1: counts&psl.jpg
|
|
| Attachment 2: 5dcounts.jpg
|
|
|
83
|
Thu Nov 8 11:40:21 2007 |
steve | Update | PEM | particle counts are up |
I turned up the psl HEPA filter to 100%
This 4 days plot shows why |
| Attachment 1: pslhepaon.jpg
|
|
|
114
|
Mon Nov 19 14:19:25 2007 |
steve | Omnistructure | PEM | jackhammer |
The construction personal successfully jackhemmered a fence around the "Drever's parking slot"
There is no parking space available close by |
| Attachment 1: jackhammer.jpg
|
|
| Attachment 2: jackhammer2.jpg
|
|
|
151
|
Fri Nov 30 20:17:26 2007 |
Andrey | Configuration | PEM | Accelerometers and alum.plates for them |
All 6 accelerometers which were located near the ITMX are turned off and disconnected from the power cords.
Actually these accelerometers are now in the office area on the electronics bench (to the left from Steve Vass' place).
I made today 4 new aluminum mounting plates for the accelerometers (I drilled holes and made threads in them). On Monday I will buy short screws and install accelerometers on these new mounting plates. These mounting plates will be screwed directly into the metallic frame which is firmly cemented to the ground. Before yesterday accelerometers were mounted on top of blue stack towers, not on the ground directly, so we hope that new measurements of the ground noise will be more realistic.
The 4 mounting plates are on the same desk -> on the electronics bench (to the left from Steve Vass' place). Please do not displace them.
Attached is a drawing of the aluminum mountain plate. |
| Attachment 1: Scheme_Aluminum_Piece-inches.pdf
|
|
|
152
|
Fri Nov 30 21:27:24 2007 |
rana | DAQ | PEM | weather / stacis / c1pem1 |
I was trying to add some Seis BLRMS channels to the c1pem1 processor so that we could have DMT trends.
Then I found that none of the Weather channels have been working for a year or so. I could also not
telnet into it. I tried resetting it but no luck. There was no entry in the Wiki for it so I added
a place holder.
Have the weather channels ever worked? Do we have those sensors? I think I've never actually looked
for this. Seems like a fine ugrad job. |
|
161
|
Mon Dec 3 19:44:58 2007 |
Accelerometers on new mounts | Configuration | PEM | Andrey |
I (Andrey) continued today working with new accelerometer mounting. (see entry #151 about my Friday work).
I bought screws/washers and attached those mounts with accelerometers to metallic frames which are firmly cemented to the floor.
One such mount with three accelerometers (in X-, Y-, Z-directions) is installed near the ITMX (in the previous location, but NOT on top of the unused stack as before Friday), the other mount with three accelerometers in three orthogonal directions is installed near ETMX in the east end of the room (this set of accelerometers was installed between MC and BS before Friday). I uncoiled the cables, put them into the cable tray towards the ETMX, and hooked-up the three accelerometers near ETMX in the east end of the room.
Now all six accelerometers are hooked-up (that is, connected to power supply board with cables).
We decided with Steve Vass to put red cones (similar to those that are on highways in the road construction zones) in order to prevent people from bumping into accelerometers. Please use caution when walking along the X-arm.
I took several pictures of the new accelerometer setup. Picture "DSC_0194.JPG" shows the mount with accelerometers near the the ITMX and the beamsplitter chamber,
picture "DSC_0195.JPG" is the "zoomed-in" view of the same accelerometers, while picture "DSC_0196.JPG" shows the mount with accelerometers near ETMX in the east end of the room.
Many thanks to Mr. Steve Vass for his thorough explanation/showing me how to drill the metal and put threads in the holes. |
| Attachment 1: DSC_0194.JPG
|
|
| Attachment 2: DSC_0195.JPG
|
|
| Attachment 3: DSC_0196.JPG
|
|
|
172
|
Wed Dec 5 23:19:03 2007 |
Andrey | Configuration | PEM | Accelerometers are turned on |
All accelerometers have been turned on, as Alan asked during Wednesday meeting.
Typical power spectra and coherence plots are attached below.
"East" in the name means that the previous location of accelerometrs was to the east from "Beamsplitter" (the location for "east" accelerometers was not changed, actually, it is still near ITMX), while "west" means that previously accelerometers were to the west from the BS, but now their new location is near the ETMX.
I will change the names of the channels tomorrow (Thursday) when someone (Tobin?) will show to me how to do it.
P.S. (addition made on Dec. 19th, 2007, by Andrey) I intended to change the names of accelerometers the next day, Thursday Dec. 06,
but I did not do it that day (did not understand how to do it), then I fell ill, and eventually
I changed the names of accelerometers on December 19th, see entry to ELOG #204) |
| Attachment 1: Power_Sp_and_Coh_XY-EAST.pdf
|
|
| Attachment 2: Coherence-ZX_East.pdf
|
|
| Attachment 3: Coherence-ZY_East.pdf
|
|
| Attachment 4: Power_Sp_WEST.pdf
|
|
| Attachment 5: Coherence-ZX_West.pdf
|
|
| Attachment 6: Coherence-XY_West.pdf
|
|
| Attachment 7: Coherence-YZ_West.pdf
|
|
|
189
|
Wed Dec 12 22:24:48 2007 |
tobin | Frogs | PEM | weather station |
I poked at the weather station briefly this evening.
* There's almost nothing in the elog about it.
* It exists. It is located on the North wall, just north of the beam splitter.
* It seems to be displaying reasonable data for the indoors, but nothing for the outdoor sensors.
* c1pem didn't seem to be starting up (couldn't telnet into it, etc). I altered its startup file and reset it several times, and eventually it came to life.
* the weather station has a serial cable that goes all the way to c1pem. I plugged it in.
* however, the Weather.st program complains "NO COMM"--it gets no data from the weather station
* The next thing to do is to plug in a laptop to that serial cable and see if the weather station can be convinced to talk. |
|
200
|
Wed Dec 19 11:31:01 2007 |
steve | Omnistructure | PEM | aircond filter maintenance |
Jeff is working on all air condiontion units of the 40m lab
This we do every six months. |
| Attachment 1: acfilters6m.jpg
|
|
|
204
|
Wed Dec 19 20:28:27 2007 |
Andrey | DAQ | PEM | Names for all 6 accelerometers have been changed |
I eventually changed the names for all 6 accelerometers (see my ELOG entry # 172 from Dec. 05 about my intent to do that).
I removed the word "BS" from their names,
and I changed the word combination "ACC_BS_EAST" in the old name for "ACC_ITMX" in the new name;
as well "ACC_BS_WEST" is now replaced by "ACC_ETMX".
(the reasoning behind such a change should become clear from my ELOG entry #172).
New accelerometer names are:
(note: there are no spaces (nowhere!) in the names of accelerometers, but ELOG replaces ": P" written without a space by a strange symbol  )
C1 : PEM - ACC _ ETMX _ X ;
C1 : PEM - ACC _ ETMX _ Y ;
C1 : PEM - ACC _ ETMX _ Z ;
C1 : PEM - ACC _ ITMX _ X ;
C1 : PEM - ACC _ ITMX _ Y ;
C1 : PEM - ACC _ ITMX _ Z .
One can find them in "C1 : PEM - ACC" in Dataviewer.
|
|
255
|
Wed Jan 23 11:41:06 2008 |
steve | Update | PEM | sulfur smell in 40m |
Led - acid car batteries were overcharged in the machine shop next door
and sulfuric acid smell is coming over to the ifo room.
Entry room 103 is specially bad. |
|
258
|
Thu Jan 24 11:52:56 2008 |
steve | Update | PEM | the mud is cleaned up & MOPA shutter is opened |
Safety glasses are required again!
I have just opened the mopa shutter.
One janitor came to help with muddy floor.
Rack 1x1 toward ITMX chamber and the south wall of these area
were completely covered by mud. I wiped the floor of bottom of the rack
with towels. The cables were lifted and still should be wiped.
The bottom of LSC rack got less water, only on the west side.
We are ready to bring up the computers.
Thanks to ALL with the clean up, including Alan Rice
who was really helpful.
|
|
264
|
Fri Jan 25 09:22:21 2008 |
steve | Update | PEM | burned toast award goes to |
DYM for collaborating with the enemy.
In order to minimize the number of ants visiting our lab we have to take out side the lab
all food left over and organic waste. If you are eating here do not expect someone else
to clean up after you.
Thanks for your cooperation. |
|
283
|
Mon Jan 28 19:35:55 2008 |
rana | Summary | PEM | Accelerometer and Seismometer Coherences |
The attached PDF shows that there is some strange behavior at low frequencies.
From the plot it looks like to me that the Wilcoxon accelerometers (which are supposed to have good response down to 0.05 Hz) are not displaying real seismic motion below 0.3 Hz. Because the coherence length for seismic waves at those frequencies should be 100's of meters we should expect that the accelerometers would have good coherence (>0.8) down there. Instead, my guess is that its all air currents, temperature, or electronics noise. These sensors are not reliable indicators for the microseism.
The Ranger seismometer, however, seems to work fine down to just below the microseism. The Ranger is mounted down around the X end and pointing in the z-direction. The coherence I plotted between it and EX_Z is larger than any other acc/seis pair (as expected).
JM and I discussed what could be done; if we get a SURF student who's into building stuff we can ask them to make a styrofoam hut for the Wilcoxons to see if that helps anything. JM also asked what the point of all this is.
IF we want to do good Adaptive Noise subtraction then we need sensors which can sense the motion which disturbs the mirrors and they need to sense it with a good SNR to get a good subtraction ratio. If the styrofoam thing doesn't work, we should probably look into getting a Guralp 3-axis seismometer for the corner area and just move the accelerometers down to the ends. The sites have Guralp CMG-40T units (~ 8k$). I think we should check out the CMG-3T or the CMG-3ESP.
Does anyone know someone in the Geo depts that we can borrow one from? |
| Attachment 1: Acc.pdf
|
|
|
295
|
Sun Feb 3 05:02:41 2008 |
rana | Update | PEM | Seism 4 day |
|
| Attachment 1: Screenshot.png
|
|
|
298
|
Tue Feb 5 17:39:05 2008 |
jweiner | Configuration | PEM | PEM-AS_MIC taken down from AS table, will put in PSL enclosure soon |
| I took down the microphone that Andrey hung above the AS table his first week in lab. I want to hang the microphone above the PMC to check the effect of acoustic noise on the PMC. The cables were a little more tangled than I thought so I've only taken the microphone down and haven't hung it back up yet, but on Thursday I'll have enough time to carefully put it up inside the PSL and see what I can find out about acoustic noise inside the PSL. I think the microphone should be sensitive enough for the frequencies we're interested in, and I'll hopefully find out if it's sufficient once I put it up in the PSL. The microphone cable and microphone are on top of the PSL for now. |
|
299
|
Wed Feb 6 09:17:31 2008 |
steve | Update | PEM | IST building construction continoues |
| The bulldosers at work |
| Attachment 1: seismic8d.jpg
|
|
| Attachment 2: seisioo.jpg
|
|
|
303
|
Fri Feb 8 17:55:53 2008 |
josh | Configuration | PEM | PEM-AS_MIC now in PSL enclosure |
I have moved the microphone to the PSL enclosure, hanging near the south (Y) side from a support rod for the overhanging storage area so that it's reasonably close to the PMC. I've fastened it in many places using cable ties to make sure that it won't fall.
Alberto helped me solder together a female BNC-female 3.5 mm stereo adapter so that I can use the DAQ to output through BNC to PC speakers. My plan is do sweep sine output through PC speakers to find the transfer function of sound from outside the enclosure to inside the enclosure and by moving the microphone more centrally over the PSL table, check if there are any strong resonances. Hopefully I can use this technique at other places around the interferometer or measure the effect of installing acoustic foam. |
|
308
|
Mon Feb 11 14:24:19 2008 |
steve | Update | PEM | more earthquakes |
ITMX and ITMY sus damping restored after Baja earthquake 5.1 mag at 10:29 this morning.
The ground preparation for The ITS building is almost finished.
Activity is winding down, however the Baja California_ Mexico earthquake zone
"Guadala Victoria" started acting up on Friday.
http://earthquake.usgs.gov/eqcenter/recenteqsus/Maps/special/California_Nevada_eqs.php |
| Attachment 1: eqfeb11.jpg
|
|
|
324
|
Tue Feb 19 18:28:41 2008 |
John | Update | PEM | More seismic in Baja California |
Steve spotted more activity from the same quake.
Reset watchdog on ETMY. |
|
330
|
Fri Feb 22 02:51:20 2008 |
Andrey | Update | PEM | Accelerometer ITMX seems to be broken |
As people probably know,
I am trying (for a long time) to create a computational program that calculates the evolution of accelerometer time-domain data through stacks and pendulum transfer functions to test masses, and calculate the RMS of differential arm lenght spectrum.
I noticed on Tuesday that time-domain signals from the two accelerometers (one is near ETMX, the other one is near ITMX) seem to have different amplitudes of fluctuations around the mean value. I suspected that this is the main reason why I cannot get the awaited result of minimum of RMS for equal values of Q-factors for ETMX and ITMX suspensions (because we subtract two very different numbers, so we cannot get anything close to zero). I took amplitude spectra of the accelerometer data (dttfft2), and they look very differently for ETMX and ITMX accelerometers. I believe that spectrum of ETMX accelerometer represents seismic noise, but accelerometer ITMX seems to provide us with irrelevant and wrong data. No peaks, just almost monotoneous decreasing curve, and 10 times smaller amplitude. Therefore, ITMX seems to be broken.
I will try tomorrow to clap my hands, shout, yell, near the broken accelerometer to confirm that the accelerometer is broken (more precisely, that either accelerometer itself is broken,
or cable connections, or DAQ channel, but something is wrong). Now it is very late, and I am going home.
See attached figures: time-scale is 10^(-1), 10^0, 10^1, 10^2 Hz. |
| Attachment 1: Accelerom-EYMX-Feb22.jpg
|
|
| Attachment 2: Accelerom-ITMX-Feb22.jpg
|
|
|
336
|
Fri Feb 22 15:16:33 2008 |
Andrey | Update | PEM | ITMX Accelerometer is NOT broken |
As I wrote in message 330, there was a bad signal from ITMX accelerometer. I have found the reason: the BNC-cable which goes from the black board with switches for accelerometer gain (1,10,100) towards DAQ-tower was completely disconnected from that black board with gain-switches. The end of the long BNC-cable was on the floor. Therefore, it was totally impossible to see any accelerometer signal. The cable that I am writing about should transport the signal from ITMX_X accelerometer.
Now all the BNC-connections seem to be in good shape, and spectra of accelerometers near ITMX and ETMX , both of them are in x-directions, are very much similar. |
| Attachment 1: Accelerom-ITMX-Feb23.jpg
|
|
| Attachment 2: Accelerom-ETMX-Feb23.jpg
|
|
|
363
|
Fri Mar 7 00:47:54 2008 |
rana | Configuration | PEM | Ranger SS-1 |
Yesterday evening around 7:30 PM, I changed the Ranger seismometer from a
vertical to a horizontal seismometer. To do this I followed the instructions
in the manual.
1) Lock it down.
2) Turn it sideways. Use the leveling screws to center the bubble level.
3) Carefully loosen the hanger rod and release slowly the tension to allow
the mass to recenter.
4) Look through the little viewhole next to the rod to make the white lines
line up. This means the mass is centered.
5) Look at the output on a scope. It should be freely moving with a ~1 sec.
period.
The attached plot shows the before and after spectra. |
| Attachment 1: ss1.pdf
|
|
|
368
|
Tue Mar 11 23:14:01 2008 |
rana | Configuration | PEM | Accelerometer and Seismometer movements |
Steve and Matt moved the accelerometers and seismometers today.
The accelerometers are now placed around the MC and the seismometer is in-between MC1 & MC2.
We have changed the names of the acc channels to reflect whether they are close to MC1/MC3
or MC2. We tested the accelerometer to channel name mapping by switching gains at the wilcoxon
breakout box and also by tapping. It seems now that the previous setup near the ITMX/ETMX had
some few channels mislabeled which would have given some confusing results.
Alex, Jay, and Rolf came over today and installed, then de-installed some of the hardware for
sending the PEM channels over to the C1ASS machine where the adaptive filter front end will go.
Everything should be back to the way it was...hopefully, the guys will modify the ADCU PEM
code to send the signals to the new FE over the reflective memory net and then send them to the
MCL inputs of the suspensions. So the first incarnation should use the accelerometers and seismometer
to drive MC1 and/or MC3. |
| Attachment 1: Acc.pdf
|
|
|
369
|
Wed Mar 12 00:36:52 2008 |
rana | Configuration | PEM | Accelerometer and Seismometer movements |
I used the MISO FIR Wiener matlab code to see how well we might do in principle.
The attached 3 page PDF file shows the MC_L control signal (force on MC2) and the residual
after subtracting off the accelerometer and seismometer using a 32 Hz sample rate and
512 taps (page 1), 1024 taps (page 2), and 2048 taps (page 3). As Matt smarmily points out,
there's not a lot to win by going beyond 512; maybe a factor of sqrt(2) for a factor of 4
tap number. |
| Attachment 1: finished.pdf
|
|
|
370
|
Wed Mar 12 00:40:35 2008 |
rana | Configuration | PEM | Accelerometer and Seismometer movements |
| Same as above but with 2048 taps and a 128 Hz sample rate. Does much better at the 16 Hz bounce mode. |
| Attachment 1: mc2048-128.pdf
|
|
|
371
|
Wed Mar 12 00:47:26 2008 |
rana | Configuration | PEM | Accelerometer and Seismometer movements |
| And this is a cool snapshot showing how this operation used 16 cores on menkar ! |
| Attachment 1: Screenshot.png
|
|
|
384
|
Mon Mar 17 18:30:48 2008 |
mevans | Configuration | PEM | Adaptive Filtering |
| It seems that adaptive filtering can achieve results similar to those of the MISO FIR Wiener (entry 369). The adaptive code simulates real-time operation, but uses the same data used by Rana for the Wiener filter. I ran the adaptive filter over the data 100 times to ensure that it was well trained... maybe too well. |
| Attachment 1: mcacc_adaptive.png
|
|
|
391
|
Fri Mar 21 23:15:11 2008 |
rana | Configuration | PEM | Ranger SS-1: New Setup |
The Ranger seismometer has been in a bad state. Its output had been sent into a SR560 without any termination.
The seismometer is, internally, just a mass on a flexure with a magnet and a pickup coil for readout.
The damping of the system depends on the resistor hooked up across the coil. With the SR560 this is
the 1 Meg input impedance of it and so the mass is undamped.
I installed a 4300 Ohm resistor in there which seems to nearly critically damp it. However, this will not
allow us to reach the ultimate quantum noise limited performance. We will have to analyze the thermal, voltage,
and current noise to get that.
I then also increased the gain from 10 to 100 on the SR560. This should now make the front end noise of the
seismometer/SR560 close to equal to the noise of the PEM ADC. |
|
418
|
Tue Apr 8 09:08:54 2008 |
steve | Configuration | PEM | weather station disconnected |
We can not leave cables in the walkways and go on vacation.
I disconnected the weather station from the DAQ
Every Tuesday is janitor day in the 40m.
We have to give him free space for doing a good job.
Burned toast award goes to Andrey. |
| Attachment 1: toast.jpg
|
|
|
420
|
Wed Apr 16 09:47:35 2008 |
Andrey | Summary | PEM | Weather Station |
The weather station is functional again.
The long ethernet Cat5 cable connecting 'WeatherLink' and processor 'c1pem1' was repaired yesterday, namely the RJ45 connector was replaced,
and information about weather conditions is now again continuously being transferred from the 'Weather Monitor' to the control UNIX computers. We can see this information in 'c0Checklist.adl' screen and in Dataviewer.
Below are the two sets of trends for the temperature, wind speed and direction, pressure and the amount of precipitation.
The upper set of trends ("Attachment 1") is "Full Data" in Dataviewer for the 3 hours from 6.30AM till 9.30AM this morning,
and the lower set of trends ("Attachment 2") is "Minute Trend" in Dataviewer for 15 hours from 6.30PM yesterday till 9.30AM this morning.
I also updated the wiki-40 page describing the Weather Station and added to there a description of the process of attaching the RJ45 connector to the end of ethernet Cat5 cable. To access the wiki-40 page about the "weather station" you should go from the main page to "PEM" section and click on "Weather Station". |
| Attachment 1: Weather-FullData_3hrs.png
|
|
| Attachment 2: Weather_Trend_15hrs.png
|
|
|
424
|
Thu Apr 17 20:17:37 2008 |
Andrey | Update | PEM | Two issues with our weather station |
I encountered two difficulties working with the "Weather Station".
(1) It turns out that there is no indication for "outside humidity" on the "weather monitor" (a small black box located on the north wall of the interferometer). I realized that "outside humidity" is absent in our system when I tried to see the Dataviewer trend and real-time value from the channel "C1: PEM-weather-outsideHumid". It shows impossible number 128%.
It follows from the "Davis" technical documentation that the outside sensor can be of two types: either "External Temperature Sensor" or "External Temperature/Humidity Sensor". I suspect (I do not know for sure) that we have the first type of sensor "external temperature only" and therefore we in principle cannot have information about outside humidity. I propose to Steve to climb to the roof on Friday to resolve this uncertainty looking at the sensor.
(2) I wanted to change the units of pressure from "Pascal" (force/area) to other units, "mbar" for example. For this purpose I need to edit the file "Weather.st" in the directory /cvs/cds/caltech/target/c1pem1 (this file is run on the VME processor "c1pem1"). Unfortunately, when I try to open the file with emacs, I get the message that the file exists but protected from modifications. I do not know how to unblock the file "Weather.st". I need some help with that.
I thought that switching-off the processor "c1pem1" could resolve the issue, so I switched-off the whole crate where the processor "c1pem1" is installed for about 5 minutes, turning the metallic key. As it did not make any difference for the accessibility of the file "Weather.st", I switched-on the crate after 5 minutes. There are other processors besides "c1pem1", so they were turned-off for several minutes earlier today.
Also, I created a new MEDM screen which has information about weather only, a smaller version of the "C0Checklist.adl" MEDM screen. Both screens are now located under the most top-left button "Checklist" of the main MEDM screen. |
|
427
|
Fri Apr 18 16:48:13 2008 |
Andrey | Update | PEM | Rain collector of weather station |
Today the rain collector of our weather station was cleaned. As a result, we checked that the rain indication on the weather monitor and on the MEDM screens is alive and working properly. I am adding some details about the roof sensors to the wiki-40 page about the weather station. See especially the link "More description of the roof sensors and their interaction with UNIX computers" from the main Weather Station page in wiki-40.
Pictures of the rain collector before (dirty, the opening is fully clogged with dust and dirt) and after (clean opening in the bottom of the bowl) the cleaning are attached. |
| Attachment 1: DSC_0520--before.JPG
|
|
| Attachment 2: DSC_0537--after.JPG
|
|
|
445
|
Thu Apr 24 23:27:48 2008 |
rana | Update | PEM | acoustic noise in MC_F |
I looked at the coherence between the Microphone in the PSL (PEM-AS_MIC) and the MC_F channel.
We want to use a microphone to do Wiener/Adaptive noise cancellation on the MC and so we need to
have a coherence of more than ~0.1 in order for that to have any useful effect.
The attached plot shows the spectrum and coherence with and without the HEPA turned up. As you can
see, the HEPA noise is just barely noticeable in this microphone. 
We will need to get something with at least 20 dB more sensitivity.:P |
|
448
|
Fri Apr 25 13:20:04 2008 |
Andrey | Update | PEM | Microphone test |
In response to Rana's request, I tested the microphone (if it is alive or not) by clapping my hands and speaking aloud nearby.
The microphone is alive, see the attached "Full Data" for 5 minutes from Dataviewer. |
| Attachment 1: Microphone.png
|
|
|
452
|
Sat Apr 26 01:45:38 2008 |
Andrey | Summary | PEM | Weather Station enhancement |
Two more things concerning weather monitoring have been done during this week.
1) A Dataviewer template was created, so that it allows to see "real-time" information from weather channels immediately, without adding many channels "manually".
If one wants to use this template,
open Dataviewer -> "File" -> "Restore Settings", /cvs/cds/caltech/users/Templates/Dataviewer_Templates/Weather.xml.
2) I wrote a couple of Matlab scripts that allow to read data (minute trends) from the Dataviewer channels over some time in the past, save the received data in mat-files, and plot those minute-trends. Thus, one can get plots that are very much similar to what one can see in Dataviewer. These two Matlab files are located in the directory
"/cvs/cds/caltech/users/weather_station". File "WeatherReading.m" allows reading from the weather channels (paths to mDV directory must be configured before using my script), file "WeatherTrends.m" allows plotting of those minute trends.
Unfortunately, hardware problems arise very often if we want to read for a somewhat long time in the past, so until now I have not succeeded in getting trends for more than 20 minutes. As an example, see the attached png-file with the 20-minutes trends of data from Thursday evening.
3) So far I did not have success in learning how to recalculate pressure from Pascals to mbars in EPICS (although I tried google-search).
4) I am making every effort in recent weeks not to put any personal or non-scientific information into elog, but this message could be important for all of us, so I cannot resist:
a shark in the Pacific Ocean has killed a swimmer near San-Diego (I saw this in russian news and then made a quick google-search).
http://latimesblogs.latimes.com/lanow/2008/04/this-just-in-fa.html |
| Attachment 1: Matlab_Weather_Trends.png
|
|
|
460
|
Tue Apr 29 21:30:49 2008 |
Andrey | Update | PEM | In the process of renaming channels for Weather Station |
I startted renaming channels for the weather station, and I will continue this tomorrow, on Wednesday.
I have restarted 'c1pem1' several times and reconfigured "C0DCU1" on the framebuilder MEDM screen.
Framebuilder now does not work. |
|
461
|
Wed Apr 30 20:48:58 2008 |
Andrey | Summary | PEM | New Weather Channels |
I created the new channels for the weather station, all letters are capital ones. They are of the form "C1 : PEM-WS_PARAMETER" where "PARAMETER" is temperature, pressure, wind,... characteristics (names are self-obvious).
These new weather channels are indicated on the "Weather Checklist" MEDM screen. Also, units of pressure were changed from Pascal to torr and mbars.
The new weather channels are also visible in Dataviewer. I updated the template, and as an example of Dataviewer data I attach the following 5-hour trends of weather parameters from 3.30PM to 8.30PM on April 30th. |
| Attachment 1: April30-5hours.png
|
|
|
480
|
Thu May 15 14:39:33 2008 |
Caryn | Summary | PEM | filtering mode cleaner with mic |
| Tried filtering for mode cleaner data(C1:IOO-MC_L) using a siso-firwiener filter and microphone data(C1:PEM-AS_MIC) for noise input. The noise reduction in mode cleaner data using the microphone-filter is comparable to the noise reduction when an accelerometer(C1:PEM-ACC_MC1_X) filter is used. See attached graphs. |
| Attachment 1: MC_L_with_PEM-AS_MIC_filter.pdf
|
|
| Attachment 2: MC_L_with_PEM-ACC_MC1_X_filter.pdf
|
|
|
518
|
Wed Jun 4 16:25:06 2008 |
Caryn | Summary | PEM | microphone moved |
The microphone 'C1:PEM-AS_MIC' has been moved right a bit. This change didn't seem to have much effect on filtering the 'C1:IOO-MC_L' signal, at least not compared to how the filter changes with time. Also used microphone data to filter MC_L data using firwiener filter/levinson. The N(order) and sample rate were varied to see how the filter changed. Attached are graphs of the max(rms(noise_estimate)) vs N or IR for varying srate. Note that filtered_signal=signal-noise_estimate. So, the larger the noise_estimate, the more the filter subtracts from the signal.
Green-filtered signal
blue-noise estimate
red-MC_L signal
note decreasing sample rate is more effective than increasing N (higher N takes more time to compute)
note sample rate doesn't change the max(rms(noise_estimate)) very much if impulse response time remains constant
note the 64hz, N=7000 (impulse response about 110s) filter is a better filter than the 512Hz, N=7000(impulse response about 14s) |
| Attachment 1: 1_MC_L.pdf
|
|
|
522
|
Fri Jun 6 11:19:13 2008 |
Caryn | Summary | PEM | Filtering MC_L and MC_F with PEM:ACC and microphone |
Tried to filter MC_L and MC_F with acc/seis data and microphone data using wiener filter (levinson)
-Used get_mic_data.m and miso_filter_lev.m to make SISO filter for 2 minutes of IOO-MC_F data. Used PEM-AS_MIC signal as noise input data. Filters calculated at initial time were applied to later data in 1 hour intervals.
-microphone filter did not seem to filter MC_F very well in high frequency range using this filtering procedure.
-residual is larger than est (see MC_F pdf)
-Used do_all_time_lev.m to make graph of max(rms(residual)) to N(order) for different times.(note for each N, filter was calculated for initial time and then applied to data at other times).
-relation of max(rms(residual)) to N(order) is time sensitive (note-on graph, time interval is 1hour) (see MC_F pdf)
-Presumably, max(rms(residual)) should decrease as N increases and increase as time increases since the filter probably becomes worse with time. I think the reason this isn't always true in this case is that the max(rms(residual)) corresponds to a peak (possibly a 60Hz multiple) and the wiener filter isn't filtering out that peak very well.
-Used get_z_data.m and miso_filter_lev.m to make MISO filter for 2 minutes of IOO-MC_L used the following signals as noise input data
PEM-ACC_MC1_X
PEM-ACC_MC2_X
PEM-ACC_MC1_Y
PEM-ACC_MC2_Y
PEM-ACC_MC1_Z
PEM-ACC_MC2_Z
PEM-SEIS_MC1_Y
-Filter was applied to later data in 2hour intervals.
-Used do_all_time_lev.m to make graph of max(rms(residual)) to N(order) for different times.(note for each N, filter was calculated for initial time and then applied to data at other times).
-acc/seis filter seemed to filter MC_L OK for 128,256,512Hz srates. 64 Hz wasn't ok for certain N's after a period of time.
-residual is smaller than est for srates not 64Hz (see MC_L pdf)
-residual is larger than est for 64Hz at N=1448 for later times (see MC_L pdf)
-relation of max(rms(residual)) to N is not as time sensitive for higher sample rates (note-on graph, time interval is 2hours) (see MC_L pdf). Perhaps the levinson 64Hz sample rate filter doesn't do as well as time passes for these signals. When the filter didn't do well, the max(rms(residual)) seemed to increase with N.
-For 512Hz sample rate filter the max(rms(residual)) decreased with time. If the max(rms(residual)) were an indication of filter performance, it would mean that the 512Hz filter calculated at the initial time was performing better later as hours passed by! Perhaps max(rms(residual)) isn't always great at indicating filter performance.
Programming notes
-I had to modify values in do_all_time_lev.m to get the program to loop over the srates,N's,times I wanted
-do_all_time_lev.m is not as clean as do_all_lev.m
-for making the plots do_all_lev.m (which isn't really a procedure and is messy) has some examples of how to plot things from do_all_time_lev.m. |
| Attachment 1: MC_F.pdf
|
|
| Attachment 2: MC_L.pdf
|
|
| Attachment 3: miso_filter_lev.m
|
function [s] = miso_filter_lev(N,srate,rat,z)
%MISO_FILTER_LEV(N,srate,z) uses miso_firlev to get levinson
% FIR Wiener filter of order N-1, using impulse response of
% N/srate. z is a structure gotten from the get_data function.
% z(end) is the signal which is filtered using z(i) for all i.
% 'rat' is the fraction of z which will be put into filter
% funtion. The data from z is downsampled using srate and
% detrended. Let rat=1. I don't have that part working yet.
... 107 more lines ...
|
| Attachment 4: get_mic_data.m
|
function[z,t0,duration]=get_mic_data(t,d_t,d)
%get_mic_data gets data for'C1:IOO-MC_F', 'C1:PEM-AS_MIC,
% Example: z = get_mic_data('now',120,60)
% start time is 't- d_t' so d_t should be given in seconds. t should be given
% as a number like 893714452. d is duration in seconds. get_mic_data saves
% data to a file in current directory named 'temp_mic'. You will be asked to
% save file as 'mic_(start_time)_(duration)'.
duration = d;
... 32 more lines ...
|
| Attachment 5: do_all_time_lev.m
|
function[r] = do_all_time_lev(n,t0,int,duration,N,srate,rat,order,time,MC_L,MC_F,sample_rate)
%do_all_time explores how filter performance changes with time, sample rate,
%and order of filter. Outputs data,noise estimate, structure of max
%rms error and other info. It uses get_data, miso_filter_lev, and miso_filter_int and retrives
%MC_Ldata or MC_Fdata for multiple times, calculates a miso_filter for initial-time data
%file, applies filter to the other data files, and keeps track of the...
%max(rms(residual)) for each filter. n+1 is number of data files. int is time interval between
%data files, t0 is start time, duration is duration of each data file, srate
%is the sample rate for which filter is calculated, n_N is number of orders
%of the filter you want the program to calculate,int_N is interval by which N
... 215 more lines ...
|
| Attachment 6: do_all_lev.m
|
function[r] = do_all_lev(n,t0,int,duration,n_N,int_N,n_srate,int_srate,rat,MC_L,MC_F)
%do_all_lev explores how filter performance changes with time, sample rate,
%and order of filter. Outputs data,noise estimate, structure of max
%rms error and other info. It uses get_data, miso_filter_lev, and miso_filter_int and retrives
%MC_Ldata or MC_Fdata for multiple times, calculates a miso_filter for initial-time data
%file, applies filter to the other data files, and graphs the rms of the cost
%function vs time. n+1 is number of data files. int is time interval between
%data files, t0 is start time, duration is duration of each data file, srate
%is the sample rate for which filter is calculated, n_N is number of orders
%of the filter you want the program to calculate,int_N is interval by which N
... 283 more lines ...
|
| Attachment 7: do_all_plot.m
|
function[r] = do_all_plot(r,x,v)
%do_all_plot plots variables contained in r(structure from
%do_all_time_lev).Plots error(r.B.y) vs x. x can be
%'s'(srate),'N'(order),'t'(time),'p'(impulse response). v can be 's','N','t'.
%example: do_all_plot(r,'s','t') makes a plot of error vs srate for
%different times.
kk=1
err_N_srate=0
... 388 more lines ...
|
| Attachment 8: miso_filter_int.m
|
function [s] = miso_filter_int(s,y)
%miso_filter_int inputs a filter and a structure array of data sets y, applies filter to data, and
%outputs a structure with fields: ppos(signal frequ spectrum),perr(cost
%function frequ spectrum),pest(signal estimate frequency
%spectrum),f(frequency),target(signal),est_darm(noise estimate),t(time).
%data file for which filter has been calculated is s (obtained using miso_filter).
%y consists of data structures which will be filtered using
%filter from s. Then the power spectrum of the difference between signal and filtered-data is
%graphed for all the data files of y for comparison too see how well filter performs
%over time. Note if you want to create a y, take z1,z2,z3,etc. structures
... 120 more lines ...
|
|
547
|
Fri Jun 20 01:38:55 2008 |
rana | Update | PEM | 20 day Weather |
Yoichi showed me that its possible to make PNG images from PS using GS:
gs -sDEVICE=png16m -sOutputFile=foo.png bar.ps |
| Attachment 1: test.png
|
|
|
571
|
Thu Jun 26 01:10:18 2008 |
rana | Update | PEM | Alarm Handler indicates that dust level is high |
In its first useful act, the Alarm Handler started beeping indicating that the dust particle
counts for particles of diameter less than 0.5 micron had exceeded 5000 /cu. ft. Here's the
80 day trend of particles, temperature and humidity: |
| Attachment 1: Untitled.png
|
|
|
572
|
Thu Jun 26 10:56:15 2008 |
Max Jones | Update | PEM | Removed Magnetometer |
| I removed the Bartington Magnetometer on the x arm to one of the outside benches. I'll be trying to determine if and how it works today. It makes a horrible high pitched sound which is due to the fact that the battery is probably 16 yrs old. It still works with ac power though and I want to see if it is still operating correctly before I ask to buy a new battery. Sorry for the bother. |
|
581
|
Fri Jun 27 09:20:15 2008 |
steve | Frogs | PEM | dust particle count is up & alarm handler is on |
This 3 years plot show the trend of seasons.
When outside air quality goes bad ( 0.5 micron > 1 million ) the lab follows.
I will demonstrate this effect with a 4th of July fire works calibration.
Let's do not forget the construction activity next door either.
The alarm handler is busy:
It's sound level were reset to a modest level yesterday.
It would be nice to change the alarm sound so it can play Wagner:Der Ring des Nibelungen:
Das Rheingold, Die Walkure & Siegfried and Gotterdammerung
or something more appropriate than the present frog call
1, half micron count is climbing ( it's close to 16k now )
2, MZ refl signal is too high
3, MC lenght servo LO is too |
| Attachment 1: dust3y.jpg
|
|
|
640
|
Mon Jul 7 13:58:37 2008 |
Eric, josephb | DAQ | PEM | Using unused PEM channels to test camera code |
Joe and I have taken control of the EPICS channels C1:PEM-Stacis_EEEX_geo and C1:PEM-Stacis_EEEY_geo since we heard that they are no longer in use. We are currently
using them to test the ability for the Snap camera code to read and write from EPICS channels. Thus, the information being written to these channels is completely unrelated
to their names or previous use. This is only temporary; we'll create our own channels for the camera code shortly (probably within the next couple of days).
- Eric |
|
646
|
Tue Jul 8 10:20:10 2008 |
steve | Update | PEM | HEPA turned on |
It is specially important to run the PSL-HEPA filters when inside counts is peaking at 30,000 counts
There is a small label at the hepa on/off switch:
enclosure open 100%,
low noise off,
normal 60% of Variac voltage setting on the top of the enclosure
Not running the HEPAs will lower the temp fluctuations from 1.5 to o.5 degree C
at the cost of particle counts from 0 to immidiate room counts. |
| Attachment 1: hepa.jpg
|
|