static int multiWigTotalHeight(struct track *tg, enum trackVisibility vis)
/* Return total height of multiWigcontainer. */
{
char *aggregate = cartOrTdbString(cart, tg->tdb, "aggregate", NULL);
boolean overlay = isOverlayTypeAggregate(aggregate);
int totalHeight = 0;
if (overlay)
totalHeight = wigTotalHeight(tg, vis);
struct track *subtrack;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
// Logic is slightly complicated by fact we want to call the totalHeight
// method for each subtrack even if in overlay mode.
int oneHeight = subtrack->totalHeight(subtrack, vis);
if (!overlay)
{
if (totalHeight != 0)
totalHeight += 1;
totalHeight += oneHeight;
}
}
}
tg->height = totalHeight;
return totalHeight;
}
static void multiWigDraw(struct track *tg, int seqStart, int seqEnd,
struct hvGfx *hvg, int xOff, int yOff, int width,
MgFont *font, Color color, enum trackVisibility vis)
/* Draw items in multiWig container. */
{
struct track *subtrack;
char *aggregate = cartOrTdbString(cart, tg->tdb, "aggregate", NULL);
boolean overlay = isOverlayTypeAggregate(aggregate);
int y = yOff;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
int height = subtrack->totalHeight(subtrack, vis);
hvGfxSetClip(hvg, xOff, y, width, height);
if (sameString(WIG_AGGREGATE_TRANSPARENT, aggregate))
hvGfxSetWriteMode(hvg, MG_WRITE_MODE_MULTIPLY);
if (overlay)
subtrack->lineHeight = tg->lineHeight;
subtrack->drawItems(subtrack, seqStart, seqEnd, hvg, xOff, y, width, font, color, vis);
if (!overlay)
y += height + 1;
hvGfxSetWriteMode(hvg, MG_WRITE_MODE_NORMAL);
hvGfxUnclip(hvg);
}
}
char *url = trackUrl(tg->track, chromName);
//fprintf(stderr, "mapBoxHgcOrHgGene(id: %s;) in multiWigDraw\n",
// tg->track);
mapBoxHgcOrHgGene(hvg, seqStart, seqEnd, xOff, y, width, tg->height, tg->track, tg->track, NULL,
url, TRUE, NULL);
}
static int multiWigTotalHeight(struct track *tg, enum trackVisibility vis)
/* Return total height of multiWigcontainer. */
{
struct wigCartOptions *wigCart = tg->wigCartData;
int totalHeight = 0;
if (wigCart->aggregateFunction != wiggleAggregateNone)
totalHeight = wigTotalHeight(tg, vis);
struct track *subtrack;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
// Logic is slightly complicated by fact we want to call the totalHeight
// method for each subtrack even if in overlay mode.
int oneHeight = subtrack->totalHeight(subtrack, vis);
if (wigCart->aggregateFunction == wiggleAggregateNone)
{
if (totalHeight != 0)
totalHeight += 1;
totalHeight += oneHeight;
}
}
}
tg->height = totalHeight;
return totalHeight;
}
static void multiWigLeftLabels(struct track *tg, int seqStart, int seqEnd,
struct hvGfx *hvg, int xOff, int yOff, int width, int height,
boolean withCenterLabels, MgFont *font, Color color,
enum trackVisibility vis)
/* Draw left labels - by deferring to first subtrack. */
{
struct wigCartOptions *wigCart = tg->wigCartData;
if (wigCart->aggregateFunction != wiggleAggregateNone)
{
struct track *firstVisibleSubtrack = NULL;
boolean showNumbers = graphLimitsAllSame(tg->subtracks, &firstVisibleSubtrack);
struct track *subtrack = (showNumbers ? firstVisibleSubtrack : tg->subtracks);
wigLeftAxisLabels(tg, seqStart, seqEnd, hvg, xOff, yOff, width, height, withCenterLabels,
font, color, vis, tg->shortLabel, subtrack->graphUpperLimit,
subtrack->graphLowerLimit, showNumbers);
}
else
{
struct track *subtrack;
int y = yOff;
if (withCenterLabels)
y += tl.fontHeight+1;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
int height = subtrack->totalHeight(subtrack, vis);
if (vis == tvDense)
{
/* Avoid wigLeftAxisLabels here because it will repeatedly add center label
* offsets, and in dense mode we will only draw the one label. */
hvGfxTextRight(hvg, xOff, y, width - 1, height, subtrack->ixColor, font,
subtrack->shortLabel);
}
else
{
// black labels for readability
wigLeftAxisLabels(subtrack, seqStart, seqEnd, hvg, xOff, y, width, height,
withCenterLabels, font, MG_BLACK, vis, subtrack->shortLabel,
subtrack->graphUpperLimit, subtrack->graphLowerLimit, TRUE);
}
y += height+1;
}
}
}
}
static boolean graphLimitsAllSame(struct track *trackList, struct track **retFirstTrack)
/* Return TRUE if graphUpperLimit and graphLowerLimit same for all tracks. */
{
struct track *firstTrack = NULL;
struct track *track;
for (track = trackList; track != NULL; track = track->next)
{
if (isSubtrackVisible(track))
{
if (firstTrack == NULL)
*retFirstTrack = firstTrack = track;
else if (track->graphUpperLimit != firstTrack->graphUpperLimit
|| track->graphLowerLimit != firstTrack->graphLowerLimit)
return FALSE;
}
}
return firstTrack != NULL;
}
static void multiWigLeftLabels(struct track *tg, int seqStart, int seqEnd,
struct hvGfx *hvg, int xOff, int yOff, int width, int height,
boolean withCenterLabels, MgFont *font, Color color,
enum trackVisibility vis)
/* Draw left labels - by deferring to first subtrack. */
{
char *aggregate = cartOrTdbString(cart, tg->tdb, "aggregate", NULL);
boolean overlay = isOverlayTypeAggregate(aggregate);
if (overlay)
{
struct track *firstVisibleSubtrack = NULL;
boolean showNumbers = graphLimitsAllSame(tg->subtracks, &firstVisibleSubtrack);
struct track *subtrack = (showNumbers ? firstVisibleSubtrack : tg->subtracks);
wigLeftAxisLabels(tg, seqStart, seqEnd, hvg, xOff, yOff, width, height, withCenterLabels,
font, color, vis, tg->shortLabel, subtrack->graphUpperLimit,
subtrack->graphLowerLimit, showNumbers);
}
else
{
struct track *subtrack;
int y = yOff;
if (withCenterLabels)
y += tl.fontHeight+1;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
int height = subtrack->totalHeight(subtrack, vis);
wigLeftAxisLabels(subtrack, seqStart, seqEnd, hvg, xOff, y, width, height, withCenterLabels,
font, subtrack->ixColor, vis, subtrack->shortLabel, subtrack->graphUpperLimit,
subtrack->graphLowerLimit, TRUE);
y += height+1;
}
}
}
}
static void multiWigDraw(struct track *tg, int seqStart, int seqEnd,
struct hvGfx *hvg, int xOff, int yOff, int width,
MgFont *font, Color color, enum trackVisibility vis)
/* Draw items in multiWig container. */
{
struct track *subtrack;
boolean errMsgShown = FALSE;
int y = yOff;
boolean errMsgFound = FALSE;
int numTracks = 0;
// determine if any subtracks had errors and count them up
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack) )
{
numTracks++;
if (subtrack->networkErrMsg)
errMsgFound = TRUE;
}
}
if (errMsgFound)
{
Color yellow = hvGfxFindRgb(hvg, &undefinedYellowColor);
hvGfxBox(hvg, xOff, yOff, width, tg->height, yellow);
}
struct wigCartOptions *wigCart = tg->wigCartData;
struct wigGraphOutput *wgo = setUpWgo(xOff, yOff, width, tg->height, numTracks, wigCart, hvg);
// we want to the order to be the same in all the modes
if (wigCart->aggregateFunction == wiggleAggregateStacked)
slReverse(&tg->subtracks);
/* Cope with autoScale and stacked bars - we do it here rather than in the child tracks, so that
* all children can be on same scale. */
double minVal, maxVal;
/* Force load of all predraw arrays so can do calcs. Build up list, and then
* figure out max/min. No worries about multiple loading, the loaders protect
* themselves. */
struct slRef *refList = NULL;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
struct preDrawContainer *pre = subtrack->loadPreDraw(subtrack, seqStart, seqEnd, width);
if (pre != NULL) // pre maybe null if the load fails
{
preDrawWindowFunction(pre->preDraw, pre->preDrawSize, wigCart->windowingFunction,
wigCart->transformFunc, wigCart->doNegative);
preDrawSmoothing(pre->preDraw, pre->preDrawSize, wigCart->smoothingWindow);
pre->smoothingDone = TRUE;
refAdd(&refList, pre);
}
}
}
slReverse(&refList);
minMaxVals(refList, &minVal, &maxVal, wigCart->windowingFunction, wigCart->alwaysZero, wgo->yOffsets);
slFreeList(&refList);
if (!wigCart->autoScale)
{
minVal = wigCart->minY;
maxVal = wigCart->maxY;
}
/* Loop through again setting up the wigCarts of the children to have minY/maxY for
* our limits and autoScale off. */
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (!isSubtrackVisible(subtrack))
continue;
struct wigCartOptions *wigCart = subtrack->wigCartData;
wigCart->minY = minVal;
wigCart->maxY = maxVal;
wigCart->autoScale = wiggleScaleManual;
struct preDrawContainer *pre = subtrack->preDrawContainer;
if (pre != NULL) // pre maybe null if the load fails
{
pre->graphUpperLimit = maxVal;
pre->graphLowerLimit = minVal;
}
}
int numTrack = 0;
for (subtrack = tg->subtracks; subtrack != NULL; subtrack = subtrack->next)
{
if (isSubtrackVisible(subtrack))
{
if (!subtrack->networkErrMsg || !errMsgShown)
{
if (subtrack->networkErrMsg)
errMsgShown = TRUE;
wgo->numTrack = numTrack++;
subtrack->wigGraphOutput = wgo;
int height = subtrack->totalHeight(subtrack, vis);
hvGfxSetClip(hvg, xOff, y, width, height);
if (wigCart->aggregateFunction != wiggleAggregateNone)
subtrack->lineHeight = tg->lineHeight;
subtrack->drawItems(subtrack, seqStart, seqEnd, hvg, xOff, y, width, font, color, vis);
if (wigCart->aggregateFunction == wiggleAggregateNone)
{
y += height + 1;
wgo->yOff = y;
}
hvGfxUnclip(hvg);
}
}
}
if (wigCart->aggregateFunction == wiggleAggregateTransparent)
{
floatPicIntoHvg(wgo->image, xOff, yOff, hvg);
floatPicFree((struct floatPic **)&wgo->image);
}
char *url = trackUrl(tg->track, chromName);
mapBoxHgcOrHgGene(hvg, seqStart, seqEnd, xOff, yOff, width, tg->height, tg->track, tg->track, NULL,
url, TRUE, NULL);
}
