From d2b554365dd1e8e04356096882662b7291b8ff78 Mon Sep 17 00:00:00 2001
From: pinyili2 <pinyili2@illinois.edu>
Date: Fri, 4 Oct 2024 18:43:44 -0500
Subject: [PATCH] set angle to 180
---
mrdna/readers/segmentmodel_from_cadnano.py | 1 +
1 file changed, 1 insertion(+)
diff --git a/mrdna/readers/segmentmodel_from_cadnano.py b/mrdna/readers/segmentmodel_from_cadnano.py
index 2936c4d..7633620 100644
--- a/mrdna/readers/segmentmodel_from_cadnano.py
+++ b/mrdna/readers/segmentmodel_from_cadnano.py
@@ -105,6 +105,7 @@ def get_helix_angle(part, helix_id, indices, index_offset=0, fwd=True):
twist_per_base = tpr*360./bpr
# angle = eulerZ - twist_per_base*indices + 0.5*mgroove + 180
angle = eulerZ - twist_per_base*indices - 0.5*mgroove
+ angle=angle*(np.pi/180)
cos,sin = (np.cos(angle), np.sin(angle))
R = np.zeros( (len(angle),3,3) )
R[:,0,0] = cos
--
GitLab