#!/usr/bin/perl

@datasets = ("hsTTPKO-m", "hsL1KO-m", "hsL2KO-m", "hsTripleKO-m", "hsTripleKO-f", "hsTripleKOadd-m");
@conditions = ("WT", "KO");
@reps = (1,2,3,4);

%ds2color = (
"hsTTPKO-m" => "204,229,255|170,212,255|136,195,255|102,178,255", #blue
"hsL1KO-m" => "255,229,204|255,212,170|255,195,136|255,178,102", # orange
"hsL2KO-m" => "240,230,140|232,209,104|225,187,68|218,165,32",  # yellow/brown
"hsTripleKO-m" => "229,204,255|212,170,255|195,136,255|178,102,255", # purple
"hsTripleKO-f" => "255,204,255|255,170,255|255,136,255|255,102,255", # pink
"hsTripleKOadd-m" => "255,204,204|255,170,170|255,136,136|255,102,102"  # red
);
# "204,255,255|153,255,204|102,255,178|0,204,102",  #green


open(OUT, ">project1.trackDb.txt");

$supertrack_name = "st_NormDepth1";
print OUT "track $supertrack_name\n";
print OUT "superTrack on show\n";
print OUT "shortLabel\tRNAseq Proj1 Normalized Depth\n";
print OUT "longLabel\tRNAseq Proj1 Normalized Depth\n\n";
$pri = 1.0;
foreach $ds (@datasets) { foreach $cnd (@conditions) { $root = "$ds.$cnd";
  $pri+=0.1; $prio = sprintf("%.1f", $pri);
  $mwtrack = "st_NormDepth_$root"."_multiwig";
  @colors = split/\|/, $ds2color{$ds};
  print OUT "\ttrack $mwtrack\n";
  print OUT "\tparent $supertrack_name\n";
  print OUT "\tcontainer multiWig\n";
  print OUT "\tshortLabel\t$root depth\n";
  print OUT "\tlongLabel\t$root depth\n";
  print OUT "\ttype bigWig\n";
#  print OUT "\tviewLimits 0:100\n";
  print OUT "\twindowingFunction mean\n";
  print OUT "\tgridDefault on\n";
  print OUT "\tautoScale on\n";
  print OUT "\talwaysZero on\n";
  print OUT "\tconfigurable on\n";
  print OUT "\tvisibility full\n";
  print OUT "\taggregate transparentOverlay\n";
  print OUT "\tshowSubtrackColorOnUi on\n";
  print OUT "\tmaxHeightPixels 150:40:11\n";
  print OUT "\tpriority $prio\n\n";
  foreach $r (@reps) { $id = "$ds.$cnd.$r";
    $rr	= $r-1;	$rgb = $colors[$rr]; 
    print OUT "\t\ttrack st_NormDepth_$id\n";
    print OUT "\t\tbigDataUrl http://snpinfo.niehs.nih.gov/ucscview/blackshear/RNAseqHub_May2016/mm10/DATA/$id.norm.bigWig\n";
    print OUT "\t\tshortLabel\t$id depth\n";
    print OUT "\t\tlongLabel\tRNAseq depth for $id\n";
    print OUT "\t\tparent $mwtrack\n";
    print OUT "\t\ttype bigWig\n";
#    print OUT "\t\tgraphTypeDefault points\n";
    print OUT "\t\tcolor $rgb\n\n";
  }
  print OUT "\n";
} }
print OUT "\n\n\n";


close(OUT);