#!/usr/bin/perl

@datasets = ("M_fed_9-11am_8-10wk", "M_fasted_9-11am_10-12wk", "M_fed_8am_8-10wk", "M_fed_8pm_8-10wk", "F_fed_8am_12-16wk");
@conditions = ("WT", "KO");
@reps = (1,2,3,4);

%ds2color = (
"M_fed_9-11am_8-10wk" => "204,229,255|170,212,255|136,195,255|102,178,255", #blue
"M_fasted_9-11am_10-12wk" => "255,229,204|255,212,170|255,195,136|255,178,102", # orange
"M_fed_8am_8-10wk" => "229,204,255|212,170,255|195,136,255|178,102,255", # purple
"M_fed_8pm_8-10wk" => "204,255,255|153,255,204|102,255,178|0,204,102",  #green
"F_fed_8am_12-16wk" => "255,204,255|255,170,255|255,136,255|255,102,255" # pink
);
# "255,204,204|255,170,170|255,136,136|255,102,102",  # red
# "240,230,140|232,209,104|225,187,68|218,165,32",  # yellow/brown


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

$supertrack_name = "st_NormDepthPB";
print OUT "track $supertrack_name\n";
print OUT "superTrack on show\n";
print OUT "shortLabel\tRNAseq Normalized Depth\n";
print OUT "longLabel\tRNAseq 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";
    $bwfile = "DATA/$ds\_$cnd\_rep$r.depthNorm.bw";
    $rr	= $r-1;	$rgb = $colors[$rr]; 
    print OUT "\t\ttrack st_NormDepth_$id\n";
    print OUT "\t\tbigDataUrl $bwfile\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);