src/hg/hgc/bigBedClick.c 1.2

1.2 2009/02/10 22:17:39 kent
Creating bigBed type custom tracks.
Index: src/hg/hgc/bigBedClick.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/hgc/bigBedClick.c,v
retrieving revision 1.1
retrieving revision 1.2
diff -b -B -U 1000000 -r1.1 -r1.2
--- src/hg/hgc/bigBedClick.c	3 Feb 2009 21:37:40 -0000	1.1
+++ src/hg/hgc/bigBedClick.c	10 Feb 2009 22:17:39 -0000	1.2
@@ -1,81 +1,103 @@
 /* Handle details pages for wiggle tracks. */
 
 #include "common.h"
 #include "wiggle.h"
 #include "cart.h"
 #include "hgc.h"
 #include "hCommon.h"
 #include "hgColors.h"
 #include "bigBed.h"
 
 static char const rcsid[] = "$Id$";
 
 
-void genericBigBedClick(struct sqlConnection *conn, struct trackDb *tdb, 
+static void bigBedClick(char *fileName, struct trackDb *tdb, 
 		     char *item, int start, int bedSize)
 /* Handle click in generic bigBed track. */
 {
 char *chrom = cartString(cart, "c");
 
-char query[256];
-safef(query, sizeof(query), "select fileName from %s", tdb->tableName);
-char *fileName = sqlQuickString(conn, query);
-if (fileName == NULL)
-    errAbort("Missing fileName in %s table", tdb->tableName);
-
 /* Open BigWig file and get interval list. */
 struct bbiFile *bbi = bigBedFileOpen(fileName);
 struct lm *lm = lmInit(0);
 struct bigBedInterval *bbList = bigBedIntervalQuery(bbi, chrom, winStart, winEnd, lm);
-bbiFileClose(&bbi);
+
+/* Get bedSize if it's not already defined. */
+if (bedSize == 0)
+    bedSize = bbi->definedFieldCount;
 
 /* Find particular item in list - matching start, and item if possible. */
 struct bigBedInterval *bbMatch = NULL, *bb;
 for (bb = bbList; bb != NULL; bb = bb->next)
     {
     if (bb->start == start)
         {
 	if (bedSize > 3)
 	    {
 	    char *name = cloneFirstWordInLine(bb->rest);
 	    boolean match = sameString(name, item);
 	    freez(&name);
 	    if (match)
 	        {
 		bbMatch = bb;
 		break;
 		}
 	    }
 	else
 	    {
 	    bbMatch = bb;
 	    break;
 	    }
 	}
     }
 
 if (bbMatch != NULL)
     {
     char *fields[bedSize];
     char startBuf[16], endBuf[16];
     char *rest = cloneString(bbMatch->rest);
     int bbFieldCount = bigBedIntervalToRow(bbMatch, chrom, startBuf, endBuf, fields, bedSize); 
     if (bbFieldCount != bedSize)
         {
 	errAbort("Disagreement between trackDb field count (%d) and %s fieldCount (%d)", 
 		bedSize, fileName, bbFieldCount);
 	}
     struct bed *bed = bedLoadN(fields, bedSize);
     bedPrintPos(bed, bedSize, tdb);
     printf("Full bed record:<BR><PRE><TT>%s\t%u\t%u\t%s\n</TT></PRE>\n",
 	    chrom, bb->start, bb->end, rest);
     }
 else
     {
     printf("No item %s starting at %d\n", emptyForNull(item), start);
     }
 
 lmCleanup(&lm);
+bbiFileClose(&bbi);
 }
 
+void genericBigBedClick(struct sqlConnection *conn, struct trackDb *tdb, 
+		     char *item, int start, int bedSize)
+/* Handle click in generic bigBed track. */
+{
+char query[256];
+safef(query, sizeof(query), "select fileName from %s", tdb->tableName);
+char *fileName = sqlQuickString(conn, query);
+if (fileName == NULL)
+    errAbort("Missing fileName in %s table", tdb->tableName);
+bigBedClick(fileName, tdb, item, start, bedSize);
+}
 
+void bigBedCustomClick(struct trackDb *tdb)
+/* Display details for BigWig custom tracks. */
+{
+char *fileName = trackDbSetting(tdb, "dataUrl");
+char *item = cartOptionalString(cart, "i");
+int start = cartInt(cart, "o");
+uglyf("fileName=%s, item=%s, start=%d<BR>\n", fileName, item, start);
+struct bbiFile *bbi = bigBedFileOpen(fileName);
+uglyf("bbi->fieldCount = %d, ->definedFieldCount=%d<BR>\n", bbi->fieldCount, bbi->definedFieldCount);
+bigBedClick(fileName, tdb, item, start, 0);
+#ifdef SOON
+#endif
+}