no_peaks.cpp

#include "no_peaks.h"
#include <unordered_set>
#include <iostream>

class ZeroPeaksCounter {
private:
	int all_no_peaks, good_no_peaks;

	deconv_program program;
	
	std::unordered_map < int, Scan > &scan_map;

	const double EVALUE_BORDER;
public:
	ZeroPeaksCounter(deconv_program new_program, std::unordered_map < int, Scan > &new_map,
		double new_border): all_no_peaks(0), good_no_peaks(0), program(new_program),
		scan_map(new_map), EVALUE_BORDER(new_border) {};

	void operator() (Scan &scan) {
		std::unordered_map < int, Scan>::iterator prediction = scan_map.find(scan.id);
		if (scan.peaks.size() == 0 && prediction != scan_map.end()) {
			if (prediction->second.e_value <= EVALUE_BORDER) {
				good_no_peaks++;
				std::cout << program << ' ' << scan.id << ' ' << prediction->second.e_value << std::endl;
			}
			all_no_peaks++;
		}
	}

	int get_all_no_peaks() {
		return all_no_peaks;
	}

	int get_dood_no_peaks() {
		return good_no_peaks;
	}
};

void find_no_peaks_scans(std::unordered_map < int, Scan > &theoretic_scans_map,
	deconv_program program, std::string filename, double evalue_border) {
	std::cout << "Looking for scans without peaks in " << filename << std::endl;
	ZeroPeaksCounter zero_peaks(program, theoretic_scans_map, evalue_border);
	go_through_mgf(program, filename, zero_peaks);
	std::cout << "Found " << zero_peaks.get_all_no_peaks() << " scans without peaks including "
		<< zero_peaks.get_dood_no_peaks() << " good scans." << std::endl;
}

void check_no_peaks_scans(std::vector < std::string > filenames, double evalue_border) {
	for (std::string pref : filenames) {
		ScansMapCreator map_creator;
		go_through_tsv(pref + TSV_SUF, map_creator);
		std::unordered_map < int, Scan > theoretic_scans_map = map_creator.get_map();

		find_no_peaks_scans(theoretic_scans_map, MS_Deconv, pref + MSDECONV_SUF, evalue_border);
		find_no_peaks_scans(theoretic_scans_map, Thermo_Xtract, pref + XTRACT_SUF, evalue_border);
	}
	std::cout << std::endl;
}