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 247 of file miniAOD_tools.py.

References miniAOD_customizeCommon(), and miniAOD_customizeData().

247 
248 def miniAOD_customizeAllData(process):
249  miniAOD_customizeCommon(process)
250  miniAOD_customizeData(process)
251  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:247
miniAOD_tools.miniAOD_customizeData
def miniAOD_customizeData
Definition: miniAOD_tools.py:243
def miniAOD_tools.miniAOD_customizeAllMC (   process)

Definition at line 252 of file miniAOD_tools.py.

References miniAOD_customizeCommon(), and miniAOD_customizeMC().

252 
253 def miniAOD_customizeAllMC(process):
254  miniAOD_customizeCommon(process)
255  miniAOD_customizeMC(process)
256  return process
miniAOD_tools.miniAOD_customizeAllMC
def miniAOD_customizeAllMC
Definition: miniAOD_tools.py:252
miniAOD_tools.miniAOD_customizeMC
def miniAOD_customizeMC
Definition: miniAOD_tools.py:217
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  #VID Electron IDs
137  electron_ids = ['RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_CSA14_50ns_V1_cff',
138  'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_CSA14_PU20bx25_V0_cff',
139  'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV50_CSA14_25ns_cff',
140  'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV50_CSA14_startup_cff']
141  switchOnVIDElectronIdProducer(process)
142  process.egmGsfElectronIDs.physicsObjectSrc = \
143  cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
144  process.electronIDValueMapProducer.src = \
145  cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
146  process.electronIDValueMapProducer.ebReducedRecHitCollection = \
147  cms.InputTag("reducedEgamma","reducedEBRecHits")
148  process.electronIDValueMapProducer.eeReducedRecHitCollection = \
149  cms.InputTag("reducedEgamma","reducedEERecHits")
150  process.electronIDValueMapProducer.esReducedRecHitCollection = \
151  cms.InputTag("reducedEgamma","reducedESRecHits")
152  for idmod in electron_ids:
153  setupAllVIDIdsInModule(process,idmod,setupVIDElectronSelection)
154 
155  # Adding puppi jets
156  process.load('CommonTools.PileupAlgos.Puppi_cff')
157  process.load('RecoJets.JetProducers.ak4PFJetsPuppi_cfi')
158  process.ak4PFJetsPuppi.doAreaFastjet = True # even for standard ak4PFJets this is overwritten in RecoJets/Configuration/python/RecoPFJets_cff
159  #process.puppi.candName = cms.InputTag('packedPFCandidates')
160  #process.puppi.vertexName = cms.InputTag('offlineSlimmedPrimaryVertices')
161 
162  from RecoJets.JetAssociationProducers.j2tParametersVX_cfi import j2tParametersVX
163  process.ak4PFJetsPuppiTracksAssociatorAtVertex = cms.EDProducer("JetTracksAssociatorAtVertex",
164  j2tParametersVX,
165  jets = cms.InputTag("ak4PFJetsPuppi")
166  )
167  process.patJetPuppiCharge = cms.EDProducer("JetChargeProducer",
168  src = cms.InputTag("ak4PFJetsPuppiTracksAssociatorAtVertex"),
169  var = cms.string('Pt'),
170  exp = cms.double(1.0)
171  )
172 
173  addJetCollection(process, postfix = "", labelName = 'Puppi', jetSource = cms.InputTag('ak4PFJetsPuppi'),
174  jetCorrections = ('AK4PFchs', ['L2Relative', 'L3Absolute'], ''),
175  algo= 'AK', rParam = 0.4, btagDiscriminators = map(lambda x: x.value() ,process.patJets.discriminatorSources)
176  )
177 
178  process.patJetGenJetMatchPuppi.matched = 'slimmedGenJets'
179 
180  process.patJetsPuppi.userData.userFloats.src = cms.VInputTag(cms.InputTag(""))
181  process.patJetsPuppi.jetChargeSource = cms.InputTag("patJetPuppiCharge")
182 
183  process.selectedPatJetsPuppi.cut = cms.string("pt > 20")
184 
185  process.load('PhysicsTools.PatAlgos.slimming.slimmedJets_cfi')
186  process.slimmedJetsPuppi = process.slimmedJets.clone()
187  process.slimmedJetsPuppi.src = cms.InputTag("selectedPatJetsPuppi")
188  process.slimmedJetsPuppi.packedPFCandidates = cms.InputTag("packedPFCandidates")
189 
190  ## puppi met
191  process.load('RecoMET.METProducers.PFMET_cfi')
192  process.pfMetPuppi = process.pfMet.clone()
193  process.pfMetPuppi.src = cms.InputTag("puppi")
194  process.pfMetPuppi.alias = cms.string('pfMetPuppi')
195  ## type1 correction, from puppi jets
196  process.corrPfMetType1Puppi = process.corrPfMetType1.clone(
197  src = 'ak4PFJetsPuppi',
198  jetCorrLabel = 'ak4PFCHSL2L3Corrector',
199  )
200  del process.corrPfMetType1Puppi.offsetCorrLabel # no L1 for PUPPI jets
201  process.pfMetT1Puppi = process.pfMetT1.clone(
202  src = 'pfMetPuppi',
203  srcCorrections = [ cms.InputTag("corrPfMetType1Puppi","type1") ]
204  )
205 
206  from PhysicsTools.PatAlgos.tools.metTools import addMETCollection
207  addMETCollection(process, labelName='patMETPuppi', metSource='pfMetT1Puppi') # T1
208  addMETCollection(process, labelName='patPFMetPuppi', metSource='pfMetPuppi') # RAW
209 
210  process.load('PhysicsTools.PatAlgos.slimming.slimmedMETs_cfi')
211  process.slimmedMETsPuppi = process.slimmedMETs.clone()
212  process.slimmedMETsPuppi.src = cms.InputTag("patMETPuppi")
213  process.slimmedMETsPuppi.rawUncertainties = cms.InputTag("patPFMetPuppi") # only central value
214  process.slimmedMETsPuppi.type1Uncertainties = cms.InputTag("patPFMetT1") # only central value for now
215  del process.slimmedMETsPuppi.type1p2Uncertainties # not available
216 
vid_id_tools.setupAllVIDIdsInModule
def setupAllVIDIdsInModule
Definition: vid_id_tools.py:31
vid_id_tools.switchOnVIDElectronIdProducer
def switchOnVIDElectronIdProducer
Electrons.
Definition: vid_id_tools.py:45
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 243 of file miniAOD_tools.py.

