#!/usr/bin/perl -w

# usage : perl mapper.pl
#Olivier Rué - Sarah Maman - 13/09/2012

use strict;

#genome de reference
my $input1 = $ARGV[3];
#pour les options "$collapse"
my $collapse = $ARGV[2];
my $output1 = $ARGV[0];
my $output2 = $ARGV[1];
#reads
my $input_nr = @ARGV[4..$#ARGV];
my $nbFastaFiles = @ARGV; 


#pour connaître le path vers le script mapper.pl
my $MAPPER = '/usr/local/bioinfo/src/galaxy/galaxy-dist/tools/sm_mirdeep2/mapper.pl';
my $cmd = '';
my $collapse_code = '';
my $collapse_code2 = '';

#Pour recuperer la date et l'heure
my ($sec,$min,$hour,$mday,$mon,$year,$wday,$yday,$isdst) = localtime(time); 
my $ALL =  $mday."-".($mon+1)."_".$hour."h".$min."mn".$sec."sec"; 

#Ouverture des fichiers entrants
open (IN, "<$input1") or die "Cannot open $input1 !";


#TESTS
if ($collapse eq "YES"){
#$collapse_code = "-m -s @ARGV[4..$#ARGV]_collapse.fa";
$collapse_code = "-m -s ${ALL}_collapse.fa";
$collapse_code2 = "-m -s ${ALL}_collapse.fa";      
}
else {
$collapse_code = "";
$collapse_code2 = "";
}
  

#convertir le fichier entrant nr en fichier fasta

#creer le fichier de sortie au format fasta

foreach $input_nr (@ARGV[4..$#ARGV]){
		#lire les lignes du fichier nr entrant
		@ARGV = ("$input_nr");
                my $inputNR = "fileWithoutIntra_".$ALL.".fasta";
                open FILE,">$inputNR";
		while (my $lignenr = <>){
        		my $seq = substr ($lignenr, 0, length($lignenr)-3);
			print FILE ">seq", "\n", $seq , "\n"; 	
		}	
 	}
close FILE || die "Problème à la fermeture du fichier fasta sortant";

#script mapper

if ($nbFastaFiles > 5){
	#Commande mapper POUR PLUSIEURS FICHIERS FASTA
	#perl mapper.pl TRUITE/config.txt -d -c -j -k ATCTCGTATGCCGTCTTCTGCTTG -l 18 -m -p TRUITE/Oncorhynchus_mykiss_V1.0.fa -s TESTS/reads_used.fa -t TESTS/mapped.arf -v
	#creer le fichier config.txt
#commande shell pour recuperer tous les fichiers fasta et ajouter un identifiant a coté (nombre max de 3 caracteres pas géré)
#find . -name \*.fa | perl -ne 'chomp;($n) = $_ =~ /.\/(\S+)/;print $n,"\t","s".$i,"\n";$i++;' > config.txt

	my $config = "config_".$ALL.".txt";
	open FILE,">$config";

	my $cpt=1;

	#foreach my $file (@ARGV){
        #foreach my $file (@ARGV[4..$#ARGV]){
	 foreach my $file (@ARGV[3..$#ARGV]){
		my $id="s";
		if(length($cpt)<2){
			$id = $id."0".$cpt;
		}
		else{
			$id = $id."$cpt";
		}
		print FILE $file,"\t",$id,"\n"; 
		$cpt++;
	}

	#close FILE;



	#lancer la commande (sans nettoyer les adaptateurs car cela est déjà fait avec le cutadapt
	$cmd = "($MAPPER $config -d -c -j $collapse_code2 -p $input1 -t ${ALL}_mapped.arf -v) >& ./mapperPlrsFasta.log 2>&1";
	system $cmd;
print "\n *Commande - plusieurs fastas: $cmd \n";

		if (! -e "${ALL}_mapped.arf"){print STDERR "Echec a la generation du fichier mapped.arf\n";}
		else {$cmd = "(cp -a \"${ALL}_mapped.arf\" $output1) >& ./cp_mapper_output1.log 2>&1";
		system $cmd;}

		if (! -e "${ALL}_collapse.fa"){print STDERR "Echec a la generation du fichier genome.fa \n";}
		else {$cmd = "(cp -a \"${ALL}_collapse.fa\" $output2) >& ./cp_mapper_output2.log 2>&1";
		system $cmd;} 

close FILE;
                
}
else{
	#commande mapper POUR UN SEUL FASTA
	#perl mapper.pl TRUITE/s_1_uT21.fa -c -j -k TCGTATGCCGTCTTCTGCTTGT -l 18 -m -p TRUITE/Oncorhynchus_mykiss_V1.0.fa -s TESTS/ reads_collapsed.fa -t TESTS/reads_collapsed_vs_genome.arf -v
	#$cmd = "($MAPPER $input2 -c -j $collapse_code -p $input1 -t ${input2}_reads_collapsed_vs_genome.arf -v) >& ./mapper1fasta.log 2>&1";
 	#$cmd = "($MAPPER @ARGV[4..$#ARGV] -c -j $collapse_code -p $input1 -t @ARGV[4..$#ARGV]_reads_collapsed_vs_genome.arf -v) >& ./mapper1fasta.log 2>&1";
        $cmd = "($MAPPER @ARGV[4..$#ARGV] -c -j $collapse_code -p $input1 -t ${ALL}_reads_collapsed_vs_genome.arf -v) >& ./mapper1fasta.log 2>&1";
	system $cmd; 
#print "\n ** Voici la commande utilisée avec un fasta en entree: $cmd \n";


		if (! -e "${ALL}_reads_collapsed_vs_genome.arf"){print STDERR "Echec a la generation du fichier collapse.arf\n";}
		else {$cmd = "(cp -a \"${ALL}_reads_collapsed_vs_genome.arf\" $output1) >& ./cp_mapper_output1.log 2>&1";
		system $cmd;}


		if (! -e "${ALL}_collapse.fa"){print STDERR "Echec a la generation du fichier genome.fa \n";}
		else {$cmd = "(cp -a \"${ALL}_collapse.fa\" $output2) >& ./cp_mapper_output2.log 2>&1";
		system $cmd;} 

}
close( IN );
