File indexing completed on 2023-03-17 11:16:48
0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TagProbe")
0004
0005
0006
0007
0008
0009 process.load('FWCore.MessageService.MessageLogger_cfi')
0010 process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0011 process.MessageLogger.cerr.FwkReport.reportEvery = 100
0012
0013 process.source = cms.Source("PoolSource",
0014 fileNames = cms.untracked.vstring(
0015 '/store/relval/CMSSW_7_2_1/RelValZMM_13/GEN-SIM-RECO/PU50ns_PHYS14_25_V1_Phys14-v1/00000/287B9489-B85E-E411-95DF-02163E00EB3F.root'
0016 )
0017 )
0018 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(500) )
0019
0020
0021
0022
0023 process.tagMuons = cms.EDFilter("MuonRefSelector",
0024 src = cms.InputTag("muons"),
0025 cut = cms.string("isGlobalMuon && pt > 20 && abs(eta) < 2"),
0026 )
0027
0028 process.probeMuons = cms.EDFilter("MuonRefSelector",
0029 src = cms.InputTag("muons"),
0030 cut = cms.string("isTrackerMuon && pt > 10"),
0031 )
0032
0033
0034
0035 process.probesPassingCal = cms.EDFilter("MuonRefSelector",
0036 src = cms.InputTag("muons"),
0037 cut = cms.string(process.probeMuons.cut.value() + " && caloCompatibility > 0.6"),
0038 )
0039
0040
0041 process.drToNearestJet = cms.EDProducer("DeltaRNearestJetComputer",
0042 probes = cms.InputTag("muons"),
0043
0044
0045 objects = cms.InputTag("ak5CaloJets"),
0046 objectSelection = cms.InputTag("et > 20 && abs(eta) < 3 && n60 > 3 && (.05 < emEnergyFraction < .95)"),
0047 )
0048
0049
0050 process.tpPairs = cms.EDProducer("CandViewShallowCloneCombiner",
0051 decay = cms.string("tagMuons@+ probeMuons@-"),
0052 cut = cms.string("40 < mass < 200"),
0053 )
0054
0055
0056 process.muMcMatch = cms.EDFilter("MCTruthDeltaRMatcherNew",
0057 pdgId = cms.vint32(13),
0058 src = cms.InputTag("muons"),
0059 distMin = cms.double(0.3),
0060 matched = cms.InputTag("genParticles")
0061 )
0062
0063
0064 process.muonEffs = cms.EDAnalyzer("TagProbeFitTreeProducer",
0065
0066 tagProbePairs = cms.InputTag("tpPairs"),
0067 arbitration = cms.string("OneProbe"),
0068
0069 variables = cms.PSet(
0070
0071 eta = cms.string("eta"),
0072 pt = cms.string("pt"),
0073
0074 nsegm = cms.string("numberOfMatches"),
0075
0076 drj = cms.InputTag("drToNearestJet"),
0077 ),
0078
0079 flags = cms.PSet(
0080
0081 passingCal = cms.InputTag("probesPassingCal"),
0082
0083 passingGlb = cms.string("isGlobalMuon"),
0084 passingIso = cms.string("(isolationR03.hadEt+isolationR03.emEt+isolationR03.sumPt) < 0.1 * pt"),
0085 ),
0086
0087 isMC = cms.bool(True),
0088 motherPdgId = cms.vint32(22,23),
0089 makeMCUnbiasTree = cms.bool(True),
0090 checkMotherInUnbiasEff = cms.bool(True),
0091 tagMatches = cms.InputTag("muMcMatch"),
0092 probeMatches = cms.InputTag("muMcMatch"),
0093 allProbes = cms.InputTag("probeMuons"),
0094 )
0095
0096
0097
0098
0099
0100
0101 process.tagAndProbe = cms.Path(
0102 (process.tagMuons + process.probeMuons) *
0103 (process.probesPassingCal +
0104 process.drToNearestJet +
0105 process.tpPairs +
0106 process.muMcMatch) *
0107 process.muonEffs
0108 )
0109
0110 process.TFileService = cms.Service("TFileService", fileName = cms.string("testTagProbeFitTreeProducer_ZMuMu.root"))
0111
0112
0113
0114