Referenced by miniAOD_customizeAllData().

243 
244 def miniAOD_customizeData(process):
245  from PhysicsTools.PatAlgos.tools.coreTools import runOnData
246  runOnData( process, outputModules = [] )
miniAOD_tools.miniAOD_customizeData
def miniAOD_customizeData
Definition: miniAOD_tools.py:243
def miniAOD_tools.miniAOD_customizeMC (   process)

Definition at line 217 of file miniAOD_tools.py.

Referenced by miniAOD_customizeAllMC().

217 
218 def miniAOD_customizeMC(process):
219  process.muonMatch.matched = "prunedGenParticles"
220  process.electronMatch.matched = "prunedGenParticles"
221  process.electronMatch.src = cms.InputTag("reducedEgamma","reducedGedGsfElectrons")
222  process.photonMatch.matched = "prunedGenParticles"
223  process.photonMatch.src = cms.InputTag("reducedEgamma","reducedGedPhotons")
224  process.tauMatch.matched = "prunedGenParticles"
225  process.tauGenJets.GenParticles = "prunedGenParticles"
226  process.patJetPartonMatch.matched = "prunedGenParticles"
227  process.patJetPartonMatch.mcStatus = [ 3, 23 ]
228  process.patJetGenJetMatch.matched = "slimmedGenJets"
229  process.patJetGenJetMatchAK8.matched = "slimmedGenJetsAK8"
230  process.patMuons.embedGenMatch = False
231  process.patElectrons.embedGenMatch = False
232  process.patPhotons.embedGenMatch = False
233  process.patTaus.embedGenMatch = False
234  process.patJets.embedGenPartonMatch = False
235  #also jet flavour must be switched to ak4
236  process.patJetFlavourAssociation.rParam = 0.4
miniAOD_tools.miniAOD_customizeMC
def miniAOD_customizeMC
Definition: miniAOD_tools.py:217
def miniAOD_tools.miniAOD_customizeOutput (   out)

Definition at line 237 of file miniAOD_tools.py.

237 
238 def miniAOD_customizeOutput(out):
239  out.dropMetaData = cms.untracked.string('ALL')
240  out.fastCloning= cms.untracked.bool(False)
241  out.overrideInputFileSplitLevels = cms.untracked.bool(True)
242  out.compressionAlgorithm = cms.untracked.string('LZMA')
miniAOD_tools.miniAOD_customizeOutput
def miniAOD_customizeOutput
Definition: miniAOD_tools.py:237
Generated for CMSSW Reference Manual by doxygen 1.8.5