]> git.donarmstrong.com Git - mothur.git/blobdiff - screenseqscommand.h
added sequence name to error string in fastq.info. Changed np_shannon to npshannon.
[mothur.git] / screenseqscommand.h
index 8f62ae48af466148c56b2b7b53aef7240cb95761..49d992ac193c7aa42ad4492c9190bda40671fc0b 100644 (file)
@@ -16,16 +16,26 @@ class ScreenSeqsCommand : public Command {
        
 public:
        ScreenSeqsCommand(string);
-       ~ScreenSeqsCommand();
-       int execute();
-       void help();
+       ScreenSeqsCommand();
+       ~ScreenSeqsCommand() {}
+       
+       vector<string> setParameters();
+       string getCommandName()                 { return "screen.seqs";                         }
+       string getCommandCategory()             { return "Sequence Processing";         }
+       string getHelpString(); 
+       string getCitation() { return "http://www.mothur.org/wiki/Screen.seqs"; }
+       string getDescription()         { return "enables you to keep sequences that fulfill certain user defined criteria"; }
+
+       int execute(); 
+       void help() { m->mothurOut(getHelpString()); }  
+       
        
 private:
 
        struct linePair {
-               unsigned long int start;
-               unsigned long int end;
-               linePair(unsigned long int i, unsigned long int j) : start(i), end(j) {}
+               unsigned long long start;
+               unsigned long long end;
+               linePair(unsigned long long i, unsigned long long j) : start(i), end(j) {}
        };
 
        vector<int> processIDS;   //processid
@@ -34,18 +44,27 @@ private:
        int screenNameGroupFile(set<string>);
        int screenGroupFile(set<string>);
        int screenAlignReport(set<string>);
+       int screenQual(set<string>);
+       int screenTaxonomy(set<string>);
        
-       int driver(linePair*, string, string, string, string, set<string>&);
-       int createProcesses(string, string, string, string, set<string>&);
+       int driver(linePair*, string, string, string, set<string>&);
+       int createProcesses(string, string, string, set<string>&);
        
        #ifdef USE_MPI
-       int driverMPI(int, int, MPI_File&, MPI_File&, MPI_File&, MPI_File&, vector<unsigned long int>&, set<string>&);
+       int driverMPI(int, int, MPI_File&, MPI_File&, MPI_File&, vector<unsigned long long>&, set<string>&);
        #endif
 
        bool abort;
-       string fastafile, namefile, groupfile, alignreport, outputDir;
-       int startPos, endPos, maxAmbig, maxHomoP, minLength, maxLength, processors;
+       string fastafile, namefile, groupfile, alignreport, outputDir, qualfile, taxonomy;
+       int startPos, endPos, maxAmbig, maxHomoP, minLength, maxLength, processors, criteria;
        vector<string> outputNames;
+       vector<string> optimize;
+       map<string, int> nameMap;
+       int readNames();
+       
+       int getSummary(vector<unsigned long long>&);
+       int createProcessesCreateSummary(vector<int>&, vector<int>&, vector<int>&, vector<int>&, vector<int>&, string);
+       int driverCreateSummary(vector<int>&, vector<int>&, vector<int>&, vector<int>&, vector<int>&, string, linePair*);       
 };
 
 #endif