From be132b9736d7952028bc3c8ff57630586fe5847b Mon Sep 17 00:00:00 2001 From: Daniel Asarnow Date: Sat, 30 Sep 2017 15:36:31 -0700 Subject: [PATCH] Let relion compute apix if none provided --- reconstruct.py | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/reconstruct.py b/reconstruct.py index 6f583a9..9296a13 100755 --- a/reconstruct.py +++ b/reconstruct.py @@ -65,9 +65,14 @@ def do_job(star): def do_reconstruct(star, mrc, apix, sym="C1", ctf=True, relion_path="relion_reconstruct"): - com = relion_path + \ - " --angpix %f --sym %s --ctf %s --i %s --o %s" % \ - (apix, sym, str(ctf).lower(), star, mrc) + if apix is not None: + com = relion_path + \ + " --angpix %f --sym %s --ctf %s --i %s --o %s" % \ + (apix, sym, str(ctf).lower(), star, mrc) + else: + com = relion_path + \ + " --sym %s --ctf %s --i %s --o %s" % \ + (sym, str(ctf).lower(), star, mrc) #os.system(com) try: output = subprocess.check_output(shlex.split(com), stderr=subprocess.STDOUT)