def check_file(filename):
if (not os.access(filename, os.R_OK)):
print "%s is not accessible." % filename
print usage_mesg
sys.exit(1)
filename_mzXML = sys.argv[1]
filename_hit_list = sys.argv[2]
check_file(filename_mzXML)
check_file(filename_hit_list)
ms1_mz2i = dict()
mzXML = MzXML()
mzXML.parse_file(filename_mzXML)
for tmp_ms1 in mzXML.MS1_list:
scan_id = tmp_ms1.id
ms1_mz2i[scan_id] = dict()
for i in range(0, len(tmp_ms1.mz_list)):
tmp_mz = tmp_ms1.mz_list[i]
tmp_i = tmp_ms1.intensity_list[i]
ms1_mz2i[scan_id][tmp_mz] = tmp_i
ms2_precursor_i = dict()
for tmp_ms2 in mzXML.MS2_list:
scan_id = tmp_ms2.id
ms2_precursor_i[scan_id] = tmp_ms2.precursor_intensity
mz2scan = dict()
