You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
*[challenge3_2048hz.gwf](https://dcc.ligo.org/public/0187/G2300818/001/challenge3_2048hz.gwf) <-- Downsampled version of Challenge 3 file
11
11
12
-
Workshop participants may [submit solutions via thinkific](https://learn.gwosc.org) as individuals or in teams of up to 3 people.
13
-
Check for the deadline date on the Thinkific website.
12
+
Workshop participants may [submit solutions via the online course](https://learn.gwosc.org/courses/odw2026) as individuals or in teams of up to 3 people (check for the deadline date on the website).
14
13
15
14
Challenges are ordered by difficulty. Entries will be rewarded a number of
16
15
points that scales with the difficulty of the challenge.
@@ -19,78 +18,88 @@ Good luck to all!
19
18
20
19
## Challenge 1 (1 point) -- Novice
21
20
22
-
Identify a loud binary black hole signal in white, Gaussian noise.
21
+
Identify a loud binary black hole signal in white, Gaussian noise:
23
22
24
23
* Use the data file `challenge1.gwf`. The channel name is `H1:CHALLENGE1`.
25
24
* The data are white, Gaussian noise containing a simulated BBH signal.
26
25
27
26
Instructions:
28
-
1. Load the data into memory. What are the sampling rate and duration of the data?
29
-
2. Plot the data in the time-domain.
30
-
3. Plot a spectrogram (or q-transform) of the data, and try to identify the signal.
31
27
32
-
Question:
33
-
* What is the time of the merger (in seconds from the start of the segment)?
28
+
* Load the data into memory. What are the sampling rate and duration of the data?
29
+
* Plot the data in the time-domain.
30
+
* Plot a spectrogram (or Q-transform) of the data, and try to identify the signal.
31
+
32
+
Question for Challenge 1:
33
+
34
+
1. What is the time of the merger (in seconds from the start of the segment)?
34
35
35
36
## Challenge 2 (2 points) -- Rookie
36
37
37
-
Signal in colored, Gaussian noise.
38
+
Signal in colored, Gaussian noise:
38
39
39
40
* Use the data file `challenge2.gwf`, with channel name `H1:CHALLENGE2`.
40
41
* The data contain a BBH signal with m1 = m2 = 30 solar masses, spin = 0.
41
42
42
-
Question A:
43
-
44
-
* What is the approximative time of the merger? (Hint: a plot of the q-transform could help)
45
-
46
-
Next steps:
47
-
1. Generate a time-domain template waveform using approximate `SEOBNRv4_opt`.
48
-
with the same parameters as above. Plot this waveform.
49
-
2. Calculate a PSD of the data, and plot this on a log-log scale.
50
-
Use axes ranging from 20 Hz up to the Nyquist frequency.
51
-
3. Use the template waveform and PSD to calculate the SNR time series. Plot the SNR time-series.
43
+
Instructions:
52
44
53
-
Question B:
45
+
* Plot a Q-transform of the data.
46
+
* Generate a time-domain template waveform using approximate `SEOBNRv4_opt`.
47
+
with the same parameters as above. Plot this waveform.
48
+
* Calculate a PSD of the data, and plot this on a log-log scale.
49
+
Use axes ranging from 20 Hz up to the Nyquist frequency.
50
+
* Use the template waveform and PSD to calculate the SNR time series. Plot the SNR time-series.
54
51
55
-
* What is the matched filter SNR of the signal?
52
+
Questions for Challenge 2:
56
53
54
+
2. From the Q-transform, what is the approximative time of the merger?
55
+
3. What is the matched filter SNR of the signal?
57
56
58
57
## Challenge 3 (4 points) -- Intermediate
59
58
59
+
Search for a loud event in realistic data:
60
+
60
61
* Use the data file `challenge3.gwf` with channel `H1:CHALLENGE3`.
61
62
* These are real LIGO data from O2, though we've adjusted the time labels and
62
63
added some simulated signals.
63
64
* The data contain a loud simulated signal with m1 = m2 = 10 solar masses.
64
65
65
-
Questions:
66
+
Instructions:
67
+
68
+
* Use the template waveform `SEOBNRv4_opt` and PSD to calculate the SNR time series.
66
69
67
-
A. What is the merger time of this signal?
70
+
Questions for challenge 3:
68
71
69
-
B. What is the matched-filter SNR of this signal?
72
+
4. What is the merger time of the loud signal?
73
+
5. What is the matched-filter SNR of the loud signal?
70
74
71
75
**Note:** it's not needed to use bilby for this challenge.
72
76
73
77
## Challenge 4 (8 points) -- Advanced
74
78
75
-
* Use the data file `challenge3.gwf` with channels `H1:CHALLENGE3` and `L1:CHALLENGE3`.
79
+
Realistic search and parameter estimation:
80
+
81
+
* We will use the same realistic data set `challenge3.gwf` with channels `H1:CHALLENGE3` and `L1:CHALLENGE3`.
76
82
* These are real LIGO data from O2, though we've adjusted the time labels and
77
83
added some simulated signals.
78
-
* Any simulated signals have been added to both the H1 and L1 data
79
-
* All simulated signals have spin = 0 and m1 = m2, with m1 somewhere in the range 10-50 solar masses
84
+
* All simulated signals have been added to both the H1 and L1 data.
85
+
* All simulated signals have spin = 0 and m1 = m2, with m1 somewhere in the range 10-50 solar masses.
86
+
* Watch out! These are real data, and so glitches may be present.
80
87
81
88
Questions:
82
-
A. Identify as many signals as you can. Watch out! These are real data, and so glitches may be
83
-
present. Any correct detection is +1 point but any false alarms will count -1 point
89
+
90
+
6. Identify as many signals as you can. Any correct detection is +1 point but any false alarms will count -1 point
84
91
against your score. For each signal you find, list:
85
92
86
93
* The merger time
87
94
* The SNR
88
95
* Your estimate of the component masses
89
96
90
-
B. Identify as many glitches as you can. Make a spectrogram of each one.
97
+
7. Identify as many glitches as you can. Make a spectrogram of each one.
98
+
8. For the earliest event you found, use bilby to compute the posterior distribution for the mass.
99
+
Hint:
91
100
92
-
C. For the earliest event you found, use bilby to compute the posterior distribution for the mass.
93
-
Fix the spin and mass ratio to make this run faster.
101
+
- Fix the spin and mass ratio to make this run faster.
102
+
- You can also fix other angular parameters to zero, while leaving at least the chirp mass, sky position, time and distance free. This will still give an approximately correct posterior for the chirp mass.
0 commit comments