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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
#ifndef _BLASR_HDF_CMP_REF_ALIGNMENT_GROUP_HPP_
#define _BLASR_HDF_CMP_REF_ALIGNMENT_GROUP_HPP_
#include <map>
#include <string>
#include "H5Cpp.h"
#include "HDFData.hpp"
#include "HDFGroup.hpp"
#include "HDFCmpExperimentGroup.hpp"
using namespace std;
class HDFCmpRefAlignmentGroup {
public:
HDFGroup refGroup;
string refGroupName;
vector<HDFCmpExperimentGroup*> readGroups;
HDFAtom<string> annotationStringAtom;
map<string,int> experimentNameToIndex;
// A RefAlignmentGroup may contain one or more
// ExperimentGroups. The following shows a
// RefAlignmentGroup containing two ExperimentGroups.
// /ref00001/m121219_103658_42194_c000447032559900001500000112311426_s1_p0
// /ref00001/m121219_103658_42194_c000447032559900001500000112311426_s2_p0
//
// But these ExperimentGroups are not necessarily
// grouped by movie. For example, the following groups
// can be seen in deep-sorted cmp.h5
// /ref00001/rg8953-0
// /ref00001/rg2453-1
int Initialize(H5::CommonFG &group, string _refGroupName) {
refGroupName = _refGroupName;
refGroup.Initialize(group, _refGroupName);
// annotationStringAtom.Initialize(refGroup.group, "annotationString");
return 1;
}
void Create(HDFGroup parent, string refGroupNameP) {
refGroupName = refGroupNameP;
parent.AddGroup(refGroupName);
refGroup.Initialize(parent, refGroupName);
}
HDFCmpExperimentGroup *GetExperimentGroup(string readGroupName) {
//
// In contrast to initialization, only create one group.
//
map<string,int>::iterator it = experimentNameToIndex.find(readGroupName);
if (it != experimentNameToIndex.end()) {
assert(it->second < int(readGroups.size()));
return readGroups[it->second];
}
//
// Allocate the new group structure
//
int newReadGroupIndex = readGroups.size();
HDFCmpExperimentGroup* readGroupPtr = new HDFCmpExperimentGroup;
if (readGroupPtr == nullptr) {cout << "ERROR, failed to allocate memory for HDFCmpExperimentGroup!" << endl; exit(1);}
readGroups.push_back(readGroupPtr);
experimentNameToIndex[readGroupName] = newReadGroupIndex;
//
// Now add it to the cmp.h5 file.
//
if (readGroupPtr->Create(refGroup, readGroupName) == 0) {
delete readGroupPtr;
experimentNameToIndex[readGroupName] = -1;
return NULL;
}
return readGroupPtr;
}
bool ContainsExperimentGroup(string readGroupName) {
return experimentNameToIndex.find(readGroupName) != experimentNameToIndex.end();
}
HDFCmpExperimentGroup* InitializeExperimentGroup(string experimentGroupName, set<string> &includedFields) {
if (refGroup.ContainsObject(experimentGroupName)) {
HDFCmpExperimentGroup* newGroup = new HDFCmpExperimentGroup;
if (newGroup == nullptr) {cout << "ERROR, failed to allocate memory for HDFCmpExperimentGroup!" << endl; exit(1);}
if (newGroup->Initialize(refGroup, experimentGroupName, includedFields) == 0) {
cout << "ERROR, could not initialize the exp group." << endl;
exit(1);
}
experimentNameToIndex[experimentGroupName] = readGroups.size();
readGroups.push_back(newGroup);
return newGroup;
}
else {
return NULL;
}
}
HDFCmpExperimentGroup* InitializeExperimentGroup(string experimentGroupName) {
set<string> EMPTYIncludedFields;
return InitializeExperimentGroup(experimentGroupName, EMPTYIncludedFields);
}
};
#endif
|
