Skip to content

Commit

Permalink
fix savejetcollection
Browse files Browse the repository at this point in the history
  • Loading branch information
alefisico committed Jan 17, 2019
1 parent 3901344 commit 5a63ff3
Show file tree
Hide file tree
Showing 2 changed files with 21 additions and 24 deletions.
29 changes: 13 additions & 16 deletions python/jetToolbox_cff.py
Original file line number Diff line number Diff line change
Expand Up @@ -52,7 +52,6 @@ def jetToolbox( proc, jetType, jetSequence, outputFile,
addQGTagger=False, QGjetsLabel='chs',
addEnergyCorrFunc=False, ecfType = "N", ecfBeta = 1.0, ecfN3 = False,
addEnergyCorrFuncSubjets=False, ecfSubjetType = "N", ecfSubjetBeta = 1.0, ecfSubjetN3 = False,
saveJetCollection=False, # set this to true to enable creation of edm root file
verbosity=2, # 0 = no printouts, 1 = warnings only, 2 = warnings & info, 3 = warnings, info, debug
):

Expand Down Expand Up @@ -1204,21 +1203,19 @@ def jetToolbox( proc, jetType, jetSequence, outputFile,
#################################################################################
###### Adding to outputModule OR creating output file
setattr(proc, jetSequence, jetSeq)
if outputFile!='':
if hasattr(proc, outputFile): getattr(proc, outputFile).outputCommands += elemToKeep
else: setattr( proc, outputFile,
cms.OutputModule('PoolOutputModule',
fileName = cms.untracked.string('jettoolbox.root'),
outputCommands = cms.untracked.vstring( elemToKeep ) ) )

#################################################################################
##### fix to replace unschedule mode
if saveJetCollection:
task = getPatAlgosToolsTask(proc)
if hasattr(proc, 'endpath'):
getattr(proc, 'endpath').associate(task)
else:
setattr( proc, 'endpath', cms.EndPath(getattr(proc, outputFile), task) )
if hasattr(proc, outputFile): getattr(proc, outputFile).outputCommands += elemToKeep
else: setattr( proc, outputFile,
cms.OutputModule('PoolOutputModule',
fileName = cms.untracked.string('jettoolbox.root'),
outputCommands = cms.untracked.vstring( elemToKeep ) ) )

task = getPatAlgosToolsTask(proc)
if hasattr(proc, 'endpath'):
getattr(proc, 'endpath').associate(task)
else:
if outputFile=='noOutput':
setattr( proc, 'endpath', cms.EndPath(task) )
else: setattr( proc, 'endpath', cms.EndPath(getattr(proc, outputFile), task) )

#################################################################################
#### removing mc matching for data
Expand Down
16 changes: 8 additions & 8 deletions test/ClusterWithToolboxAndMakeHistos.py
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,7 @@
from JMEAnalysis.JetToolbox.jetToolbox_cff import *

# AK R=0.4 jets from CHS inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'ak4', 'ak4JetSubs', 'out',
jetToolbox( process, 'ak4', 'ak4JetSubs', 'noOutput',
PUMethod='CHS',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand All @@ -49,7 +49,7 @@
)

# AK R=0.8 jets from PF inputs with basic grooming, W tagging, and top tagging
#jetToolbox( process, 'ak8', 'ak8JetSubs', 'out',
#jetToolbox( process, 'ak8', 'ak8JetSubs', 'noOutput',
# PUMethod='Plain',
# addPruning=True, addSoftDrop=True , # add basic grooming
# addTrimming=True, addFiltering=True,
Expand All @@ -59,7 +59,7 @@
#)

# AK R=0.8 jets from CHS inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'ak8', 'ak8JetSubs', 'out',
jetToolbox( process, 'ak8', 'ak8JetSubs', 'noOutput',
PUMethod='CHS',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand All @@ -69,7 +69,7 @@
)

# AK R=0.8 from PUPPI inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'ak8', 'ak8JetSubs', 'out',
jetToolbox( process, 'ak8', 'ak8JetSubs', 'noOutput',
PUMethod='Puppi',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand All @@ -79,7 +79,7 @@
)

# CA R=0.8 jets from CHS inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'ca8', 'ca8JetSubs', 'out',
jetToolbox( process, 'ca8', 'ca8JetSubs', 'noOutput',
PUMethod='CHS',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand All @@ -89,7 +89,7 @@
)

# KT R=0.8 jets from CHS inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'kt8', 'kt8JetSubs', 'out',
jetToolbox( process, 'kt8', 'kt8JetSubs', 'noOutput',
PUMethod='CHS',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand All @@ -99,7 +99,7 @@
)

# AK R=1.2 jets from CHS inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'ak12', 'ak12JetSubs', 'out',
jetToolbox( process, 'ak12', 'ak12JetSubs', 'noOutput',
PUMethod='CHS',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand All @@ -109,7 +109,7 @@
)

# AK R=1.5 jets from CHS inputs with basic grooming, W tagging, and top tagging
jetToolbox( process, 'ak15', 'ak15JetSubs', 'out',
jNoSaveetToolbox( process, 'ak15', 'ak15JetSubs', 'noOutput',
PUMethod='CHS',
addPruning=True, addSoftDrop=True , # add basic grooming
addTrimming=True, addFiltering=True,
Expand Down

0 comments on commit 5a63ff3

Please sign in to comment.