library("tima")
# Calculate m/z for a protonated molecule
calculate_mz_from_mass(neutral_mass = 122.45, adduct_string = "[M+H]+")[1] 123.4578
# Expected: ~123.4567
# Verify round-trip calculation
mass <- 122.45
adduct <- "[M+H]+"
mz <- calculate_mz_from_mass(mass, adduct)
mass_back <- calculate_mass_of_m(mz, adduct)
all.equal(mass, mass_back) # Should be TRUE[1] TRUE