1:  2:  3:  4:  5:  6:  7:  8:  9: 10: 11: 12: 13: 14: 15: 16: 17: 18: 19: 20: 21: 22: 23: 24: 25: 26: 27: 28: 29: 30: 31: 32: 33: 34: 35: 36: 37: 38: 39: 40: 41: 42: 43: 44: 45: 46: 47: 48: 49: 50: 51: 52: 53: 54: 55: 56: 57: 58: 59: 60: 61: 62: 63: 

<?php

namespace DispatchModule\Helpers;

/**
 * RnaplotHelper handles generation of svg visualizations of predicted structures using RNAplot
 */
class RnaplotHelper {

    /**
     * Deletes all rna plot svg files - if there's more than maximum limit (see maxRnaPlotSvgStored in config)
     *
     * @param string $rnaplotTemp path to rna plot temporary dir
     * @param int $maxfiles maximum ammount of files to keep or -1 for infinity
     */
    public static function deleteRnaplotFiles($rnaplotTemp, $maxfiles) {
        if ($maxfiles != -1) {
            $files = \Nette\Utils\Finder::findFiles('*.svg')->in($rnaplotTemp);
            if (count($files) > $maxfiles) {
                foreach ($files as $key => $file) {
                    $deleteFile = new \BaseModule\Models\FileModel($key);
                    $deleteFile->deleteFile();
                }
            }
        }
    }

    /**
     * Genereates RNAplot svg file
     *
     * @param string $rnaplotTemp path to rna plot temporary dir
     * @param string $idString svg name
     * @param string $sequence whole sequence
     * @param string $prediction predicted structure in dot-parent notation
     * @param int $maxfiles max svg files to keep in temporary storage (-1 for infinity)
     * @throws \BaseModule\Exceptions\BaseException
     */
    public static function generateRnaplot($rnaplotTemp, $idString, $sequence, $prediction, $maxfiles = 200) {
        // delete all files, if there's too many saved
        self::deleteRnaplotFiles($rnaplotTemp, $maxfiles);

        $resultPath = $rnaplotTemp . $idString . '_ss.svg';
        $timestamp = microtime(true);
        $tempFilePath = $rnaplotTemp . $timestamp . '.in';

        // create input file nad put content in it (proper format for RNAplot: >file_name\n sequence\n prediction)
        $inputFile = new \BaseModule\Models\FileModel($tempFilePath);
        $inputFile->createFile();
        $inputFile->appendContent(">$idString\n" . \DispatchModule\Sequence::makeSimplePrediction($sequence) . "\n" . $prediction);

        // run RNAplot
        exec("cd $rnaplotTemp && " . RNAPLOT_PATH . " -o svg < $tempFilePath", $out);

        // delete input temporary file
//        $inputFile->deleteFile();

        if (!file_exists($resultPath)) {
            throw new \BaseModule\Exceptions\BaseException('Something went wrong in the visualization factory');
        }
    }

}