Note
Click here to download the full example code
Short transient search examples: Make data
An example to generate data with a short transient signal.
This can be run either stand-alone (will just generate SFT files and nothing else); or it is also being imported from PyFstat_example_short_transient_grid_search.py and PyFstat_example_short_transient_MCMC_search.py
14 import os
15
16 import pyfstat
17
18 outdir = os.path.join("PyFstat_example_data", "PyFstat_example_short_transient_search")
19
20 F0 = 30.0
21 F1 = -1e-10
22 F2 = 0
23 Alpha = 0.5
24 Delta = 1
25
26 tstart = 1000000000
27 duration = 2 * 86400
28
29 transient_tstart = tstart + 0.25 * duration
30 transient_duration = 0.5 * duration
31 transientWindowType = "rect"
32 tref = tstart
33
34 h0 = 1e-23
35 cosi = 0
36 psi = 0
37 phi = 0
38 sqrtSX = 1e-22
39 detectors = "H1,L1"
40
41 Tsft = 1800
42
43 if __name__ == "__main__":
44
45 transient = pyfstat.Writer(
46 label="simulated_transient_signal",
47 outdir=outdir,
48 tref=tref,
49 tstart=tstart,
50 duration=duration,
51 F0=F0,
52 F1=F1,
53 F2=F2,
54 Alpha=Alpha,
55 Delta=Delta,
56 h0=h0,
57 cosi=cosi,
58 detectors=detectors,
59 sqrtSX=sqrtSX,
60 transientStartTime=transient_tstart,
61 transientTau=transient_duration,
62 transientWindowType=transientWindowType,
63 Tsft=Tsft,
64 Band=0.1,
65 )
66 transient.make_data()
67 print("Predicted 2F from injection Writer:", transient.predict_fstat())
Total running time of the script: ( 0 minutes 0.000 seconds)