40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log, Page 240 of 339  Not logged in ELOG logo
ID Date Author Type Categoryup Subject
  949   Tue Sep 16 10:57:45 2008 YoichiConfigurationPEMParticle counter gain
Summary:
Since we reduced the integration time of the particle counter by a factor of 10, we had to add a gain of 10
to the EPICS channels C1:PEM-count_full and C1:PEM-count_half.
I asked Alex to change it and he did it. I forgot to ask him to change the gain of C1:PEM-count_half. So now only
C1:PEM-count_full has x10 gain.

Detail:
C1:PEM-count_full and C1:PEM-count_half are 'Soft Channel' records in the database (Pcount.db). The values are actually
written into the VAL fields directly by an SNL code Particle.o.
Particle.o reads data from the RS-232C port, to which the particle counter is connected. Then it parses the data and put values
into relevant EPICS channels using channel access. This means we cannot change the gain of the channels by modifying the
database file. For example, ASLO field does not have any effect when the value is directly written into the VAL field.
We had to modify the SNL code. Alex modified Particle.st and the new SNL object file is Particle_x10.o sitting in 
/cvs/cds/caltech/target/c1psl/. I modified seq.load so that c1psl loads Particle_x10.o when rebooted.
The source code for the old Particle.st can be found on lesath.ligo.caltech.edu in
/export/CDS/d/epics/apple/Caltech/40mVac/40mVacScipe/dev/src
I asked Alex to disclose the location of the source of the new code.
In order to compile the SNL code into an object file for Motorola CPU by ourselves, we have to call Dave Barker at LHO.
  950   Tue Sep 16 13:04:22 2008 YoichiConfigurationPEMC1:PSL-FSS_RMTEMP alarm level changed
At the request of Steve, I modified the HIGH value of C1:PSL-FSS_RMTEMP from 21.27 to 23.0.
The HIHI is set to 23.50.
  1057   Mon Oct 20 09:45:56 2008 steveUpdatePEMPSL HEPA on
The PSL HEPA filter was turned on.
It should be running all times.
The 0.5 micron particle count is up to 20,000 this morning.
  1089   Fri Oct 24 21:49:15 2008 JenneConfigurationPEMShort Seismometer Cable
Bad news regarding the cable that goes between the Guralp seismometer and the box that I've been working on: it's too short by about a factor of 2. Dang it. I've placed the seismometer underneath the Beam Splitter Chamber (where it needs to go), and started running the cable toward the ADC rack where box was planned to go, and as Rana guessed earlier tonight, the cable isn't nearly long enough. We have some options: the seismometer can go back into the half-height rack near the BS, SRM, PRM oplev's optical table where I think it used to be, or it can go into the rack with the Kepco high voltage power supplies and the laser's supply. The cable won't reach any farther than that.

I think that we can just add BNC extensions onto the octopus cable that Bob made for the Guralp box, so all we need to figure out after we decide on a rack is the power for the box.
  1090   Fri Oct 24 22:30:38 2008 Jenne,ranaUpdatePEMNoise from Guralp Seismometer
Attached is a Power Spectrum of the noise on the Vert1 channel of the Guralp seismometer. The noise is in the several hundreds of nV/rtHz up near 50Hz and higher, but is in the several microV/rtHz range at lower frequencies. Our high frequency noise is almost definitely below the noise of the ADC, but the lower frequencies, where we actually care, it's not as clear.

To Do list:
  • Measure the noise of the ADC - is the Guralp Box lower for all frequencies?
  • Use conversion factors to convert this measured noise into the minimum ground motion that we can measure. Is this at least a factor of 100 lower than our regular ground motion?

** UPDATE: This is actually the noise of the Guralp breakout box, not the Guralp itself. It is the noise measured on the output of the box
with the input shorted. The board is configured to have a gain of 20 (10 from the AD620 and 2x for differential drive). We also measured
directly at the AD620 output and all of this noise comes directly from that chip. If Jenne calculates that this noise is too high we would
have to find a replacement with a better low frequency floor (e.g. LT1012 or LT1007 depending on the Guralps source impedance).
Attachment 1: Vert1_Noise_24Oct2008.png
Vert1_Noise_24Oct2008.png
  1100   Wed Oct 29 12:54:28 2008 JenneUpdatePEMCalibrated Guralp Noise compared to average ground motion
Here is a calibrated noise plot of the Guralp seismometer box. This is the same noise measured on Friday, measured at TP3 (just after the first gain stage), with the inputs shorted.

The Guralp calibration is:
                TP3 noise
noise in m/s = -------------------
                 10 * 802(V/(m/s))
The 10 is from the gain of 10 between the output of the seismometer and the input of the breakout box, and the 802 V/(m/s) is from the calibration data that came with the seismometer.

