CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Functions
miniAOD_tools Namespace Reference

Functions

def miniAOD_customizeAllData
 
def miniAOD_customizeAllMC
 
def miniAOD_customizeCommon
 
def miniAOD_customizeData
 
def miniAOD_customizeMC
 
def miniAOD_customizeOutput
 

Function Documentation

def miniAOD_tools.miniAOD_customizeAllData (   process)

Definition at line 257 of file miniAOD_tools.py.

References miniAOD_customizeCommon(), and miniAOD_customizeData().

257 
258 def miniAOD_customizeAllData(process):
259  miniAOD_customizeCommon(process)
260  miniAOD_customizeData(process)
261  return process
miniAOD_tools.miniAOD_customizeCommon
def miniAOD_customizeCommon
Definition: miniAOD_tools.py:5
miniAOD_tools.miniAOD_customizeAllData
def miniAOD_customizeAllData
Definition: miniAOD_tools.py:257
miniAOD_tools.miniAOD_customizeData
def miniAOD_customizeData
Definition: miniAOD_tools.py:253
def miniAOD_tools.miniAOD_customizeAllMC (   process)

Definition at line 262 of file miniAOD_tools.py.

References miniAOD_customizeCommon(), and miniAOD_customizeMC().

262 
263 def miniAOD_customizeAllMC(process):
264  miniAOD_customizeCommon(process)
265  miniAOD_customizeMC(process)
266  return process
miniAOD_tools.miniAOD_customizeAllMC
def miniAOD_customizeAllMC
Definition: miniAOD_tools.py:262
miniAOD_tools.miniAOD_customizeMC
def miniAOD_customizeMC
Definition: miniAOD_tools.py:227
miniAOD_tools.miniAOD_customizeCommon
def miniAOD_customizeCommon
Definition: miniAOD_tools.py:5
def miniAOD_tools.miniAOD_customizeCommon (   process)

Definition at line 5 of file miniAOD_tools.py.

References applySubstructure_cff.applySubstructure(), python.multivaluedict.map(), vid_id_tools.setupAllVIDIdsInModule(), trigTools.switchOnTriggerStandAlone, and vid_id_tools.switchOnVIDElectronIdProducer().

Referenced by miniAOD_customizeAllData(), and miniAOD_customizeAllMC().

