Update 06.05.2018 -- Run 9 embedding: new matching algorithm

We recently discovered that I wasn't handling the detector-to-particle matching in the best manner.  What tipped us off was the (very small) population of jets at large 'pTdetector' but small 'pTparticle' in the response matrix:

This was a consequence of how I was handling the selection of the "best" match candidate.  The matching algorithm I was using went like so:

  1. Select a particle-level jet.
  2. For each detector-level jet compute the separation 'dR = sqrt((phiPar - phiDet)^2 + (etaPar - etaDet)^2)' between the particle- and detector-level axes, and the momentum fraction 'qT = pTdet / pTpar'
  3. If a detector-level jet satisfies 'dR < 0.3' (where 0.3 is the jet resolution parameter) and 'qT > 0.15', then consider it a match candidate.
  4. From the set of match candidates, the detector-level jet which has the smallest 'dR' and largest 'qT' is defined to be the best match candidate.
  5. Remove the particle-level jet and its best match from their respective lists, and iterate until there are no more particle-level jets.

This also explains why my jet reconstruction efficiency was a little lower than anticipated (~80%).  So we modified the algorithm to go like so:

  1. Select a particle-level jet.
  2. For each detector-level jet compute the separation 'dR = sqrt((phiPar - phiDet)^2 + (etaPar - etaDet)^2)' between the particle- and detector-level axes, and the momentum fraction 'qT = pTdet / pTpar'
  3. If a detector-level jet satisfies 'dR < 0.3' (where 0.3 is the jet resolution parameter) and 'qT > 0.15', then consider it a match candidate.
  4. From the set of match candidates, the detector-level jet which has the 'qT' closest to one is defined to be the best match candidate.
  5. Remove the particle-level jet and its best match from their respective lists, and iterate until there are no more particle-level jets.

Step 4 is implemented by requiring that the best match has the smallest 'abs(qT - 1)'.  This had a substantial impact on our reconstruction efficiency and response matrix:


The above response matrix (and efficiency) is the combined RFF and FF response matrices (and efficiencies).  Here's how the RFF and FF responses compare to each other:

And here's the new efficiency and response matrix with variable binning: