Mercurial > repos > bgruening > rxdock_rbcavity
diff rbdock.py @ 0:e99f861ba68a draft
"planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/chemicaltoolbox/rxdock commit bb19570293b920983b6856b30b42203a09543bc5"
author | bgruening |
---|---|
date | Sat, 04 Apr 2020 15:04:13 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/rbdock.py Sat Apr 04 15:04:13 2020 -0400 @@ -0,0 +1,28 @@ +import subprocess +import argparse + +def main(): + parser = argparse.ArgumentParser(description='Simple wrapper for rbdock') + parser.add_argument('-n', '--num', type=int, help='Number of docking poses to generate') + parser.add_argument('-s', '--seed', type=int, help='Random seed') + args = parser.parse_args() + + cmd = ['rbdock', '-i', 'ligands.sdf', '-r', 'receptor.prm', '-p', 'dock.prm', '-n', str(args.num), '-o', 'rdock_output'] + if args.seed != None: + cmd += ['-s', str(args.seed)] + + ps = subprocess.Popen(cmd, stdout=subprocess.PIPE) + + error_counter = 0 + for stdout_line in iter(ps.stdout.readline, ''): + if 'RBT_DOCKING_ERROR' in str(stdout_line): + error_counter += 1 + if error_counter == 10: + print(ps.stdout) + exit(23) + if ps.poll() != None: + print(ps.stdout) + exit(int(ps.poll())) + +if __name__ == "__main__": + main()