From elog 881 by Rana, in the ~1-50Hz band, the calibration of the Ranger seismometer is 488*10^6 counts/(m/s). Using DataViewer, I estimated that the nighttime ground motion measured by the Ranger is ~3500 counts, and the max daytime ground motion is ~8000 counts. This is what was used for the nighttime/daytime lines in this plot.

It seems like the noise of the Guralp box is fine just as it is, and we don't need to worry about replacing the first gain stage (differential instrumentation amp) with a lower-noise op-amp, since at even the lowest freqs, we have almost a factor of 100 at night, and better than that at higher freqs.

NOTE about the plot: the legend isn't showing the correct colors for the night and day motion - obviously the nighttime motion is the lower RED line, and the day is the higher GREEN line.

Yet another note: When I was measuring the counts on the Ranger, I forgot to subtract the mean, so these numbers are overestimating the ambient ground motion. The blue curve is correct however.
Attachment 1: GuralpVert1Noise_mPERs_Ranger.png
GuralpVert1Noise_mPERs_Ranger.png
  1102   Thu Oct 30 20:39:47 2008 carynConfigurationPEMtemperature sensor
We attached the temperature sensor box to the MC1/MC3 chamber with a C-clamp. We connected the temp sensor to a 2nd box with a short BNC. Bob set up a power cable coming from the X-end towards the MC1/MC3 chamber(Thanks, Bob!) We soldered the end of Bob's power cable to a plug and attached it to the 2nd box (The power supply enters through the 2nd box). A ~20ft BNC cable connects the output signal of the 2nd box to the tall thing by the PSL where all the signals go labeled 1Y2. Once we had everything connected, we put in the fuses for the power supply. So, now the temperature sensor is receiving power. We checked that the power supply was working (we measured +15.08V and -14.95V, and we wanted 15V and -15V so it's OK for now). Tomorrow we will modify C1IOOF.INI file and reboot the frame builder.

About sensor-
There is an LM34 (looks like a transistor) glued w/ epoxy and thermal paste to the inside of a Pomona box ~1"x"1.5"x2". The lid to the box is covered with a 1-2mm thick piece of copper and a little thermal paste is sandwiched between the Pomona lid and the copper piece. A C-clamp attaches the copper piece to the chamber. A BNC is connected to one side of the box (the side with less copper)

About power supply box-
There is a power regulator and an op-amp inside a Pomona box ~2.5"x4"x2". The power regulator is attached to the center of lid of the pomona box with a screw and washer. There's a power plug on the front of the box
Left:+15V:red wire
Center:GND:white wire
Right:-15V:black wire
There are 2 BNC connections on the sides of the box. The left BNC connection is for the output signal and the right BNC connection is for the temperature sensor (if the power plug is coming out of the box towards you).

Sensor location-
Chamber which contains MC1/MC3. On the door facing towards the Y-end. On the bottom-left side. Behind the door. Attached with a C-clamp.

Power supply box location-
Chamber which contains MC1/MC3. On some metal leg thing near the floor facing towards the Y-end. Attached with a zip-tie

Power supply-
Coming from the X-end from a tall thing with all the fuses labeled 1X1
Fuse 160:+15V:red wire
Fuse 171:GND:white wire
Fuse 172:-15V:black wire

Signal-
Going towards the PSL to the tall thing labeled 1Y1 on the rack labeled SN208
ICS-110B
J12 (which we believe corresponds to 50-51 and channel number 13650)
Temperature sensor is connected to J12 with a ~20ft BNC attached to a BNC2LEMO connector we found lying around
  1106   Sun Nov 2 21:37:22 2008 ranaUpdatePEMRanger recovery
The ranger signal has been bad since around 11 AM on Oct 25 (last Saturday). There are no elog
entries from that day, but I am quite sure that someone must have been working around the PSL
rack area.

It looks like what happened is that someone moved the chair with the monitor on it and/or the wooden
stool next to it. That put tension on the cable connecting the SR560 and the seismometer. The SR560
connector now seems loose and I think probably the cable ground wasn't connected. I swapped the
cable over to the "B" side of the SR560 and the ranger signal is now reasonable (very small offset
and normal seismic signal).

Please be careful when working around there. Everyone always says "I didn't do anything" or "it doesn't
effect anything".

We need to clean up the cabling around there in addition to running a new power cable for the RF amplifier
on the POY table.


I have also reduced its sample rate from 2048 to 512 Hz. The data are OK after 909640694.

I also increased the sample rate of AS_MIC from 2048 to 16384 Hz but that one seems to be broken
---->> the microphone seems to be either disconnected or broken.
  1113   Tue Nov 4 01:03:01 2008 ranaSummaryPEMperiodic thump noise in MC1_ACC
There seems to be a periodic thump seen by the MC1 Accelerometers as well as the surrounding optics.

