ln -s '$matrix_h5_cooler' 'matrix.$matrix_h5_cooler.ext' &&
        mkdir output &&
        hicPCA --matrix 'matrix.$matrix_h5_cooler.ext'

        
        #if $chromosome_list:
            #set $chromosomes = []
            #for $chrom in $chromosome_list:
                #silent $chromosomes.append("'%s'" % $chrom.chromosomes)
            #end for
            #set $chromosome_string = ' '.join($chromosomes)
            --chromosomes $chromosome_string
        #end if
    

        #set $eigenvectors = ' '.join([ '\'%s\'' % $eigenvector for $eigenvector in str($whichEigenvectors).split(' ') ])
       --whichEigenvectors $eigenvectors
       
        #set $eigenvectors_output = ' '.join([ '\'output/pca%s.%s\'' % ($eigenvector, $outputFormat) for $eigenvector in str($whichEigenvectors).split(' ') ])

        --outputFileName $eigenvectors_output
        --format $outputFormat

        #if $ligation_factor:
            $ligation_factor
        #end if

        #if $ignoreMaskedBins:
            $ignoreMaskedBins
        #end if

        #if $extra_track_conditional.extra_track_selection == 'gene_density':
            --extraTrack '$extra_track_conditional.extraTrack'
        #elif $extra_track_conditional.extra_track_selection == 'histon_mark':
            --extraTrack '$extra_track_conditional.extraTrack'
            --histonMarkType '$extra_track_conditional.histonMarkType'
        #end if

        #if $pearsonMatrix:
            --pearsonMatrix pearson.$matrix_h5_cooler.ext
        #end if

        #if $obsexpMatrix:
            --obsexpMatrix obsexp.$matrix_h5_cooler.ext
        #end if

        #if $pearsonMatrix:
            && mv pearson.$matrix_h5_cooler.ext pearson
        #end if

        #if $obsexpMatrix:
            && mv obsexp.$matrix_h5_cooler.ext obsexp
        #end if