def addFlashggPFCHSJets(process,
isData,
vertexIndex = 0,
#doQGTagging = True,
label ='',
debug = False):
setattr(process, 'selectedMuons' + label, cms.EDFilter("CandPtrSelector",
src = cms.InputTag("slimmedMuons"),
cut = cms.string('''abs(eta)<2.5 && pt>10. &&
(pfIsolationR04().sumChargedHadronPt+
max(0.,pfIsolationR04().sumNeutralHadronEt+
pfIsolationR04().sumPhotonEt-
0.50*pfIsolationR04().sumPUPt))/pt < 0.20 &&
(isPFMuon && (isGlobalMuon || isTrackerMuon) )''')))
setattr(process, 'selectedElectrons' + label, cms.EDFilter("CandPtrSelector",
src = cms.InputTag("slimmedElectrons"),
cut = cms.string('''abs(eta)<2.5 && pt>20. &&
gsfTrack.isAvailable() &&
gsfTrack.hitPattern().numberOfLostHits(\'MISSING_INNER_HITS\') < 2 &&
(pfIsolationVariables().sumChargedHadronPt+
max(0.,pfIsolationVariables().sumNeutralHadronEt+
pfIsolationVariables().sumPhotonEt-
0.5*pfIsolationVariables().sumPUPt))/pt < 0.15''')))
# Simple producer which just removes the Candidates which
# don't come from the legacy vertex according to the Flashgg Vertex Map
setattr(process,'flashggCHSLegacyVertexCandidates' + label,
cms.EDProducer('FlashggMultiCHSLegacyVertexCandProducer',
PFCandidatesTag = cms.InputTag('packedPFCandidates'),
DiPhotonTag = cms.InputTag('flashggDiPhotons'),
VertexCandidateMapTag = cms.InputTag("flashggVertexMapForCHS"),
VertexTag = cms.InputTag('offlineSlimmedPrimaryVertices'),
vertexIndex = cms.uint32(vertexIndex),
debug = cms.untracked.bool(debug)
)
)
setattr(process, 'pfCHSLeg' + label, cms.EDFilter("CandPtrSelector",
src = cms.InputTag('flashggCHSLegacyVertexCandidates' + label),
cut = cms.string('')))
# then remove the previously selected muons
setattr(process, 'pfNoMuonCHSLeg' + label, cms.EDProducer("CandPtrProjector",
src = cms.InputTag("pfCHSLeg" + label),
veto = cms.InputTag("selectedMuons" + label)))
# then remove the previously selected electrons
setattr(process, 'pfNoElectronsCHSLeg' + label, cms.EDProducer("CandPtrProjector",
src = cms.InputTag("pfNoMuonCHSLeg" + label),
veto = cms.InputTag("selectedElectrons" + label)))
#Import RECO jet producer for ak4 PF and GEN jet
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
setattr(process, 'ak4PFJetsCHSLeg' + label, ak4PFJets.clone ( src = 'pfNoElectronsCHSLeg' + label, doAreaFastjet = True))
if isData:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute','L2L3Residual']
else:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute']
# NOTE: this is the 74X recipe for the jet clustering
addJetCollection(
process,
postfix = label,
labelName = 'AK4PFCHSLeg',
jetSource = cms.InputTag('ak4PFJetsCHSLeg' + label),
pvSource = cms.InputTag('offlineSlimmedPrimaryVertices'),
pfCandidates = cms.InputTag('packedPFCandidates'),
svSource = cms.InputTag('slimmedSecondaryVertices'),
btagDiscriminators = [ flashggBTag ],
jetCorrections = ('AK4PFchs', JECs, 'None'),
genJetCollection = cms.InputTag('slimmedGenJets'),
genParticles = cms.InputTag('prunedGenParticles'),
# jet param
algo = 'AK', rParam = 0.4
)
#adjust PV used for Jet Corrections
#process.patJetCorrFactorsAK4PFCHSLeg.primaryVertices = "offlineSlimmedPrimaryVertices"
getattr(process, 'patJetCorrFactorsAK4PFCHSLeg' + label).primaryVertices = "offlineSlimmedPrimaryVertices"
if not hasattr(process,"QGPoolDBESSource"):
process.QGPoolDBESSource = cms.ESSource("PoolDBESSource",
CondDBSetup,
toGet = cms.VPSet(),
connect = cms.string('sqlite:QGL_'+qgDatabaseVersion+'.db')
)
process.es_prefer_qg = cms.ESPrefer('PoolDBESSource','QGPoolDBESSource')
for type in ['AK4PFchs']:#,'AK4PFchs_antib']:
process.QGPoolDBESSource.toGet.extend(cms.VPSet(cms.PSet(
record = cms.string('QGLikelihoodRcd'),
tag = cms.string('QGLikelihoodObject_'+qgDatabaseVersion+'_'+type),
label = cms.untracked.string('QGL_'+type)
)))
from RecoJets.JetProducers.QGTagger_cfi import QGTagger
setattr( process, 'QGTaggerPFCHS' + label,
QGTagger.clone( srcJets = 'patJetsAK4PFCHSLeg' + label ,jetsLabel = 'QGL_AK4PFchs', vertexIndex = cms.uint32(vertexIndex),
srcVertexCollection = 'offlineSlimmedPrimaryVertices'))
from RecoJets.JetProducers.PileupJetIDParams_cfi import full_80x_chs
flashggJets = cms.EDProducer('FlashggJetProducer',
DiPhotonTag = cms.InputTag('flashggDiPhotons'),
VertexTag = cms.InputTag('offlineSlimmedPrimaryVertices'),
JetTag = cms.InputTag('patJetsAK4PFCHSLeg' + label),
VertexCandidateMapTag = cms.InputTag("flashggVertexMapForCHS"),
qgVariablesInputTag = cms.InputTag('QGTaggerPFCHS'+label, 'qgLikelihood'),
ComputeSimpleRMS = cms.bool(True),
PileupJetIdParameters = full_80x_chs,
rho = cms.InputTag("fixedGridRhoFastjetAll"),
JetCollectionIndex = cms.uint32(vertexIndex),
Debug = cms.untracked.bool(False)
)
setattr( process, 'flashggPFCHSJets'+ label, flashggJets)
# randomize Jets
from flashgg.MicroAOD.flashggRandomizedJetProducer_cfi import flashggRandomizedJets
flashggRandomizedPFCHSJets = flashggRandomizedJets.clone()
flashggRandomizedPFCHSJets.src = "flashggPFCHSJets" + label
setattr(process.RandomNumberGeneratorService, 'flashggRandomizedPFCHSJets' + label, cms.PSet(initialSeed = cms.untracked.uint32(36423784 + int(label))))
setattr( process, 'flashggRandomizedPFCHSJets' + label, flashggRandomizedPFCHSJets )
flashggSelectedJets = cms.EDFilter("FLASHggJetSelector",
src = cms.InputTag( 'flashggRandomizedPFCHSJets' + label ),
cut = cms.string("pt > 15.")
)
setattr( process, 'flashggSelectedPFCHSJets'+label, flashggSelectedJets )
def addFlashggPuppiJets(process,
vertexIndex = 0,
doQGTagging = True,
label ='',
useLocalJEC = True,
dbfile = 'flashgg/MetaData/data/PuppiJEC/PY8_RunIISpring15DR74_bx50_MC.db',
debug = False):
from CommonTools.PileupAlgos.flashggPuppi_cff import flashggPuppi
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
# fill the puppi parameters
setattr(process, 'flashggPuppi' + label,
flashggPuppi.clone( candName = cms.InputTag('packedPFCandidates'),
vertexName = cms.InputTag('offlineSlimmedPrimaryVertices'),
diPhotonTag = cms.InputTag('flashggDiPhotons'),
VertexCandidateMapTag = cms.InputTag('flashggVertexMapForPUPPI'),
vertexIndex = cms.uint32(vertexIndex),
debug = cms.untracked.bool(debug)
)
)
setattr ( process, 'ak4PFJetsPuppi' + label,
ak4PFJets.clone ( src = cms.InputTag('flashggPuppi' + label), doAreaFastjet = True)
)
if useLocalJEC :
print ':: using a local JEC dbfile for PUPPI :',
print '\t -- ', dbfile
from flashgg.MicroAOD.flashggJetTools_cfi import loadLocalJECDBfile
loadLocalJECDBfile(process,
dbfile = os.environ['CMSSW_BASE'] + '/src/' + dbfile,
tag = 'JetCorrectorParametersCollection_PY8_RunIISpring15DR74_bx50_MC_AK4PUPPI',
label = 'AK4PFPuppi')
# do jet clustering
addJetCollection(
process,
postfix = label,
labelName = 'AK4PUPPI',
jetSource = cms.InputTag('ak4PFJetsPuppi' + label),
pvSource = cms.InputTag('offlineSlimmedPrimaryVertices'),
pfCandidates = cms.InputTag('packedPFCandidates'),
svSource = cms.InputTag('slimmedSecondaryVertices'),
btagDiscriminators = [ flashggBTag ],
jetCorrections = ('AK4PFPuppi',['L1FastJet', 'L2Relative', 'L3Absolute'], 'None'),
genJetCollection = cms.InputTag('slimmedGenJets'),
genParticles = cms.InputTag('prunedGenParticles'),
# jet param
algo = 'AK', rParam = 0.4
)
getattr(process, 'patJetCorrFactorsAK4PUPPI' + label).primaryVertices = "offlineSlimmedPrimaryVertices"
setattr( process,'flashggPUPPIJets'+ label,
cms.EDProducer('FlashggJetProducer',
DiPhotonTag = cms.InputTag('flashggDiPhotons'),
VertexTag = cms.InputTag('offlineSlimmedPrimaryVertices'),
JetTag = cms.InputTag('patJetsAK4PUPPI' + label),
VertexCandidateMapTag = cms.InputTag("flashggVertexMapForPUPPI"),
UsePuppi = cms.untracked.bool(True),
ComputeSimpleRMS = cms.bool(True)
# PileupJetIdParameters = cms.PSet(pu_jetid)
))
# randomize Jets
from flashgg.MicroAOD.flashggRandomizedJetProducer_cfi import flashggRandomizedJets
flashggRandomizedPUPPIJets = flashggRandomizedJets.clone()
flashggRandomizedPUPPIJets.src = "flashggPUPPIJets" + label
setattr(process.RandomNumberGeneratorService, 'flashggRandomizedPUPPIJets' + label, cms.PSet(initialSeed = cms.untracked.uint32(36421523 + int(label))))
setattr( process, 'flashggRandomizedPUPPIJets' + label, flashggRandomizedPUPPIJets )
setattr( process, 'selectedFlashggPUPPIJets'+ label,
cms.EDFilter("FLASHggJetSelector",
src = cms.InputTag( 'flashggRandomizedPUPPIJets'+ label ),
cut = cms.string("pt > 15.")
))
def addFlashggPuppiJets(
process,
vertexIndex=0,
doQGTagging=True,
label='',
useLocalJEC=True,
dbfile='flashgg/MetaData/data/PuppiJEC/PY8_RunIISpring15DR74_bx50_MC.db',
debug=False):
from CommonTools.PileupAlgos.flashggPuppi_cff import flashggPuppi
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
# fill the puppi parameters
setattr(
process, 'flashggPuppi' + label,
flashggPuppi.clone(
candName=cms.InputTag('packedPFCandidates'),
vertexName=cms.InputTag('offlineSlimmedPrimaryVertices'),
diPhotonTag=cms.InputTag('flashggDiPhotons'),
VertexCandidateMapTag=cms.InputTag('flashggVertexMapForPUPPI'),
vertexIndex=cms.uint32(vertexIndex),
debug=cms.untracked.bool(debug)))
setattr(
process, 'ak4PFJetsPuppi' + label,
ak4PFJets.clone(src=cms.InputTag('flashggPuppi' + label),
doAreaFastjet=True))
if useLocalJEC:
print ':: using a local JEC dbfile for PUPPI :',
print '\t -- ', dbfile
from flashgg.MicroAOD.flashggJetTools_cfi import loadLocalJECDBfile
loadLocalJECDBfile(
process,
dbfile=os.environ['CMSSW_BASE'] + '/src/' + dbfile,
tag=
'JetCorrectorParametersCollection_PY8_RunIISpring15DR74_bx50_MC_AK4PUPPI',
label='AK4PFPuppi')
# do jet clustering
addJetCollection(
process,
postfix=label,
labelName='AK4PUPPI',
jetSource=cms.InputTag('ak4PFJetsPuppi' + label),
pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'),
pfCandidates=cms.InputTag('packedPFCandidates'),
svSource=cms.InputTag('slimmedSecondaryVertices'),
btagDiscriminators=[flashggBTag],
jetCorrections=('AK4PFPuppi',
['L1FastJet', 'L2Relative', 'L3Absolute'], 'None'),
genJetCollection=cms.InputTag('slimmedGenJets'),
genParticles=cms.InputTag('prunedGenParticles'),
# jet param
algo='AK',
rParam=0.4)
getattr(process, 'patJetCorrFactorsAK4PUPPI' +
label).primaryVertices = "offlineSlimmedPrimaryVertices"
setattr(
process,
'flashggPUPPIJets' + label,
cms.EDProducer(
'FlashggJetProducer',
DiPhotonTag=cms.InputTag('flashggDiPhotons'),
VertexTag=cms.InputTag('offlineSlimmedPrimaryVertices'),
JetTag=cms.InputTag('patJetsAK4PUPPI' + label),
VertexCandidateMapTag=cms.InputTag("flashggVertexMapForPUPPI"),
UsePuppi=cms.untracked.bool(True),
ComputeSimpleRMS=cms.bool(True),
ComputeRegVars=cms.bool(False)
# PileupJetIdParameters = cms.PSet(pu_jetid)
))
# randomize Jets
from flashgg.MicroAOD.flashggRandomizedJetProducer_cfi import flashggRandomizedJets
flashggRandomizedPUPPIJets = flashggRandomizedJets.clone()
flashggRandomizedPUPPIJets.src = "flashggPUPPIJets" + label
setattr(process.RandomNumberGeneratorService,
'flashggRandomizedPUPPIJets' + label,
cms.PSet(initialSeed=cms.untracked.uint32(36421523 + int(label))))
setattr(process, 'flashggRandomizedPUPPIJets' + label,
flashggRandomizedPUPPIJets)
setattr(
process, 'selectedFlashggPUPPIJets' + label,
cms.EDFilter("FLASHggJetSelector",
src=cms.InputTag('flashggRandomizedPUPPIJets' + label),
cut=cms.string("pt > 15.")))
def addFlashggPFCHSJets(
process,
isData,
vertexIndex=0,
#doQGTagging = True,
label='',
debug=False):
setattr(
process, 'selectedMuons' + label,
cms.EDFilter("CandPtrSelector",
src=cms.InputTag("slimmedMuons"),
cut=cms.string('''abs(eta)<2.5 && pt>10. &&
(pfIsolationR04().sumChargedHadronPt+
max(0.,pfIsolationR04().sumNeutralHadronEt+
pfIsolationR04().sumPhotonEt-
0.50*pfIsolationR04().sumPUPt))/pt < 0.20 &&
(isPFMuon && (isGlobalMuon || isTrackerMuon) )'''
)))
setattr(
process, 'selectedElectrons' + label,
cms.EDFilter("CandPtrSelector",
src=cms.InputTag("slimmedElectrons"),
cut=cms.string('''abs(eta)<2.5 && pt>20. &&
gsfTrack.isAvailable() &&
gsfTrack.hitPattern().numberOfLostHits(\'MISSING_INNER_HITS\') < 2 &&
(pfIsolationVariables().sumChargedHadronPt+
max(0.,pfIsolationVariables().sumNeutralHadronEt+
pfIsolationVariables().sumPhotonEt-
0.5*pfIsolationVariables().sumPUPt))/pt < 0.15'''
)))
# Simple producer which just removes the Candidates which
# don't come from the legacy vertex according to the Flashgg Vertex Map
setattr(
process, 'flashggCHSLegacyVertexCandidates' + label,
cms.EDProducer(
'FlashggMultiCHSLegacyVertexCandProducer',
PFCandidatesTag=cms.InputTag('packedPFCandidates'),
DiPhotonTag=cms.InputTag('flashggDiPhotons'),
VertexCandidateMapTag=cms.InputTag("flashggVertexMapForCHS"),
VertexTag=cms.InputTag('offlineSlimmedPrimaryVertices'),
vertexIndex=cms.uint32(vertexIndex),
debug=cms.untracked.bool(debug)))
setattr(
process, 'pfCHSLeg' + label,
cms.EDFilter("CandPtrSelector",
src=cms.InputTag('flashggCHSLegacyVertexCandidates' +
label),
cut=cms.string('')))
# then remove the previously selected muons
setattr(
process, 'pfNoMuonCHSLeg' + label,
cms.EDProducer("CandPtrProjector",
src=cms.InputTag("pfCHSLeg" + label),
veto=cms.InputTag("selectedMuons" + label)))
# then remove the previously selected electrons
setattr(
process, 'pfNoElectronsCHSLeg' + label,
cms.EDProducer("CandPtrProjector",
src=cms.InputTag("pfNoMuonCHSLeg" + label),
veto=cms.InputTag("selectedElectrons" + label)))
#Import RECO jet producer for ak4 PF and GEN jet
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
setattr(
process, 'ak4PFJetsCHSLeg' + label,
ak4PFJets.clone(src='pfNoElectronsCHSLeg' + label, doAreaFastjet=True))
if isData:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual']
else:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute']
# NOTE: this is the 74X recipe for the jet clustering
addJetCollection(
process,
postfix=label,
labelName='AK4PFCHSLeg',
jetSource=cms.InputTag('ak4PFJetsCHSLeg' + label),
pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'),
pfCandidates=cms.InputTag('packedPFCandidates'),
svSource=cms.InputTag('slimmedSecondaryVertices'),
btagDiscriminators=[flashggBTag],
jetCorrections=('AK4PFchs', JECs, 'None'),
genJetCollection=cms.InputTag('slimmedGenJets'),
genParticles=cms.InputTag('prunedGenParticles'),
# jet param
algo='AK',
rParam=0.4,
btagInfos=['pfImpactParameterTagInfos',
'pfSecondaryVertexTagInfos'] #Extra btagging info
)
#Recalculate btagging info
getattr(process, 'patJetsAK4PFCHSLeg' + label).addTagInfos = True
#adjust PV used for Jet Corrections
#process.patJetCorrFactorsAK4PFCHSLeg.primaryVertices = "offlineSlimmedPrimaryVertices"
getattr(process, 'patJetCorrFactorsAK4PFCHSLeg' +
label).primaryVertices = "offlineSlimmedPrimaryVertices"
if not hasattr(process, "QGPoolDBESSource"):
process.QGPoolDBESSource = cms.ESSource(
"PoolDBESSource",
CondDBSetup,
toGet=cms.VPSet(),
connect=cms.string('sqlite:QGL_' + qgDatabaseVersion + '.db'))
process.es_prefer_qg = cms.ESPrefer('PoolDBESSource',
'QGPoolDBESSource')
for type in ['AK4PFchs']: #,'AK4PFchs_antib']:
process.QGPoolDBESSource.toGet.extend(
cms.VPSet(
cms.PSet(record=cms.string('QGLikelihoodRcd'),
tag=cms.string('QGLikelihoodObject_' +
qgDatabaseVersion + '_' + type),
label=cms.untracked.string('QGL_' + type))))
from RecoJets.JetProducers.QGTagger_cfi import QGTagger
setattr(
process, 'QGTaggerPFCHS' + label,
QGTagger.clone(srcJets='patJetsAK4PFCHSLeg' + label,
jetsLabel='QGL_AK4PFchs',
vertexIndex=cms.uint32(vertexIndex),
srcVertexCollection='offlineSlimmedPrimaryVertices'))
from RecoJets.JetProducers.PileupJetIDParams_cfi import full_80x_chs
flashggJets = cms.EDProducer(
'FlashggJetProducer',
DiPhotonTag=cms.InputTag('flashggDiPhotons'),
VertexTag=cms.InputTag('offlineSlimmedPrimaryVertices'),
JetTag=cms.InputTag('patJetsAK4PFCHSLeg' + label),
VertexCandidateMapTag=cms.InputTag("flashggVertexMapForCHS"),
qgVariablesInputTag=cms.InputTag('QGTaggerPFCHS' + label,
'qgLikelihood'),
ComputeSimpleRMS=cms.bool(True),
ComputeRegVars=cms.bool(True),
PileupJetIdParameters=full_80x_chs,
rho=cms.InputTag("fixedGridRhoFastjetAll"),
JetCollectionIndex=cms.uint32(vertexIndex),
Debug=cms.untracked.bool(False))
setattr(process, 'flashggPFCHSJets' + label, flashggJets)
# randomize Jets
from flashgg.MicroAOD.flashggRandomizedJetProducer_cfi import flashggRandomizedJets
flashggRandomizedPFCHSJets = flashggRandomizedJets.clone()
flashggRandomizedPFCHSJets.src = "flashggPFCHSJets" + label
setattr(process.RandomNumberGeneratorService,
'flashggRandomizedPFCHSJets' + label,
cms.PSet(initialSeed=cms.untracked.uint32(36423784 + int(label))))
setattr(process, 'flashggRandomizedPFCHSJets' + label,
flashggRandomizedPFCHSJets)
flashggSelectedJets = cms.EDFilter(
"FLASHggJetSelector",
src=cms.InputTag('flashggRandomizedPFCHSJets' + label),
cut=cms.string("pt > 15."))
setattr(process, 'flashggSelectedPFCHSJets' + label, flashggSelectedJets)
def addFlashggPFCHSJets(process,
isData,
vertexIndex=0,
useZeroVertex=False,
label='',
debug=False):
#if useZeroVertex:
setattr(
process, 'pfCHSLeg' + label,
cms.EDFilter("CandPtrSelector",
src=cms.InputTag("packedPFCandidates"),
cut=cms.string("fromPV(%i)" % vertexIndex)))
#else:
# setattr(process, 'selectedMuons' + label, cms.EDFilter("CandPtrSelector",
# src = cms.InputTag("slimmedMuons"),
# cut = cms.string('''abs(eta)<2.5 && pt>10. &&
# (pfIsolationR04().sumChargedHadronPt+
# max(0.,pfIsolationR04().sumNeutralHadronEt+
# pfIsolationR04().sumPhotonEt-
# 0.50*pfIsolationR04().sumPUPt))/pt < 0.20 &&
# (isPFMuon && (isGlobalMuon || isTrackerMuon) )''')))
# setattr(process, 'selectedElectrons' + label, cms.EDFilter("CandPtrSelector",
# src = cms.InputTag("slimmedElectrons"),
# cut = cms.string('''abs(eta)<2.5 && pt>20. &&
# gsfTrack.isAvailable() &&
# gsfTrack.hitPattern().numberOfLostHits(\'MISSING_INNER_HITS\') < 2 &&
# (pfIsolationVariables().sumChargedHadronPt+
# max(0.,pfIsolationVariables().sumNeutralHadronEt+
# pfIsolationVariables().sumPhotonEt-
# 0.5*pfIsolationVariables().sumPUPt))/pt < 0.15''')))
# # Simple producer which just removes the Candidates which
# # don't come from the legacy vertex according to the Flashgg Vertex Map
# setattr(process,'flashggCHSLegacyVertexCandidates' + label,
# cms.EDProducer('FlashggMultiCHSLegacyVertexCandProducer',
# PFCandidatesTag = cms.InputTag('packedPFCandidates'),
# DiPhotonTag = cms.InputTag('flashggDiPhotons'),
# VertexCandidateMapTag = cms.InputTag("flashggVertexMapForCHS"),
# VertexTag = cms.InputTag('offlineSlimmedPrimaryVertices'),
# vertexIndex = cms.uint32(vertexIndex),
# debug = cms.untracked.bool(debug)
# )
# )
# setattr(process, 'pfCHSLeg' + label, cms.EDFilter("CandPtrSelector",
# src = cms.InputTag('flashggCHSLegacyVertexCandidates' + label),
# cut = cms.string('')))
# # then remove the previously selected muons
# setattr(process, 'pfNoMuonCHSLeg' + label, cms.EDProducer("CandPtrProjector",
# src = cms.InputTag("pfCHSLeg" + label),
# veto = cms.InputTag("selectedMuons" + label)))
# # then remove the previously selected electrons
# setattr(process, 'pfNoElectronsCHSLeg' + label, cms.EDProducer("CandPtrProjector",
# src = cms.InputTag("pfNoMuonCHSLeg" + label),
# veto = cms.InputTag("selectedElectrons" + label)))
#Import RECO jet producer for ak4 PF and GEN jet
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
#if useZeroVertex:
setattr(process, 'ak4PFJetsCHSLeg' + label,
ak4PFJets.clone(src='pfCHSLeg' + label, doAreaFastjet=True))
#else:
# setattr(process, 'ak4PFJetsCHSLeg' + label, ak4PFJets.clone ( src = 'pfNoElectronsCHSLeg' + label, doAreaFastjet = True))
if isData:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual']
else:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute']
# NOTE: this is the 74X recipe for the jet clustering
addJetCollection(
process,
postfix=label,
labelName='AK4PFCHSLeg',
jetSource=cms.InputTag('ak4PFJetsCHSLeg' + label),
pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'),
pfCandidates=cms.InputTag('packedPFCandidates'),
svSource=cms.InputTag('slimmedSecondaryVertices'),
elSource=cms.InputTag("slimmedElectrons"),
muSource=cms.InputTag("slimmedMuons"),
runIVF=True,
btagDiscriminators=[
flashggBTag, flashggCMVABTag, flashggDeepCSVb, flashggDeepCSVbb,
flashggDeepCSVc, flashggDeepCSVudsg
],
jetCorrections=('AK4PFchs', JECs, 'None'),
genJetCollection=cms.InputTag('slimmedGenJets'),
genParticles=cms.InputTag('prunedGenParticles'),
# jet param
algo='AK',
rParam=0.4)
#Recalculate btagging info
getattr(process, 'patJetsAK4PFCHSLeg' + label).addTagInfos = True
#adjust PV used for Jet Corrections
getattr(process, 'patJetCorrFactorsAK4PFCHSLeg' +
label).primaryVertices = "offlineSlimmedPrimaryVertices"
# Add MiniAOD jet with Deep flavour for matching with flashgg jets based on 0th vertex
from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection
updateJetCollection(
process,
jetSource=cms.InputTag("slimmedJets"),
pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'),
svSource=cms.InputTag('slimmedSecondaryVertices'),
btagDiscriminators=[
flashggDeepFlavourb,
flashggDeepFlavourbb,
flashggDeepFlavourc,
flashggDeepFlavouruds,
flashggDeepFlavourlepb,
flashggDeepFlavourg,
],
jetCorrections=('AK4PFchs', JECs, 'None'),
)
miniaodJets = "selectedUpdatedPatJets"
# Add new gq tagger database
process.QGPoolDBESSource = cms.ESSource(
"PoolDBESSource",
CondDBSetup,
toGet=cms.VPSet(),
connect=cms.string(
'sqlite_fip:flashggPlugins/flashggAnalysisNtuplizer/data/QGL_' +
qgDatabaseVersion + '.db'))
process.es_prefer_qg = cms.ESPrefer('PoolDBESSource', 'QGPoolDBESSource')
process.QGPoolDBESSource.toGet.extend(
cms.VPSet(
cms.PSet(record=cms.string('QGLikelihoodRcd'),
tag=cms.string('QGLikelihoodObject_v1_AK4'),
label=cms.untracked.string('QGL_AK4PFchs'))))
from RecoJets.JetProducers.QGTagger_cfi import QGTagger
setattr(
process, 'QGTaggerPFCHS' + label,
QGTagger.clone(srcJets='patJetsAK4PFCHSLeg' + label,
jetsLabel='QGL_AK4PFchs',
vertexIndex=cms.uint32(vertexIndex),
srcVertexCollection='offlineSlimmedPrimaryVertices'))
from RecoJets.JetProducers.PileupJetIDParams_cfi import full_81x_chs
flashggJets = cms.EDProducer(
'FlashggJetProducer',
DiPhotonTag=cms.InputTag('flashggDiPhotons'),
VertexTag=cms.InputTag('offlineSlimmedPrimaryVertices'),
JetTag=cms.InputTag('patJetsAK4PFCHSLeg' + label),
VertexCandidateMapTag=cms.InputTag("flashggVertexMapForCHS"),
qgVariablesInputTag=cms.InputTag('QGTaggerPFCHS' + label,
'qgLikelihood'),
ComputeSimpleRMS=cms.bool(True),
PileupJetIdParameters=full_81x_chs,
rho=cms.InputTag("fixedGridRhoFastjetAll"),
JetCollectionIndex=cms.uint32(vertexIndex),
Debug=cms.untracked.bool(debug),
DoPuJetID=cms.bool(True),
ComputeRegVars=cms.bool(True),
MinPtForEneSum=cms.double(0.),
MaxEtaForEneSum=cms.double(2.5),
NJetsForEneSum=cms.uint32(0),
MiniAodJetTag=cms.InputTag(miniaodJets))
setattr(process, 'flashggPFCHSJets' + label, flashggJets)
# randomize Jets
from flashgg.MicroAOD.flashggRandomizedJetProducer_cfi import flashggRandomizedJets
flashggRandomizedPFCHSJets = flashggRandomizedJets.clone()
flashggRandomizedPFCHSJets.src = "flashggPFCHSJets" + label
setattr(process.RandomNumberGeneratorService,
'flashggRandomizedPFCHSJets' + label,
cms.PSet(initialSeed=cms.untracked.uint32(36423784 + int(label))))
setattr(process, 'flashggRandomizedPFCHSJets' + label,
flashggRandomizedPFCHSJets)
flashggSelectedJets = cms.EDFilter(
"FLASHggJetSelector",
src=cms.InputTag('flashggRandomizedPFCHSJets' + label),
cut=cms.string("pt > 15."))
setattr(process, 'flashggSelectedPFCHSJets' + label, flashggSelectedJets)
def addFlashggPFCHSJets(
process,
DeepJet,
isData,
vertexIndex=0,
#doQGTagging = True,
label='',
debug=False):
setattr(
process, 'pfCHSLeg' + label,
cms.EDFilter("CandPtrSelector",
src=cms.InputTag("packedPFCandidates"),
cut=cms.string("fromPV(%i)" % vertexIndex)))
#Import RECO jet producer for ak4 PF and GEN jet
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
setattr(process, 'ak4PFJetsCHSLeg' + label,
ak4PFJets.clone(src='pfCHSLeg' + label, doAreaFastjet=True))
if isData:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual']
else:
JECs = ['L1FastJet', 'L2Relative', 'L3Absolute']
# NOTE: this is the 74X recipe for the jet clustering - has been confirmed to give reasonable results in 10X
addJetCollection(
process,
postfix=label,
labelName='AK4PFCHSLeg',
jetSource=cms.InputTag('ak4PFJetsCHSLeg' + label),
pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'),
pfCandidates=cms.InputTag('packedPFCandidates'),
svSource=cms.InputTag('slimmedSecondaryVertices'),
elSource=cms.InputTag("slimmedElectrons"),
muSource=cms.InputTag("slimmedMuons"),
runIVF=True,
btagDiscriminators=[
flashggBTag, flashggCMVABTag, flashggDeepCSVb, flashggDeepCSVbb,
flashggDeepCSVc, flashggDeepCSVudsg
],
jetCorrections=('AK4PFchs', JECs, 'None'),
genJetCollection=cms.InputTag('slimmedGenJets'),
genParticles=cms.InputTag('prunedGenParticles'),
# jet param
algo='AK',
rParam=0.4
#btagInfos = ['pfImpactParameterTagInfos', 'pfSecondaryVertexTagInfos', 'pfInclusiveSecondaryVertexFinderTagInfos', 'softPFMuonsTagInfos', 'softPFElectronsTagInfos'] #Extra Btagging Info
)
#Recalculate btagging info
getattr(process, 'patJetsAK4PFCHSLeg' + label).addTagInfos = True
#adjust PV used for Jet Corrections
#process.patJetCorrFactorsAK4PFCHSLeg.primaryVertices = "offlineSlimmedPrimaryVertices"
getattr(process, 'patJetCorrFactorsAK4PFCHSLeg' +
label).primaryVertices = "offlineSlimmedPrimaryVertices"
if DeepJet == "rerun":
from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection
updateJetCollection(
process,
jetSource=cms.InputTag("slimmedJets"),
pvSource=cms.InputTag('offlineSlimmedPrimaryVertices'),
svSource=cms.InputTag('slimmedSecondaryVertices'),
btagDiscriminators=[
flashggDeepFlavourb,
flashggDeepFlavourbb,
flashggDeepFlavourc,
flashggDeepFlavouruds,
flashggDeepFlavourlepb,
flashggDeepFlavourg,
],
jetCorrections=('AK4PFchs', JECs, 'None'),
)
miniaodJets = "selectedUpdatedPatJets"
else:
miniaodJets = "slimmedJets"
if not hasattr(process, "QGPoolDBESSource"):
process.QGPoolDBESSource = cms.ESSource(
"PoolDBESSource",
#CondDB,
toGet=cms.VPSet(),
connect=cms.string('sqlite:QGL_' + qgDatabaseVersion + '.db'))
process.es_prefer_qg = cms.ESPrefer('PoolDBESSource',
'QGPoolDBESSource')
for type in ['AK4PFchs']: #,'AK4PFchs_antib']:
process.QGPoolDBESSource.toGet.extend(
cms.VPSet(
cms.PSet(record=cms.string('QGLikelihoodRcd'),
tag=cms.string('QGLikelihoodObject_v1_AK4'),
label=cms.untracked.string('QGL_' + type))))
from RecoJets.JetProducers.QGTagger_cfi import QGTagger
setattr(
process, 'QGTaggerPFCHS' + label,
QGTagger.clone(srcJets='patJetsAK4PFCHSLeg' + label,
jetsLabel='QGL_AK4PFchs',
vertexIndex=cms.uint32(vertexIndex),
srcVertexCollection='offlineSlimmedPrimaryVertices'))
from RecoJets.JetProducers.PileupJetIDParams_cfi import full_81x_chs
pujidparam = full_81x_chs
flashggJets = cms.EDProducer(
'FlashggJetProducer',
DiPhotonTag=cms.InputTag('flashggDiPhotons'),
VertexTag=cms.InputTag('offlineSlimmedPrimaryVertices'),
JetTag=cms.InputTag('patJetsAK4PFCHSLeg' + label),
VertexCandidateMapTag=cms.InputTag("flashggVertexMapForCHS"),
qgVariablesInputTag=cms.InputTag('QGTaggerPFCHS' + label,
'qgLikelihood'),
ComputeSimpleRMS=cms.bool(True),
PileupJetIdParameters=pujidparam,
rho=cms.InputTag("fixedGridRhoFastjetAll"),
JetCollectionIndex=cms.uint32(vertexIndex),
Debug=cms.untracked.bool(debug),
DoPuJetID=cms.bool(True),
ComputeRegVars=cms.bool(True),
MinPtForEneSum=cms.double(0.),
MaxEtaForEneSum=cms.double(2.5),
NJetsForEneSum=cms.uint32(0),
MiniAodJetTag=cms.InputTag(miniaodJets))
setattr(process, 'flashggPFCHSJets' + label, flashggJets)
# randomize Jets
from flashgg.MicroAOD.flashggRandomizedJetProducer_cfi import flashggRandomizedJets
flashggRandomizedPFCHSJets = flashggRandomizedJets.clone()
flashggRandomizedPFCHSJets.src = "flashggPFCHSJets" + label
setattr(process.RandomNumberGeneratorService,
'flashggRandomizedPFCHSJets' + label,
cms.PSet(initialSeed=cms.untracked.uint32(36423784 + int(label))))
setattr(process, 'flashggRandomizedPFCHSJets' + label,
flashggRandomizedPFCHSJets)
flashggSelectedJets = cms.EDFilter(
"FLASHggJetSelector",
src=cms.InputTag('flashggRandomizedPFCHSJets' + label),
cut=cms.string("pt > 15."))
setattr(process, 'flashggSelectedPFCHSJets' + label, flashggSelectedJets)