The first 5 hour minute-trend plot shows the periodic thumping as well as the one large saturating event which ruins the
Wiener noise subtraction.

The second plot is a 30 minute second-trend zoom in.
Attachment 1: Untitled.png
Untitled.png
Attachment 2: Untitled2.png
Untitled2.png
  1120   Fri Nov 7 08:08:00 2008 steveUpdatePEMAC turned off in IFO room
All tree AC units in IFO room 104 switched off at 7:45am
I'm looking for the periodic thump signal in the Vertex area
noted in entry 1113 of PEM-ACC_MC1
  1121   Fri Nov 7 10:52:57 2008 steveUpdatePEMAC is back on
The 20 minutes and 6 sec thumps are not related to the 40m ac units
This period is precisely right on so it is not likely mechanical as a source.
The time and frequency domain of this signal at attachment 2&3

First I had the chilling water cut off and later I turned off the fans
as 5 hrs temp plot shows
Attachment 1: acf.jpg
acf.jpg
Attachment 2: acc_mc1.jpg
acc_mc1.jpg
Attachment 3: test.jpg
test.jpg
  1122   Fri Nov 7 15:58:10 2008 ranaUpdatePEMAC is back on
I'll bet Steve a dollar that it is mechanical. The attached PDF compares all of the accelerometers from right now.

You can see that the RMS in MC2 is way bigger than MC1.

In the second PDF file you can see the time series. I had to play around a lot with DTT to get it to work. The DTT/Foton
combo on Allegra is not stable, so make sure your work early and often.

In the plots shown, I am bandpassing the time series from 600-700 Hz. I found that doing so allowed the burp in MC1 to remain
large and reduce the extraneous fuzz in MC2. As you can see there is no such noise in MC2.

So its a noise around 600-700 Hz that comes on quickly and then shuts off after several seconds. Its also very periodic in that
it comes on around every 20 minutes. Steve also tells me (although he refuses to put in the elog) that it started up around
August 20th (?). I feel like someone in the 40m lab ought to be able to guess what this is at this point.

Please convince Steve to elog his findings about when the noise started.

If one goes out there and stands next to it when the trend predicts its happening it becomes clear what it is.
Attachment 1: mc-acc.pdf
mc-acc.pdf
Attachment 2: mc-acc-quad.pdf
mc-acc-quad.pdf
  1127   Mon Nov 10 16:44:14 2008 steveUpdatePEMparticle counter at MC2
The particle counter was moved to MC2 table temporarily.
It is clamped to the table at the south west corner.
It is pumping for 20s and waiting for 2 minutes now.
  1129   Mon Nov 10 19:39:40 2008 steveUpdatePEMthump noise from particle counter
The particle counter is back on the IOOC location on a piece of FOAM
It needs this isolation, so when the pump is running, it's not shaking things.
The counter was counting for 6 sec and it was on holding for 20 mins.
Now I set the counter for 20 sec so it is easy to recognise it's signal and it holds for 2min only.
This will set the alarm handler in action.

Atm1: 40 mins plot
PEM-ACC_MC2_x,y,z up to 13 mins: pcounter at MC2 table, clamped, counting for 20s and holds for 2 mins
PEM-ACC_MC2_x,y,z from 13 to 26 mins: pcounter at MC2 table, not clamped, seated on 2" foam, counting 20s and holds for 2 mins
PEM-ACC_MC1_x,y,z from 26 to 40 mins: pcounter at MC1_IOOC location, not clamped, seated on 2" foam, counting 20s and holds for 2 mins

Rana won the bet



Quote:
I'll bet Steve a dollar that it is mechanical. The attached PDF compares all of the accelerometers from right now.

You can see that the RMS in MC2 is way bigger than MC1.

In the second PDF file you can see the time series. I had to play around a lot with DTT to get it to work. The DTT/Foton
combo on Allegra is not stable, so make sure your work early and often.

In the plots shown, I am bandpassing the time series from 600-700 Hz. I found that doing so allowed the burp in MC1 to remain
large and reduce the extraneous fuzz in MC2. As you can see there is no such noise in MC2.

So its a noise around 600-700 Hz that comes on quickly and then shuts off after several seconds. Its also very periodic in that
it comes on around every 20 minutes. Steve also tells me (although he refuses to put in the elog) that it started up around
August 20th (?). I feel like someone in the 40m lab ought to be able to guess what this is at this point.

Please convince Steve to elog his findings about when the noise started.

If one goes out there and stands next to it when the trend predicts its happening it becomes clear what it is.
Attachment 1: partcfm.jpg
partcfm.jpg
  1131   Wed Nov 12 11:36:13 2008 JenneUpdatePEMGuralp Breakout Box is ~50% hooked up
