From a999c0cfdcb4d7dcf0782d31ae0d36b1f597fc86 Mon Sep 17 00:00:00 2001
From: Yifan Wang <yifan.wang@aei.mpg.de>
Date: Mon, 31 May 2021 16:47:53 +0200
Subject: [PATCH] add working waveform for parity violation
---
.gitignore | 2 ++
birefringence-waveform/__init__.py | 1 +
.../birefringence_waveform.py | 2 +-
birefringence-waveform/waveform.py | 17 -----------------
.../plot_waveform.ipynb | 0
.../printwfname.py | 0
6 files changed, 4 insertions(+), 18 deletions(-)
create mode 100644 .gitignore
create mode 100644 birefringence-waveform/__init__.py
delete mode 100644 birefringence-waveform/waveform.py
rename {birefringence-waveform => utils}/plot_waveform.ipynb (100%)
rename {birefringence-waveform => utils}/printwfname.py (100%)
diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..eeb514c
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1,2 @@
+__pycache__
+.ipynb_checkpoints
diff --git a/birefringence-waveform/__init__.py b/birefringence-waveform/__init__.py
new file mode 100644
index 0000000..287ff9f
--- /dev/null
+++ b/birefringence-waveform/__init__.py
@@ -0,0 +1 @@
+import .birefringence_waveform
diff --git a/birefringence-waveform/birefringence_waveform.py b/birefringence-waveform/birefringence_waveform.py
index b9d6e92..bc2bd7d 100644
--- a/birefringence-waveform/birefringence_waveform.py
+++ b/birefringence-waveform/birefringence_waveform.py
@@ -18,4 +18,4 @@ def IMRPhenomParity(parity_Amu, **kwds):
return hp_parity, hc_parity
-pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)
\ No newline at end of file
+pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)
diff --git a/birefringence-waveform/waveform.py b/birefringence-waveform/waveform.py
deleted file mode 100644
index 41ec79b..0000000
--- a/birefringence-waveform/waveform.py
+++ /dev/null
@@ -1,17 +0,0 @@
-from pycbc.waveform import get_fd_waveform
-import numpy as np
-import pycbc.conversions
-import lal
-
-def IMRPhenomParity(parity_Amu, **kwds):
-
- if 'approximant' in kwds:
- kwds.pop("approximant")
- hp, hc = get_fd_waveform(approximant="IMRPhenomXPHM", **kwds)
-
- temp = parity_Amu * 1e9 * lal.QE_SI * lal.PI * lal.PI / lal.H0_SI
-
- hp_parity = hp + hc * temp * hp.sample_frequencies **2
- hc_parity = hc - hp * temp * hp.sample_frequencies **2
-
- return hp_parity, hc_parity
diff --git a/birefringence-waveform/plot_waveform.ipynb b/utils/plot_waveform.ipynb
similarity index 100%
rename from birefringence-waveform/plot_waveform.ipynb
rename to utils/plot_waveform.ipynb
diff --git a/birefringence-waveform/printwfname.py b/utils/printwfname.py
similarity index 100%
rename from birefringence-waveform/printwfname.py
rename to utils/printwfname.py
--
GitLab