Ejemplo n.º 1
0
if numberOfFilesToProcess > 0:
    process.source.fileNames = process.source.fileNames[:numberOfFilesToProcess]

addPuppi = True

# Adding jet collection
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.GlobalTag.globaltag = 'MCRUN2_74_V9::All'
# process.GlobalTag.globaltag = 'auto:run2_mc'

process.load("Configuration.StandardSequences.Geometry_cff")
process.load("Configuration.StandardSequences.MagneticField_38T_cff")

from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection
process.tauPreSelection = tauPreSelection.clone()

process.tauCount = cms.EDFilter(
    "CandViewCountFilter",
    src=cms.InputTag("tauPreSelection"),
    minNumber=cms.uint32(1),
)

if addPuppi:
    process.load('CommonTools/PileupAlgos/Puppi_cff')

    process.puppi.candName = cms.InputTag('packedPFCandidates')
    process.puppi.vertexName = cms.InputTag('offlineSlimmedPrimaryVertices')

    process.packedPFCandidatesWoMuon = cms.EDFilter("CandPtrSelector", src=cms.InputTag("packedPFCandidates"), cut=cms.string("fromPV>=2 && abs(pdgId)!=13 "))
    process.particleFlowNoMuonPUPPI = process.puppi.clone()
import FWCore.ParameterSet.Config as cms

from CMGTools.H2TauTau.objects.cmgTauMu_cfi import cmgTauMu
from CMGTools.H2TauTau.skims.cmgTauMuSel_cfi import cmgTauMuSel

from CMGTools.H2TauTau.objects.cmgTauMuCor_cfi import cmgTauMuCor
from CMGTools.H2TauTau.objects.tauMuSVFit_cfi import tauMuSVFit

from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection
from CMGTools.H2TauTau.objects.muCuts_cff import muonPreSelection

from CMGTools.H2TauTau.skims.skim_cff import tauMuFullSelSkimSequence, tauMuFullSelCount

# tau pre-selection
tauPreSelectionTauMu = tauPreSelection.clone()
muonPreSelectionTauMu = muonPreSelection.clone()

# Correct tau pt (after MVA MET according to current baseline)
cmgTauMuCor = cmgTauMuCor.clone()

# This selector goes after the tau pt correction
cmgTauMuTauPtSel = cms.EDFilter("PATCompositeCandidateSelector",
                                src=cms.InputTag("cmgTauMuCor"),
                                cut=cms.string("daughter(0).pt()>18."))

cmgTauMuTauPtSel = cmgTauMuTauPtSel.clone()

# SVFit
cmgTauMuCorSVFitPreSel = tauMuSVFit.clone()

# If you want to apply some extra selection after SVFit, do it here
Ejemplo n.º 3
0
import FWCore.ParameterSet.Config as cms

from CMGTools.H2TauTau.skims.cmgDiTauSel_cfi import cmgDiTauSel

from CMGTools.H2TauTau.objects.cmgDiTau_cfi import cmgDiTau
from CMGTools.H2TauTau.objects.cmgDiTauCor_cfi import cmgDiTauCor
from CMGTools.H2TauTau.objects.diTauSVFit_cfi import diTauSVFit
from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection

from RecoMET.METPUSubtraction.mvaPFMET_cff import pfMVAMEt

# tau pre-selection
tauPreSelectionDiTau = tauPreSelection.clone(
    #cut = 'pt > 40. && abs(eta) < 2.5 && tauID("decayModeFinding") > 0.5')
    cut='pt > 40. && abs(eta) < 2.5 && tauID("decayModeFindingNewDMs") > 0.5'
)  # RIC: new DM. Probably we'd want to save both with an OR and decide later, useful for studying new tauID

# 2012 preselection:
# cut = 'leg1().pt()>40. && leg2().pt()>40. && leg1().tauID("byCombinedIsolationDeltaBetaCorrRaw3Hits") < 10. &&  leg2().tauID("byCombinedIsolationDeltaBetaCorrRaw3Hits") < 10.',

# mva MET

