A statistical veto method employing a back-coupling consistency - - PowerPoint PPT Presentation

a statistical veto method employing a back coupling
SMART_READER_LITE
LIVE PREVIEW

A statistical veto method employing a back-coupling consistency - - PowerPoint PPT Presentation

A statistical veto method employing a back-coupling consistency check Stefan Hild, P. Ajith and M. Hewitsion (AEI Hannover) LIGO-G060641-00-Z Stefan Hild 1 GWDAW11, Potsdam, December 2006 Standard statistical veto Noise couples into


slide-1
SLIDE 1

Stefan Hild 1 GWDAW11, Potsdam, December 2006

A statistical veto method employing a back-coupling consistency check

Stefan Hild, P. Ajith and M. Hewitsion (AEI Hannover)

LIGO-G060641-00-Z

slide-2
SLIDE 2

Stefan Hild 2 GWDAW11, Potsdam, December 2006

Standard statistical veto

  • Noise couples into both: H and X
  • Events in H are partly correlated

with events in X.

  • Veto condition: Events in H and

X occure at the same time If there is any GW-signal in X => high false veto rate Standard statistical veto works fine

  • nly for GW-free veto channels, like

microphones or magnetometers

slide-3
SLIDE 3

Stefan Hild 3 GWDAW11, Potsdam, December 2006

Veto channels containing traces of GW-signal

Unfortunately many promissing veto channels may contain traces of GW-signal, for example Interferometer signals (light powers, control signals, ...) Two populations of coincident events:

  • Events originating from noise (we want to veto)
  • GW-like events coupling back to X (we DON‘T want to veto)
slide-4
SLIDE 4

Stefan Hild 4 GWDAW11, Potsdam, December 2006

Seperate two populations by ampli- tude ratio of the coincicent events

If event X(j) originates from the event H(i) their amplitude ratio has to correspond to the transfer function for back-coupling: In order to get a safe veto method we have to compare amplitude ratio of the two coincident events with the back-coupling transfer function: If H(i) is not vetoed If H(i) gets vetoed !

slide-5
SLIDE 5

Stefan Hild 5 GWDAW11, Potsdam, December 2006

Real world scenario

In reality we have to allow for some inaccuracies:

  • Error in the amplitude estimation of the two

events

  • Error in back-coupling transfer function

(measurement, non stationarity) Allow for overall error VETO CONDITION Two coincident events H(i) and X(j) are vetoed in the case that the amplitude ratio matches one of these requirements:

slide-6
SLIDE 6

Stefan Hild 6 GWDAW11, Potsdam, December 2006

Dust falling through main output beam

1 2 3 4 5 6 7 8 200 400 600 800 1000 1200 1400 1600 1800 2000

Time from 2006-05-09 14:59:46 (831222000) (h) Frequency (Hz)

1719 events from DER_DATA_H 916 events from LSC_MID_VIS 1245 LSC_MID_VIS events coinc with DER_DATA_H 1 2 3 4 5 6 7 8 200 400 600 800 1000 1200 1400 1600 1800 2000

Time from 2006-06-28 22:59:46 (835570800) (h) Frequency (Hz)

1054 events from DER_DATA_H 102 events from LSC_MID_VIS 49 LSC_MID_VIS events coinc with DER_DATA_H

high dust concentration (broken AC) low dust concentration

Time coincidence window = 10ms Time coincidence window = 10ms

When dust is falling through the main output beam, coincidence glitches are induced to H and PDC.

slide-7
SLIDE 7

Stefan Hild 7 GWDAW11, Potsdam, December 2006

PDC contains traces of GW-signal

What is PDC ? It is the DC light from the main dark port photo detector. It contains traces of GW-signal. Hardware injections of sinusoidal signals show coherence of 1.

slide-8
SLIDE 8

Stefan Hild 8 GWDAW11, Potsdam, December 2006

Determine back-coupling transfer function

10

2

10

3

10

14

10

15

Frequency [Hz] Amplitude ratio

αback

Injecting differential arm length noise (to mimic the effect of a GW) and then measure transfer function from H to PDC ?

slide-9
SLIDE 9

Stefan Hild 9 GWDAW11, Potsdam, December 2006

Sine-Gaussian hardware injections

10

2

10

3

10

14

10

15

Frequency [Hz] Amplitude ratio

αback

aX/aH from hardware injections

277 injections detected in H => 14 Injections also detected in PDC The injections found in PDC match the back-coupling transfer function. Injecting sine-Gaussians into differential arm length servo.

slide-10
SLIDE 10

Stefan Hild 10 GWDAW11, Potsdam, December 2006

Determine overall error

Need to determine !! 1. Back-coupling TF was measured to vary less than +/-50% over months. 2. Maximum error in amplitude estimation of mHACR using 3 sigma gives 60% for events of SNR = 4 (sine-Gaussian injections into Gaussian noise) 1. For the real data we will allow for 200% error in amplitude estimation.

slide-11
SLIDE 11

Stefan Hild 11 GWDAW11, Potsdam, December 2006

Application of a statistical veto employing a back-coupling consistency check

Application to two data sets of GEO S5 data:

  • Data Set 1: Full September 2006 (low dust concentration)
  • Data Set 2: 8 hours from May 2006 (high dust concentration)

Final set of three veto conditions: Time coincidence Frequency coincidence Amplitude cut (checking that the ratio is not consistent with back-coupling)

slide-12
SLIDE 12

Stefan Hild 12 GWDAW11, Potsdam, December 2006

Data set 1

10

2

10

3

10

14

10

15

10

16

10

17

10

18

Frequency [Hz] Amplitude ratio Data set 1: Full September 2006 aX/aH Used amplitude cut

αback

aX/aH form hardware injections

slide-13
SLIDE 13

Stefan Hild 13 GWDAW11, Potsdam, December 2006

Data set 2

10

2

10

3

10

14

10

15

10

16

10

17

10

18

Frequency [Hz] Amplitude ratio Data set 2: 8 hours from May 2006 aX/aH Used amplitude cut

αback

aX/aH from hardware injections

slide-14
SLIDE 14

Stefan Hild 14 GWDAW11, Potsdam, December 2006

Summary of the Veto Performance

  • 80
  • 60
  • 40
  • 20
20 40 60 80 10 10 1 10 2 10 3 10 4 Vetoed events Time shift [sec]
  • 80
  • 60
  • 40
  • 20
20 40 60 80 10
  • 1
10 10 1 10 2 10 3 Vetoed events Time shift [sec]

Data set 1: Full September 2006 Data set 2: 8 hours of May 2006

slide-15
SLIDE 15

Stefan Hild 15 GWDAW11, Potsdam, December 2006

Summary

  • We developed a method for safe statistical vetoes

using interferometer channels (potentially containing traces of GW-signal).

  • This method employs an additional back-coupling

consistency check.

  • Application to GEO S5 data showed a good

performance.

  • The method is generally applicable.
slide-16
SLIDE 16

Stefan Hild 16 GWDAW11, Potsdam, December 2006

E N D

slide-17
SLIDE 17

Stefan Hild 17 GWDAW11, Potsdam, December 2006

Full Veto pipeline used for Data Set 1

slide-18
SLIDE 18

Stefan Hild 18 GWDAW11, Potsdam, December 2006

Example from GEO600: Mains monitor

Application of a single co- incidence window for time: Application of a multi coincidence window for time (6ms) and frequency: Efficiency to Background ratio (Significance) improved !