Match reconstructed tracks and MC truth

The relation between reconstructed and MC particles is not necessarily one-to-one. Therefore, a loop has to be constructed:

        ITK1 = ... (any given MC or reconstructed track number)
        DO 10 I = 1, KNMTCH(ITK1)
          IF (KSMTCH (ITK1,I) .LE. (min. required shared hits))  GO TO 10
          ITK2 = KMTCH (ITK1,I)
  C       ...
     10 CONTINUE

If ITK1 is a reconstructed track then ITK2 is a matching MC track. If ITK1 is a MC track then ITK2 is a matching reconstructed track.

KNMTCH (ITK)

Number of matching candidates for track ITK.

KMTCH (ITK,I)

Track number of Ith matching particle.

KSMTCH (ITK,I)

Number of shared hits between MC and reconstructed track.

KBESTM (ITK)

Best match to track ITK.

Remarks:

• The match is performed on the basis of shared hits in the TPC and ITC.
• A much better matching is provided for charged tracks since ALPHA 122 by a set of new routines: MCMATCH, JULMATCH, VDHMATCH which are described in Section 12.13 on p. .
• A specific matching is available for Energy Flow objects (charged and neutral); it is described in Section 11.1.5 on p. .