#
# Additional comment (BA, 04/20/06)
# the important lines for MET correction
# after EMJet removal have been highlighted
# in blue (flag to change from default definition),
# red (EM objects used for EMJet removal),
# or green (jets selected after EMJet removal)
#
# !/usr/bin/env cafe
#
# Top MET Treatment
######################################################

+SingleTopMET.Run: Group(MetCorrection)
+SingleTopMET.Run: METSelector(stMet)

+TTbarLepJetsMET.Run: Group(MetCorrection)
+TTbarLepJetsMET.Run: METSelector(ttMet)

+DiElectronMET.Run: ReComputeMET(diem_met) METSelector(diem_alibimet)
+DiMuonMET.Run: ReComputeMET(dimu_met) METSelector(dimu_metcut)

+MetCorrection.Run: ElectronSelector(em_for_met)
+MetCorrection.Run: EMJetMatching(jet_for_met)
+MetCorrection.Run: MuonSelector(muon_for_met)
+MetCorrection.Run: ReComputeMET(MetCorr)

######################################################


# MET correction
######################################################

MetCorr.From: Met
MetCorr.To: CorrMet
MetCorr.Tree: SelectedObjects
MetCorr.electronBranch: MetEM
MetCorr.muonBranch: MetMuon
MetCorr.jetConeSize: 0.5
MetCorr.jetBranch: MetJCCB
MetCorr.useDefaultGoodJetDef: false
#MetCorr.Debug: true

# single top MET selection
######################################################

stMet.From: CorrMet
stMet.minMET: 15
stMet.maxMET: 200
stMet.jetAlgoType: JCCB
stMet.jesMuCorrected: false
stMet.MuCorrected: true

# default value of MCsmeared is true
# set it to false if the electrons, muons are not over-smeared.
stMet.MCsmeared: false

# ttbar MET selection
######################################################

ttMet.From: CorrMet
ttMet.minMET: 15
ttMet.jetAlgoType: JCCB
ttMet.jesMuCorrected: false
ttMet.MuCorrected: true

# default value of MCsmeared is true
# set it to false if the electrons, muons are not over-smeared.
ttMet.MCsmeared: false

############ MET for ttbar -> 2e+2b+2nu (just propagate new JES) ##############

diem_met.From: Met
diem_met.To: corrMet
diem_met.Tree: SelectObjects
diem_met.jetBranch: CorrJCCB
diem_met.jetConeSize: 0.5
diem_met.MuCorrected: true

diem_alibimet.From: corrMet
diem_alibimet.minMET: 0.
diem_alibimet.maxMET: 0.

############ MET for ttbar -> 2mu+2b+2nu ##############

dimu_met.From: Met
dimu_met.To: corrMet
dimu_met.Tree: SelectObjects
dimu_met.jetBranch: CorrJCCB
dimu_met.jetConeSize: 0.5
dimu_met.muonBranch: TightMuon
diem_met.MuCorrected: true

########### Need this for the triangle cuts!

dimu_metcut.From: corrMet
dimu_metcut.minMET: 45.

######################################################