Update the django app to report "raw reads" for a sequence lane.
[htsworkflow.git] / htsworkflow / frontend / samples / views.py
index d6647d1dba38311733d9f3d07530df9f77442a6b..86c15117e8c2440d3ff7d088053a054f2047c8ff 100644 (file)
@@ -9,7 +9,8 @@ from htsworkflow.frontend import settings
 from htsworkflow.util import makebed
 from htsworkflow.util import opener
 
-from django.http import HttpResponse
+from django.http import HttpResponse, HttpResponseRedirect
+from django.shortcuts import render_to_response
 from django.template import RequestContext
 from django.template.loader import get_template
 
@@ -62,70 +63,41 @@ def library_to_flowcells(request, lib_id):
     """
     Display information about all the flowcells a library has been run on.
     """
-    t = get_template("samples/library_detail.html")
     
     try:
       lib = Library.objects.get(library_id=lib_id)
     except:
       return HttpResponse("Library %s does not exist" % (lib_id))
-    
-    output = []
-    
-    output.append('<b>Library ID:</b> %s' % (lib.library_id))
-    output.append('<b>Name:</b> %s' % (lib.library_name))
-    output.append('<b>Species:</b> %s' % (lib.library_species.scientific_name))
-    output.append('')
-    
-    output.append('<b>FLOWCELL - LANE:</b>')
-    
-    output.extend([ '%s - Lane 1 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 1)) for fc in lib.lane_1_library.all() ])
-    output.extend([ '%s - Lane 2 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 2)) for fc in lib.lane_2_library.all() ])
-    output.extend([ '%s - Lane 3 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 3)) for fc in lib.lane_3_library.all() ])
-    output.extend([ '%s - Lane 4 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 4)) for fc in lib.lane_4_library.all() ])
-    output.extend([ '%s - Lane 5 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 5)) for fc in lib.lane_5_library.all() ])
-    output.extend([ '%s - Lane 6 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 6)) for fc in lib.lane_6_library.all() ])
-    output.extend([ '%s - Lane 7 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 7)) for fc in lib.lane_7_library.all() ])
-    output.extend([ '%s - Lane 8 %s' % (fc.flowcell_id, _files(fc.flowcell_id, 8)) for fc in lib.lane_8_library.all() ])
-    
-    record_count = lib.lane_1_library.count() + \
-                    lib.lane_2_library.count() + \
-                    lib.lane_3_library.count() + \
-                    lib.lane_4_library.count() + \
-                    lib.lane_5_library.count() + \
-                    lib.lane_6_library.count() + \
-                    lib.lane_7_library.count() + \
-                    lib.lane_8_library.count()
-    
+   
     flowcell_list = []
-    flowcell_list.extend([ (fc.flowcell_id, 1) for fc in lib.lane_1_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 2) for fc in lib.lane_2_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 3) for fc in lib.lane_3_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 4) for fc in lib.lane_4_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 5) for fc in lib.lane_5_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 6) for fc in lib.lane_6_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 7) for fc in lib.lane_7_library.all() ])
-    flowcell_list.extend([ (fc.flowcell_id, 8) for fc in lib.lane_8_library.all() ])
+    interesting_flowcells = {} # aka flowcells we're looking at
+    for lane in LANE_LIST:
+        lane_library = getattr(lib, 'lane_%d_library' % (lane,))
+        for fc in lane_library.all():
+            flowcell_id, id = parse_flowcell_id(fc.flowcell_id)
+            if flowcell_id not in interesting_flowcells:
+                interesting_flowcells[flowcell_id] = get_flowcell_result_dict(flowcell_id)
+            flowcell_list.append((fc.flowcell_id, lane))
+
     flowcell_list.sort()
     
     lane_summary_list = []
+    eland_results = []
     for fc, lane in flowcell_list:
         lane_summary, err_list = _summary_stats(fc, lane)
         
+        eland_results.extend(_make_eland_results(fc, lane, interesting_flowcells))
         lane_summary_list.extend(lane_summary)
-    
-        for err in err_list:    
-            output.append(err)
-   
-    output.append('<br />')
-    output.append(t.render(RequestContext(request, {'lane_summary_list': lane_summary_list})))
-    output.append('<br />')
-    
-    if record_count == 0:
-        output.append("None Found")
-    
-    return HttpResponse('<br />\n'.join(output))
 
 
+    return render_to_response(
+        'samples/library_detail.html',
+        {'lib': lib,
+         'eland_results': eland_results,
+         'lane_summary_list': lane_summary_list,
+        },
+        context_instance = RequestContext(request))
+
 def summaryhtm_fc_cnm(request, fc_id, cnm):
     """
     returns a Summary.htm file if it exists.
@@ -170,9 +142,12 @@ def result_fc_cnm_eland_lane(request, fc_id, cnm, lane):
     
     filepath = erd[lane]
     
-    f = opener.autoopen(filepath, 'r')
+    #f = opener.autoopen(filepath, 'r')
+    # return HttpResponse(f, mimetype="application/x-elandresult")
+
+    f = open(filepath, 'r')
+    return HttpResponse(f, mimetype='application/x-bzip2')
     
-    return HttpResponse(f, mimetype="application/x-elandresult")
 
 
 def bedfile_fc_cnm_eland_lane_ucsc(request, fc_id, cnm, lane):
@@ -247,7 +222,10 @@ def _summary_stats(flowcell_id, lane_id):
                 eland_summary.flowcell_id = flowcell_id
                 eland_summary.clusters = gerald_summary[end][lane_id].cluster
                 eland_summary.cycle_width = cycle_width
-                eland_summary.summarized_reads = runfolder.summarize_mapped_reads(eland_summary.genome_map, eland_summary.mapped_reads)
+               if hasattr(eland_summary, 'genome_map'):
+                    eland_summary.summarized_reads = runfolder.summarize_mapped_reads( 
+                                                       eland_summary.genome_map, 
+                                                       eland_summary.mapped_reads)
 
                 # grab some more information out of the flowcell db
                 flowcell = FlowCell.objects.get(flowcell_id=fc_id)
@@ -335,7 +313,45 @@ def get_eland_result_type(pathname):
         return 'result'
     else:
         return 'unknown'
-    
+
+def _make_eland_results(flowcell_id, lane, interesting_flowcells):
+
+    cur_fc = interesting_flowcells.get(flowcell_id, None)
+    if cur_fc is None:
+      return []
+
+    results = []
+    for cycle in cur_fc.keys():
+        result_path = cur_fc[cycle]['eland_results'].get(lane, None)
+        result_link = make_result_link(flowcell_id, cycle, lane, result_path)
+        results.append({'flowcell_id': flowcell_id,
+                        'cycle': cycle, 
+                        'lane': lane, 
+                        'summary_url': make_summary_url(flowcell_id, cycle),
+                        'result_url': result_link[0],
+                        'result_label': result_link[1],
+                        'bed_url': result_link[2],
+        })
+    return results
+
+def make_summary_url(flowcell_id, cycle_name):
+    url = '/results/%s/%s/summary/' % (flowcell_id, cycle_name)
+    return url
+
+def make_result_link(flowcell_id, cycle_name, lane, eland_result_path):
+    if eland_result_path is None:
+        return ("", "", "")
+
+    result_type = get_eland_result_type(eland_result_path)
+    result_url = '/results/%s/%s/eland_result/%s' % (flowcell_id, cycle_name, lane)
+    result_label = 'eland %s' % (result_type,)
+    bed_url = None
+    if result_type == 'result':
+       bed_url_pattern = '/results/%s/%s/bedfile/%s'
+       bed_url = bed_url_pattern % (flowcell_id, cycle_name, lane)
+    
+    return (result_url, result_label, bed_url)
+
 def _files(flowcell_id, lane):
     """
     Sets up available files for download
@@ -370,4 +386,8 @@ def _files(flowcell_id, lane):
         return ''
     
     return '(' + '|'.join(output) + ')'
-            
+
+def library_id_to_admin_url(request, lib_id):
+    lib = Library.objects.get(library_id=lib_id)
+    return HttpResponseRedirect('/admin/samples/library/%s' % (lib.id,))
+