Fix #2582: polymer persistence length example fix

Added name to AUTHOR file
MDAnalysis · Mar 7, 2020 · 6226217 · 6226217
1 parent 1f4e492
commit 6226217
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 5 deletions.
diff --git a/package/AUTHORS b/package/AUTHORS
@@ -133,6 +133,7 @@ Chronological list of authors
   - Hao Tian
   - Hugo MacDermott-Opeskin
   - Anshul Angaria
+  - Shubham Sharma
 
 
 

diff --git a/package/CHANGELOG b/package/CHANGELOG
@@ -58,6 +58,7 @@ Fixes
     argument. The directives parsed into bonds, angles, impropers, and dihedrals now
     match TPRParser. (PR #2408)
   * Added parmed to setup.py
+  * Fixed example docs for polymer persistence length (#2582)
 
 Enhancements
   * New analysis.hole2 module for HOLE interfacing. Contains a function (hole) 

diff --git a/package/MDAnalysis/analysis/polymer.py b/package/MDAnalysis/analysis/polymer.py
@@ -48,16 +48,16 @@
   # so we can select the chains by using the .fragments attribute
   chains = u.atoms.fragments
   # select only the backbone atoms for each chain
-  backbones = [chain.select_atoms('not name O* H') for chain in chains]
+  backbones = [chain.select_atoms('not name O* H*') for chain in chains]
   # sort the chains, removing any non-backbone atoms
   sorted_backbones = [polymer.sort_backbone(bb) for bb in backbones]
-  lp = polymer.PersistenceLength(sorted_backbones)
+  persistence_length = polymer.PersistenceLength(sorted_backbones)
   # Run the analysis, this will average over all polymer chains
   # and all timesteps in trajectory
-  lp = lp.run()
-  print('The persistence length is: {}'.format(lp.pl))
+  persistence_length = persistence_length.run()
+  print('The persistence length is: {}'.format(persistence_length.lp))
   # always check the visualisation of this:
-  lp.plot()
+  persistence_length.plot()
 
 """
 from __future__ import division, absolute_import