The Guralp box is about halfway hooked up now. The seismometer is under the BSC, and the long cable from the seismometer to the breakout box is connected to "Guralps 1 Input" on the front panel. This corresponds to the set of 3 channels that Caryn stuffed with the new fancy-pants resistors few weeks ago. (When we finally get the other Guralp back from the company, we'll have to stuff the next set of 3 channels).

The Breakout Box is on the very top of 1Y1, sitting on top of the black power supplies. This should be fine, but it's pretty toasty hot up there, so if we find that there are problems with running the box at higher-than-room-temperature, step 1 will be to find a new spot for the box. (I'm not at this time anticipating a problem, but you never know....) Steve put a little foot between the Guralp box and the power supply to get some air circulation.

The ADC Octopus cable that Bob made is connected, and going up through the top of the rack. I am now going on a BNC cable hunt to extend this cable over to the PEM ADC. The PEM ADC is in 1Y7, so I'll need some medium-long BNC cable to get there.

The power cable is also ready to be connected to the rack's +/- 15VDC. I'll talk to Bob about getting this done.

Next step: pick some channels on the PEM ADC, and create them in the .ini files
  1141   Mon Nov 17 16:59:22 2008 JenneUpdatePEMSeismometer hooked up, reading channels on DataViewer
Alberto, Jenne

The Guralp Seismometer is (finally) hooked up to the PEM ADCU. Alberto helped me make channels in the c1pem1 .ini file, which correspond to:

Guralp1 VERT = channel 9 on PEM ADCU = C1:PEM-SEIS_MC1_VERT
Guralp1 NS = channel 10 on PEM ADCU = C1:PEM-SEIS_MC1_NS
Guralp1 EW = channel 11 on PEM ADCU = C1:PEM-SEIS_MC1_EW

We also renamed the Ranger seismometer's channel to C1:PEM-SEIS_MC2_Y from C1:PEM-SEIS_MC1_Y, since tomorrow I'll move the Ranger Seismometer to be underneath MC2's chamber (it's currently sitting somewhere in the middle of the Mode Cleaner).


I can see the VERT and NS channels with dataviewer, but EW looks dead. I need to figure out if this is a bad cable thing, or if the ADC channel is no good, or if something in the box on that channel is no good. All 3 channels were tested and working after all the soldering was completed by Caryn, but something may have come undone while putting the box into its new place in the top of 1Y1. (In dataviewer, it looks like the EW channel is just floating, and not connected to anything.)

Plan of Attack:
* figure out why EW looks dead on Dataviewer
* redo Rana's static Wiener filter analysis, now that we have 2 seismometers (1 Ranger and 1 Guralp)
* work on adaptive Wiener filtering with the Guralp
  1152   Fri Nov 21 16:52:48 2008 JenneUpdatePEMGuralp seismometer's Channel Problems are solved
PROBLEM noticed earlier this week: It looked like one of the seismometer channels (VERT-1) wasn't working, no matter how I put which channel into which input of the PEM ADCU. Watching the channel on Dataviewer, it looked like the ADC was measuring VERT-1 to be zero (actually measuring zero, not digital noise-type zero). I had checked the ADC by putting in a sine wave with a function generator, and saw on Dataviewer the wave I expected, so I knew that I had the correct channel, and that the channel was good.

SOLUTION: This afternoon I took the box out of the rack and opened it up. As soon as I opened it, I saw that I had left something inside the box which was causing the problem. Back when we were measuring the noise of the box, to ensure that it is lower than the ADC's noise, Rana and I had shorted the test points on the input of the VERT-1 channel with a little piece of wire. It turns out that I had closed up the box without remembering to remove the wire.

CONCLUSION of the story: I took out the piece of wire, and now all three seismometer channels (VERT-1, N/S-1, E/W-1) all work, and all detect me jumping around near the BSC. Since the seismometer breakout box reads a differential measurement, and since the input test points were connected, it was indeed measuring zero. Zero equals zero is all well and good, but it's even better now that it's measuring actual seismic motion.
  1153   Fri Nov 21 17:27:47 2008 JenneUpdatePEMGuralp noise measurement
Here is the data from the Guralp Seismometer for the past day or so, before I fixed the VERT-1 channel. The NS and EW show what's going on in the world, and VERT is measuring essentially the noise of the box, through the ADC, in counts.
Attachment 1: guralp_vert_shorted.jpg
guralp_vert_shorted.jpg
  1154   Fri Nov 21 19:47:26 2008 ranaUpdatePEMGuralp noise measurement
and here's the spectra with them connected - from the coherences, it looks like it needs to be rotated by 90 deg.

I'll next rename the channels to fix this so that we get good seismic data over the weekend with the MC.
Attachment 1: a.png
a.png
  1156   Fri Nov 21 21:20:24 2008 ranaUpdatePEMGuralp noise measurement