5 
6 def miniAOD_customizeCommon(process):
7  process.patMuons.isoDeposits = cms.PSet()
8  process.patElectrons.isoDeposits = cms.PSet()
9  process.patTaus.isoDeposits = cms.PSet()
10  process.patPhotons.isoDeposits = cms.PSet()
11  #
12  process.patMuons.embedTrack = True # used for IDs
13  process.patMuons.embedCombinedMuon = True # used for IDs
14  process.patMuons.embedMuonBestTrack = True # used for IDs
15  process.patMuons.embedStandAloneMuon = True # maybe?
16  process.patMuons.embedPickyMuon = False # no, use best track
17  process.patMuons.embedTpfmsMuon = False # no, use best track
18  process.patMuons.embedDytMuon = False # no, use best track
19  #
20  # disable embedding of electron and photon associated objects already stored by the ReducedEGProducer
21  process.patElectrons.embedGsfElectronCore = False ## process.patElectrons.embed in AOD externally stored gsf electron core
22  process.patElectrons.embedSuperCluster = False ## process.patElectrons.embed in AOD externally stored supercluster
23  process.patElectrons.embedPflowSuperCluster = False ## process.patElectrons.embed in AOD externally stored supercluster
24  process.patElectrons.embedSeedCluster = False ## process.patElectrons.embed in AOD externally stored the electron's seedcluster
25  process.patElectrons.embedBasicClusters = False ## process.patElectrons.embed in AOD externally stored the electron's basic clusters
26  process.patElectrons.embedPreshowerClusters = False ## process.patElectrons.embed in AOD externally stored the electron's preshower clusters
27  process.patElectrons.embedPflowBasicClusters = False ## process.patElectrons.embed in AOD externally stored the electron's pflow basic clusters
28  process.patElectrons.embedPflowPreshowerClusters = False ## process.patElectrons.embed in AOD externally stored the electron's pflow preshower clusters
29  process.patElectrons.embedRecHits = False ## process.patElectrons.embed in AOD externally stored the RecHits - can be called from the PATElectronProducer
30  process.patElectrons.electronSource = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
31  process.patElectrons.electronIDSources = cms.PSet(
32  # configure many IDs as InputTag <someName> = <someTag> you
33  # can comment out those you don't want to save some disk space
34  eidRobustLoose = cms.InputTag("reducedEgamma","eidRobustLoose"),
35  eidRobustTight = cms.InputTag("reducedEgamma","eidRobustTight"),
36  eidLoose = cms.InputTag("reducedEgamma","eidLoose"),
37  eidTight = cms.InputTag("reducedEgamma","eidTight"),
38  eidRobustHighEnergy = cms.InputTag("reducedEgamma","eidRobustHighEnergy"),
39  )
40  process.patElectrons.addPFClusterIso = cms.bool(True)
41  process.patElectrons.ecalPFClusterIsoMap = cms.InputTag("reducedEgamma", "eleEcalPFClusIso")
42  process.patElectrons.hcalPFClusterIsoMap = cms.InputTag("reducedEgamma", "eleHcalPFClusIso")
43 
44  process.elPFIsoDepositChargedPAT.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
45  process.elPFIsoDepositChargedAllPAT.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
46  process.elPFIsoDepositNeutralPAT.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
47  process.elPFIsoDepositGammaPAT.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
48  process.elPFIsoDepositPUPAT.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
49  #
50  process.patPhotons.embedSuperCluster = False ## whether to process.patPhotons.embed in AOD externally stored supercluster
51  process.patPhotons.embedSeedCluster = False ## process.patPhotons.embed in AOD externally stored the photon's seedcluster
52  process.patPhotons.embedBasicClusters = False ## process.patPhotons.embed in AOD externally stored the photon's basic clusters
53  process.patPhotons.embedPreshowerClusters = False ## process.patPhotons.embed in AOD externally stored the photon's preshower clusters
54  process.patPhotons.embedRecHits = False ## process.patPhotons.embed in AOD externally stored the RecHits - can be called from the PATPhotonProducer
55  process.patPhotons.addPFClusterIso = cms.bool(True)
56  process.patPhotons.ecalPFClusterIsoMap = cms.InputTag("reducedEgamma", "phoEcalPFClusIso")
57  process.patPhotons.hcalPFClusterIsoMap = cms.InputTag("reducedEgamma", "phoHcalPFClusIso")
58  process.patPhotons.photonSource = cms.InputTag("reducedEgamma","reducedGedPhotons")
59  process.patPhotons.electronSource = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
60  process.patPhotons.photonIDSources = cms.PSet(
61  PhotonCutBasedIDLoose = cms.InputTag('reducedEgamma',
62  'PhotonCutBasedIDLoose'),
63  PhotonCutBasedIDTight = cms.InputTag('reducedEgamma',
64  'PhotonCutBasedIDTight')
65  )
66 
67  process.phPFIsoDepositChargedPAT.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
68  process.phPFIsoDepositChargedAllPAT.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
69  process.phPFIsoDepositNeutralPAT.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
70  process.phPFIsoDepositGammaPAT.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
71  process.phPFIsoDepositPUPAT.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
72  #
73  process.selectedPatJets.cut = cms.string("pt > 10")
74  process.selectedPatMuons.cut = cms.string("pt > 5 || isPFMuon || (pt > 3 && (isGlobalMuon || isStandAloneMuon || numberOfMatches > 0 || muonID('RPCMuLoose')))")
75  process.selectedPatElectrons.cut = cms.string("")
76  process.selectedPatTaus.cut = cms.string("pt > 18. && tauID('decayModeFindingNewDMs')> 0.5")
77  process.selectedPatPhotons.cut = cms.string("")
78 
79  from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection
80 
81  from PhysicsTools.PatAlgos.slimming.applySubstructure_cff import applySubstructure
82  applySubstructure( process )
83 
84 
85  #
86  from PhysicsTools.PatAlgos.tools.trigTools import switchOnTriggerStandAlone
87  switchOnTriggerStandAlone( process, outputModule = '' )
88  process.patTrigger.packTriggerPathNames = cms.bool(True)
89  #
90  # apply type I/type I + II PFMEt corrections to pat::MET object
91  # and estimate systematic uncertainties on MET
92  # FIXME: are we 100% sure this should still be PF and not PFchs?
93  from PhysicsTools.PatUtils.tools.runType1PFMEtUncertainties import runType1PFMEtUncertainties
94  addJetCollection(process, postfix = "ForMetUnc", labelName = 'AK4PF', jetSource = cms.InputTag('ak4PFJets'), jetCorrections = ('AK4PF', ['L1FastJet', 'L2Relative', 'L3Absolute'], ''))
95  process.patJetsAK4PFForMetUnc.getJetMCFlavour = False
96  runType1PFMEtUncertainties(process,
97  addToPatDefaultSequence=False,
98  jetCollectionUnskimmed="patJetsAK4PFForMetUnc",
99  jetCollection="selectedPatJetsAK4PFForMetUnc",
100  electronCollection="selectedPatElectrons",
101  muonCollection="selectedPatMuons",
102  tauCollection="selectedPatTaus",
103  makeType1p2corrPFMEt=True,
104  doSmearJets=False,
105  outputModule=None)
106 
107  from PhysicsTools.PatAlgos.tools.metTools import addMETCollection
108  addMETCollection(process,
109  labelName = "patCaloMet",
110  metSource = "caloMetM"
111  )
112 
113 
114  #keep this after all addJetCollections otherwise it will attempt computing them also for stuf with no taginfos
115  #Some useful BTAG vars
116  process.patJets.userData.userFunctions = cms.vstring(
117  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).p4.M):(0)',
118  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).numberOfSourceCandidatePtrs):(0)',
119  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").flightDistance(0).value):(0)',
120  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").flightDistance(0).significance):(0)',
121  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).p4.x):(0)',
122  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).p4.y):(0)',
123  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).p4.z):(0)',
124  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).vertex.x):(0)',
125  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).vertex.y):(0)',
126  '?(tagInfoCandSecondaryVertex("pfSecondaryVertex").nVertices()>0)?(tagInfoCandSecondaryVertex("pfSecondaryVertex").secondaryVertex(0).vertex.z):(0)',
127  )
128  process.patJets.userData.userFunctionLabels = cms.vstring('vtxMass','vtxNtracks','vtx3DVal','vtx3DSig','vtxPx','vtxPy','vtxPz','vtxPosX','vtxPosY','vtxPosZ')
129  process.patJets.tagInfoSources = cms.VInputTag(cms.InputTag("pfSecondaryVertexTagInfos"))
130  process.patJets.addTagInfos = cms.bool(True)
131  #
132  ## PU JetID
133  process.load("PhysicsTools.PatAlgos.slimming.pileupJetId_cfi")
134  process.patJets.userData.userFloats.src = [ cms.InputTag("pileupJetId:fullDiscriminant"), ]
135 
136  ## CaloJets
137  process.caloJetMap = cms.EDProducer("RecoJetDeltaRValueMapProducer",
138  src = process.patJets.jetSource,
139  matched = cms.InputTag("ak4CaloJets"),
140  distMax = cms.double(0.4),
141  values = cms.vstring('pt','emEnergyFraction'),
142  valueLabels = cms.vstring('pt','emEnergyFraction'),
143  lazyParser = cms.bool(True) )
144  process.patJets.userData.userFloats.src += [ cms.InputTag("caloJetMap:pt"), cms.InputTag("caloJetMap:emEnergyFraction") ]
145 
146  #VID Electron IDs
147  electron_ids = ['RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_CSA14_50ns_V1_cff',
148  'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_CSA14_PU20bx25_V0_cff',
149  'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV50_CSA14_25ns_cff',
150  'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV50_CSA14_startup_cff']
151  switchOnVIDElectronIdProducer(process, DataFormat.MiniAOD)
152  process.egmGsfElectronIDs.physicsObjectSrc = \
153  cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
154  process.electronIDValueMapProducer.src = \
155  cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
156  process.electronIDValueMapProducer.ebReducedRecHitCollection = \
157  cms.InputTag("reducedEgamma","reducedEBRecHits")
158  process.electronIDValueMapProducer.eeReducedRecHitCollection = \
159  cms.InputTag("reducedEgamma","reducedEERecHits")
160  process.electronIDValueMapProducer.esReducedRecHitCollection = \
161  cms.InputTag("reducedEgamma","reducedESRecHits")
162  for idmod in electron_ids:
163  setupAllVIDIdsInModule(process,idmod,setupVIDElectronSelection)
164 
165  # Adding puppi jets
166  process.load('CommonTools.PileupAlgos.Puppi_cff')
167  process.load('RecoJets.JetProducers.ak4PFJetsPuppi_cfi')
168  process.ak4PFJetsPuppi.doAreaFastjet = True # even for standard ak4PFJets this is overwritten in RecoJets/Configuration/python/RecoPFJets_cff
169  #process.puppi.candName = cms.InputTag('packedPFCandidates')
170  #process.puppi.vertexName = cms.InputTag('offlineSlimmedPrimaryVertices')
171 
172  from RecoJets.JetAssociationProducers.j2tParametersVX_cfi import j2tParametersVX
173  process.ak4PFJetsPuppiTracksAssociatorAtVertex = cms.EDProducer("JetTracksAssociatorAtVertex",
174  j2tParametersVX,
175  jets = cms.InputTag("ak4PFJetsPuppi")
176  )
177  process.patJetPuppiCharge = cms.EDProducer("JetChargeProducer",
178  src = cms.InputTag("ak4PFJetsPuppiTracksAssociatorAtVertex"),
179  var = cms.string('Pt'),
180  exp = cms.double(1.0)
181  )
182 
183  addJetCollection(process, postfix = "", labelName = 'Puppi', jetSource = cms.InputTag('ak4PFJetsPuppi'),
184  jetCorrections = ('AK4PFchs', ['L2Relative', 'L3Absolute'], ''),
185  algo= 'AK', rParam = 0.4, btagDiscriminators = map(lambda x: x.value() ,process.patJets.discriminatorSources)
186  )
187 
188  process.patJetGenJetMatchPuppi.matched = 'slimmedGenJets'
189 
190  process.patJetsPuppi.userData.userFloats.src = cms.VInputTag(cms.InputTag(""))
191  process.patJetsPuppi.jetChargeSource = cms.InputTag("patJetPuppiCharge")
192 
193  process.selectedPatJetsPuppi.cut = cms.string("pt > 20")
194 
195  process.load('PhysicsTools.PatAlgos.slimming.slimmedJets_cfi')
196  process.slimmedJetsPuppi = process.slimmedJets.clone()
197  process.slimmedJetsPuppi.src = cms.InputTag("selectedPatJetsPuppi")
198  process.slimmedJetsPuppi.packedPFCandidates = cms.InputTag("packedPFCandidates")
199 
200  ## puppi met
201  process.load('RecoMET.METProducers.PFMET_cfi')
202  process.pfMetPuppi = process.pfMet.clone()
203  process.pfMetPuppi.src = cms.InputTag("puppi")
204  process.pfMetPuppi.alias = cms.string('pfMetPuppi')
205  ## type1 correction, from puppi jets
206  process.corrPfMetType1Puppi = process.corrPfMetType1.clone(
207  src = 'ak4PFJetsPuppi',
208  jetCorrLabel = 'ak4PFCHSL2L3Corrector',
209  )
210  del process.corrPfMetType1Puppi.offsetCorrLabel # no L1 for PUPPI jets
211  process.pfMetT1Puppi = process.pfMetT1.clone(
212  src = 'pfMetPuppi',
213  srcCorrections = [ cms.InputTag("corrPfMetType1Puppi","type1") ]
214  )
215 
216  from PhysicsTools.PatAlgos.tools.metTools import addMETCollection
217  addMETCollection(process, labelName='patMETPuppi', metSource='pfMetT1Puppi') # T1
218  addMETCollection(process, labelName='patPFMetPuppi', metSource='pfMetPuppi') # RAW
219 
220  process.load('PhysicsTools.PatAlgos.slimming.slimmedMETs_cfi')
221  process.slimmedMETsPuppi = process.slimmedMETs.clone()
222  process.slimmedMETsPuppi.src = cms.InputTag("patMETPuppi")
223  process.slimmedMETsPuppi.rawUncertainties = cms.InputTag("patPFMetPuppi") # only central value
224  process.slimmedMETsPuppi.type1Uncertainties = cms.InputTag("patPFMetT1") # only central value for now
225  del process.slimmedMETsPuppi.type1p2Uncertainties # not available
226 
vid_id_tools.setupAllVIDIdsInModule
def setupAllVIDIdsInModule
Definition: vid_id_tools.py:38
vid_id_tools.switchOnVIDElectronIdProducer
def switchOnVIDElectronIdProducer
Definition: vid_id_tools.py:57
applySubstructure_cff.applySubstructure
def applySubstructure
Definition: applySubstructure_cff.py:3
miniAOD_tools.miniAOD_customizeCommon
def miniAOD_customizeCommon
Definition: miniAOD_tools.py:5
j2tParametersVX_cfi
trigTools.switchOnTriggerStandAlone
tuple switchOnTriggerStandAlone
Definition: trigTools.py:276
def miniAOD_tools.miniAOD_customizeData (   process)