mvaMETDiTau = pfMVAMEt.clone()

mvaMETDiTau.srcPFCandidates = cms.InputTag("packedPFCandidates")
mvaMETDiTau.srcVertices = cms.InputTag("offlineSlimmedPrimaryVertices")
mvaMETDiTau.srcLeptons = cms.VInputTag(
    cms.InputTag("tauPreSelectionDiTau", "", ""),
    cms.InputTag("tauPreSelectionDiTau", "", ""),
)
mvaMETDiTau.permuteLeptons = cms.bool(True)
import FWCore.ParameterSet.Config as cms

from CMGTools.H2TauTau.objects.cmgTauEle_cfi import cmgTauEle
from CMGTools.H2TauTau.skims.cmgTauEleSel_cfi import cmgTauEleSel

from CMGTools.H2TauTau.objects.cmgTauEleCor_cfi import cmgTauEleCor
from CMGTools.H2TauTau.objects.tauEleSVFit_cfi import tauEleSVFit

from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection
from CMGTools.H2TauTau.objects.eleCuts_cff import electronPreSelection

from CMGTools.H2TauTau.skims.skim_cff import tauEleFullSelSkimSequence, tauEleFullSelCount

# tau pre-selection
tauPreSelectionTauEle = tauPreSelection.clone()
electronPreSelectionTauEle = electronPreSelection.clone()

# Correct tau pt (after MVA MET according to current baseline)
cmgTauEleCor = cmgTauEleCor.clone()

# This selector goes after the tau pt correction
cmgTauEleTauPtSel = cms.EDFilter("PATCompositeCandidateSelector",
                                 src=cms.InputTag("cmgTauEleCor"),
                                 cut=cms.string("daughter(0).pt()>18."))

cmgTauEleTauPtSel = cmgTauEleTauPtSel.clone()

# SVFit
cmgTauEleCorSVFitPreSel = tauEleSVFit.clone()

# If you want to apply some extra selection after SVFit, do it here
import FWCore.ParameterSet.Config as cms

from CMGTools.H2TauTau.objects.cmgTauEle_cfi import cmgTauEle
from CMGTools.H2TauTau.skims.cmgTauEleSel_cfi import cmgTauEleSel

from CMGTools.H2TauTau.objects.cmgTauEleCor_cfi import cmgTauEleCor 
from CMGTools.H2TauTau.objects.tauEleSVFit_cfi import tauEleSVFit 

from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection
from CMGTools.H2TauTau.objects.eleCuts_cff import electronPreSelection

from RecoMET.METPUSubtraction.mvaPFMET_cff import pfMVAMEt

# tau pre-selection
tauPreSelectionTauEle = tauPreSelection.clone()
electronPreSelectionTauEle = electronPreSelection.clone()

# mva MET
mvaMETTauEle = cms.EDProducer('PFMETProducerMVATauTau', 
                             **pfMVAMEt.parameters_())

mvaMETTauEle.srcPFCandidates = cms.InputTag("packedPFCandidates")
mvaMETTauEle.srcVertices = cms.InputTag("offlineSlimmedPrimaryVertices")
mvaMETTauEle.srcLeptons = cms.VInputTag(
  cms.InputTag("tauPreSelectionTauEle", "", ""),
  cms.InputTag("electronPreSelectionTauEle", "", ""),
  )
mvaMETTauEle.permuteLeptons = cms.bool(True)


Ejemplo n.º 6
0
    process.source.fileNames = process.source.fileNames[:
                                                        numberOfFilesToProcess]

addPuppi = True

