| 1 |
import operator,sys,os
|
| 2 |
from FWCore.ParameterSet import Config as cms
|
| 3 |
from PhysicsTools.PatAlgos.tools import coreTools,pfTools
|
| 4 |
|
| 5 |
def tags(stuff) :
|
| 6 |
return ( cms.InputTag(stuff) if type(stuff)!=list else
|
| 7 |
cms.VInputTag([tags(item) for item in stuff]) )
|
| 8 |
|
| 9 |
class TopRefPF2PAT(object) :
|
| 10 |
'''Implement the Top Reference configuration of PF2PAT.
|
| 11 |
|
| 12 |
https://twiki.cern.ch/twiki/bin/viewauth/CMS/TWikiTopRefEventSel
|
| 13 |
Modify the PF selections to match the PAT selections, for consistent TopProjection cross-cleaning.
|
| 14 |
'''
|
| 15 |
|
| 16 |
def __init__(self, process, options) :
|
| 17 |
self.stdout = sys.stdout
|
| 18 |
sys.stdout = open(os.devnull, 'w')
|
| 19 |
print >>self.stdout, "usePF2PAT() output suppressed. (%s)"%__file__
|
| 20 |
|
| 21 |
self.before = set(dir(process))
|
| 22 |
process.load( "PhysicsTools.PatAlgos.patSequences_cff" )
|
| 23 |
pfTools.usePF2PAT( process,
|
| 24 |
runPF2PAT = True,
|
| 25 |
postfix = options.postfix,
|
| 26 |
runOnMC = not options.isData,
|
| 27 |
pvCollection = tags( 'goodOfflinePrimaryVertices' ),
|
| 28 |
typeIMetCorrections = True,
|
| 29 |
jetAlgo = 'AK5',
|
| 30 |
jetCorrections = ('AK5PFchs', ['L1FastJet','L2Relative','L3Absolute','L2L3Residual'][:None if options.isData else -1] )
|
| 31 |
)
|
| 32 |
|
| 33 |
getattr( process, 'pfPileUp'+options.postfix).checkClosestZVertex = False # https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookJetEnergyCorrections#JetEnCorPFnoPU2012
|
| 34 |
getattr( process, 'pfPileUpIso'+options.postfix).checkClosestZVertex = True
|
| 35 |
getattr( process, 'patJetCorrFactors'+options.postfix).flavorType = 'T' # top events, not dijets, but only relevant for L5Flavor,L7Parton
|
| 36 |
#getattr( process, 'patJetCorrFactors'+options.postfix).levels.append('L5Flavor') # apparently not available in GlobalTags
|
| 37 |
if options.isData: coreTools.runOnData( process, names = [ 'PFAll' ], postfix = options.postfix )
|
| 38 |
coreTools.removeSpecificPATObjects( process, names = [ 'Photons', 'Taus' ], postfix = options.postfix )
|
| 39 |
|
| 40 |
self.isoValues = {'el':0.15,'mu':0.20}
|
| 41 |
self.eleID = 'mvaTrigV0'
|
| 42 |
self.minEleID = 0.
|
| 43 |
self.isoEA = 'elPFIsoValueEA03'
|
| 44 |
self.doElectronEA = True
|
| 45 |
self.dBFactorEl = -1.0 if self.doElectronEA else -0.5
|
| 46 |
self.dBFactorMu = -0.5
|
| 47 |
self.cuts = {'el': ['abs(eta)<2.5',
|
| 48 |
'pt>20.',
|
| 49 |
'gsfTrackRef.isNonnull',
|
| 50 |
'gsfTrackRef.trackerExpectedHitsInner.numberOfLostHits<2',
|
| 51 |
'(chargedHadronIso+max(0.,neutralHadronIso+photonIso%+.1f*%s))/pt < %.2f'%(self.dBFactorEl, 'userIsolation("User1Iso")' if self.doElectronEA else 'puChargedHadronIso',self.isoValues['el']),
|
| 52 |
'electronID("%s") > %.2f'%(self.eleID,self.minEleID),
|
| 53 |
],
|
| 54 |
'mu' :['abs(eta)<2.5',
|
| 55 |
'pt>10.',
|
| 56 |
'(chargedHadronIso+neutralHadronIso+photonIso%+.2f*puChargedHadronIso)/pt < %.2f'%(self.dBFactorMu, self.isoValues['mu']),
|
| 57 |
'(isPFMuon && (isGlobalMuon || isTrackerMuon) )',
|
| 58 |
],
|
| 59 |
'jet' : ['abs(eta)<2.5', # Careful! these jet cuts affect the typeI met corrections
|
| 60 |
'pt > 15.',
|
| 61 |
# PF jet ID:
|
| 62 |
'numberOfDaughters > 1',
|
| 63 |
'neutralHadronEnergyFraction < 0.99',
|
| 64 |
'neutralEmEnergyFraction < 0.99',
|
| 65 |
'(chargedEmEnergyFraction < 0.99 || abs(eta) >= 2.4)',
|
| 66 |
'(chargedHadronEnergyFraction > 0. || abs(eta) >= 2.4)',
|
| 67 |
'(chargedMultiplicity > 0 || abs(eta) >= 2.4)']
|
| 68 |
}
|
| 69 |
|
| 70 |
self.process = process
|
| 71 |
self.options = options
|
| 72 |
self.fix = options.postfix
|
| 73 |
self.patSeq = getattr(process, 'patPF2PATSequence' + self.fix)
|
| 74 |
self.btags = self.options.btags
|
| 75 |
self.taginfos = ["impactParameter","secondaryVertex"]
|
| 76 |
|
| 77 |
self.configTopProjections()
|
| 78 |
self.configPfLepton('el', iso3 = True)
|
| 79 |
self.configPfLepton('mu', iso3 = False)
|
| 80 |
self.configPatMuons()
|
| 81 |
self.configPatElectrons()
|
| 82 |
self.configPatJets()
|
| 83 |
self.configLeptonFilter()
|
| 84 |
self.removeCruft()
|
| 85 |
sys.stdout = self.stdout
|
| 86 |
|
| 87 |
def attr(self, item) : return getattr(self.process, item)
|
| 88 |
def show(self, item) : print >> (sys.stdout if self.options.quiet else self.stdout), item, '=', self.attr(item).dumpPython() if hasattr(self.process, item) else 'Not Found.'
|
| 89 |
|
| 90 |
def configTopProjections(self) :
|
| 91 |
'''Enable TopProjections except for Taus (do not remove tau cands from jet collection).'''
|
| 92 |
|
| 93 |
for key,val in {'PileUp':1,'Muon':1,'Electron':1,'Jet':1,'Tau':0}.items() :
|
| 94 |
self.attr('pfNo'+key+self.fix).enable = bool(val)
|
| 95 |
self.attr('pfNoElectron').verbose = True
|
| 96 |
return
|
| 97 |
|
| 98 |
def configPfLepton(self, lep, iso3) :
|
| 99 |
full = {'el':'Electrons','mu':'Muons'}[lep] + self.fix
|
| 100 |
iso = self.attr('pfIsolated'+full)
|
| 101 |
iso.isolationCut = self.isoValues[lep]
|
| 102 |
iso.doDeltaBetaCorrection = True
|
| 103 |
iso.deltaBetaFactor = getattr(self, 'dBFactor' + full[:2])
|
| 104 |
sel = self.attr('pfSelected'+full)
|
| 105 |
sel.cut = ' && '.join(self.cuts[lep][:-2])
|
| 106 |
isoEA = getattr(self.process, self.isoEA)
|
| 107 |
print >>self.stdout, ""
|
| 108 |
if lep == 'el' :
|
| 109 |
idName = 'pfIdentifiedElectrons'+self.fix
|
| 110 |
id = cms.EDFilter("ElectronIDPFCandidateSelector",
|
| 111 |
src = cms.InputTag('pfElectronsFromVertex'+self.fix),
|
| 112 |
recoGsfElectrons = cms.InputTag("gsfElectrons"),
|
| 113 |
electronIdMap = cms.InputTag(self.eleID),
|
| 114 |
electronIdCut = cms.double(self.minEleID))
|
| 115 |
setattr(self.process, idName, id )
|
| 116 |
self.patSeq.replace( sel, id*sel*isoEA)
|
| 117 |
sel.src = idName
|
| 118 |
if not iso3 : isoEA.EffectiveAreaType = "kEleGammaAndNeutralHadronIso04"
|
| 119 |
self.show(idName)
|
| 120 |
|
| 121 |
if iso3:
|
| 122 |
pf = {'el':'PFId','mu':''}[lep] + self.fix
|
| 123 |
for item in ['pf','pfIsolated'] :
|
| 124 |
col = self.attr( item+full )
|
| 125 |
col.deltaBetaIsolationValueMap = tags( lep+'PFIsoValuePU03'+pf )
|
| 126 |
col.isolationValueMapsCharged = tags( [lep+'PFIsoValueCharged03'+pf] )
|
| 127 |
col.isolationValueMapsNeutral = tags( [ lep+'PFIsoValueNeutral03'+pf, lep+'PFIsoValueGamma03'+pf] )
|
| 128 |
val = self.attr( 'pat'+full ).isolationValues
|
| 129 |
val.pfNeutralHadrons = tags( lep+'PFIsoValueNeutral03' + pf )
|
| 130 |
val.pfChargedAll = tags( lep+'PFIsoValueChargedAll03' + pf )
|
| 131 |
val.pfPUChargedHadrons = tags( lep+'PFIsoValuePU03' + pf )
|
| 132 |
val.pfPhotons = tags( lep+'PFIsoValueGamma03' + pf )
|
| 133 |
val.pfChargedHadrons = tags( lep+'PFIsoValueCharged03' + pf )
|
| 134 |
|
| 135 |
if lep == 'el' and self.doElectronEA :
|
| 136 |
self.attr( 'pfIsolated' + full).deltaBetaIsolationValueMap = self.isoEA
|
| 137 |
|
| 138 |
self.show('pfSelected'+full)
|
| 139 |
self.show('pfIsolated'+full)
|
| 140 |
return
|
| 141 |
|
| 142 |
def configPatMuons(self) :
|
| 143 |
for mod,attr,val in [('patMuons','usePV',False), # use beam spot rather than PV, which is necessary for 'dB' cut
|
| 144 |
('patMuons','embedTrack',True), # embedded track needed for muon ID cuts
|
| 145 |
('selectedPatMuons','cut', ' && '.join(self.cuts['mu'])),
|
| 146 |
] : setattr( self.attr(mod+self.fix), attr, val )
|
| 147 |
self.show('selectedPatMuons'+self.fix)
|
| 148 |
return
|
| 149 |
|
| 150 |
def configPatElectrons(self) :
|
| 151 |
electronIDSources = cms.PSet(**dict([(i,tags(i)) for i in [self.eleID]]))
|
| 152 |
for mod,attr,val in [('patElectrons','electronIDSources', electronIDSources),
|
| 153 |
('selectedPatElectrons','cut', ' && '.join(c for c in self.cuts['el'] if 'gsfTrackRef' not in c)),
|
| 154 |
] : setattr( self.attr( mod+self.fix), attr, val )
|
| 155 |
self.attr('patElectrons'+self.fix).isolationValues.user = tags( [self.isoEA] )
|
| 156 |
|
| 157 |
self.show('selectedPatElectrons'+self.fix)
|
| 158 |
return
|
| 159 |
|
| 160 |
def configPatJets(self) :
|
| 161 |
for mod,attr,val in [('patJets','discriminatorSources',[tags(tag+'BJetTagsAOD'+self.fix) for tag in self.btags]),
|
| 162 |
('patJets','tagInfoSources',[tags(tag+'TagInfosAOD'+self.fix) for tag in self.taginfos]),
|
| 163 |
('patJets','addTagInfos',True),
|
| 164 |
('selectedPatJets','cut', ' && '.join(self.cuts['jet'])),
|
| 165 |
] : setattr( self.attr(mod+self.fix), attr, val )
|
| 166 |
self.show('selectedPatJets'+self.fix)
|
| 167 |
return
|
| 168 |
|
| 169 |
def configLeptonFilter(self) :
|
| 170 |
if not self.options.requireLepton : return
|
| 171 |
muons20Name = 'pfIsolatedMuons20'+self.fix
|
| 172 |
leptonsName = 'pfLeptons'+self.fix
|
| 173 |
requireLeptonName = 'requireLepton'+self.fix
|
| 174 |
|
| 175 |
muons20 = cms.EDFilter("GenericPFCandidateSelector", src = tags('pfIsolatedMuons'+self.fix), cut = cms.string('pt>20 && abs(eta) < 2.4'))
|
| 176 |
leptons = cms.EDProducer("CandViewMerger", src = tags(['pfIsolated%s%s'%(lep,self.fix) for lep in ['Electrons','Muons20']]))
|
| 177 |
requireLepton = cms.EDFilter("CandViewCountFilter", src = tags(leptonsName), minNumber = cms.uint32(1) )
|
| 178 |
|
| 179 |
setattr(self.process, muons20Name, muons20)
|
| 180 |
setattr(self.process, leptonsName, leptons)
|
| 181 |
setattr(self.process, requireLeptonName, requireLepton)
|
| 182 |
|
| 183 |
jets = self.attr('pfJets'+self.fix)
|
| 184 |
self.patSeq.replace(jets, muons20*leptons*requireLepton*jets)
|
| 185 |
|
| 186 |
self.show(muons20Name)
|
| 187 |
self.show(leptonsName)
|
| 188 |
self.show(requireLeptonName)
|
| 189 |
return
|
| 190 |
|
| 191 |
def removeCruft(self) :
|
| 192 |
nothanks = [mod for mod in set(str(self.patSeq).split('+'))
|
| 193 |
if ( any( keyword in mod for keyword in ['count','Legacy','pfJetsPiZeros','ak7','soft','iterativeCone',
|
| 194 |
'tauIsoDeposit','tauGenJet','hpsPFTau','tauMatch','pfTau','hpsSelection',
|
| 195 |
'photonMatch','phPFIso','pfIsolatedPhotons','pfCandsNotInJet','pfCandMETcorr',
|
| 196 |
'Negative','ToVertex','particleFlowDisplacedVertex',
|
| 197 |
'ype2Corr','ype0','patType1p2CorrectedPFMet','atCandidateSummaryTR','atPFParticles'] ) or
|
| 198 |
( 'JetTagsAOD' in mod and not any( (btag+'B') in mod for btag in self.btags ) ) or
|
| 199 |
( 'TagInfosAOD' in mod and not any( (info+'TagInfosAOD') in mod for info in self.taginfos ) ) or
|
| 200 |
( 'elPFIsoValue' in mod and ('04' in mod or 'NoPFId' in mod ) ) or
|
| 201 |
( 'muPFIsoValue' in mod and ('03' in mod))) ]
|
| 202 |
for mod in nothanks : self.patSeq.remove( self.attr(mod) )
|
| 203 |
|
| 204 |
nothanks = ['ic5','kt4','kt6','JPT','ak7','ak5Calo']
|
| 205 |
now = set(dir(self.process))
|
| 206 |
for item in now - self.before - set(str(self.patSeq).split('+')) :
|
| 207 |
if any(part in item for part in nothanks) :
|
| 208 |
delattr(self.process,item)
|
| 209 |
return
|