Definition at line 253 of file miniAOD_tools.py.

Referenced by miniAOD_customizeAllData().

253 
254 def miniAOD_customizeData(process):
255  from PhysicsTools.PatAlgos.tools.coreTools import runOnData
256  runOnData( process, outputModules = [] )
miniAOD_tools.miniAOD_customizeData
def miniAOD_customizeData
Definition: miniAOD_tools.py:253
def miniAOD_tools.miniAOD_customizeMC (   process)

Definition at line 227 of file miniAOD_tools.py.

Referenced by miniAOD_customizeAllMC().

227 
228 def miniAOD_customizeMC(process):
229  process.muonMatch.matched = "prunedGenParticles"
230  process.electronMatch.matched = "prunedGenParticles"
231  process.electronMatch.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
232  process.photonMatch.matched = "prunedGenParticles"
233  process.photonMatch.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
234  process.tauMatch.matched = "prunedGenParticles"
235  process.tauGenJets.GenParticles = "prunedGenParticles"
236  process.patJetPartonMatch.matched = "prunedGenParticles"
237  process.patJetPartonMatch.mcStatus = [ 3, 23 ]
238  process.patJetGenJetMatch.matched = "slimmedGenJets"
239  process.patJetGenJetMatchAK8.matched = "slimmedGenJetsAK8"
240  process.patMuons.embedGenMatch = False
241  process.patElectrons.embedGenMatch = False
242  process.patPhotons.embedGenMatch = False
243  process.patTaus.embedGenMatch = False
244  process.patJets.embedGenPartonMatch = False
245  #also jet flavour must be switched to ak4
246  process.patJetFlavourAssociation.rParam = 0.4
miniAOD_tools.miniAOD_customizeMC
def miniAOD_customizeMC
Definition: miniAOD_tools.py:227
def miniAOD_tools.miniAOD_customizeOutput (   out)

Definition at line 247 of file miniAOD_tools.py.

247 
248 def miniAOD_customizeOutput(out):
249  out.dropMetaData = cms.untracked.string('ALL')
250  out.fastCloning= cms.untracked.bool(False)
251  out.overrideInputFileSplitLevels = cms.untracked.bool(True)
252  out.compressionAlgorithm = cms.untracked.string('LZMA')
miniAOD_tools.miniAOD_customizeOutput
def miniAOD_customizeOutput
Definition: miniAOD_tools.py:247
Generated for CMSSW Reference Manual by doxygen 1.8.5