#!/usr/bin/perl

@A_proj = ("EMC_Mitochondrial_Disease", "EMC_Bone", "ULaval_Obesity.IA", "ULaval_Obesity.SC");
@A_dir  = ("forward", "reverse");

foreach $proj (@A_proj) {

	foreach $dir (@A_dir) {

	        $txt = qq {
track $proj.$dir
type bigWig
container multiWig
shortLabel $proj.$dir
longLabel $proj.$dir
visibility full
aggregate stacked
autoScale on
showSubtrackColorOnUi on
maxHeightPixels 100:48:8
priority 1
};
		print $txt;
	}

	$cnt = 0;

	$val = 150;

	$proj2 = $proj;
	if ($proj =~ /ULaval/) { $proj2 = "ULaval_Obesity"; }

	my @A_files;
	@A_files = </mnt/www/webroot/hg19/rna_seq/v_2/$proj2/hg19/wiggle_tracks/*forward.bigwig>;

	foreach $file (@A_files) {

		($sample) = $file =~ /wiggle_tracks\/(.*)\.forward\.bigwig/;
		next if ($sample !~ /IA/ && $proj eq "ULaval_Obesity.IA");
		next if ($sample !~ /SC/ && $proj eq "ULaval_Obesity.SC");
		
		foreach $dir (@A_dir) {
			if ($dir eq "forward") { $color = "0,0,$val"; }
			else                   { $color = "$val,0,0"; }

			$txt = qq{
track $sample.$dir
bigDataUrl http://pastinenlab.genome.mcgill.ca/hg19/rna_seq/v_2/$proj2/hg19/wiggle_tracks/$sample.$dir.bigwig
shortLabel $sample.$dir
longLabel $sample.$dir
parent $proj.$dir
type bigWig
autoScale on
color $color
};
			print "$txt\n";
			$val += 25;
		}
		last if ($cnt > 5);
	}
}