This is the spectra and time series of the Guralp channels along with the Ranger (MC2). Looks like we could reduce the gain
on the ranger. The Guralp channels run into ADC noise around 40 Hz (which is OK). We'll have to look at the weekday trends
to see if they saturate.
Attachment 1: a.png
a.png
  1165   Mon Dec 1 15:09:27 2008 robUpdatePEMhalf-micron particle count is alarming
  1167   Tue Dec 2 19:18:10 2008 ranaSummaryPEMRanger SS-1
In entry http://dziban.ligo.caltech.edu:40/40m/881 and a follow up from Jenne I put in the Ranger calibration.
Since then, we've reduced the SR560 gain from 200 to 100 so the calibration factor is now:

1e-9 (m/s)/count and then 2 poles at 0 Hz, and a Q~1 zero pair at 1 Hz.
in DTT:
G = 1e-9
p = 0, 0
z = 0.7 0.7
  1168   Tue Dec 2 19:51:32 2008 ranaUpdatePEMhalf-micron particle count is alarming
The 0.5 micron dust monitor count is now pretty high (36000). I wandered around the lab to see if there was anything
nasty going on but I didn't see or smell anything in particular. Since today Alberto was sitting around where the
dust monitor is while aligning the PSL beam, we should blame him. Its either garlic, cologne, or time to bathe.

The 400 day hour trend shows that while the counts are not so unusual, the 40m is dirtier than it was last year.
Attachment 1: Untitled.png
Untitled.png
Attachment 2: dust.png
dust.png
  1171   Wed Dec 3 19:21:09 2008 ranaConfigurationPEMRanger move
I looked at the Ranger signals. Somehow it has a relative transfer function of 'f' between it and the Guralp.
      Ranger
i.e.  ------ ~ f
      Guralp

which is strange since according to their manuals, they should both be giving us a voltage output which is proportional
to velocity. I checked that the Ranger only has a load resistor and then an SR560 low pass at 300 Hz. Jenne assures
me that the Guralp breakout box shouldn't have any poles either (to be double checked). Its a mystery.

We made sure that the SR560 now is DC coupled, G = 100, & 1-pole low pass at 300 Hz. I moved it over next to the Guralp
(went through the mass recentering procedure after forgetting to lock it before moving). It is behaving as it was
before.

Attached is a 2 page PDF of the comparisons. The 'MC1' channels are Guralp and 'MC2' is Ranger.

The second attachment compares our seismometers (in counts) with the LHO Guralp seismometers. There's no high frequency
rolloff there like what we see here so I bet a dollar that there's a pole in the Guralp box somewhere.
Attachment 1: c.pdf
c.pdf c.pdf
Attachment 2: wsnb.pdf
wsnb.pdf
  1172   Wed Dec 3 20:10:09 2008 Jenne, RanaUpdatePEMComparing Wiener subtraction with different seismometers
Attached is a plot of MC_L, and then the residual MC_L after static Wiener filtering, using different combinations of our accelerometers and seismometers.

This is the same type of plot that Rana has included in the past few weeks, using Wiener filters calculated with c1wino.m

This data is from GPS 912312914, duration = 7200 sec, sometime during the night last night.

Unfortunately, it doesn't look like adding the Guralp seismometer to the Accelerometers and the Ranger did much, especially at low frequencies (all sensors = black curve). We'll have to investigate why this is true, and what we can do to get some low-frequency subtraction going on.

In the legend, "Residuals Accels, Guralp, Ranger" implies that the residual has been calculated using all of the sensors listed.
Attachment 1: Dec032008_c1wino_seisCombos.png
Dec032008_c1wino_seisCombos.png
  1173   Wed Dec 3 20:36:07 2008 Jenne, RanaUpdatePEMComparing Wiener subtraction with different seismometers
The Ranger has now been moved over to sit underneath the MC2 tank (it was previously close to the PSL rack). It
is still pointed in the +Y direction (towards ETMY, aka south).

New spectra attached - looks like the coherence is still there between the Guralp and the Ranger which are now
seperated by the MC length (~12 m). At LLO, I have witnessed a coherence of less than 0.3 above 1 Hz for these
distances. Curious.

