source: trunk/grails-app/domain/nl/tno/metagenomics/AssaySample.groovy @ 2

package nl.tno.metagenomics

/**
 * Represents a samples that is used in an assay.
 * 
 * @author Robert Horlings (robert@isdat.nl)
 *
 */
class AssaySample {
        // To be computed in run time
        private long _numSequences = -1;
        private float _averageQuality = -1.0;
        
        Integer numUniqueSequences      // Number of unique sequences / OTUs. Is only available after preprocessing

        String oligoNumber              // Oligonumber used to identify the sample
        String tagSequence              // Tag originally used to identify the sample

        static belongsTo  = [ assay: Assay, sample: Sample ]
        static hasMany    = [ sequenceData: SequenceData ]

        static constraints = {
                numUniqueSequences(nullable: true)
                oligoNumber(nullable: true)
                tagSequence(nullable: true)
        }

        static mapping = {
                columns {
                        numSequences index:'numsequences_idx'
                }
        }

        /**
         * Returns the number of files in the system, belonging to this 
         * assay-sample combination. 
         * 
         * @return
         */
        public int numFiles() {
                if( !sequenceData )
                        return 0

                int numFiles = 0;
                sequenceData.each { numFiles += it.numFiles() }

                return numFiles;
        }

        /**
         * Returns the number of sequences in the files on the system, belonging to this
         * assay-sample combination.
         *
         * @return
         */
        public long numSequences() {
                if( _numSequences > -1 )
                        return _numSequences;
                        
                if( !sequenceData )
                        return 0

                long numSequences = 0;
                sequenceData.each { numSequences += it.numSequences }

                // Save as cache
                _numSequences = numSequences;

                return numSequences;
        }

        /**
         * Returns the average quality of the sequences in the files on the system, 
         * belonging to this assay-sample combination.
         *
         * @return
         */
        public float averageQuality() {
                if( _averageQuality > -1 )
                        return _averageQuality;

                if( !sequenceData )
                        return 0.0

                int numSequences = 0;
                float averageQuality = 0.0;

                sequenceData.each {
                        numSequences += it.numSequences
                        averageQuality = averageQuality + ( it.averageQuality - averageQuality ) / numSequences * it.numSequences;
                }

                // Save as cache
                _averageQuality = averageQuality;
                
                return averageQuality;
        }
}