annotate BEDTools-Version-2.14.3/src/annotateBed/annotateBed.h @ 0:dfcd8b6c1bda

Uploaded
author aaronquinlan
date Thu, 03 Nov 2011 10:25:04 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
1 /*****************************************************************************
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
2 annotateBed.h
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
3
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
4 (c) 2009 - Aaron Quinlan
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
5 Hall Laboratory
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
6 Department of Biochemistry and Molecular Genetics
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
7 University of Virginia
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
8 aaronquinlan@gmail.com
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
9
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
10 Licenced under the GNU General Public License 2.0 license.
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
11 ******************************************************************************/
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
12 #ifndef ANNOTATEBED_H
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
13 #define ANNOTATEBED_H
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
14
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
15 #include "bedFile.h"
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
16 #include <vector>
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
17 #include <algorithm>
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
18 #include <iostream>
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
19 #include <iomanip>
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
20 #include <fstream>
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
21 #include <stdlib.h>
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
22
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
23 using namespace std;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
24
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
25 //************************************************
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
26 // Class methods and elements
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
27 //************************************************
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
28 class BedAnnotate {
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
29
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
30 public:
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
31
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
32 // constructor
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
33 BedAnnotate(const string &mainFile, const vector<string> &annoFileNames,
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
34 const vector<string> &annoTitles, bool sameStrand, bool diffStrand, bool reportCounts, bool reportBoth);
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
35
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
36 // destructor
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
37 ~BedAnnotate(void);
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
38
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
39 // annotate the master file with all of the annotation files.
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
40 void AnnotateBed();
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
41
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
42 private:
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
43
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
44 // input files.
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
45 string _mainFile;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
46 vector<string> _annoFileNames;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
47 vector<string> _annoTitles;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
48
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
49 // instance of a bed file class.
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
50 BedFile *_bed;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
51 vector<BedFile*> _annoFiles;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
52
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
53 // do we care about strandedness when counting coverage?
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
54 bool _sameStrand;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
55 bool _diffStrand;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
56
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
57 bool _reportCounts;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
58 bool _reportBoth;
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
59
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
60 // private function for reporting coverage information
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
61 void ReportAnnotations();
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
62
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
63 void OpenAnnoFiles();
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
64
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
65 void CloseAnnoFiles();
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
66
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
67 void PrintHeader();
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
68
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
69 void InitializeMainFile();
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
70 };
dfcd8b6c1bda Uploaded
aaronquinlan
parents:
diff changeset
71 #endif /* ANNOTATEBED_H */