L960019-00-F describes measurements done at SLAC on seismic coherence. The iLIGO LSC PDD
(http://www.ligo.caltech.edu/docs/T/T970122-00.pdf) discusses in sec 4.2 how this was incorporated into the LSC design.

When we get our next Guralp, it will be interesting to move them around and determine what the cross-spectrum
is between different points in the lab during typical times.

In the second attachment, I have plotted the square of the quantity used in the LSC PDD (S_xy) which I think is
what we now plot in DTT as 'Coherence'.

The third attachment shows the coherences among the TM SUSPOS_INs. I've turned off the oplev servos for this but
the OSEM damping is still on. Its not quite the same as the theory, but we could probably measure/tweak the
seismic velocity and then get better agreement.
Attachment 1: d.pdf
d.pdf
Attachment 2: sco.png
sco.png
Attachment 3: fly.pdf
fly.pdf
  1174   Thu Dec 4 13:49:39 2008 JenneUpdatePEMMore of: Comparing Wiener subtraction using different sensors
Here is another version of the same type plot I put in the elog yesterday. This plot is looking at the 7200 seconds after 04Dec2008 08:45:00 UTC. This time was last night, when there was no crazy seismic activity, and well after the Ranger seismometer was moved to its new place under MC2.

This plot includes all possible combinations of the accelerometers, Guralp seismometer and Ranger seismometer (taking all 6 accelerometers as a set, and all 3 Guralp channels as a set). It is good to see that for the set of traces which do not include the accelerometers - brown, dark green and light blue - the subtraction at higher frequencies isn't all that great. Thus, the accelerometers are doing their job, and work well with the Wiener subtraction.

Still under investigation is why we don't see a whole lot of improvement at low frequency.
Attachment 1: Dec042008_c1wino_seisCombos.png
Dec042008_c1wino_seisCombos.png
  1219   Fri Jan 9 09:14:17 2009 steveUpdatePEMifo is recovered after eq
There is no obvious damage from last night earth quake.
All sus dampings were turned on, MC locked and the arms locked right on
Attachment 1: eq4.5sb.jpg
eq4.5sb.jpg
  1337   Wed Feb 25 11:48:02 2009 JenneUpdatePEMWiener filtering update - work on filtering some S5 DARM_CTRL data

Quick update on my wiener filtering status:

Joe has been helping me get on the GRID, so I now have  a grid certificate, and accounts on most/all of the clusters.

Joe also helped me get menkar to get S5 data so that I can do wiener filtering to the back-data. 

 

I've been running the wiener filtering algorithm, and right now, it doesn't do anything to improve the DARM_CTRL data.  I am confident that this is because something is funky in the wiener filtering algorithm somewhere.  The indicator of this is that the wiener filtering calculation takes the same amount of time (~95 seconds) to calculate a filter for 64 seconds of data as for 1 hour of data (both for N = 2000 taps). 

 

For reference, attached are my plots for the wiener filtering result for (1) 64 seconds of S5 data, and for (2) 3600 seconds of S5 data.

These plots were made using H1:DARM_CTRL as the signal to minimize, with 4 seismometers as the witness channels (EX_SEISX, EY_SEISY, LVEA_SEISX, LVEA_SEISY)

 

I'm working on figuring out what's going on with the filtering algorithm, and why it does work for C1:MC_L minimization, but does not work for H1:DARM_CTRL minimization.

 

 

 

Attachment 1: h1_DARM_64s_4seis_25Feb09.png
h1_DARM_64s_4seis_25Feb09.png
Attachment 2: h1_DARM_3600s_4seis_25Feb09.png
h1_DARM_3600s_4seis_25Feb09.png
  1347   Tue Mar 3 08:44:31 2009 steveUpdatePEMair cond. maintenance today
IFO room 104 air conditions will be shut down for maintenance today.
This should be finished by noon.
The temperature and particle count variation can be more than usual.
  1411   Fri Mar 20 11:01:02 2009 steveUpdatePEMparticle counts are high

The outside particle counts for 0.5 micron are 3 million this morning at 9am. Low clouds, foggy condition with low inversion layer.

This makes the 40m lab 30-50K

I just turned on the HEPA filter at the PSL enclosure.

Please, leave it on high

 

Attachment 1: particles32d.jpg
particles32d.jpg
  1447   Tue Mar 31 09:42:32 2009 steveUpdatePEMETMY sus damping restored again

The Caltech gasoline storage tank is being upgraded.

They are jack hammering and digging with bulldozer 50 yards south of  ETMY

  1453   Fri Apr 3 14:52:38 2009 JenneOmnistructurePEMGuralp is finally back!

After many, many "it'll be there in 2 weeks" from the Guralp people, our seismometer is finally back!

I have it plugged into the Guralp breakout box's Channel 1xyz (so I have unplugged the other Guralp).  Both of the Guralp's are currently sitting under the MC1/MC3 chamber.

Before we can have both Guralps up and running, I need to stuff the next 3 channels of the breakout box (back in the fall, I only had Caryn do 1x, 1y, 1z, and now I need 2x, 2y and 2z done with the fancy low-noise resistors), so all the gains match between the 2 sets of channels.

I'm leaving the new Guralp plugged in so we can see how it behaves for the next couple days, until I take out the breakout box for stuffing.

  1455   Mon Apr 6 19:09:15 2009 JenneUpdatePEMOld Guralp is hooked back up to the ADC

Old Guralp is hooked back up, the new one is sitting next to it, disconnected for now.

  1528   Tue Apr 28 12:55:57 2009 CarynDAQPEMUnplugged Guralp channels

For the purpose of testing out the temperature sensors, I stole the PEM-SEIS_MC1X,Y,Z channels.

I unplugged Guralp NS1b, Guralp Vert1b, Guralp EW1b cables from the PEM ADCU(#10,#11,#12) near 1Y7 and put temp sensors in their place (temporarily).

  1540   Sat May 2 16:34:31 2009 carynDAQPEMGuralp channels plugged back in

I plugged the Guralp cables back into the PEM ADCU

       Guralp NS1b ---> #11

       Guralp Vert1b --->#10

       Guralp EW1b --->#12

  1546   Tue May 5 09:22:46 2009 carynUpdatePEMzeros

For several of the channels on the PEM ADCU, zeros are occuring at the same time. Does anyone know why that might happen or how to fix it?

Attachment 1: zerotest2.png
zerotest2.png
Attachment 2: zerotest.png
zerotest.png
  1566   Fri May 8 16:03:31 2009 JenneUpdatePEMUpdate on Jenne's Filtering Stuff

To include the plots that I've been working on in some form other than on my computer, here they are:

First is the big surface plot of all the amplitude spectra, taken in 10min intervals on one month of S5 data. The times when the IFO is unlocked are represented by vertical black stripes (white was way too distracting).  For the paper, I need to recreate this plot, with traces only at selected times (once or twice a week) so that it's not so overwhelmingly large.  But it's pretty cool to look at as-is.

Second is the same information, encoded in a pseudo-BLRMS.  (Pseudo on the RMS part - I don't ever actually take the RMS of the spectra, although perhaps I should).  I've split the data from the surface plot into bands (The same set of bands that we use for the DMF stuff, since those seem like reasonable seismic bands), and integrated under the spectra for each band, at each time.  i.e. one power spectra gives me 5 data points for the BLRMS - one in each band.  This lets us see how good the filter is doing at different times.

At the lower frequencies, after ~25 days, the floor starts to pick up.  So perhaps that's about the end of how long we can use a given Wiener filter for.  Maybe we have to recalculate them about every 3 weeks.  That wouldn't be tragic. 

I don't really know what the crazy big peak in the 0.1-0.3Hz plot is (it's the big yellow blob in the surface plot).  It is there for ~2 days, and it seems awfully symmetric about it's local peak.  I have not yet correlated my peaks to high-seismic times in the H1 elog.  Clearly that's on the immediate todo list. 

Also perhaps on the todo list is to indicate in some way (analagous to the black stripes in the surface plot) times when the data in the band-limited plot is just extrapolated, connecting the dots between 2 valid data points.

 

A few other thoughts:  The time chosen for the training of the filter for these plots is 6:40pm-7:40pm PDT on Sept 9, 2007 (which was a Sunday night).  I need to try training the filter on a more seismically-active time, to see if that helps reduce the diurnal oscillations at high frequency.  If that doesn't do it, then perhaps having a "weekday filter" and an "offpeak" filter would be a good idea.  I'll have to investigate.

Attachment 1: H1S5OneMonthWienerCompBLACK.png
H1S5OneMonthWienerCompBLACK.png
Attachment 2: H1S5BandLimitedTimePlot.png
H1S5BandLimitedTimePlot.png
  1571   Sun May 10 13:34:32 2009 carynUpdatePEMUnplugged Guralp channels

I unplugged Guralp EW1b and Guralp Vert1b and plugged in temp sensors temporarily. Guralp NS1b is still plugged in.

  1597   Mon May 18 01:54:35 2009 ranaUpdatePEMUnplugged Guralp channels
To see if Caryn's data dropouts were happening, I looked at a trend of all of our temperature channels. Looks OK now.

Although you can't see it because I zoomed in, there's a ~24 hour relaxation happening before Caryn's sensors equilibrate.
I guess that's the insulating action of the cooler? We need a picture of the cooler in the elog for posterity.
Attachment 1: Untitled.png
Untitled.png
  1599   Mon May 18 10:06:56 2009 carynUpdatePEMTemp sensor

Quote:
To see if Caryn's data dropouts were happening, I looked at a trend of all of our temperature channels. Looks OK now.

Although you can't see it because I zoomed in, there's a ~24 hour relaxation happening before Caryn's sensors equilibrate.
I guess that's the insulating action of the cooler? We need a picture of the cooler in the elog for posterity.[/quote


Dropouts can't been seen with a minute trend, only a second trend. No big deal, but they are still occurring. See plot below.

The 24hr relaxation period is due to the cooler and some metal blocks that were cooled in the freezer and then put in the cooler to see if the relationship between the temp sensors changed with temperature. The relationship is not linear, which probably means there is some non-linearity in each temperature sensor's relationship to temperature. So, when calibrating them with Bob's temp sensor, more than 2 data points need to be collected.

Picture of cooler for posterity is attached
Attachment 1: datadropout.png
datadropout.png
Attachment 2: coolerpic1.jpg
coolerpic1.jpg
Attachment 3: coolerpic2.jpg
coolerpic2.jpg
  1600   Mon May 18 15:31:11 2009 ranaUpdatePEMTemp sensor

Quote:
Picture of cooler for posterity is attached


I'm puzzled as to why the minute trend doesn't pick this up; its clearly there in the full data.

Looks like its several samples too. Can someone please reboot this DCU and see if the problem goes away?
Attachment 1: Untitled.png
Untitled.png
  1606   Tue May 19 15:54:29 2009 JenneUpdatePEMMore Plots for the S5 H1:DARM Wiener Filtering....

Even more plots for the Wiener filtering!

We have a set of spectrograms, which show (in color) the amplitude spectrum, at various times during a one month stretch of time, during S5. Each vertical data-'stripe' is 10min long.

We also have a set of band-limited plots, which take the spectra at each time, and integrate under it, for different frequency bands.

Each set of plots has the following 3 plots:  The raw DARM spectrum, a ratio of residual/raw, and the residuals, normalized to the first one (on which the wiener filter was trained).

The residuals are the DARM spectrum, after subtracting the Wiener-filtered seismometer witness data.

 

From the ratio plots, it looks like the wiener filter is pretty much equally effective at the time on which the filter was trained, as one month later.  Static filters may be okey-dokey for a long period of time with for the seismic stuff.

Attachment 1: H1darmCompSpecgramRAW.png
H1darmCompSpecgramRAW.png
Attachment 2: H1darmCompSpecgramRATIO.png
H1darmCompSpecgramRATIO.png
Attachment 3: H1darmCompSpecgramRESIDUALS.png
H1darmCompSpecgramRESIDUALS.png
Attachment 4: H1darmCompWienerRAW.png
H1darmCompWienerRAW.png
Attachment 5: H1darmCompWienerRATIO.png
H1darmCompWienerRATIO.png
Attachment 6: H1darmCompWienerRESIDUALS.png
H1darmCompWienerRESIDUALS.png
  1612   Wed May 20 09:55:18 2009 steveUpdatePEMoplev servos turned off

All oplevs servos turned off to protect our suspentions from vibration due to drilling and pounding in CES high bay area.

This activity will be done from 10 am till 3 pm today.

 

Meanwhile our IFO-air conditions are turned off for maintenance.

Their performance of 6 months is shown on plot.

Attachment 1: ACpart6m.jpg
ACpart6m.jpg
  1624   Mon May 25 21:31:47 2009 carynUpdatePEMplugged in Guralp channels

Guralp Vert1b and Guralp EW1b are plugged back in to PEM ADCU #10 and #12 respectively. Guralp NS1b remains plugged in. So,  PEM-SEIS_MC1_X,Y,Z should now corrsp to seismometer as before.

  1647   Wed Jun 3 11:28:01 2009 carynUpdatePEMUnplugged Guralp channels
  1648   Wed Jun 3 12:31:13 2009 carynUpdatePEMplugged in guralp channels
  1653   Thu Jun 4 23:39:23 2009 peteUpdatePEM5 days, 20 days of accelerometers

Looks like yesterday was particularly noisy.  It's unclear to me why diurnal variation much more visible in MC1_Y, and why the floor wanders.

 

The first plot shows 5 days.  The second plot shows 20 days.

Attachment 1: acc_5day.png
acc_5day.png
Attachment 2: acc_20days.png
acc_20days.png
  1665   Wed Jun 10 09:06:12 2009 steveUpdatePEMparticle counts and turbulance

I moved the mobile HEPA filter from ITMX's north door to ITMX-ISCT and covered it up with a merostate tent to accommodate the aluminum foil particle measurement on June 5

It lowered the 40m baseline counts by about a factor of 3 of 0.5 micron and a factor of 2 of 1.0 micron.

The HEPA filter is sweeping the floor and blowing the particles upwards. The MET ONE counter is on the top of the IOOC looking south at ~75 degrees upward.

 

Attachment 1: parturbulnc.jpg
parturbulnc.jpg
  1692   Tue Jun 23 23:14:36 2009 ClaraConfigurationPEMAccelerometers relocated

Both accelerometers have been moved in an attempt to optimize their positions. The MC1 accelerometer was moved from one green bar to the other (I don't know what to call them) at the base of the MC1 and MC3 chambers. That area is pretty tight, as there is an optical table right there, and I did my best to be careful, but if you suspect something has been knocked loose, you might check in that area. The MC2 accelerometer was moved from the horizontal bar down to the metal table on which the MC2 chamber rests.

ELOG V3.1.3-