chimeraslayer.h

   1 #ifndef CHIMERASLAYER_H
   2 #define CHIMERASLAYER_H
   3
   4 /*
   5  *  chimeraslayer.h
   6  *  Mothur
   7  *
   8  *  Created by westcott on 9/25/09.
   9  *  Copyright 2009 Schloss Lab. All rights reserved.
  10  *
  11  */
  12
  13
  14 #include "chimera.h"
  15 #include "maligner.h"
  16 #include "slayer.h"
  17
  18
  19
  20 //***********************************************************************/
  21 //This class was modeled after the chimeraSlayer written by the Broad Institute
  22 /***********************************************************************/
  23
  24
  25 class ChimeraSlayer : public Chimera {
  26
  27         public:
  28                 ChimeraSlayer(string, string, string, int, int, int, int, float, int, int, int, int, int, int, int, int, bool);
  29                 ChimeraSlayer(string, string, string, string, string, int, int, int, int, float, int, int, int, int, int, int, int, int, bool);
  30
  31                 ~ChimeraSlayer();
  32
  33                 int getChimeras(Sequence*);
  34                 int print(ostream&, ostream&);
  35                 void printHeader(ostream&);
  36                 int doPrep();
  37
  38                 #ifdef USE_MPI
  39                 int print(MPI_File&, MPI_File&);
  40                 #endif
  41
  42         private:
  43                 Sequence* querySeq;
  44                 DeCalculator* decalc;
  45                 map<int, int>  spotMap;
  46                 Database* databaseRight;
  47                 Database* databaseLeft;
  48                 map<string, vector<string> > nameMapRank;  //sequence name to rank so you can construct a template of the abundant sequences if the user uses itself as template
  49
  50                 vector<data_struct>  chimeraResults;
  51                 string chimeraFlags, searchMethod, fastafile, includeAbunds;
  52                 bool realign;
  53                 int window, numWanted, kmerSize, match, misMatch, minSim, minCov, minBS, minSNP, parents, iters, increment;
  54                 float divR;
  55
  56                 void printBlock(data_struct, string, ostream&);
  57                 string getBlock(data_struct, string);
  58                 int readNameFile(string);
  59                 vector<Sequence*> getTemplate(Sequence*);
  60
  61 };
  62
  63 /************************************************************************/
  64
  65 #endif
  66
  67