/*
 *                    BioJava development code
 *
 * This code may be freely distributed and modified under the
 * terms of the GNU Lesser General Public Licence.  This should
 * be distributed with the code.  If you do not have a copy,
 * see:
 *
 *      http://www.gnu.org/copyleft/lesser.html
 *
 * Copyright for this code is held jointly by the individual
 * authors.  These should be listed in @author doc comments.
 *
 * For more information on the BioJava project and its aims,
 * or to join the biojava-l mailing list, visit the home page
 * at:
 *
 *      http://www.biojava.org/
 *
 * Created on June 29, 2010
 * Author: Mark Chapman
 */

package org.biojava.nbio.alignment;

import org.biojava.nbio.core.alignment.matrices.SubstitutionMatrixHelper;
import org.biojava.nbio.alignment.template.GapPenalty;
import org.biojava.nbio.alignment.template.PairwiseSequenceScorer;
import org.biojava.nbio.core.alignment.template.SubstitutionMatrix;
import org.biojava.nbio.core.exceptions.CompoundNotFoundException;
import org.biojava.nbio.core.sequence.ProteinSequence;
import org.biojava.nbio.core.sequence.compound.AminoAcidCompound;
import org.junit.Before;
import org.junit.Test;

import static org.junit.Assert.assertEquals;
import static org.junit.Assert.assertNotNull;

public class FractionalIdentityScorerTest {

	private static final double PRECISION = 0.00000001;

	private ProteinSequence query, target;
	private GapPenalty gaps;
	private SubstitutionMatrix<AminoAcidCompound> blosum62;
	private NeedlemanWunsch<ProteinSequence, AminoAcidCompound> alignment, self;
	private PairwiseSequenceScorer<ProteinSequence, AminoAcidCompound> scorer1, scorer2;

	@Before
	public void setup() throws CompoundNotFoundException {
		query = new ProteinSequence("ARND");
		target = new ProteinSequence("RDG");
		gaps = new SimpleGapPenalty((short) 10, (short) 1);
		blosum62 = SubstitutionMatrixHelper.getBlosum62();
		alignment = new NeedlemanWunsch<ProteinSequence, AminoAcidCompound>(query, target, gaps, blosum62);
		self = new NeedlemanWunsch<ProteinSequence, AminoAcidCompound>(query, query, gaps, blosum62);
		scorer1 = new FractionalIdentityScorer<ProteinSequence, AminoAcidCompound>(alignment);
		scorer2 = new FractionalIdentityScorer<ProteinSequence, AminoAcidCompound>(self);
	}

	@Test
	public void testFractionalIdentityScorerPairwiseSequenceAlignerOfSC() {
		assertNotNull(new FractionalIdentityScorer<ProteinSequence, AminoAcidCompound>(alignment));
	}

	@Test
	public void testFractionalIdentityScorerSequencePairOfSC() {
		assertNotNull(new FractionalIdentityScorer<ProteinSequence, AminoAcidCompound>(alignment.getPair()));
	}

	@Test
	public void testGetQuery() {
		assertEquals(scorer1.getQuery(), query);
		assertEquals(scorer2.getQuery(), query);
	}

	@Test
	public void testGetTarget() {
		assertEquals(scorer1.getTarget(), target);
		assertEquals(scorer2.getTarget(), query);
	}

	@Test
	public void testGetMaxScore() {
		assertEquals(scorer1.getMaxScore(), 4, PRECISION);
		assertEquals(scorer2.getMaxScore(), 4, PRECISION);
	}

	@Test
	public void testGetMinScore() {
		assertEquals(scorer1.getMinScore(), 0, PRECISION);
		assertEquals(scorer2.getMinScore(), 0, PRECISION);
	}

	@Test
	public void testGetScore() {
		assertEquals(scorer1.getScore(), 1, PRECISION);
		assertEquals(scorer2.getScore(), 4, PRECISION);
	}

}
