meName = 'DiTauCandidateQuantities/PzetaDiff',
        title = "P_{#zeta} - 1.5*P_{#zeta}^{vis} (final Event sample)",
        xAxis = 'GeV',
        name = "finalSamplePlots_#PROCESSNAME#_PzetaDiff"
    ),
]

# Make plots for each analysis chain
for analyzer, cut in analyzers:
    # The basic draw job configurator
    drawFakeRateHist_background = drawFakeRateHistConfigurator(
        template = cms.PSet(
            xAxis = cms.string('unlabeled'),
            #yAxis = cms.string('numEntries_linear'),
            yAxis = cms.string('numEntries_log'),
            legend = cms.string('regular'),
            labels = cms.vstring('mcNormScale')
        ),
        dqmDirectories = dqmDirectories[analyzer],
        legendEntries = legendEntries,
        frTypes = fr_types
    )
    # Add all our of our fake background sources
    for source in fake_sources:
        drawFakeRateHist_background.addProcess(
            source, cms.string('/harvested/%s' % source))

    drawFakeRateHist_background.addPlots(
        afterCut = cut,
        plots = bgEstFakeRatePlots
    )
                       name="bgEstFakeRatePlots_#PROCESSNAME#_mtMuonMET"),
    drawJobConfigEntry(meName='DiTauCandidateQuantities/VisMass',
                       title="M_{vis}(Muon + Tau) (final Event sample)",
                       xAxis='Mass',
                       name="bgEstFakeRatePlots_#PROCESSNAME#_mVisible")
]

drawFakeRateHistConfiguratorZtoMuTau_background = drawFakeRateHistConfigurator(
    template=cms.PSet(
        xAxis=cms.string('unlabeled'),
        #yAxis = cms.string('numEntries_linear'),
        yAxis=cms.string('numEntries_log'),
        legend=cms.string('regular'),
        labels=cms.vstring('mcNormScale')),
    dqmDirectories=dqmDirectories,
    legendEntries=legendEntries,
    frTypes=[
        "tauIdDiscr",
        "frQCDmuEnriched",
        "frQCDdiJetLeadJet",
        "frQCDdiJetSecondLeadJet",
        "frWplusJets",
        ##"frGammaPlusJets",
        "frSysUncertainty"
    ])

drawFakeRateHistConfiguratorZtoMuTau_background.addProcess(
    "Zmumu", processZtoMuTau_Zmumu.config_dqmHistPlotter.dqmDirectory)
drawFakeRateHistConfiguratorZtoMuTau_background.addProcess(
    "WplusJets", processZtoMuTau_WplusJets.config_dqmHistPlotter.dqmDirectory)
drawFakeRateHistConfiguratorZtoMuTau_background.addProcess(
    "TTplusJets",
        title = "M_{vis}(Muon + Tau) (final Event sample)",
        xAxis = 'Mass',
        name = "bgEstFakeRatePlots_#PROCESSNAME#_mVisible"
    )
]

drawFakeRateHistConfiguratorZtoMuTau_background = drawFakeRateHistConfigurator(
    template = cms.PSet(
        xAxis = cms.string('unlabeled'),
        #yAxis = cms.string('numEntries_linear'),
        yAxis = cms.string('numEntries_log'),
        legend = cms.string('regular'),
        labels = cms.vstring('mcNormScale')
    ),
    dqmDirectories = dqmDirectories,
    legendEntries = legendEntries,
    frTypes = [
        "tauIdDiscr",
        "frQCDmuEnriched",
        "frQCDdiJetLeadJet",
        "frQCDdiJetSecondLeadJet",
        "frWplusJets",
        ##"frGammaPlusJets",
        "frSysUncertainty"
    ]
)

drawFakeRateHistConfiguratorZtoMuTau_background.addProcess(
    "Zmumu", processZtoMuTau_Zmumu.config_dqmHistPlotter.dqmDirectory)
drawFakeRateHistConfiguratorZtoMuTau_background.addProcess("WplusJets", processZtoMuTau_WplusJets.config_dqmHistPlotter.dqmDirectory)
drawFakeRateHistConfiguratorZtoMuTau_background.addProcess("TTplusJets", processZtoMuTau_TTplusJets.config_dqmHistPlotter.dqmDirectory)
drawFakeRateHistConfiguratorZtoMuTau_background.addProcess("QCD", cms.string('harvested/qcdSum'))