X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=flowdata.cpp;h=1fe7d7faf1e31ee486890980c4731ddebfdca8fb;hb=96dbe925073caefaed6e6db85659c144a806aeb1;hp=eb6f9bfec3a88b05bf2dcd03a751b01b320430d0;hpb=73fb1fe84d565ea9a0ff793de990546b177310c8;p=mothur.git diff --git a/flowdata.cpp b/flowdata.cpp index eb6f9bf..1fe7d7f 100644 --- a/flowdata.cpp +++ b/flowdata.cpp @@ -19,14 +19,14 @@ FlowData::~FlowData(){ /* do nothing */ } //********************************************************************************************************************** -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 = ""; } @@ -42,17 +42,16 @@ FlowData::FlowData(int numFlows, float signal, float noise, int maxHomoP) : bool FlowData::getNext(ifstream& flowFile){ try { - - string lengthString; - string flowString; - - flowFile >> seqName >> endFlow; - for(int i=0;i> flowData[i]; } - - updateEndFlow(); - translateFlow(); - - m->gobble(flowFile); + flowFile >> seqName >> endFlow; + if (seqName.length() != 0) { + //cout << "in Flowdata " + seqName << endl; + for(int i=0;i> flowData[i]; } + //cout << "in Flowdata read " << seqName + " done" << endl; + updateEndFlow(); + translateFlow(); + m->gobble(flowFile); + }else{ m->mothurOut("Error in reading your flowfile, at position " + toString(flowFile.tellg()) + ". Blank name."); m->mothurOutEndLine(); } + if(flowFile){ return 1; } else { return 0; } } @@ -138,7 +137,8 @@ void FlowData::capFlows(int mF){ try{ maxFlows = mF; - if(endFlow > maxFlows){ endFlow = maxFlows; } + if(endFlow > maxFlows){ endFlow = maxFlows; } + translateFlow(); } catch(exception& e) { @@ -212,3 +212,16 @@ void FlowData::printFlows(ofstream& outFlowFile, string scrapCode){ } //********************************************************************************************************************** + +string FlowData::getName(){ + + try{ + return seqName; + } + catch(exception& e) { + m->errorOut(e, "FlowData", "getName"); + exit(1); + } +} + +//**********************************************************************************************************************