DJR plots & Matching/Merging

Hello everyone,

I am trying to understand how to correctly set xqcut and qcut, so I am following the W+jet production example in this tutorial [https://indico.cern.ch/event/555228/sessions/203428/attachments/1315471/1970459/tutorial-CMSandATLAS-2016.pdf] by Marco Zaro. There, around slide 146, the discussion about matching & merging begins for the example mentioned above.

By running the following MG5 lines,
===========
import model sm-ckm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
generate p p > w+ @0
add process p p > w+ j @1
add process p p > w+ j j @2
output WPROD
launch WPROD
shower = Pythia8
detector = OFF
analysis = OFF
set use_syst False
set iseed 1234
set ickkw 1
set pjt 0
set xqcut 30
done

launch WPROD
set ickkw 1
set pjt 0
set xqcut 10
done

launch WPROD
set ickkw 1
set pjt 0
set xqcut 60
done

launch WPROD
set ickkw 1
set pjt 0
set xqcut 100
done

launch WPROD
set ickkw 1
set pjt 0
set xqcut 30
set pythia8_card JetMatching:qCut 42
done

launch WPROD
set ickkw 1
set pjt 0
set xqcut 30
set pythia8_card JetMatching:qCut 30
done

launch WPROD
set ickkw 1
set pjt 0
set xqcut 30
set pythia8_card JetMatching:qCut 55
done
===========

I expected to compare my DJR plots with those on the slides. Surprisingly, mine don't exhibit the same features (discontinuities and kinks) as those on M. Zaro's slides, so I don't know what I am doing wrong and how to correct it. I have put some of them here [https://www.dropbox.com/t/p5kliRXe4sF47zpg].

My ultimate goal is to use these types of plots for a multijet background study, and I've learned that the proper choice of xqcut & qcut can be inferred from the smoothness of the DJR plots.

Many thanks, Sebastián