#!/usr/bin/perl

@datasets = ("TTPKO-m", "DoubleKO-m", "TripleKO-m");
@conditions = ("WT", "KO");
@reps = (1,2,3,4);

%ds2color = (
"TTPKO-m" => "204,229,255|170,212,255|136,195,255|102,178,255", #blue
"DoubleKO-m" => "204,255,255|153,255,204|102,255,178|0,204,102",  #green
"TripleKO-m" => "229,204,255|212,170,255|195,136,255|178,102,255" # purple
);


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

$supertrack_name = "st_NormDepth2";
print OUT "track $supertrack_name\n";
print OUT "superTrack on show\n";
print OUT "shortLabel\tRNAseq Proj2 Normalized Depth\n";
print OUT "longLabel\tRNAseq Proj2 Normalized Depth\n\n";
$pri = 3.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);