aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorVolker Hoffmann <volker@cheleb.net>2016-11-26 22:29:58 +0100
committerVolker Hoffmann <volker@cheleb.net>2016-11-26 22:30:05 +0100
commit76071bb79b6c08229b1fc2925438f70269ad7756 (patch)
tree817f64234ddaddb9d3b1ddc4211ad5efea1d0027
parentb361f1fbd272ce1f2a8d647e0f99cc290828e5bb (diff)
moved script to cut moons to G3 repoHEADmaster
-rw-r--r--IC_Helpers/modic_cutmoons.py42
1 files changed, 0 insertions, 42 deletions
diff --git a/IC_Helpers/modic_cutmoons.py b/IC_Helpers/modic_cutmoons.py
deleted file mode 100644
index 5e5e5b7..0000000
--- a/IC_Helpers/modic_cutmoons.py
+++ /dev/null
@@ -1,42 +0,0 @@
-"""
-Remove moons from stdin (requires default Genga column ordering).
-
-Call Signature:
-$ python modic_cutmoons.py \
- --minpid 5 \
- --ce 95000 < $output_to_change > $changed_output
-"""
-
-import numpy as np
-import argparse
-import sys
-
-# Parse arguments
-parser = argparse.ArgumentParser()
-parser.add_argument('--ce', type=int, required=True, \
- help="Number of close encounters required for removal")
-parser.add_argument('--minpid', type=int, default=1, \
- help="Minimum PID that is cut.")
-args = parser.parse_args()
-
-# Read lines from stdin
-lines_in = sys.stdin.read().rstrip("\n").split("\n")
-
-# Remove particles with too many CEs
-lines_out = []; cce = 0
-for line in lines_in:
- pid = int(line.strip().split()[1])
- cnt = int(line.strip().split()[19])
- if (cnt < args.ce) or (pid < args.minpid):
- lines_out.append(line)
- else:
- sys.stderr.write("Removed Particle %i (%i CEs) \n" % \
- ( int(line.strip().split()[1]), int(line.strip().split()[19]) ) )
- cce += 1
-
-# Removed Stats
-sys.stderr.write("Removed %i Particles\n" % cce)
-
-# Dump to stdout
-for line in lines_out:
- print line