@@ -70,17 +70,19 @@ process(path; df, out, mode, εt, εm, fmt_target, fmts) = begin
7070 s = s .& ((df. start .- εt) .≤ ms. retention_time .≤ (df. stop .+ εt))
7171 return map (r -> UniMZ. Ion (r. mz, r. z), eachrow (df[s, :]))
7272 end
73- @info " MS2 preprocessing using isotopic pattern..."
74- M2 = @showprogress map (M2) do ms
75- ps = map (ms. peaks) do p
76- if any (z -> UniMZ. query_ε (ms. peaks, p. mz - Δ / z, εm) |> ! isempty, 1 : 3 )
77- return UniMZ. Peak[]
73+ if false
74+ @info " MS2 preprocessing using isotopic pattern..."
75+ M2 = @showprogress map (M2) do ms
76+ ps = map (ms. peaks) do p
77+ if any (z -> UniMZ. query_ε (ms. peaks, p. mz - Δ / z, εm) |> ! isempty, 1 : 3 )
78+ return UniMZ. Peak[]
79+ end
80+ zs = filter (z -> UniMZ. query_ε (ms. peaks, p. mz + Δ / z, εm) |> ! isempty, 1 : 3 )
81+ zs = isempty (zs) ? [1 ] : zs
82+ return map (z -> UniMZ. Peak (UniMZ. mz_to_mh (p. mz, z), p. inten), zs)
7883 end
79- zs = filter (z -> UniMZ. query_ε (ms. peaks, p. mz + Δ / z, εm) |> ! isempty, 1 : 3 )
80- zs = isempty (zs) ? [1 ] : zs
81- return map (z -> UniMZ. Peak (UniMZ. mz_to_mh (p. mz, z), p. inten), zs)
84+ return UniMZ. fork (ms; peaks= reduce (vcat, ps)|> sort)
8285 end
83- return UniMZ. fork (ms; peaks= reduce (vcat, ps)|> sort)
8486 end
8587 @info " MS2 preprocessing using xic pattern..."
8688 S = @showprogress map (eachindex (M2), M2, I) do idx_ms, ms, ions
0 commit comments