Method calculateCopyNumberProbs_and_exomeLength

src/GenomeCopyNumber.cpp:1759–2120

void GenomeCopyNumber::calculateCopyNumberProbs_and_exomeLength(int breakPointType)
    {
	//estimationOfExomeSize_ = 0;
	vector<ChrCopyNumber>::iterator it;
	unsigned long long count = 0;
    int endsSize=0;
	string NormalBAF,NormalBAF_XY;
    float normalXYploidy=ploidy_*0.5; //will only use it when the genome is male
    estimationOfGenomeSize_=0;
    if (hasBAF_)
        {
            NormalBAF = getNormalBAFforPloidy(ploidy_);
            NormalBAF_XY=getXYBAFforPloidy(ploidy_);//will only use it when the genome is male
        }

    CNVs_.clear();
    copyNumberProbs_.clear();

	for ( it=chrCopyNumber_.begin() ; it != chrCopyNumber_.end(); it++ )
        {
		float previousLevel = NA;
		int lengthOfPreviousLevel = 0;
		float nextLevel = NA;
		int lengthOfNextLevel = 0;
		it->setIsSmoothed(true);
		string chr = it->getChromosome();
		cout << "..Annotation of CNVs for "<< chr <<"\n";
		//will use these three variables to collect CNVs:
		int start = 0;
		int end = NA;
		int cnumber = NA;
		string BAFSym = "-";
        string BAFprev = "-";
        string BAFnext = "-";
        string lBAF = "-";

        float lUncertainty = NA;
        float BAFUncertainty = NA;

        float normalLevel=1;
        if (sex_.compare("XY")==0 && (it->getChromosome().find("X")!=string::npos || it->getChromosome().find("Y")!=string::npos)) {
            normalLevel = 0.5;
        }
		for (int i = 0; i<(int)it->getMedianValues().size(); i++) {
			float level = it->getLevelAt(i, ploidy_);
			int fragmentLength =  it->getFragmentLengths()[i];
            if (hasBAF_) {
                BAFSym = it->getBAFsymbPerFrg(i);
                BAFUncertainty = it->getEstimatedBAFuncertaintyAtBin(i);
//                cout << "..Control: read "<< BAFSym << "\n";
            }

			if (level == NA)
                {
				int nextIndex = it->nextNoNAIndex(i, ploidy_,0);
				BAFUncertainty = NA;
				if (nextIndex == NA) {
					nextLevel = NA;