File: OpenAssemblerChooser.h

package info (click to toggle)
ray 2.3.1-9
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 6,008 kB
  • sloc: cpp: 49,973; sh: 339; makefile: 281; python: 168
file content (60 lines) | stat: -rw-r--r-- 1,823 bytes parent folder | download | duplicates (5)
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
/*
 	Ray
    Copyright (C) 2010, 2011, 2012 Sébastien Boisvert

	http://DeNovoAssembler.SourceForge.Net/

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, version 3 of the License.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You have received a copy of the GNU General Public License
    along with this program (gpl-3.0.txt).  
	see <http://www.gnu.org/licenses/>

*/

#ifndef _OpenAssemblerChooser
#define _OpenAssemblerChooser

#include "Chooser.h" // for IMPOSSIBLE_CHOICE
#include "NovaEngine.h"

/**
 * de Bruijn heuristic to choose extension direction in a graph, described in paper 
 * 	
 * 	Ray: simultaneous assembly of reads from a mix of high-throughput sequencing technologies.
 * 	Sébastien Boisvert, François Laviolette, and Jacques Corbeil.
 * 	Journal of Computational Biology (Mary Ann Liebert, Inc. publishers).
 * 	November 2010, 17(11): 1519-1533.
 * 	doi:10.1089/cmb.2009.0238
 * 	http://dx.doi.org/doi:10.1089/cmb.2009.0238
 * \author Sébastien Boisvert
 */
class OpenAssemblerChooser{
	NovaEngine m_novaEngine;

	double m_singleEndMultiplicator;
	double m_pairedEndMultiplicator;
	void updateMultiplicators();

	int getWinner(vector<set<int> >*battleVictories,int choices);
	
	/**
 * choose where to go based on coverage
 */
	void chooseWithCoverage(ExtensionData*ed,int minCoverage,vector<set<int> >*battleVictories);
public:
	int choose(ExtensionData*m_ed,Chooser*m_c,int m_minimumCoverage,
	Parameters*parameters);

	void constructor();

};

#endif