X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=chimeraslayer.h;h=7bf663afeadeb3517daf4f9bfa25135fbf872cab;hp=14a22beed74f0d0fe6dbea9e5d255117d3187008;hb=cf9987b67aa49777a4c91c2d21f96e58bf17aa82;hpb=9c23307c583d4e8595f75278c13e708788f2f058 diff --git a/chimeraslayer.h b/chimeraslayer.h index 14a22be..7bf663a 100644 --- a/chimeraslayer.h +++ b/chimeraslayer.h @@ -15,36 +15,60 @@ #include "maligner.h" #include "slayer.h" -/***********************************************************************/ + + +//***********************************************************************/ //This class was modeled after the chimeraSlayer written by the Broad Institute /***********************************************************************/ - class ChimeraSlayer : public Chimera { public: - ChimeraSlayer(string, bool, string); + ChimeraSlayer(string, string, bool, string, int, int, int, int, float, int, int, int, int, int, int, int, int, bool, string, int); + ChimeraSlayer(string, string, bool, map&, string, int, int, int, int, float, int, int, int, int, int, int, int, int, bool, string, int); + ChimeraSlayer(string, string, bool, map&, string, int, int, int, int, float, int, int, int, int, int, int, int, int, bool, string, int, bool); + ~ChimeraSlayer(); int getChimeras(Sequence*); - void print(ostream&, ostream&); + Sequence print(ostream&, ostream&); + Sequence print(ostream&, ostream&, data_results, data_results); void printHeader(ostream&); - void doPrep(); + int doPrep(); + int getNumNoParents() { return numNoParents; } + data_results getResults() { return printResults; } + + #ifdef USE_MPI + Sequence print(MPI_File&, MPI_File&); + Sequence print(MPI_File&, MPI_File&, data_results, data_results, bool&); + #endif private: - Sequence* querySeq; - DeCalculator* decalc; - Maligner* maligner; - Slayer* slayer; - map spotMap; + Sequence querySeq; + Sequence trimQuery; + DeCalculator decalc; Database* databaseRight; Database* databaseLeft; - + map priority; //for template=self, seqname, seqAligned, abundance + set chimericSeqs; //for template=self, so we don't add chimeric sequences to the userTemplate set + int numNoParents, threadID; + vector chimeraResults; - string chimeraFlags, searchMethod, fastafile; - bool realign; + data_results printResults; + string chimeraFlags, searchMethod, fastafile, blastlocation; + bool realign, trimChimera; + int window, numWanted, kmerSize, match, misMatch, minSim, minCov, minBS, minSNP, parents, iters, increment; + float divR; - void printBlock(data_struct, ostream&); + void printBlock(data_struct, string, ostream&); + void printBlock(data_results, data_results, bool, bool, string, ostream&); + string getBlock(data_struct, string); + string getBlock(data_results, data_results, bool, bool, string); + //int readNameFile(string); + vector getTemplate(Sequence, vector&); + vector getRefSeqs(Sequence, vector&, vector&); + vector getBlastSeqs(Sequence, vector&, int); + vector getKmerSeqs(Sequence, vector&, int); };