# Adding jet collection
process.load(
    'Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.GlobalTag.globaltag = 'MCRUN2_74_V9::All'
# process.GlobalTag.globaltag = 'auto:run2_mc'

process.load("Configuration.StandardSequences.Geometry_cff")
process.load("Configuration.StandardSequences.MagneticField_38T_cff")

from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection
process.tauPreSelection = tauPreSelection.clone()

process.tauCount = cms.EDFilter(
    "CandViewCountFilter",
    src=cms.InputTag("tauPreSelection"),
    minNumber=cms.uint32(1),
)

if addPuppi:
    process.load('CommonTools/PileupAlgos/Puppi_cff')

    process.puppi.candName = cms.InputTag('packedPFCandidates')
    process.puppi.vertexName = cms.InputTag('offlineSlimmedPrimaryVertices')

    process.packedPFCandidatesWoMuon = cms.EDFilter(
        "CandPtrSelector",
import FWCore.ParameterSet.Config as cms

from CMGTools.H2TauTau.skims.cmgDiTauSel_cfi   import cmgDiTauSel

from CMGTools.H2TauTau.objects.cmgDiTau_cfi    import cmgDiTau
from CMGTools.H2TauTau.objects.cmgDiTauCor_cfi import cmgDiTauCor 
from CMGTools.H2TauTau.objects.diTauSVFit_cfi  import diTauSVFit 
from CMGTools.H2TauTau.objects.tauCuts_cff     import tauPreSelection

from RecoMET.METPUSubtraction.mvaPFMET_cff     import pfMVAMEt

# tau pre-selection
tauPreSelectionDiTau = tauPreSelection.clone(
  #cut = 'pt > 40. && abs(eta) < 2.5 && tauID("decayModeFinding") > 0.5')
  cut = 'pt > 40. && abs(eta) < 2.5 && tauID("decayModeFindingNewDMs") > 0.5') # RIC: new DM. Probably we'd want to save both with an OR and decide later, useful for studying new tauID

# 2012 preselection:
# cut = 'leg1().pt()>40. && leg2().pt()>40. && leg1().tauID("byCombinedIsolationDeltaBetaCorrRaw3Hits") < 10. &&  leg2().tauID("byCombinedIsolationDeltaBetaCorrRaw3Hits") < 10.',

# mva MET

mvaMETDiTau = cms.EDProducer('PFMETProducerMVATauTau', 
                             **pfMVAMEt.parameters_())

mvaMETDiTau.srcPFCandidates = cms.InputTag("packedPFCandidates")
mvaMETDiTau.srcVertices = cms.InputTag("offlineSlimmedPrimaryVertices")
mvaMETDiTau.srcLeptons = cms.VInputTag(
  cms.InputTag("tauPreSelectionDiTau", "", ""),
  cms.InputTag("tauPreSelectionDiTau", "", ""),
  )
mvaMETDiTau.permuteLeptons = cms.bool(True)
import FWCore.ParameterSet.Config as cms

from CMGTools.H2TauTau.objects.cmgTauMu_cfi import cmgTauMu
from CMGTools.H2TauTau.skims.cmgTauMuSel_cfi import cmgTauMuSel

from CMGTools.H2TauTau.objects.cmgTauMuCor_cfi import cmgTauMuCor 
from CMGTools.H2TauTau.objects.tauMuSVFit_cfi import tauMuSVFit 

from CMGTools.H2TauTau.objects.tauCuts_cff import tauPreSelection
from CMGTools.H2TauTau.objects.muCuts_cff import muonPreSelection

from RecoMET.METPUSubtraction.mvaPFMET_cff import pfMVAMEt

# tau pre-selection
tauPreSelectionTauMu = tauPreSelection.clone()
muonPreSelectionTauMu = muonPreSelection.clone()

# mva MET
mvaMETTauMu = cms.EDProducer('PFMETProducerMVATauTau', 
                             **pfMVAMEt.parameters_())#pfMVAMEt.clone()

mvaMETTauMu.srcPFCandidates = cms.InputTag("packedPFCandidates")
mvaMETTauMu.srcVertices = cms.InputTag("offlineSlimmedPrimaryVertices")
mvaMETTauMu.srcLeptons = cms.VInputTag(
  cms.InputTag("tauPreSelectionTauMu", "", ""),
  cms.InputTag("muonPreSelectionTauMu", "", ""),
  )
mvaMETTauMu.permuteLeptons = cms.bool(True)


# Correct tau pt (after MVA MET according to current baseline)