Skip to content
GitLab
Menu
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
lscsoft
bilby
Commits
44809982
Commit
44809982
authored
May 10, 2018
by
Colm Talbot
Browse files
Merge branch 'master' of git.ligo.org:Monash/tupak
parents
f3afa08b
ef6e8151
Pipeline
#18653
passed with stages
in 2 minutes and 35 seconds
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
tutorials/create_your_own_source_model.py
0 → 100644
View file @
44809982
#!/bin/python
"""
"""
from
__future__
import
division
,
print_function
import
tupak
import
numpy
as
np
tupak
.
utils
.
setup_logger
()
def
sine_gaussian
(
f
,
A
,
f0
,
tau
,
phi0
,
geocent_time
,
ra
,
dec
,
psi
):
arg
=
-
(
np
.
pi
*
tau
*
(
f
-
f0
))
**
2
+
1j
*
phi0
plus
=
np
.
sqrt
(
np
.
pi
)
*
A
*
tau
*
np
.
exp
(
arg
)
/
2.
cross
=
plus
*
np
.
exp
(
1j
*
np
.
pi
/
2
)
return
{
'plus'
:
plus
,
'cross'
:
cross
}
outdir
=
'outdir'
label
=
'GW150914_sine_gaussian'
time_of_event
=
1126259462.422
H1
=
tupak
.
detector
.
get_interferometer
(
'H1'
,
time_of_event
,
version
=
1
,
outdir
=
outdir
)
L1
=
tupak
.
detector
.
get_interferometer
(
'L1'
,
time_of_event
,
version
=
1
,
outdir
=
outdir
)
interferometers
=
[
H1
,
L1
]
prior
=
dict
()
prior
[
'A'
]
=
tupak
.
prior
.
Uniform
(
0
,
1e-20
,
'A'
)
prior
[
'f0'
]
=
tupak
.
prior
.
Uniform
(
0
,
10
,
'f'
)
prior
[
'tau'
]
=
tupak
.
prior
.
Uniform
(
0
,
10
,
'tau'
)
prior
[
'geocent_time'
]
=
tupak
.
prior
.
Uniform
(
time_of_event
-
0.1
,
time_of_event
+
0.1
,
'geocent_time'
)
prior
[
'phi0'
]
=
0
#tupak.prior.Uniform(0, 2*np.pi, 'phi')
prior
[
'ra'
]
=
0
prior
[
'dec'
]
=
0
prior
[
'psi'
]
=
0
waveform_generator
=
tupak
.
waveform_generator
.
WaveformGenerator
(
sine_gaussian
,
H1
.
sampling_frequency
,
H1
.
duration
)
likelihood
=
tupak
.
likelihood
.
Likelihood
(
interferometers
,
waveform_generator
)
result
=
tupak
.
sampler
.
run_sampler
(
likelihood
,
prior
,
sampler
=
'pymultinest'
,
outdir
=
outdir
,
label
=
label
,
resume
=
False
)
result
.
plot_walks
()
result
.
plot_distributions
()
result
.
plot_corner
()
print
(
result
)
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment