Use calojet methods from event.
Cut on 4th calojet pt instead of pf jet.
Turn on cleaning filters.
geo2/3ddist methods.
min/max_geo2ddist.
Get rid of mfvSelectedVerticesMedium (didn't I already?).
Fiducial cut: max_geo2ddist = 2.5.
One place to get analysis-wide constants from in C++ or python.
Use the AnalysisConstants header/module instead of hardcoding e.g. int_lumi everywhere.