"質量分析の方法であって、飛行時間領域内に入るようにイオンを加速する電極と、前記飛行時間領域を前記イオンが通過した後にイオンを検出するように構成されるイオン検出器と、を備える飛行時間型質量分析器と準備する工程と、前記イオン検出器から出力される第1の信号をデジタル化して、第1のデジタル化信号を生成する工程と、前記第1のデジタル化信号のデコンボリューションを行なって、前記第1のデジタル化信号に関して1つ以上の第1のイオン到達時間と1つ以上の第1のイオン到達強度とを求める工程と、前記イオン検出器から出力される第2の信号をデジタル化して、第2のデジタル化信号を生成する工程と、前記第2のデジタル化信号のデコンボリューションを行なって、前記第2のデジタル化信号に関して1つ以上の第2のイオン到達時間と1つ以上の第2のイオン到達強度とを求める工程と、前記イオン検出器から出力される第3及びさらなる信号をデジタル化して、第3及びさらなるデジタル化信号を生成する工程と、前記第3及びさらなるデジタル化信号のデコンボリューションを行なって、前記第3及びさらなるデジタル化信号に関して1つ以上の第3及びさらなるイオン到達時間と1つ以上の第3及びさらなるイオン到達強度とを求める工程と、前記1つ以上の第1のイオン到達時間と前記1つ以上の第2のイオン到達時間と前記1つ以上の第3及びさらなるイオン到達時間とを組み合わせると共に、前記1つ以上の第1のイオン到達強度と前記1つ以上の第2のイオン到達強度と前記1つ以上の第3及びさらなるイオン到達強度とを組み合わせて、複合イオン到達時間-強度スペクトルを生成する工程と、を備える質量分析の方法。請求項1に記載の質量分析の方法であって、前記イオン検出器から出力される前記第1の信号をデジタル化する工程と、前記イオン検出器から出力される前記第2の信号をデジタル化する工程と、前記イオン検出器から出力される前記第3及びさらなる信号をデジタル化する工程とが、アナログ-デジタル変換器を用いて、前記第1の信号と、前記第2の信号と、前記第3及びさらなる信号とをデジタル化する工程を備える、質量分析の方法。請求項1又は2のいずれか一項に記載の質量分析の方法であって、前記第1のデジタル化信号のデコンボリューションを行なう工程と、前記第2のデジタル化信号のデコンボリューションを行なう工程と、前記第3及びさらなるデジタル化信号のデコンボリューションを行なう工程とが、i前記イオン検出器に到達して前記イオン検出器により検出されるイオンの特徴を示す点広がり関数を求める工程、又は、ii前記イオン検出器に到達して前記イオン検出器により検出されるイオンの特徴を示す所定の点広がり関数を用いる工程、のいずれかの工程を備える、質量分析の方法。請求項1から3までのいずれか一項に記載の質量分析の方法であって、i前記第1のデジタル化信号のデコンボリューションを行なう工程が、前記イオン検出器に到達して前記イオン検出器により検出されるイオンの特徴を示す点広がり関数の逆数を用いて、前記第1のデジタル化信号のコンボリューションを行なう工程を備える、及び、ii前記第2のデジタル化信号のデコンボリューションを行なう工程が、前記イオン検出器に到達して前記イオン検出器により検出されるイオンの特徴を示す点広がり関数の逆数を用いて、前記第2のデジタル化信号のコンボリューションを行なう工程を備える、及び、iii前記第3及びさらなるデジタル化信号のデコンボリューションを行なう工程が、前記イオン検出器に到達して前記イオン検出器により検出されるイオンの特徴を示す点広がり関数の逆数を用いて、前記第3及びさらなるデジタル化信号のコンボリューションを行なう工程を備える、質量分析の方法。請求項1から4までのいずれか一項に記載の質量分析の方法であって、i前記第1のデジタル化信号のデコンボリューションを行なう工程が、各イオン到達が既知の点広がり関数によって表わされる応答を生じる場合に前記第1のデジタル化信号に最も良く適合するイオン到達時間の分布を求める工程を備える、及び、ii前記第2のデジタル化信号のデコンボリューションを行なう工程が、各イオン到達が既知の点広がり関数によって表わされる応答を生じる場合に前記第2のデジタル化信号に最も良く適合するイオン到達時間の分布を求める工程を備える、及び、iii前記第3及びさらなるデジタル化信号のデコンボリューションを行なう工程が、各イオン到達が既知の点広がり関数によって表わされる応答を生じる場合に前記第3及びさらなるデジタル化信号に最も良く適合するイオン到達時間の分布を求める工程を備える、質量分析の方法。請求項1から5までのいずれか一項に記載の質量分析の方法であって、前記第1のデジタル化信号と、前記第2のデジタル化信号と、前記第3及びさらなるデジタル化信号とに関して前記イオン到達時間と前記イオン到達強度とを求める工程が、高速デコンボリューションアルゴリズムを用いる工程を備える、質量分析の方法。請求項6に記載の質量分析の方法であって、前記高速デコンボリューションアルゴリズムが、i修正CLEANアルゴリズムとii最大エントロピー法とiii高速フーリエ変換とiv非負の最小二乗法とからなる群から選択される、質量分析の方法。請求項6又は7のいずれか一項に記載の質量分析の方法であって、前記高速デコンボリューションアルゴリズムが、既知の線幅と、前記イオン検出器により生成された後に個々のイオン到達に応じてデジタル化される信号の特徴を示す形状とを用いる、質量分析の方法。請求項1から8までのいずれか一項に記載の質量分析の方法であって、さらに、得られたイオンの到達時間T0を第1の到達時間Tnと第2の到達時間Tn1とに変換する工程であって、nがT0に最も近いデジタル化時間ビンを表わす工程と、得られたイオンの強度S0を第1の強度Snと第2の強度Sn1とで表わす工程と、を備え、である、質量分析の方法。請求項1から9までのいずれか一項に記載の質量分析の方法であって、前記第1のデジタル化信号と前記第2のデジタル化信号と前記第3及びさらなるデジタル化信号とのデコンボリューションを行なう工程が、前記第1のデジタル化信号と前記第2のデジタル化信号と前記第3及びさらなるデジタル化信号との後処理により行なわれる、質量分析の方法。請求項1から9までのいずれかに記載の質量分析の方法であって、前記第1のデジタル化信号と前記第2のデジタル化信号と前記第3及びさらなるデジタル化信号とのデコンボリューションを行なう工程が、フィールドプログラマブルゲートアレイFPGA:FieldProgrammableGateArray又はグラフィカルプロセッサユニットGPA:GraphicalProcessorUnitを用いてリアルタイムで行なわれる、質量分析の方法。請求項1から11までのいずれか一項に記載の質量分析の方法であって、さらに、i前記第1の信号をデジタル化する工程及び/又は前記第1のデジタル化信号のデコンボリューションを行なう工程の前に、前記飛行時間領域内に入るように第1のイオン群を加速する工程、及び/又は、ii前記第2の信号をデジタル化する工程及び/又は前記第2のデジタル化信号のデコンボリューションを行なう工程の前に、前記飛行時間領域内に入るように第2のイオン群を加速する工程、及び/又は、iii前記第3の信号をデジタル化する工程及び/又は前記第3のデジタル化信号のデコンボリューションを行なう工程の前に、前記飛行時間領域内に入るように第3のイオン群を加速する工程、を備える質量分析の方法。質量分析計であって、飛行時間領域内に入るようにイオンを加速する電極と、前記飛行時間領域を前記イオンが通過した後にイオンを検出するように構成されるイオン検出器と、を備える飛行時間型質量分析器と、制御システムであって、i前記イオン検出器から出力される第1の信号をデジタル化して、第1のデジタル化信号を生成し、ii前記第1のデジタル化信号のデコンボリューションを行なって、前記第1のデジタル化信号に関して1つ以上の第1のイオン到達時間と1つ以上の第1のイオン到達強度とを求め、iii前記イオン検出器から出力される第2の信号をデジタル化して、第2のデジタル化信号を生成し、iv前記第2のデジタル化信号のデコンボリューションを行なって、前記第2のデジタル化信号に関して1つ以上の第2のイオン到達時間と1つ以上の第2のイオン到達強度とを求め、v前記イオン検出器から出力される第3及びさらなる信号をデジタル化して、第3及びさらなるデジタル化信号を生成し、vi前記第3及びさらなるデジタル化信号のデコンボリューションを行なって、前記第3及びさらなるデジタル化信号に関して1つ以上の第3及びさらなるイオン到達時間と1つ以上の第3及びさらなるイオン到達強度とを求め、vii前記1つ以上の第1のイオン到達時間と前記1つ以上の第2のイオン到達時間と前記1つ以上の第3及びさらなるイオン到達時間とを組み合わせると共に、前記1つ以上の第1のイオン到達強度と前記1つ以上の第2のイオン到達強度と前記1つ以上の第3及びさらなるイオン到達強度とを組み合わせて、複合イオン到達時間-強度スペクトルを生成する、ように配置及び構成される制御システムと、を備える質量分析計。請求項13に記載の質量分析計であって、前記制御システムが、i前記第1の信号をデジタル化する前及び/又は前記第1のデジタル化信号のデコンボリューションを行なう前に、前記飛行時間領域内に入るように第1のイオン群を加速する、及び/又は、ii前記第2の信号をデジタル化する前及び/又は前記第2のデジタル化信号のデコンボリューションを行なう前に、前記飛行時間領域内に入るように第2のイオン群を加速する、及び/又は、iii前記第3の信号をデジタル化する前及び/又は前記第3のデジタル化信号のデコンボリューションを行なう前に、前記飛行時間領域内に入るように第3のイオン群を加速する、ように配置及び構成される、質量分析計。質量分析の方法であって、飛行時間領域内に入るようにイオンを加速する電極と、前記飛行時間領域を前記イオンが通過した後にイオンを検出するように構成されるイオン検出器と、を備える飛行時間型質量分析器と準備する工程と、アナログ-デジタル変換器を用いて前記イオン検出器から出力される第1の信号をデジタル化して、第1のデジタル化信号を生成する工程と、前記第1のデジタル化信号のデコンボリューションを行なって、前記第1のデジタル化信号に関して1つ以上の第1のイオン到達時間と1つ以上の第1のイオン到達強度とを求める工程であって、前記第1のデジタル化信号のデコンボリューションを行なう工程が、各イオン到達が既知の点広がり関数により表わされる応答を生じる場合に前記第1のデジタル化信号に最も良く適合するイオン到達時間の分布を求める工程を備える、工程と、アナログ-デジタル変換器を用いて前記イオン検出器から出力される第2の信号をデジタル化して、第2のデジタル化信号を生成する工程と、前記第2のデジタル化信号のデコンボリューションを行なって、前記第2のデジタル化信号に関して1つ以上の第2のイオン到達時間と1つ以上の第2のイオン到達強度とを求める工程であって、前記第2のデジタル化信号のデコンボリューションを行なう工程が、各イオン到達が既知の点広がり関数により表わされる応答を生じる場合に前記第2のデジタル化信号に最も良く適合するイオン到達時間の分布を求める工程を備える、工程と、アナログ-デジタル変換器を用いて前記イオン検出器から出力される第3及びさらなる信号をデジタル化して、第3及びさらなるデジタル化信号を生成する工程と、前記第3及びさらなるデジタル化信号のデコンボリューションを行なって、前記第3及びさらなるデジタル化信号に関して1つ以上の第3及びさらなるイオン到達時間と1つ以上の第3及びさらなるイオン到達強度とを求める工程であって、前記第3及びさらなるデジタル化信号のデコンボリューションを行なう工程が、各イオン到達が既知の点広がり関数により表わされる応答を生じる場合に前記第3及びさらなるデジタル化信号に最も良く適合するイオン到達時間の分布を求める工程を備える、工程と、前記1つ以上の第1のイオン到達時間と前記1つ以上の第2のイオン到達時間と前記1つ以上の第3及びさらなるイオン到達時間とを組み合わせると共に、前記1つ以上の第1のイオン到達強度と前記1つ以上の第2のイオン到達強度と前記1つ以上の第3及びさらなるイオン到達強度とを組み合わせて、複合イオン到達時間-強度スペクトルを生成する工程と、を備える質量分析の方法。請求項15に記載の質量分析の方法であって、さらに、i前記第1の信号をデジタル化する工程及び/又は前記第1のデジタル化信号のデコンボリューションを行なう工程の前に、前記飛行時間領域内に入るように第1のイオン群を加速する工程、及び/又は、ii前記第2の信号をデジタル化する工程及び/又は前記第2のデジタル化信号のデコンボリューションを行なう工程の前に、前記飛行時間領域内に入るように第2のイオン群を加速する工程、及び/又は、iii前記第3の信号をデジタル化する工程及び/又は前記第3のデジタル化信号のデコンボリューションを行なう工程の前に、前記飛行時間領域内に入るように第3のイオン群を加速する工程、を備える質量分析の方法。質量分析計であって、飛行時間領域内に入るようにイオンを加速する電極と、前記飛行時間領域を前記イオンが通過した後にイオンを検出するように構成されるイオン検出器と、を備える飛行時間型質量分析器と、制御システムであって、iアナログ-デジタル変換器を用いて前記イオン検出器から出力される第1の信号をデジタル化して、第1のデジタル化信号を生成し、ii前記第1のデジタル化信号のデコンボリューションを行なって、前記第1のデジタル化信号に関して1つ以上の第1のイオン到達時間と1つ以上の第1のイオン到達強度とを求め、ここで、前記制御システムが、各イオン到達が既知の点広がり関数により表わされる応答を生じる場合に前記第1のデジタル化信号に最も良く適合するイオン到達時間の分布を求めるように配置及び構成され、iiiアナログ-デジタル変換器を用いて前記イオン検出器から出力される第2の信号をデジタル化して、第2のデジタル化信号を生成し、iv前記第2のデジタル化信号のデコンボリューションを行なって、前記第2のデジタル化信号に関して1つ以上の第2のイオン到達時間と1つ以上の第2のイオン到達強度とを求め、ここで、前記制御システムが、各イオン到達が既知の点広がり関数により表わされる応答を生じる場合に前記第2のデジタル化信号に最も良く適合するイオン到達時間の分布を求めるように配置及び構成され、vアナログ-デジタル変換器を用いて前記イオン検出器から出力される第3及びさらなる信号をデジタル化して、第3及びさらなるデジタル化信号を生成し、vi前記第3及びさらなるデジタル化信号のデコンボリューションを行なって、前記第3及びさらなるデジタル化信号に関して1つ以上の第3及びさらなるイオン到達時間と1つ以上の第3及びさらなるイオン到達強度とを求め、ここで、前記制御システムが、各イオン到達が既知の点広がり関数により表わされる応答を生じる場合に前記第3及びさらなるデジタル化信号に最も良く適合するイオン到達時間の分布を求めるように配置及び構成され、vii前記1つ以上の第1のイオン到達時間と前記1つ以上の第2のイオン到達時間と前記1つ以上の第3及びさらなるイオン到達時間とを組み合わせると共に、前記1つ以上の第1のイオン到達強度と前記1つ以上の第2のイオン到達強度と前記1つ以上の第3及びさらなるイオン到達強度とを組み合わせて、複合イオン到達時間-強度スペクトルを生成する、ように配置及び構成される制御システムと、を備える質量分析計。請求項17に記載の質量分析計であって、前記制御システムが、さらに、i前記第1の信号をデジタル化する前及び/又は前記第1のデジタル化信号のデコンボリューションを行なう前に、前記飛行時間領域内に入るように第1のイオン群を加速する、及び/又は、ii前記第2の信号をデジタル化する前及び/又は前記第2のデジタル化信号のデコンボリューションを行なう前に、前記飛行時間領域内に入るように第2のイオン群を加速する、及び/又は、iii前記第3の信号をデジタル化する前及び/又は前記第3のデジタル化信号のデコンボリューションを行なう前に、前記飛行時間領域内に入るように第3のイオン群を加速する、ように配置及び構成される、質量分析計。質量分析の方法であって、飛行時間領域内に入るようにイオンを加速する電極と、前記飛行時間領域を前記イオンが通過した後にイオンを検出するように構成されるイオン検出器と、を備える飛行時間型質量分析器と準備する工程と、i前記飛行時間領域内に入るようにイオン群を加速する工程と、ii前記イオン検出器から出力される信号をデジタル化して、デジタル化信号を生成する工程と、を1回以上繰り返す工程と、前記デジタル化信号を組み合わせて、第1の複合デジタル化信号を形成する工程と、前記第1の複合デジタル化信号のデコンボリューションを行なって、前記第1の複合デジタル化信号に関して1つ以上の第1のイオン到達時間と1つ以上の第1のイオン到達強度とを求める工程と、iii前記飛行時間領域内に入るようにイオン群を加速する工程と、iv前記イオン検出器から出力される信号をデジタル化して、デジタル化信号を生成する工程と、を1回以上繰り返す工程と、前記デジタル化信号を組み合わせて、第2の複合デジタル化信号を形成する工程と、前記第2の複合デジタル化信号のデコンボリューションを行なって、前記第2の複合デジタル化信号に関して1つ以上の第2のイオン到達時間と1つ以上の第2のイオン到達強度とを求める工程と、前記1つ以上の第1のイオン到達時間と前記1つ以上の第2のイオン到達時間とを組み合わせると共に、前記1つ以上の第1のイオン到達強度と前記1つ以上の第2のイオン到達強度とを組み合わせて、複合イオン到達時間-強度スペクトルを生成する工程と、を備える質量分析の方法。質量分析計であって、飛行時間領域内に入るようにイオンを加速する電極と、前記飛行時間領域を前記イオンが通過した後にイオンを検出するように構成されるイオン検出器と、を備える飛行時間型質量分析器と、制御システムであって、i前記飛行時間領域内に入るようにイオン群を加速する工程と、ii前記イオン検出器から出力される信号をデジタル化して、デジタル化信号を生成する工程と、を1回以上繰り返し、前記デジタル化信号を組み合わせて、第1の複合デジタル化信号を形成し、前記第1の複合デジタル化信号のデコンボリューションを行なって、前記第1の複合デジタル化信号に関して1つ以上の第1のイオン到達時間と1つ以上の第1のイオン到達強度とを求め、iii前記飛行時間領域内に入るようにイオン群を加速する工程と、iv前記イオン検出器から出力される信号をデジタル化して、デジタル化信号を生成する工程と、を1回以上繰り返し、前記デジタル化信号を組み合わせて、第2の複合デジタル化信号を形成し、前記第2の複合デジタル化信号のデコンボリューションを行なって、前記第2の複合デジタル化信号に関して1つ以上の第2のイオン到達時間と1つ以上の第2のイオン到達強度とを求め、前記1つ以上の第1のイオン到達時間と前記1つ以上の第2のイオン到達時間とを組み合わせると共に、前記1つ以上の第1のイオン到達強度と前記1つ以上の第2のイオン到達強度とを組み合わせて、複合イオン到達時間-強度スペクトルを生成する、ように構成及び配置される制御システムと、を備える質量分析計。本発明の好適な実施形態を説明する。好適な実施例では、イオン検出器を備える飛行時間型質量分析器を準備する。各飛行時間分析においてイオン検出器からの出力をアナログ-デジタル変換器ADC:AnalogtoDigitalConverterデジタル化することが望ましい。好適な実施形態に置いて、各飛行時間スペクトルにデコンボリューションアルゴリズムを適用する。このデコンボリューションアルゴリズムは、整数計算のみを行なうように構成されている。デコンボリューションの方法は、以下で詳細に説明するように、この環境で重なっている信号源を取り扱うようにさらに拡張されたものでもよい。一つの実施形態において、デューティ比に損失を与えることなく個々の飛行時間スペクトル上でデコンボリューションを実行可能な高速フィールドプログラマブルゲートアレイFPGA:FieldProgrammableGateArrayのアーキテクチャを用いるようにしてもよい。アナログ-デジタル変換器ADC:AnalogtoDigitalConverterにより生成されるデジタル化信号の解析には、好適な実施形態において用いられる整数計算が特に適している。好適な実施形態の態様を説明するために、コンボリューション畳み込みによりリアルマップfxをデータ領域gxに変換する空間不変量点広がり関数PSF:pointspreadfunctionを仮定する。点広がり関数は、イオン検出器の応答を平均強度の単一イオン到達に変換する理想化プロファイルを表わす。リアルマップfxは個々のイオンの実際の到達時間を示し、データ領域gxは最終的に記録される飛行時間スペクトルを示す。イオン検出器から出力されるアナログ信号をデジタル化する際に、有限グリッド上で測定結果が得られると考えることができる。このグリッドの粗さはアナログ-デジタル変換器のデジタル化速度に依存する。信号もノイズの影響を受ける。上述の式3で示す変換を反転させる代わりに、好適な実施形態ではfを推定する。説明を容易にするために、リアルマップfx及びデータ領域gxが同じグリッド上でデジタル化されると仮定する。記録データgiは、測定値yiにノイズによる誤差を含むものである。ノイズは、独立ガウス分布であると仮定すると、一様に、以下の分散式で表わすことができる。あるいは、以下の行列-ベクトル式で表わすことができる。fに関する法線を求めることにより式6を最小化することができる。始点fから徐々に値を増加させて、χ2等を減少させるインクリメント増加量Δfを選択することにより、これを実施するようにしてもよい。ベクトルfは、イオン到達時間をデジタル化計算を行なう。点広がり関数は平均高さのイオン検出器からの電圧パルスであり、yは同じグリッド上でデジタル化された検出器電圧波形測定値である。図1にデジタル化点広がり関数pxの例を示す。この関数は、値2、6、11、14、15、14、11、6、2を有し、整数計算で閾値t=2262112142152/2=469を与える。マップを以下のようにインクリメントすることにより生成されるχ2のデクリメント減少量は以下のように表わされる。fの自然インクリメント増加は、所定の時間指数jにおいて単一のイオン到達を加算するものである。この結果、以下の式が設定される。すなわち、j番目の要素のみがゼロではない。単一指数jを選択して、イオン計数を1インクリメントする増加することにより以下の式が得られる。ここで、ここで、rは、不明確な残差ベクトルを示す。Δχ2に関する式の第1項は、不明確なデータと二重に不明確なマップとの間の差が最大になるような時間指数、すなわち、rの最大値における時間指数を選択することにより、χ2の最大デクリメント減少量が得られることを示す。すなわち、不明確なデータと二重に不明確なマップとの間の差が自己コンボリューションにより点広がり関数のピーク値の半分未満になるとインクリメント増加を停止する自然停止基準が提示される。実際には、許容閾値よりも大きな、単一繰り返しにおける不明確な残差ベクトルrのすべての極大値において、イオン計数をインクリメント増加させることができる。好適な実施形態において、用いられる修正CLEANアルゴリズムは、以下の工程を含むものでもよい。1.いずれの場所においてもfをゼロに初期化する。2.n=1~N繰り返しに関して、不明確な残差を計算する。3.よりも大きなrの各極大値rjに関して、以下の式を設定する。上述の処理は、ある程度分離された多数の点源の位置及び強度を見つけるのに特に適している。閾値を調節することにより、非ゼロのバックグラウンドレベルに対応可能である。ここで、bはバックグラウンドレベルを示す。ただし、各プッシュ毎の質量スペクトルピークにおける何十ものイオンのイオン到達速度に関して、上述の処理の効果が充分に得られるのに十分な程度にイオン到達を分離できない場合もある。電圧パルスが重なっている場合には、生成される極大値がイオン到達時間に対応しない可能性があるという問題が生じる。このような場合、選択される最初の極大値は、マップをインクリメント後に測定される後続の極大値よりも誤差が大きい可能性がある。特に好適な実施形態において、以下の説明でCLEANER処理と称する工程を含むように、上述した修正CLEAN処理をさらに修正するようにしてもよい。CLEANER処理は、以下の工程を含むものでもよい。1.いずれの場所においてもfをゼロに初期化する。2.n=1~N繰り返しに関して、不明確な残差を計算する。3.で与えられるエロージョンの確率を有する各に関して、以下の式を設定する。4.よりも大きなrの各極大値rjに関して、以下の式を設定する。繰り返し数nが増加するにつれて、エロージョンの確率qnは直線的に低下する。異なるプッシュに関して取得されたデータに対応する多数のデータセットが利用可能であれば、Hogbom1974年に記載されるループゲインγの漸増として、エロージョンの確率qの低下を観測することができる。実際には、真のイオン到達位置に関する不確実性が高い場合には、低いγ値が用いられる。好適な実施形態の様々な態様を説明するために、エレクトロスプレー型イオン源を介して、直交加速飛行時間型質量分析計に牛インスリンのサンプルを注入した。M5H5イオンにより生成されるイオン信号が、イオン検出器に入力され、3GHzのデジタル化速度で、8ビットのアナログ-デジタル変換器を用いて記録された。926個の飛行時間スペクトルが記録され、上述した好適なCLEANER処理を128回繰り返すことにより、各飛行時間スペクトルのデコンボリューションを行なった。各飛行時間スペクトルに関して求めたイオン到達位置を最終的なスペクトルにまとめた。図2に単一の飛行時間スペクトルを示す。このスペクトルでは、2つの単一イオン到達がはっきりと分かる。イオンは、牛インスリンのM5H5イオンの同位体クラスタ由来のものである。図2に示す飛行時間スペクトルの検討及び個々のイオン到達を含む他のスペクトルの検討から、イオン到達の特徴的な形状を表わす点広がり関数を誘導するようにしてもよい。この例における点広がり関数は、図3に示すように、強度値1、2、5、17、23、16、6、2、2、4、3、2、1を有する。この例では、単一イオンプロファイルは非対称であり、立下りエッジの後で有意のサテライト又はリンギングピークを有する。サテライトは、検出器エレクトロニクスにおけるインピーダンスの不一致により生じ、非常に高速の単一イオン応答では大なり小なり見られる現象である。図4Aに、同じデータセットから得られた飛行時間スペクトルNo.449を示す。この場合、イオン検出器に複数のイオンが到達した。図4Aに示す飛行時間スペクトルにおいて、ピーク1は、単一イオン到達から予測される単一応答よりも高さ及び幅が大きい。したがって、このピークは、狭い時間ウィンドウの間に到達した複数のイオン信号が重なって生成されている可能性が高い。図4Bに,好適な実施形態に従って計算した場合の到達時間位置を示す。図4Bに示すように、ピーク1は、各々が図3に示す点広がり関数を有する複数のイオン到達値に割り当てられる。これに対して、当業者には周知のように、有限インパルス応答フィルタに基づく処理等のピーク検出処理では、信号の質量中心又は頂点に対応する信号の単一飛行時間値のみを検出することができる。図4Aのピーク1により表わされる単一のイオンピークを4つのピークに分解することにより、短時間の間に7個のイオン到達が生じていることが示される。このことから、周知の方法と比較して、本発明の好適な実施形態の効果が明らかである。図5Aに、926個の飛行時間スペクトルをまとめて、閾値バックグラウンド減算法を適用して得られた飛行時間スペクトルを示す。牛インスリンの5イオンの同位体エンベロープが明確に観察される。ただし、図2に示すような各単一イオン到達に伴う非対称性により、最終的なスペクトルの各同位体ピークにおいても対応する非対称性が明確に観察される。図5Bに、好適な実施形態に従う処理後の図5Aのデータを示す。図5Aと比較して、ピークの対称性が大幅に向上していることが明らかである。これは、より好ましいピーク形状、さらに、より明確な谷の分離につながる。デコンボリューション処理に用いられる点広がり関数を検出システムの特徴的なイオンプロファイルに一致させることにより、最終的なデータにおけるアーチファクト不自然な結果やテーリングを抑制することが可能になる。これらの質的向上に加えて、質量分解能も増大する。これは、好適な実施形態において、図5Aに明示されるイオン到達プロファイルからピーク幅に対する影響を効果的に取り除くことができるためである。この例では、比較データを与えるために取得したデータの後処理を行なったが、好適な実施形態において、フィールドプログラマブルゲートアレイFPGA:FieldProgrammableGateArray又はグラフィカルプロセッサユニットGPA:GraphicalProcessorUnitアーキテクチャを用いてリアルタイムで処理を行なうことが望ましい。上述した方法では、デジタル化ビン幅の側及び-側に半分の精度でイオン到達時間を求めることが望ましい。ただし、他の実施形態において、入力信号のデジタル化精度の半分未満の精度でイオン到達時間を求めるようにしてもよい。この場合には、データと比べて点広がり関数のサンプリングを増大させればよい、及び/又は、デコンボリューションの前に補間によりデータのサンプリングを増大させればよい。あるいは、1つのデジタル化ビンにおいて許容閾値を超える不明確な残差における応答の極大値を記録する代わりに、不明確な残差の頂点を補間することによって、又は、信号の重み付け質量中心を算出することによって、より精度よく極大値を記録するようにしてもよい。イオン到達時間が高精度で得られる場合、デコンボリューション後の個々の飛行時間スペクトルを組み合わせる際に、元のデジタル化データのグリッド間隔よりも細かいグリッド間隔を用いるようにしてもよい。これにより、元のデータと比べて、最終的に得られる質量スペクトルのデジタル化速度が明らかに高くなる。さらに、イオン到達時間が高精度で得られる場合、得られたイオン到達時間T0を第1の到達時間Tnと第2の到達時間Tn1とに変換することにより、また、イオンの測定強度S0を、第1の強度Snと第2の強度Sn1とで表わすことにより、この精度を最終的なデータでも維持することができる。ここで、nはT0に最も近いデジタル化時間ビンである。以上、本発明をその好適な実施形態を参照して詳述したが、当業者には自明のことであるが、特許請求の範囲に記載される本発明の要旨を逸脱しない範囲において、形態や詳細において、種々の変形や変更が可能である。"