src/hg/lib/hgFind.c 1.224

1.224 2010/02/10 22:54:25 angie
Added hook for specifying a search type: CGI var (not cart) singleSearch=knownCanonical (and maybe other types someday). This will support stable links to hgTracks views of gene symbols short-circuited to knownCanonical locations.
Index: src/hg/lib/hgFind.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/lib/hgFind.c,v
retrieving revision 1.223
retrieving revision 1.224
diff -b -B -U 4 -r1.223 -r1.224
--- src/hg/lib/hgFind.c	9 Feb 2010 23:27:21 -0000	1.223
+++ src/hg/lib/hgFind.c	10 Feb 2010 22:54:25 -0000	1.224
@@ -2962,8 +2962,59 @@
     hgPosFreeList(&(hgp->tableList->posList->next));
     }
 }
 
+static boolean searchKnownCanonical(char *db, char *term, struct hgPositions *hgp)
+/* Look for term in kgXref.geneSymbol, and if found, put knownCanonical coords and 
+ * knownGene.name in hgp. */
+{
+boolean foundIt = FALSE;
+struct sqlConnection *conn = hAllocConn(db);
+if (sqlTableExists(conn, "knownGene") && sqlTableExists(conn, "knownCanonical") &&
+    sqlTableExists(conn, "kgXref"))
+    {
+    char query[512];
+    safef(query, sizeof(query), "select knownCanonical.chrom,chromStart,chromEnd,name "
+	  "from knownCanonical,knownGene,kgXref where kgXref.geneSymbol = '%s' "
+	  "and kgXref.kgId = knownGene.name and knownGene.name = knownCanonical.transcript;",
+	  term);
+    struct sqlResult *sr = sqlGetResult(conn, query);
+    char **row;
+    if ((row = sqlNextRow(sr)) != NULL)
+	{
+	singlePos(hgp, "UCSC Genes", term, "knownGene", row[3], row[3],
+		  cloneString(row[0]), atoi(row[1]), atoi(row[2]));
+	foundIt = TRUE;
+	}
+    sqlFreeResult(&sr);
+    }
+hFreeConn(&conn);
+return foundIt;
+}
+
+static boolean singleSearch(char *db, char *term, struct cart *cart, struct hgPositions *hgp)
+/* If a search type is specified in the CGI line (not cart), perform that search. 
+ * If the search is successful, fill in hgp as a single-pos result and return TRUE. */
+{
+char *search = cgiOptionalString("singleSearch");
+if (search == NULL)
+    return FALSE;
+
+cartRemove(cart, "singleSearch");
+boolean foundIt = FALSE;
+if (sameString(search, "knownCanonical"))
+    foundIt = searchKnownCanonical(db, term, hgp);
+else
+    warn("Unrecognized singleSearch=%s in URL", search);
+if (foundIt)
+    {
+    fixSinglePos(hgp);
+    if (cart != NULL)
+        cartSetString(cart, "hgFind.matches", hgp->tableList->posList->browserName);
+    }
+return foundIt;
+}
+
 struct hgPositions *hgPositionsFind(char *db, char *term, char *extraCgi,
 	char *hgAppNameIn, struct cart *cart, boolean multiTerm)
 /* Return table of positions that match term or NULL if none such. */
 {
@@ -2990,8 +3041,11 @@
 if (extraCgi == NULL)
     extraCgi = "";
 hgp->extraCgi = cloneString(extraCgi);
 
+if (singleSearch(db, term, cart, hgp))
+    return hgp;
+
 chimpSpecialChrom(hgp, &term, hgAppNameIn);
 
 /* Allow any search term to end with a :Start-End range -- also support stuff 
  * pasted in from BED (chrom start end) or SQL query (chrom | start | end).