//**********************************************************************************************************************
-FlowData::FlowData(int numFlows, float signal, float noise, int maxHomoP) :
- numFlows(numFlows), signalIntensity(signal), noiseIntensity(noise), maxHomoP(maxHomoP){
+FlowData::FlowData(int numFlows, float signal, float noise, int maxHomoP, string baseFlow) :
+ numFlows(numFlows), signalIntensity(signal), noiseIntensity(noise), maxHomoP(maxHomoP), baseFlow(baseFlow){
try {
m = MothurOut::getInstance();
flowData.assign(numFlows, 0);
- baseFlow = "TACG";
+// baseFlow = "TACG";
seqName = "";
locationString = "";
}
bool FlowData::getNext(ifstream& flowFile){
try {
-
- string lengthString;
- string flowString;
-
- flowFile >> seqName >> endFlow;
- for(int i=0;i<numFlows;i++) { flowFile >> flowData[i]; }
-
- updateEndFlow();
+ flowFile >> seqName >> endFlow;
+ //cout << "in Flowdata " + seqName << endl;
+ for(int i=0;i<numFlows;i++) { flowFile >> flowData[i]; }
+ //cout << "in Flowdata read " << seqName + " done" << endl;
+ updateEndFlow();
translateFlow();
m->gobble(flowFile);
void FlowData::updateEndFlow(){
try{
- int currLength = 0;
+ //int currLength = 0;
float maxIntensity = (float) maxHomoP + 0.49;
int deadSpot = 0;
try{
maxFlows = mF;
- if(endFlow > maxFlows){ endFlow = maxFlows; }
+ if(endFlow > maxFlows){ endFlow = maxFlows; }
+ translateFlow();
}
catch(exception& e) {
}
//**********************************************************************************************************************
+
+string FlowData::getName(){
+
+ try{
+ return seqName;
+ }
+ catch(exception& e) {
+ m->errorOut(e, "FlowData", "getName");
+ exit(1);
+ }
+}
+
+//**********************************************************************************************************************