annotate galaxy-tools/biobank/utils/delete_flowcell_results.py @ 12:46f08bb8dd68 draft default tip

Uploaded
author ric
date Wed, 28 Sep 2016 04:59:02 -0400
parents 43be74e62bfe
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
43be74e62bfe Uploaded
ric
parents:
diff changeset
1 import argparse, sys
43be74e62bfe Uploaded
ric
parents:
diff changeset
2
43be74e62bfe Uploaded
ric
parents:
diff changeset
3 from bl.vl.kb import KnowledgeBase as KB
43be74e62bfe Uploaded
ric
parents:
diff changeset
4 from bl.vl.utils import get_logger, LOG_LEVELS
43be74e62bfe Uploaded
ric
parents:
diff changeset
5 import bl.vl.utils.ome_utils as vlu
43be74e62bfe Uploaded
ric
parents:
diff changeset
6
43be74e62bfe Uploaded
ric
parents:
diff changeset
7
43be74e62bfe Uploaded
ric
parents:
diff changeset
8 def make_parser():
43be74e62bfe Uploaded
ric
parents:
diff changeset
9 parser = argparse.ArgumentParser(description='remove datasamples connected to a specific sample of a flowcell')
43be74e62bfe Uploaded
ric
parents:
diff changeset
10 parser.add_argument('--logfile', type=str, help='log file (default=stderr)')
43be74e62bfe Uploaded
ric
parents:
diff changeset
11 parser.add_argument('--loglevel', type=str, choices=LOG_LEVELS,
43be74e62bfe Uploaded
ric
parents:
diff changeset
12 help='logging level (default=INFO)', default='INFO')
43be74e62bfe Uploaded
ric
parents:
diff changeset
13 parser.add_argument('-H', '--host', type=str, help='omero server hostname')
43be74e62bfe Uploaded
ric
parents:
diff changeset
14 parser.add_argument('-U', '--user', type=str, help='omero user')
43be74e62bfe Uploaded
ric
parents:
diff changeset
15 parser.add_argument('-P', '--passwd', type=str, help='omero password')
43be74e62bfe Uploaded
ric
parents:
diff changeset
16 parser.add_argument('--flowcell-label', type=str, required=True,
43be74e62bfe Uploaded
ric
parents:
diff changeset
17 help='flowcell label')
43be74e62bfe Uploaded
ric
parents:
diff changeset
18 parser.add_argument('--sample-label', type=str, required=True,
43be74e62bfe Uploaded
ric
parents:
diff changeset
19 help='sample label')
43be74e62bfe Uploaded
ric
parents:
diff changeset
20 parser.add_argument('--dry-run', action='store_true',
43be74e62bfe Uploaded
ric
parents:
diff changeset
21 help='run a simulation, don\'t delete anything')
43be74e62bfe Uploaded
ric
parents:
diff changeset
22 return parser
43be74e62bfe Uploaded
ric
parents:
diff changeset
23
43be74e62bfe Uploaded
ric
parents:
diff changeset
24
43be74e62bfe Uploaded
ric
parents:
diff changeset
25 def get_flowcell_samples_map(flowcell, kb, logger):
43be74e62bfe Uploaded
ric
parents:
diff changeset
26 fc = kb.get_by_label(kb.FlowCell, flowcell)
43be74e62bfe Uploaded
ric
parents:
diff changeset
27 if not fc:
43be74e62bfe Uploaded
ric
parents:
diff changeset
28 logger.info('No flowcell with label %s', flowcell)
43be74e62bfe Uploaded
ric
parents:
diff changeset
29 sys.exit(0)
43be74e62bfe Uploaded
ric
parents:
diff changeset
30 logger.info('Loading data for flowcell %s', flowcell)
43be74e62bfe Uploaded
ric
parents:
diff changeset
31 dsamples = kb.dt.get_connected(fc, aklass = kb.SeqDataSample)
43be74e62bfe Uploaded
ric
parents:
diff changeset
32 dsamples_map = {}
43be74e62bfe Uploaded
ric
parents:
diff changeset
33 for ds in dsamples:
43be74e62bfe Uploaded
ric
parents:
diff changeset
34 if ds.sample:
43be74e62bfe Uploaded
ric
parents:
diff changeset
35 dsamples_map.setdefault(ds.sample.label, []).append(ds)
43be74e62bfe Uploaded
ric
parents:
diff changeset
36 else:
43be74e62bfe Uploaded
ric
parents:
diff changeset
37 dsamples_map.setdefault('NO_SAMPLE', []).append(ds)
43be74e62bfe Uploaded
ric
parents:
diff changeset
38 return dsamples_map
43be74e62bfe Uploaded
ric
parents:
diff changeset
39
43be74e62bfe Uploaded
ric
parents:
diff changeset
40
43be74e62bfe Uploaded
ric
parents:
diff changeset
41 def print_report(dsamples_map, sample_label, kb, logger):
43be74e62bfe Uploaded
ric
parents:
diff changeset
42 dsamples = dsamples_map.get(sample_label)
43be74e62bfe Uploaded
ric
parents:
diff changeset
43 if not dsamples:
43be74e62bfe Uploaded
ric
parents:
diff changeset
44 logger.info('No sample with label %s is related to the flowcell', sample_label)
43be74e62bfe Uploaded
ric
parents:
diff changeset
45 sys.exit(0)
43be74e62bfe Uploaded
ric
parents:
diff changeset
46 for ds in dsamples:
43be74e62bfe Uploaded
ric
parents:
diff changeset
47 dobjs = kb.get_data_objects(ds)
43be74e62bfe Uploaded
ric
parents:
diff changeset
48 logger.info('## data sample: %s', ds.label)
43be74e62bfe Uploaded
ric
parents:
diff changeset
49 for dob in dobjs:
43be74e62bfe Uploaded
ric
parents:
diff changeset
50 logger.info('### data object: %s --- mimetype: %s', dob.path, dob.mimetype)
43be74e62bfe Uploaded
ric
parents:
diff changeset
51
43be74e62bfe Uploaded
ric
parents:
diff changeset
52
43be74e62bfe Uploaded
ric
parents:
diff changeset
53 def delete(dsamples_map, sample_label, kb, logger):
43be74e62bfe Uploaded
ric
parents:
diff changeset
54 for ds in dsamples_map[sample_label]:
43be74e62bfe Uploaded
ric
parents:
diff changeset
55 # this is a hack specific for the automator workflow
43be74e62bfe Uploaded
ric
parents:
diff changeset
56 if not ds.label.startswith('stage1'):
43be74e62bfe Uploaded
ric
parents:
diff changeset
57 logger.info('Deleting data for %s', ds.label)
43be74e62bfe Uploaded
ric
parents:
diff changeset
58 dobjs = kb.get_data_objects(ds)
43be74e62bfe Uploaded
ric
parents:
diff changeset
59 for d in dobjs:
43be74e62bfe Uploaded
ric
parents:
diff changeset
60 kb.delete(d)
43be74e62bfe Uploaded
ric
parents:
diff changeset
61 a = ds.action
43be74e62bfe Uploaded
ric
parents:
diff changeset
62 kb.delete(ds)
43be74e62bfe Uploaded
ric
parents:
diff changeset
63 try:
43be74e62bfe Uploaded
ric
parents:
diff changeset
64 kb.delete(a)
43be74e62bfe Uploaded
ric
parents:
diff changeset
65 except:
43be74e62bfe Uploaded
ric
parents:
diff changeset
66 pass
43be74e62bfe Uploaded
ric
parents:
diff changeset
67
43be74e62bfe Uploaded
ric
parents:
diff changeset
68
43be74e62bfe Uploaded
ric
parents:
diff changeset
69 def main(argv):
43be74e62bfe Uploaded
ric
parents:
diff changeset
70 parser = make_parser()
43be74e62bfe Uploaded
ric
parents:
diff changeset
71 args = parser.parse_args(argv)
43be74e62bfe Uploaded
ric
parents:
diff changeset
72
43be74e62bfe Uploaded
ric
parents:
diff changeset
73 logger = get_logger('delete_flowcell_results', level=args.loglevel,
43be74e62bfe Uploaded
ric
parents:
diff changeset
74 filename=args.logfile)
43be74e62bfe Uploaded
ric
parents:
diff changeset
75
43be74e62bfe Uploaded
ric
parents:
diff changeset
76 try:
43be74e62bfe Uploaded
ric
parents:
diff changeset
77 host = args.host or vlu.ome_host()
43be74e62bfe Uploaded
ric
parents:
diff changeset
78 user = args.user or vlu.ome_user()
43be74e62bfe Uploaded
ric
parents:
diff changeset
79 passwd = args.passwd or vlu.ome_passwd()
43be74e62bfe Uploaded
ric
parents:
diff changeset
80 except ValueError, ve:
43be74e62bfe Uploaded
ric
parents:
diff changeset
81 logger.critical(ve)
43be74e62bfe Uploaded
ric
parents:
diff changeset
82 sys.exit(ve)
43be74e62bfe Uploaded
ric
parents:
diff changeset
83
43be74e62bfe Uploaded
ric
parents:
diff changeset
84 kb = KB(driver='omero')(host, user, passwd)
43be74e62bfe Uploaded
ric
parents:
diff changeset
85
43be74e62bfe Uploaded
ric
parents:
diff changeset
86 dsamples_map = get_flowcell_samples_map(args.flowcell_label, kb, logger)
43be74e62bfe Uploaded
ric
parents:
diff changeset
87 print_report(dsamples_map, args.sample_label, kb, logger)
43be74e62bfe Uploaded
ric
parents:
diff changeset
88 if not args.dry_run:
43be74e62bfe Uploaded
ric
parents:
diff changeset
89 delete(dsamples_map, args.sample_label, kb. logger)
43be74e62bfe Uploaded
ric
parents:
diff changeset
90 pass
43be74e62bfe Uploaded
ric
parents:
diff changeset
91 else:
43be74e62bfe Uploaded
ric
parents:
diff changeset
92 logger.debug('SIMULATION, exit now')
43be74e62bfe Uploaded
ric
parents:
diff changeset
93 logger.info('Job completed')
43be74e62bfe Uploaded
ric
parents:
diff changeset
94
43be74e62bfe Uploaded
ric
parents:
diff changeset
95
43be74e62bfe Uploaded
ric
parents:
diff changeset
96 if __name__ == '__main__':
43be74e62bfe Uploaded
ric
parents:
diff changeset
97 main(sys.argv[1:])