mvaMETSequence = cms.Sequence( goodPVFilter +
mvaMETDiTau +
recoilCorMETDiTau +
mvaBaseMETDiTau
# # +
# # cmgDiTauMVAPreSel
)
cmgDiTauCorPreSel = cmgDiTauCor.clone()
cmgDiTauCorPreSel.cfg.metCollection = 'mvaBaseMETDiTau'
cmgDiTauCorPreSel.cfg.diObjectCollection = 'cmgDiTauPreSel'
# SVFit
cmgDiTauCorSVFitPreSel = diTauSVFit.clone()
cmgDiTauCorSVFitPreSel.diTauSrc = 'cmgDiTauCorPreSel'
cmgDiTauCorSVFitPreSel.metsigSrc = 'mvaMETDiTau'
cmgDiTauCorSVFitFullSel = cmgDiTauSel.clone( src = 'cmgDiTauCorSVFitPreSel',
cut = 'leg1().pt()>45 && leg2().pt()>45',
)
diTauCorSVFitSequence = cms.Sequence( #
mvaMETSequence +
cmgDiTauCorPreSel +
cmgDiTauCorSVFitPreSel +
cmgDiTauCorSVFitFullSel
)
diTauSequence = cms.Sequence( diTauStdSequence +
diTauCorSVFitSequence
mvaMETDiTau.srcVertices = cms.InputTag("offlineSlimmedPrimaryVertices")
mvaMETDiTau.srcLeptons = cms.VInputTag(
cms.InputTag("tauPreSelectionDiTau", "", ""),
cms.InputTag("tauPreSelectionDiTau", "", ""),
)
mvaMETDiTau.permuteLeptons = cms.bool(True)
# correct TauES (after MVA MET according to current baseline)
cmgDiTauCor = cmgDiTauCor.clone()
# this selector goes after the TauES correction
cmgDiTauTauPtSel = cms.EDFilter(
"PATCompositeCandidateSelector",
src=cms.InputTag("cmgDiTauCor"),
cut=cms.string("daughter(0).pt()>45. && daughter(1).pt()>45."))
# recoil correction ----------------------------------------------------
# JAN: We don't know yet if we need this in 2015; re-include if necessary
# sequence -------------------------------------------------------------
diTauMVAMetSequence = cms.Sequence(mvaMETDiTau)
# SVFit ----------------------------------------------------------------
cmgDiTauCorSVFitPreSel = diTauSVFit.clone()
cmgDiTauCorSVFitFullSel = cmgDiTauSel.clone()
diTauSequence = cms.Sequence(tauPreSelectionDiTau + diTauMVAMetSequence +
cmgDiTau + cmgDiTauCor + cmgDiTauTauPtSel +
cmgDiTauCorSVFitPreSel + cmgDiTauCorSVFitFullSel)
mvaMETDiTau +
recoilCorMETDiTau
# +
# mvaBaseMETDiTau
# # +
# # cmgDiTauMVAPreSel
)
cmgDiTauCorPreSVFitSel = cmgDiTauSel.clone( src = 'recoilCorMETDiTau',
cut = 'leg1().pt()>45 && leg2().pt()>45',
)
# SVFit
cmgDiTauSVFit = diTauSVFit.clone()
cmgDiTauSVFit.diTauSrc = 'cmgDiTauCorPreSVFitSel'
# cmgDiTauSVFit.metsigSrc = 'mvaMETDiTau'
cmgDiTauCorSVFitFullSel = cmgDiTauSel.clone( src = 'cmgDiTauSVFit',
cut = '',
)
diTauCorSVFitSequence = cms.Sequence( #
mvaMETSequence +
# cmgDiTauCorPreSel +
cmgDiTauCorPreSVFitSel +
cmgDiTauSVFit +
cmgDiTauCorSVFitFullSel
)
diTauSequence = cms.Sequence( diTauStdSequence +