generator = cms.EDFilter(
"Pythia8GeneratorFilter",
pythiaHepMCVerbosity=cms.untracked.bool(False),
maxEventsToPrint=cms.untracked.int32(0),
pythiaPylistVerbosity=cms.untracked.int32(0),
filterEfficiency=cms.untracked.double(1.38e-3),
crossSection=cms.untracked.double(540000000.),
comEnergy=cms.double(14000.0),
ExternalDecays=cms.PSet(EvtGen130=cms.untracked.PSet(
decay_table=cms.string(
'GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
particle_property_file=cms.FileInPath(
'GeneratorInterface/EvtGenInterface/data/evt.pdl'),
user_decay_embedded=cms.vstring(
'#', 'Alias MyBs B_s0', 'Alias Myanti-Bs anti-B_s0',
'ChargeConj Myanti-Bs MyBs', '#', 'Decay MyBs',
'1.000 e+ e- PHOTOS PHSP;', 'Enddecay',
'#', 'Decay Myanti-Bs',
'1.000 e+ e- PHOTOS PHSP;', 'Enddecay',
'End'),
list_forced_decays=cms.vstring('MyBs', 'Myanti-Bs'),
operates_on_particles=cms.vint32(),
),
parameterSets=cms.vstring('EvtGen130')),
PythiaParameters=cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters=cms.vstring(
"SoftQCD:nonDiffractive = on",
'PTFilter:filter = on', # this turn on the filter
'PTFilter:quarkToFilter = 5', # PDG id of q quark
'PTFilter:scaleToFilter = 1.0'),
parameterSets=cms.vstring(
'pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters',
)))
process.generator = cms.EDFilter(
"Pythia8GeneratorFilter",
PythiaParameters=cms.PSet(
parameterSets=cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings', 'processParameters'),
processParameters=cms.vstring(
'SUSY:all = off', 'SUSY:gg2gluinogluino = on',
'SUSY:qqbar2gluinogluino = on', '1000024:isResonance = false',
'1000024:oneChannel = 1 1.0 100 1000022 211',
'1000024:tau0 = 100000.0', 'ParticleDecays:tau0Max = 1000000.0'),
pythia8CP5Settings=cms.vstring(
'Tune:pp 14', 'Tune:ee 7',
'MultipartonInteractions:ecmPow=0.03344', 'PDF:pSet=20',
'MultipartonInteractions:bProfile=2',
'MultipartonInteractions:pT0Ref=1.41',
'MultipartonInteractions:coreRadius=0.7634',
'MultipartonInteractions:coreFraction=0.63',
'ColourReconnection:range=5.176', 'SigmaTotal:zeroAXB=off',
'SpaceShower:alphaSorder=2', 'SpaceShower:alphaSvalue=0.118',
'SigmaProcess:alphaSvalue=0.118', 'SigmaProcess:alphaSorder=2',
'MultipartonInteractions:alphaSvalue=0.118',
'MultipartonInteractions:alphaSorder=2',
'TimeShower:alphaSorder=2', 'TimeShower:alphaSvalue=0.118'),
pythia8CommonSettings=cms.vstring(
'Tune:preferLHAPDF = 2', 'Main:timesAllowErrors = 10000',
'Check:epTolErr = 0.01', 'Beams:setProductionScalesFromLHEF = off',
'SLHA:keepSM = on', 'SLHA:minMassSM = 1000.',
'ParticleDecays:limitTau0 = on', 'ParticleDecays:tau0Max = 10',
'ParticleDecays:allowPhotonRadiation = on')),
SLHATableForPythia8=cms.string(
'\n# ISAJET SUSY parameters in SUSY Les Houches Accord 2 format\n# Created by ISALHA 2.0 Last revision: C. Balazs 21 Apr 2009\nBlock SPINFO # Program information\n 1 ISASUGRA from ISAJET # Spectrum Calculator\n 2 7.80 29-OCT-2009 12:50:36 # Version number\nBlock MODSEL # Model selection\n 1 3 # Minimal anomaly mediated (AMSB) model\nBlock SMINPUTS # Standard Model inputs\n 1 1.27836266E+02 # alpha_em^(-1)\n 2 1.16570000E-05 # G_Fermi\n 3 1.17200002E-01 # alpha_s(M_Z)\n 4 9.11699982E+01 # m_{Z}(pole)\n 5 4.19999981E+00 # m_{b}(m_{b})\n 6 1.73070007E+02 # m_{top}(pole)\n 7 1.77699995E+00 # m_{tau}(pole)\nBlock MINPAR # SUSY breaking input parameters\n 1 1.50000000E+03 # m_0\n 2 2.82870000E+05 # m_{3/2}\n 3 5.00000000E+00 # tan(beta)\n 4 1.00000000E+00 # sign(mu)\nBlock EXTPAR # Non-universal SUSY breaking parameters\n 0 1.00218262E+16 # Input scale\nBlock MASS # Scalar and gaugino mass spectrum\n# PDG code mass particle\n 24 8.04229965E+01 # W^+\n 25 1.17473404E+02 # h^0\n 35 4.63466455E+03 # H^0\n 36 4.60420850E+03 # A^0\n 37 4.62050830E+03 # H^+\n 1000001 5.26110645E+03 # dnl\n 1000002 5.26050488E+03 # upl\n 1000003 5.26110645E+03 # stl\n 1000004 5.26050488E+03 # chl\n 1000005 4.59925488E+03 # b1\n 1000006 3.83768652E+03 # t1\n 1000011 1.01695599E+03 # el-\n 1000012 9.82228577E+02 # nuel\n 1000013 1.01695599E+03 # mul-\n 1000014 9.82228577E+02 # numl\n 1000015 8.14143616E+02 # tau1\n 1000016 9.60563965E+02 # nutl\n 1000021 1400 # glss\n 1000022 7.99861450E+02 # z1ss\n 1000023 2.61068994E+03 # z2ss\n 1000024 8.00035156E+02 # w1ss\n 1000025 -4.40436084E+03 # z3ss\n 1000035 4.40551318E+03 # z4ss\n 1000037 4.41120557E+03 # w2ss\n 2000001 5.34743115E+03 # dnr\n 2000002 5.29265430E+03 # upr\n 2000003 5.34743115E+03 # str\n 2000004 5.29265430E+03 # chr\n 2000005 5.30680615E+03 # b2\n 2000006 4.64177100E+03 # t2\n 2000011 8.16208069E+02 # er-\n 2000013 8.16208069E+02 # mur-\n 2000015 9.79908020E+02 # tau2\nBlock ALPHA # Effective Higgs mixing parameter\n -1.97611898E-01 # alpha\nBlock STOPMIX # stop mixing matrix\n 1 1 7.59116933E-02 # O_{11}\n 1 2 -9.97114539E-01 # O_{12}\n 2 1 9.97114539E-01 # O_{21}\n 2 2 7.59116933E-02 # O_{22}\nBlock SBOTMIX # sbottom mixing matrix\n 1 1 9.99985218E-01 # O_{11}\n 1 2 5.43311611E-03 # O_{12}\n 2 1 -5.43311611E-03 # O_{21}\n 2 2 9.99985218E-01 # O_{22}\nBlock STAUMIX # stau mixing matrix\n 1 1 1.09243631E-01 # O_{11}\n 1 2 9.94014978E-01 # O_{12}\n 2 1 -9.94014978E-01 # O_{21}\n 2 2 1.09243631E-01 # O_{22}\nBlock NMIX # neutralino mixing matrix\n 1 1 -7.94264197E-04 #\n 1 2 9.99831140E-01 #\n 1 3 -1.76664572E-02 #\n 1 4 5.00912219E-03 #\n 2 1 9.99844968E-01 #\n 2 2 1.10197510E-03 #\n 2 3 1.46711469E-02 #\n 2 4 -9.67472792E-03 #\n 3 1 -3.54152918E-03 #\n 3 2 8.94684903E-03 #\n 3 3 7.07005918E-01 #\n 3 4 7.07142174E-01 #\n 4 1 1.72303095E-02 #\n 4 2 -1.60176214E-02 #\n 4 3 -7.06834614E-01 #\n 4 4 7.06987619E-01 #\nBlock UMIX # chargino U mixing matrix\n 1 1 -9.99664187E-01 # U_{11}\n 1 2 2.59140767E-02 # U_{12}\n 2 1 -2.59140767E-02 # U_{21}\n 2 2 -9.99664187E-01 # U_{22}\nBlock VMIX # chargino V mixing matrix\n 1 1 -9.99951184E-01 # V_{11}\n 1 2 9.88030620E-03 # V_{12}\n 2 1 -9.88030620E-03 # V_{21}\n 2 2 -9.99951184E-01 # V_{22}\nBlock GAUGE Q= 4.02882178E+03 #\n 1 3.57524961E-01 # g`\n 2 6.52378559E-01 # g_2\n 3 1.21928000E+00 # g_3\nBlock YU Q= 4.02882178E+03 #\n 3 3 8.36005747E-01 # y_t\nBlock YD Q= 4.02882178E+03 #\n 3 3 6.48209378E-02 # y_b\nBlock YE Q= 4.02882178E+03 #\n 3 3 5.15556112E-02 # y_tau\nBlock HMIX Q= 4.02882178E+03 # Higgs mixing parameters\n 1 4.40950830E+03 # mu(Q)\n 2 5.00000000E+00 # tan(beta)(M_GUT)\n 3 2.51802856E+02 # Higgs vev at Q\n 4 2.11987360E+07 # m_A^2(Q)\nBlock MSOFT Q= 4.02882178E+03 # DRbar SUSY breaking parameters\n 1 2.64877222E+03 # M_1(Q)\n 2 7.60767517E+02 # M_2(Q)\n 3 -5.11175781E+03 # M_3(Q)\n 31 9.90100281E+02 # MeL(Q)\n 32 9.90100281E+02 # MmuL(Q)\n 33 9.69131042E+02 # MtauL(Q)\n 34 8.43792542E+02 # MeR(Q)\n 35 8.43792542E+02 # MmuR(Q)\n 36 7.78837646E+02 # MtauR(Q)\n 41 5.00258447E+03 # MqL1(Q)\n 42 5.00258447E+03 # MqL2(Q)\n 43 4.39563135E+03 # MqL3(Q)\n 44 5.03392529E+03 # MuR(Q)\n 45 5.03392529E+03 # McR(Q)\n 46 3.69262158E+03 # MtR(Q)\n 47 5.08750146E+03 # MdR(Q)\n 48 5.08750146E+03 # MsR(Q)\n 49 5.11095166E+03 # MbR(Q)\nBlock AU Q= 4.02882178E+03 #\n 1 1 4.39446729E+03 # A_u\n 2 2 4.39446729E+03 # A_c\n 3 3 4.39446729E+03 # A_t\nBlock AD Q= 4.02882178E+03 #\n 1 1 1.04963750E+04 # A_d\n 2 2 1.04963750E+04 # A_s\n 3 3 1.04963750E+04 # A_b\nBlock AE Q= 4.02882178E+03 #\n 1 1 2.95507666E+03 # A_e\n 2 2 2.95507666E+03 # A_mu\n 3 3 2.95507666E+03 # A_tau\n#\n#\n#\n# =================\n# |The decay table|\n# =================\n#\n# PDG Width\nDECAY 1000021 5.50675438E+00 # gluino decay\n# BR NDA ID1 ID2 ID3\n 2.50000000E-01 3 1 -1 1000022\n 2.50000000E-01 3 2 -2 1000022\n 2.50000000E-01 3 1 -2 1000024\n 2.50000000E-01 3 -1 2 -1000024\n#\n# PDG Width\nDECAY 1000024 1.97326979e-18 # chargino decay\n#\n'
),
comEnergy=cms.double(13000.0),
crossSection=cms.untracked.double(0.0284),
filterEfficiency=cms.untracked.double(-1),
hscpFlavor=cms.untracked.string('stau'),
massPoint=cms.untracked.int32(800),
maxEventsToPrint=cms.untracked.int32(0),
particleFile=cms.untracked.string(
'DisappTrks/SignalMC/data/geant4/geant4_AMSB_chargino_800GeV_ctau10000cm.slha'
),
processFile=cms.untracked.string(
'SimG4Core/CustomPhysics/data/RhadronProcessList.txt'),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(0),
slhaFile=cms.untracked.string(''),
useregge=cms.bool(False))
'hltBLifetimeRegionalCtfWithMaterialTracksbbPhiL1FastJetFastPV',
'hltCtfActivityWithMaterialTracks'),
debugModules=cms.untracked.vstring(),
suppressError=cms.untracked.vstring('hltOnlineBeamSpot',
'hltL3MuonCandidates',
'hltL3TkTracksFromL2OIState',
'hltPFJetCtfWithMaterialTracks',
'hltL3TkTracksFromL2IOHit',
'hltL3TkTracksFromL2OIHit'))
process.hltGetConditions = cms.EDAnalyzer("EventSetupRecordDataGetter",
toGet=cms.VPSet(),
verbose=cms.untracked.bool(False))
process.hltGetRaw = cms.EDAnalyzer(
"HLTGetRaw", RawDataCollection=cms.InputTag("rawDataCollector"))
process.hltBoolFalse = cms.EDFilter("HLTBool", result=cms.bool(False))
process.hltTriggerType = cms.EDFilter("HLTTriggerTypeFilter",
SelectedTriggerType=cms.int32(1))
process.hltGtStage2Digis = cms.EDProducer(
"L1TRawToDigi",
lenSlinkTrailer=cms.untracked.int32(8),
lenAMC13Header=cms.untracked.int32(8),
CTP7=cms.untracked.bool(False),
lenAMC13Trailer=cms.untracked.int32(8),
Setup=cms.string("stage2::GTSetup"),
MinFeds=cms.uint32(0),
InputLabel=cms.InputTag("rawDataCollector"),
lenSlinkHeader=cms.untracked.int32(8),
MTF7=cms.untracked.bool(False),
FWId=cms.uint32(0),
TMTCheck=cms.bool(True),
def configure(process, options):
"""Apply configuration to a process."""
# create the main module path
process.add_path('path')
# enable the JSON filter (if given)
if options.jsonFilterFile:
process.enable_json_lumi_filter(options.jsonFilterFile)
# == configure CMSSW modules ==========================================
# -- Jets (default from miniAOD) --------------------------------------
if options.withPATCollections:
# just write out miniAOD jets
process.add_output_commands(
'keep patJets_slimmedJets_*_*',
'keep patJets_slimmedJetsAK8_*_*',
)
# -- Jets (from miniAOD, but with possibly new JECs from GT) ----------
from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection
updateJetCollection(
process,
jetSource = cms.InputTag('slimmedJets'),
labelName = 'UpdatedJEC',
jetCorrections = ('AK4PFchs', cms.vstring(['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual']), 'None') # Update: Safe to always add 'L2L3Residual' as MC contains dummy L2L3Residual corrections (always set to 1)
)
updateJetCollection(
process,
jetSource = cms.InputTag('slimmedJetsAK8'),
labelName = 'UpdatedJECAK8',
jetCorrections = ('AK8PFchs', cms.vstring(['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual']), 'None') # Update: Safe to always add 'L2L3Residual' as MC contains dummy L2L3Residual corrections (always set to 1)
)
process.jecSequence = cms.Sequence(process.patJetCorrFactorsUpdatedJEC * process.updatedPatJetsUpdatedJEC)
process.jecSequenceAK8 = cms.Sequence(process.patJetCorrFactorsUpdatedJECAK8 * process.updatedPatJetsUpdatedJECAK8)
process.path *= process.jecSequence
process.path *= process.jecSequenceAK8
if options.withPATCollections:
process.add_output_commands(
'keep patJets_updatedPatJetsUpdatedJEC_*_*',
'keep patJets_updatedPatJetsUpdatedJECAK8_*_*',
)
# -- Jets (reclustered with jet toolbox) ------------------------------
from Karma.Common.Sequences.jetToolbox_cff import addJetToolboxSequences
# create reclustering sequences
jet_collection_names = []
# AK4CHS jets (include pileupJetID)
jet_collection_names += addJetToolboxSequences(
process, isData=options.isData,
min_jet_pt=15,
jet_algorithm_specs=('ak4',),
pu_subtraction_methods=('CHS',),
do_pu_jet_id=True
)
# AK8CHS jets (no pileupJetID available)
jet_collection_names += addJetToolboxSequences(
process, isData=options.isData,
min_jet_pt=15,
jet_algorithm_specs=('ak8',),
pu_subtraction_methods=('CHS',),
do_pu_jet_id=False
)
# AK4Puppi and AK8Puppi jets
jet_collection_names += addJetToolboxSequences(
process, isData=options.isData,
min_jet_pt=15,
jet_algorithm_specs=('ak4', 'ak8',),
pu_subtraction_methods=('Puppi',),
do_pu_jet_id=False
)
# put reclustering sequences on path
for _jet_collection_name in jet_collection_names:
process.path *= getattr(process, _jet_collection_name)
## write out reclustered jets
#process.add_output_commands('keep patJets_{}_*_*'.format(_jet_collection_name))
# -- Jet ID (precomputed and embedded as userInts) -------------------
for _jet_collection_name in jet_collection_names:
_id_producer_name = "{}IDValueMap".format(_jet_collection_name)
_enriched_jet_collection_name = "{}WithJetIDUserData".format(_jet_collection_name)
# produce the jet id value map
setattr(
process,
_id_producer_name,
cms.EDProducer("PatJetIDValueMapProducer",
filterParams = cms.PSet(
version = cms.string('WINTER16'),
quality = cms.string('TIGHTLEPVETO'),
),
src = cms.InputTag(_jet_collection_name)
)
)
# embed jet id information in pat::Jet itprocess
setattr(
process,
_enriched_jet_collection_name,
cms.EDProducer("PATJetUserDataEmbedder",
src = cms.InputTag(_jet_collection_name),
userInts = cms.PSet(
jetIdWinter16TightLepVeto = cms.InputTag(_id_producer_name),
),
)
)
# add modules to path
process.path *= getattr(process, _id_producer_name)
process.path *= getattr(process, _enriched_jet_collection_name)
# write out ID-enriched jet collection
if options.withPATCollections:
process.add_output_commands(
'keep patJets_{}_*_*'.format(_enriched_jet_collection_name)
)
# -- MET --------------------------------------------------------------
from PhysicsTools.PatUtils.tools.runMETCorrectionsAndUncertainties import runMetCorAndUncFromMiniAOD
# run this to keep MET Type-I correction up-to-date with currently applied JECs
runMetCorAndUncFromMiniAOD(
process,
isData=True,
)
process.path *= process.fullPatMetSequence
if options.withPATCollections:
process.add_output_commands(
'keep patMETs_slimmedMETs_*_*',
)
# -- Electrons --------------------------------------------------------
# just write out miniAOD electrons
if options.withPATCollections:
process.add_output_commands(
"keep patElectrons_slimmedElectrons_*_*",
)
# Note: electron scale/smearing correction information is contained in
# the following userFloats: 'ecalEnergyPreCorr' and 'ecalEnergyPostCorr'
# Note: electron ID information is stored as "pseudo-userData" in
# PAT::Electrons (not in the inherited PAT::Object userData variables)
# and can be accessed using PAT::Electrons::electronID() with the
# corresponding tag (e.g. 'cutBasedElectronID-Summer16-80X-V1-loose')
# -- Muons ------------------------------------------------------------
# just write out miniAOD muons
if options.withPATCollections:
process.add_output_commands(
"keep patMuons_slimmedMuons_*_*",
)
# -- Primary Vertices -------------------------------------------------
from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector
# "good" primary vertices
process.add_module(
'goodOfflinePrimaryVertices',
cms.EDFilter(
'PrimaryVertexObjectFilter',
src = cms.InputTag("offlineSlimmedPrimaryVertices"),
filterParams = pvSelector.clone(
maxZ = 24.0
), # ndof >= 4, rho <= 2
),
on_path='path',
write_out=False,
)
# == END configure CMSSW modules ======================================
# == configure Karma modules ==========================================
# -- preliminary checks
if options.useHLTFilter and not options.hltRegexes:
raise ValueError(
"Option 'useHLTFilter' is true, but 'hltRegexes' "
"is empty: no events would be written out. Aborting!")
elif not options.hltRegexes:
print("[karmaSkim] WARNING: Option 'hltRegexes' is empty:"
"no trigger information will be written out!")
# -- General Event Information ----------------------------------------
from Karma.Skimming.EventProducer_cfi import karmaEventProducer
process.add_module(
'karmaEvents',
karmaEventProducer(isData=options.isData).clone(
goodPrimaryVerticesSrc = cms.InputTag("goodOfflinePrimaryVertices"),
hltProcessName = cms.string("HLT"),
# interesting trigger paths must match one of these regexes:
hltRegexes = cms.vstring(*options.hltRegexes),
#hltRegexes = cms.vstring("HLT_(AK8)?PFJet[0-9]+_v[0-9]+", "HLT_DiPFJetAve[0-9]+_v[0-9]+"),
metFiltersSrc = cms.InputTag("TriggerResults", "", options.metFiltersProcess),
),
on_path='path',
write_out=True,
)
# filter out event if no interesting HLT path fired (if requested)
if options.useHLTFilter:
process.add_module(
'karmaEventHLTFilter',
cms.EDFilter("EventHLTFilter",
cms.PSet(
karmaEventSrc = cms.InputTag("karmaEvents")
)
),
on_path='path',
write_out=False, # don't write out the TriggerResults object
)
# -- MC-specific event information ------------------------------------
if not options.isData:
from Karma.Skimming.GeneratorQCDInfoProducer_cfi import karmaGeneratorQCDInfoProducer
process.add_module(
'karmaGeneratorQCDInfos',
karmaGeneratorQCDInfoProducer.clone(
genEventInfoProductSrc = cms.InputTag("generator"),
),
on_path='path',
write_out=True,
)
# -- Trigger Objects --------------------------------------------------
from Karma.Skimming.TriggerObjectCollectionProducer_cfi import karmaTriggerObjectCollectionProducer
process.add_module(
'karmaTriggerObjects',
karmaTriggerObjectCollectionProducer.clone(
karmaRunSrc = cms.InputTag("karmaEvents"),
),
on_path='path',
write_out=True,
)
# -- Gen-Particles (MC only) ------------------------------------------
if not options.isData:
from Karma.Skimming.GenParticleCollectionProducer_cfi import karmaGenParticleCollectionProducer
process.add_module(
'karmaGenParticles',
karmaGenParticleCollectionProducer.clone(
inputCollection = cms.InputTag("prunedGenParticles"),
),
on_path='path',
write_out=True,
)
# -- MET --------------------------------------------------------------
from Karma.Skimming.METCollectionProducer_cfi import karmaMETCollectionProducer
process.add_module(
'karmaMETs',
karmaMETCollectionProducer.clone(
inputCollection = cms.InputTag("slimmedMETs"),
),
on_path='path',
write_out=True,
)
# -- MET correction levels (as edm::ValueMaps) ------------------------
if options.withMETCorrectionLevels:
from Karma.Skimming.METCorrectedLVValueMapProducer_cfi import karmaMETCorrectedLVValueMapProducer
process.add_module(
'karmaMETCorrectedLVs',
karmaMETCorrectedLVValueMapProducer.clone(
inputCollection = cms.InputTag("karmaMETs"),
associationSpec = cms.VPSet(
# uncorrected MET
cms.PSet(
name = cms.string("Raw"),
transientMapKey = cms.string("corP4Raw"),
),
# uncorrected MET (from CHS candidates)
cms.PSet(
name = cms.string("RawCHS"),
transientMapKey = cms.string("corP4RawCHS"),
),
# uncorrected MET
cms.PSet(
name = cms.string("Type1"),
transientMapKey = cms.string("corP4Type1"),
),
)
),
on_path='path',
write_out=True,
)
from Karma.Skimming.METCorrectedSumEtValueMapProducer_cfi import karmaMETCorrectedSumEtValueMapProducer
process.add_module(
'karmaMETCorrectedSumEts',
karmaMETCorrectedSumEtValueMapProducer.clone(
inputCollection = cms.InputTag("karmaMETs"),
associationSpec = cms.VPSet(
# uncorrected MET
cms.PSet(
name = cms.string("Raw"),
transientMapKey = cms.string("corSumEtRaw"),
),
# uncorrected MET (from CHS candidates)
cms.PSet(
name = cms.string("RawCHS"),
transientMapKey = cms.string("corSumEtRawCHS"),
),
# uncorrected MET
cms.PSet(
name = cms.string("Type1"),
transientMapKey = cms.string("corSumEtType1"),
),
)
),
on_path='path',
write_out=True,
)
# -- Jets -------------------------------------------------------------
from Karma.Skimming.JetCollectionProducer_cfi import karmaJets
from Karma.Skimming.JetCorrectedLVValueMapProducer_cfi import karmaJetCorrectedLVValueMapProducer, karmaJetCorrectedLVValueMapProducerForPuppi
from Karma.Skimming.JetIdValueMapProducers_cfi import karmaJetIdValueMapProducer, karmaJetPileupIdValueMapProducer, karmaJetPileupIdDiscriminantValueMapProducer
# create "karma::Jet" collections from pat::Jets
for _jet_collection_name in jet_collection_names:
# add karma modules for producing the skimmed jet collections
_module_name = "karma{}{}".format(_jet_collection_name[0].upper(), _jet_collection_name[1:])
process.add_module(
_module_name,
karmaJets.clone(
inputCollection = cms.InputTag("{}WithJetIDUserData".format(_jet_collection_name)),
),
on_path='path',
write_out=True,
)
# write out jet ID information to transients (used to fill value maps)
_t = getattr(process, _module_name).transientInformationSpec
_t.fromUserIntAsBool = cms.PSet(
jetIdWinter16TightLepVeto = cms.string("jetIdWinter16TightLepVeto"),
)
if 'AK4PFCHS' in _jet_collection_name:
_t.fromUserFloat = cms.PSet(
pileupJetId = cms.string("AK4PFCHSpileupJetIdEvaluator:fullDiscriminant"),
)
_t.fromUserInt = cms.PSet(
pileupJetId = cms.string("AK4PFCHSpileupJetIdEvaluator:fullId"),
)
# add karma module for producing the Jet ID value map
_valuemap_module_name = "karma{}{}JetIds".format(_jet_collection_name[0].upper(), _jet_collection_name[1:])
process.add_module(
_valuemap_module_name,
karmaJetIdValueMapProducer.clone(
inputCollection = cms.InputTag(_module_name),
),
on_path='path',
write_out=True,
)
# add karma modules for producing the pileup jet ID value maps (AK4CHS-only)
if 'AK4PFCHS' in _jet_collection_name:
_valuemap_module_name = "karma{}{}JetPileupIds".format(_jet_collection_name[0].upper(), _jet_collection_name[1:])
process.add_module(
_valuemap_module_name,
karmaJetPileupIdValueMapProducer.clone(
inputCollection = cms.InputTag(_module_name),
),
on_path='path',
write_out=True,
)
_valuemap_module_name = "karma{}{}JetPileupIdDiscriminants".format(_jet_collection_name[0].upper(), _jet_collection_name[1:])
process.add_module(
_valuemap_module_name,
karmaJetPileupIdDiscriminantValueMapProducer.clone(
inputCollection = cms.InputTag(_module_name),
),
on_path='path',
write_out=True,
)
if 'Puppi' in _jet_collection_name:
_valuemap_producer = karmaJetCorrectedLVValueMapProducerForPuppi
else:
_valuemap_producer = karmaJetCorrectedLVValueMapProducer
# add karma modules for producing the correction level value maps
_valuemap_module_name = "karma{}{}JECs".format(_jet_collection_name[0].upper(), _jet_collection_name[1:])
process.add_module(
_valuemap_module_name,
_valuemap_producer.clone(
inputCollection = cms.InputTag(_module_name),
),
on_path='path',
write_out=True,
)
# -- Gen-Jets ---------------------------------------------------------
if not options.isData:
from Karma.Skimming.GenJetCollectionProducer_cfi import karmaGenJetsAK4, karmaGenJetsAK8
process.add_module(
'karmaGenJetsAK4',
karmaGenJetsAK4.clone(
inputCollection = cms.InputTag("ak4GenJetsNoNu"),
),
on_path='path',
write_out=True,
)
process.add_module(
'karmaGenJetsAK8',
karmaGenJetsAK8.clone(
inputCollection = cms.InputTag("ak8GenJetsNoNu"),
),
on_path='path',
write_out=True,
)
# -- Electrons --------------------------------------------------------
from Karma.Skimming.ElectronCollectionProducer_cfi import karmaElectronCollectionProducer
process.add_module(
'karmaElectrons',
karmaElectronCollectionProducer.clone(),
on_path='path',
write_out=True,
)
# -- Electron IDs -----------------------------------------------------
from Karma.Skimming.ElectronIdValueMapProducer_cfi import karmaElectronIdValueMapProducer
process.add_module(
'karmaElectronIds',
karmaElectronIdValueMapProducer.clone(
inputCollection = cms.InputTag("karmaElectrons")
),
on_path='path',
write_out=True
)
# -- Muons ------------------------------------------------------------
from Karma.Skimming.MuonCollectionProducer_cfi import karmaMuonCollectionProducer
process.add_module(
'karmaMuons',
karmaMuonCollectionProducer.clone(),
on_path='path',
write_out=True,
)
# -- Primary Vertices -------------------------------------------------
from Karma.Skimming.VertexCollectionProducer_cfi import karmaVertexCollectionProducer
process.add_module(
'karmaVertices',
karmaVertexCollectionProducer.clone(),
on_path='path',
write_out=True,
)
# == END configure Karma modules ======================================
# selective writeout based on path decisions
process.enable_selective_writeout('path')
# just in case we need it
return process
passthrough = cms.string("isGlobalMuon && numberOfMatches >= 2"),
fractionOfSharedSegments = cms.double(0.499),
)
# Kalman Muon Calibrations
process.calibratedMuons = cms.EDProducer("KalmanMuonCalibrationsProducer",
muonsCollection = cms.InputTag("boostedMuons"),
isMC = cms.bool(False),
isSync = cms.bool(False),
useRochester = cms.untracked.bool(True),
year = cms.untracked.int32(2017)
)
process.selectedElectrons = cms.EDFilter("PATElectronSelector",
src = cms.InputTag("slimmedElectrons"),
cut = cms.string("pt > 5 && abs(eta)<2.5")
)
process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
calibratedPatElectrons = cms.PSet(
#initialSeed = cms.untracked.uint32(SEED), # for HPC
initialSeed = cms.untracked.uint32(123456), # for crab
engineName = cms.untracked.string('TRandom3')
)
)
from RecoEgamma.EgammaTools.EgammaPostRecoTools import setupEgammaPostRecoSeq
setupEgammaPostRecoSeq(process,
runEnergyCorrections=False,
runVID=True,
phoIDModules=['RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_94X_V2_cff'],
generator = cms.EDFilter(
'Pythia6GeneratorFilter',
comEnergy=cms.double(8000.0),
crossSection=cms.untracked.double(7.168112e+05),
filterEfficiency=cms.untracked.double(1),
maxEventsToPrint=cms.untracked.int32(0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(0),
PythiaParameters=cms.PSet(
pythiaUESettingsBlock,
processParameters=cms.vstring(
'MSEL = 0 ! user defined processes',
'MSUB(15) = 1 ! ff -> Z0 f',
'MSUB(30) = 1 ! ff -> Z0 g',
'MDME(174,1) = 0 ! Z decay into d dbar',
'MDME(175,1) = 0 ! Z decay into u ubar',
'MDME(176,1) = 0 ! Z decay into s sbar',
'MDME(177,1) = 0 ! Z decay into c cbar',
'MDME(178,1) = 0 ! Z decay into b bbar',
'MDME(179,1) = 0 ! Z decay into t tbar',
'MDME(182,1) = 0 ! Z decay into e- e+',
'MDME(183,1) = 0 ! Z decay into nu_e nu_ebar',
'MDME(184,1) = 1 ! Z decay into mu- mu+',
'MDME(185,1) = 0 ! Z decay into nu_mu nu_mubar',
'MDME(186,1) = 0 ! Z decay into tau- tau+',
'MDME(187,1) = 0 ! Z decay into nu_tau nu_taubar',
'CKIN(3) = 15 ! minimum pt hat for hard interactions',
'CKIN(4) = 3000 ! maximum pt hat for hard interactions',
'MSTP(142) = 2 ! Turns on the PYWEVT Pt reweighting routine',
),
CSAParameters=cms.vstring(
'CSAMODE = 7 ! towards a flat QCD spectrum',
'PTPOWER = 4.5 ! reweighting of the pt spectrum',
),
parameterSets=cms.vstring(
'pythiaUESettings',
'processParameters',
'CSAParameters',
)))
import FWCore.ParameterSet.Config as cms
import HLTrigger.HLTfilters.hltHighLevel_cfi
processName = "HLT"
#REDIGI311X
HLTFilterMuOrEleMC = cms.EDFilter(
'SingleTopTriggers',
HLTriggerResults=cms.InputTag("TriggerResults", "", processName),
isMC=cms.untracked.bool(False),
triggerList=cms.vstring(
"HLT_Ele22_SW_TighterCaloIdIsol_L1R_v2",
"HLT_Ele22_SW_TighterCaloIdIsol_L1R_v3",
"HLT_Ele17_SW_TighterEleIdIsol_L1R_v3",
"HLT_IsoMu17_v4",
),
channel=cms.untracked.int32(1), #Useless now
)
HLTFilterMuOrEle = cms.EDFilter(
'SingleTopTriggers',
HLTriggerResults=cms.InputTag("TriggerResults", "", processName),
isMC=cms.untracked.bool(False),
triggerList=cms.vstring(
"HLT_IsoMu17_v5",
"HLT_IsoMu17_v6",
"HLT_Ele27_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v1",
"HLT_Ele27_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v2",
"HLT_Ele27_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v3",
),
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
generator = cms.EDFilter("Pythia8GeneratorFilter",
comEnergy=cms.double(13000.0),
crossSection=cms.untracked.double(0.4255),
filterEfficiency=cms.untracked.double(1),
maxEventsToPrint=cms.untracked.int32(0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(1),
PythiaParameters=cms.PSet(
pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters=cms.vstring(
'ExcitedFermion:dg2dStar = on',
'ExcitedFermion:ug2uStar = on',
'ExcitedFermion:Lambda = 4000',
'4000001:m0 = 4000',
'4000001:onMode = off',
'4000001:onIfMatch = 2 24',
'4000002:m0 = 4000',
'4000002:onMode = off',
'4000002:onIfMatch = 1 24',
),
parameterSets=cms.vstring(
'pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters')))
ProductionFilterSequence = cms.Sequence(generator)
DECAY 6 1.37127534E+00 # Gamma(top)
1.00000000E+00 2 5 24 # BR(top -> bottom W)
"""
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
generator = cms.EDFilter(
"Pythia8GeneratorFilter",
pythiaPylistVerbosity=cms.untracked.int32(1),
filterEfficiency=cms.untracked.double(1),
pythiaHepMCVerbosity=cms.untracked.bool(False),
SLHATableForPythia8=cms.string('%s' % SLHA_TABLE),
comEnergy=cms.double(COM_ENERGY),
crossSection=cms.untracked.double(CROSS_SECTION),
maxEventsToPrint=cms.untracked.int32(1),
PythiaParameters=cms.PSet(
pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters=cms.vstring('Higgs:useBSM = on', PROCESS,
'SLHA:allowUserOverride = off',
'SLHA:minMassSM = 100.',
'PhaseSpace:mHatMin = 56.0'),
parameterSets=cms.vstring('pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters')))
ProductionFilterSequence = cms.Sequence(generator)
eventAutoFlushCompressedSize=cms.untracked.int32(20971520),
fileName=cms.untracked.string('file:output.root'),
outputCommands=process.RAWSIMEventContent.outputCommands,
splitLevel=cms.untracked.int32(0))
# Additional output definition
# Other statements
process.XMLFromDBSource.label = cms.string("Extended")
process.genstepfilter.triggerConditions = cms.vstring("generation_step")
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag,
'102X_upgrade2018_realistic_v11', '')
process.gencount = cms.EDFilter("CandViewCountFilter",
minNumber=cms.uint32(1),
src=cms.InputTag("genfilter"))
process.generator = cms.EDFilter(
"Pythia8HadronizerFilter",
PythiaParameters=cms.PSet(
parameterSets=cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'pythia8PSweightsSettings',
'pythia8PowhegEmissionVetoSettings',
'processParameters'),
processParameters=cms.vstring(
'POWHEG:nFinal = 1', 'ParticleDecays:limitTau0= off', '25:m0 =125',
'25:addChannel = 1 1.0 102 4900101 -4900101', '25:0:onMode=0',
'25:1:onMode=0', '25:2:onMode=0', '25:3:onMode=0', '25:4:onMode=0',
'25:5:onMode=0', '25:6:onMode=0', '25:7:onMode=0', '25:8:onMode=0',
secondaryFileNames=cms.untracked.vstring('file:testRunMerge0.root',
'file:testRunMerge1.root',
'file:testRunMerge2.root',
'file:testRunMerge3.root',
'file:testRunMerge4.root',
'file:testRunMerge5.root'),
duplicateCheckMode=cms.untracked.string('checkEachRealDataFile'),
noEventSort=cms.untracked.bool(False))
process.testGetterOfProducts = cms.EDFilter(
"TestGetterOfProducts",
processName=cms.string('PROD'),
expectedInputTagLabels=cms.vstring(
'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J', 'K', 'L',
'aliasForThingToBeDropped2', 'dependsOnThingToBeDropped1', 'm1', 'm2',
'm3', 'makeThingToBeDropped', 'makeThingToBeDropped1',
'makeThingToBeDropped2', 'thingWithMergeProducer', 'tryNoPut'),
expectedLabelsAfterGet=cms.vstring('A', 'B', 'C', 'D', 'E', 'F', 'G', 'H',
'I', 'J', 'K', 'L',
'aliasForThingToBeDropped2',
'dependsOnThingToBeDropped1', 'm1',
'm2', 'm3', 'makeThingToBeDropped',
'thingWithMergeProducer'))
process.testGetterOfProductsA = cms.EDAnalyzer(
"TestGetterOfProductsA",
processName=cms.string('PROD'),
branchType=cms.int32(0),
expectedInputTagLabels=cms.vstring(
'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J', 'K', 'L',
'aliasForThingToBeDropped2', 'dependsOnThingToBeDropped1', 'm1', 'm2',
'm3', 'makeThingToBeDropped', 'makeThingToBeDropped1',
generator = cms.EDFilter(
'Pythia6GeneratorFilter',
comEnergy=cms.double(7000.0),
crossSection=cms.untracked.double(1.0),
filterEfficiency=cms.untracked.double(1.0),
maxEventsToPrint=cms.untracked.int32(0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(0),
ExternalDecays=cms.PSet(
Tauola=cms.untracked.PSet(
TauolaPolar, # tau polarisation switch on
TauolaDefaultInputCards #set default TAUOLA input card
),
parameterSets=cms.vstring('Tauola')),
UseExternalGenerators=cms.untracked.bool(True),
PythiaParameters=cms.PSet(
pythiaUESettingsBlock,
processParameters=cms.vstring(
'MSEL=0 ! User defined processes',
'MSUB(186)= 1 ! gg->QQbarH (MSSM)',
'KFPR(186,2)= 5 ! Q = b Registered by [email protected]',
# MSSM settings
'IMSS(1)= 1 ! MSSM with parameters',
'RMSS(5)= 30. ! tan beta',
'RMSS(19)= 800. ! m_A',
'RMSS(1)= 100. ! M1',
'RMSS(2)= 200. ! M2',
'RMSS(3)= 800. ! Mg',
'RMSS(4)= 200. ! mu',
'RMSS(6)= 1000. ! MS',
'RMSS(7)= 1000. ! MS',
'RMSS(8)= 1000. ! MS',
'RMSS(9)= 1000. ! MS',
'RMSS(10)= 1000. ! MS',
'RMSS(11)= 1000. ! MS',
'RMSS(12)= 1000. ! MS',
'RMSS(13)= 1000. ! MS',
'RMSS(14)= 1000. ! MS',
'RMSS(15)= 2000. ! Ab',
'RMSS(16)= 2000. ! At',
'RMSS(17)= 2000. ! Atau',
# Switch off / on desirable channels for A
'MDME(420,1)=0 ! Higgs(H) decay into d dbar',
'MDME(421,1)=0 ! Higgs(H) decay into u ubar',
'MDME(422,1)=0 ! Higgs(H) decay into s sbar',
'MDME(423,1)=0 ! Higgs(H) decay into c cbar',
'MDME(424,1)=1 ! Higgs(H) decay into b bbar',
'MDME(425,1)=0 ! Higgs(H) decay into t tbar',
'MDME(426,1)=0 ! Higgs(H) decay into b bbar',
'MDME(427,1)=0 ! Higgs(H) decay into t tbar',
'MDME(428,1)=0 ! Higgs(H) decay into e- e+',
'MDME(429,1)=0 ! Higgs(H) decay into mu- mu+',
'MDME(430,1)=0 ! Higgs(H) decay into tau- tau+',
'MDME(431,1)=0 ! Higgs(H) decay into tau- tau+',
'MDME(432,1)=0 ! Higgs(H) decay into g g',
'MDME(433,1)=0 ! Higgs(H) decay into gamma gamma',
'MDME(434,1)=0 ! Higgs(H) decay into gamma Z0',
'MDME(435,1)=0 ! Higgs(H) decay into Z0 Z0',
'MDME(436,1)=0 ! Higgs(H) decay into W+ W-',
'MDME(437,1)=0 ! Higgs(H) decay into Z0 h0',
'MDME(438,1)=0 ! Higgs(H) decay into h0 h0',
'MDME(439,1)=0 ! Higgs(H) decay into W+ H-',
'MDME(440,1)=0 ! Higgs(H) decay into H+ W-',
'MDME(441,1)=0 ! Higgs(H) decay into Z0 A0',
'MDME(442,1)=0 ! Higgs(H) decay into h0 A0',
'MDME(443,1)=0 ! Higgs(H) decay into A0 A0',
'MDME(444,1)=0 ! Higgs(H) decay into ~chi_10 ~chi_10',
'MDME(445,1)=0 ! Higgs(H) decay into ~chi_20 ~chi_10',
'MDME(446,1)=0 ! Higgs(H) decay into ~chi_20 ~chi_20',
'MDME(447,1)=0 ! Higgs(H) decay into ~chi_30 ~chi_10',
'MDME(448,1)=0 ! Higgs(H) decay into ~chi_30 ~chi_20',
'MDME(449,1)=0 ! Higgs(H) decay into ~chi_30 ~chi_30',
'MDME(450,1)=0 ! Higgs(H) decay into ~chi_40 ~chi_10',
'MDME(451,1)=0 ! Higgs(H) decay into ~chi_40 ~chi_20',
'MDME(452,1)=0 ! Higgs(H) decay into ~chi_40 ~chi_30',
'MDME(453,1)=0 ! Higgs(H) decay into ~chi_40 ~chi_40',
'MDME(454,1)=0 ! Higgs(H) decay into ~chi_1+ ~chi_1-',
'MDME(455,1)=0 ! Higgs(H) decay into ~chi_1+ ~chi_2-',
'MDME(456,1)=0 ! Higgs(H) decay into ~chi_2+ ~chi_1-',
'MDME(457,1)=0 ! Higgs(H) decay into ~chi_2+ ~chi_2-',
'MDME(458,1)=0 ! Higgs(H) decay into ~d_L ~d_Lbar',
'MDME(459,1)=0 ! Higgs(H) decay into ~d_R ~d_Rbar',
'MDME(460,1)=0 ! Higgs(H) decay into ~d_L ~d_Rbar',
'MDME(461,1)=0 ! Higgs(H) decay into ~d_Lbar ~d_R',
'MDME(462,1)=0 ! Higgs(H) decay into ~u_L ~u_Lbar',
'MDME(463,1)=0 ! Higgs(H) decay into ~u_R ~u_Rbar',
'MDME(464,1)=0 ! Higgs(H) decay into ~u_L ~u_Rbar',
'MDME(465,1)=0 ! Higgs(H) decay into ~u_Lbar ~u_R',
'MDME(466,1)=0 ! Higgs(H) decay into ~s_L ~s_Lbar',
'MDME(467,1)=0 ! Higgs(H) decay into ~s_R ~s_Rbar',
'MDME(468,1)=0 ! Higgs(H) decay into ~s_L ~s_Rbar',
'MDME(469,1)=0 ! Higgs(H) decay into ~s_Lbar ~s_R',
'MDME(470,1)=0 ! Higgs(H) decay into ~c_L ~c_Lbar',
'MDME(471,1)=0 ! Higgs(H) decay into ~c_R ~c_Rbar',
'MDME(472,1)=0 ! Higgs(H) decay into ~c_L ~c_Rbar',
'MDME(473,1)=0 ! Higgs(H) decay into ~c_Lbar ~c_R',
'MDME(474,1)=0 ! Higgs(H) decay into ~b_1 ~b_1bar',
'MDME(475,1)=0 ! Higgs(H) decay into ~b_2 ~b_2bar',
'MDME(476,1)=0 ! Higgs(H) decay into ~b_1 ~b_2bar',
'MDME(477,1)=0 ! Higgs(H) decay into ~b_1bar ~b_2',
'MDME(478,1)=0 ! Higgs(H) decay into ~t_1 ~t_1bar',
'MDME(479,1)=0 ! Higgs(H) decay into ~t_2 ~t_2bar',
'MDME(480,1)=0 ! Higgs(H) decay into ~t_1 ~t_2bar',
'MDME(481,1)=0 ! Higgs(H) decay into ~t_1bar ~t_2',
'MDME(482,1)=0 ! Higgs(H) decay into ~e_L- ~e_L+',
'MDME(483,1)=0 ! Higgs(H) decay into ~e_R- ~e_R+',
'MDME(484,1)=0 ! Higgs(H) decay into ~e_L- ~e_R+',
'MDME(485,1)=0 ! Higgs(H) decay into ~e_L+ ~e_R-',
'MDME(486,1)=0 ! Higgs(H) decay into ~nu_eL ~nu_eLbar',
'MDME(487,1)=0 ! Higgs(H) decay into ~nu_eR ~nu_eRbar',
'MDME(488,1)=0 ! Higgs(H) decay into ~nu_eL ~nu_eRbar',
'MDME(489,1)=0 ! Higgs(H) decay into ~nu_eLbar ~nu_eR',
'MDME(490,1)=0 ! Higgs(H) decay into ~mu_L- ~mu_L+',
'MDME(491,1)=0 ! Higgs(H) decay into ~mu_R- ~mu_R+',
'MDME(492,1)=0 ! Higgs(H) decay into ~mu_L- ~mu_R+',
'MDME(493,1)=0 ! Higgs(H) decay into ~mu_L+ ~mu_R-',
'MDME(494,1)=0 ! Higgs(H) decay into ~nu_muL ~nu_muLbar',
'MDME(495,1)=0 ! Higgs(H) decay into ~nu_muR ~nu_muRbar',
'MDME(496,1)=0 ! Higgs(H) decay into ~nu_muL ~nu_muRbar',
'MDME(497,1)=0 ! Higgs(H) decay into ~nu_muLbar ~nu_muR',
'MDME(498,1)=0 ! Higgs(H) decay into ~tau_1- ~tau_1+',
'MDME(499,1)=0 ! Higgs(H) decay into ~tau_1- ~tau_1+',
'MDME(500,1)=0 ! Higgs(H) decay into ~tau_1- ~tau_1+',
'MDME(501,1)=0 ! Higgs(H) decay into ~tau_1- ~tau_1+',
'MDME(502,1)=0 ! Higgs(H) decay into ~tau_2- ~tau_2+'),
# This is a vector of ParameterSet names to be read, in this order
parameterSets=cms.vstring('pythiaUESettings', 'processParameters')))
# Z2star tune with pT-ordered showers
from Configuration.Generator.PythiaUEZ2starSettings_cfi import *
from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
generator = cms.EDFilter(
"Pythia6HadronizerFilter",
pythiaPylistVerbosity=cms.untracked.int32(1),
filterEfficiency=cms.untracked.double(1.0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
comEnergy=cms.double(8000.0),
# crossSection = cms.untracked.double(1.0),
maxEventsToPrint=cms.untracked.int32(0),
ExternalDecays=cms.PSet(Tauola=cms.untracked.PSet(TauolaPolar,
TauolaDefaultInputCards),
parameterSets=cms.vstring('Tauola')),
PythiaParameters=cms.PSet(
pythiaUESettingsBlock,
processParameters=cms.vstring(
'MSEL=0 ! User defined processes',
'PMAS(5,1)=4.75 ! b quark mass',
'PMAS(6,1)=172.5 ! t quark mass',
'MSTJ(1)=1 ! Fragmentation/hadronization on or off',
'MSTP(61)=1 ! Parton showering on or off', 'PARP(64)=0.25',
'PARP(72)=0.5'),
parameterSets=cms.vstring('pythiaUESettings', 'processParameters')))
configurationMetadata = cms.untracked.PSet(
version=cms.untracked.string('$Revision: 1.1 $'),
name=cms.untracked.string(
'$Source: /cvs/CMSSW/CMSSW/Configuration/GenProduction/python/EightTeV/POWHEG_PYTHIA6_top_tauola_cff.py,v $'
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
generator = cms.EDFilter("Pythia8GeneratorFilter",
comEnergy=cms.double(13000.0),
crossSection=cms.untracked.double(0.013),
filterEfficiency=cms.untracked.double(1),
maxEventsToPrint=cms.untracked.int32(1),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(1),
PythiaParameters=cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters=cms.vstring(
'NewGaugeBoson:ffbar2gmZZprime = on',
'Zprime:gmZmode = 3',
'32:m0 = 3000',
'32:onMode = off',
'32:onIfAny = 11 13 15',
),
parameterSets=cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters')))
ProductionFilterSequence = cms.Sequence(generator)
process.load('HLTrigger/HLTfilters/hltHighLevel_cfi')
process.hltHighLevel.throw = cms.bool(False)
process.hltHighLevel.HLTPaths = cms.vstring()
#process.pfMVAMEt.isTestSample = cms.bool(False)
from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
from RecoMET.METProducers.PFMET_cfi import pfMet
process.load("CondCore.DBCommon.CondDBCommon_cfi")
from CondCore.DBCommon.CondDBSetup_cfi import *
## load Puppi stuff
process.load('CommonTools/PileupAlgos/Puppi_cff')
process.load('BaconProd/Ntupler/myPUPPICorrections_cff')
process.pfCandNoLep = cms.EDFilter(
"CandPtrSelector",
src=cms.InputTag("particleFlow"),
cut=cms.string("abs(pdgId) != 13 && abs(pdgId) != 11 && abs(pdgId) != 15"))
process.pfCandLep = cms.EDFilter(
"CandPtrSelector",
src=cms.InputTag("particleFlow"),
cut=cms.string("abs(pdgId) == 13 || abs(pdgId) == 11 || abs(pdgId) == 15"))
process.puppinolep = process.puppi.clone()
process.puppinolep.candName = 'pfCandNoLep'
process.puppimetinput = cms.EDProducer("CandViewMerger",
src=cms.VInputTag(
"pfCandLep", "puppinolep"))
from PhysicsTools.SelectorUtils.tools.vid_id_tools import *
my_id_modules = [
'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_PHYS14_PU20bx25_V2_cff'
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, '94X_mc2017_realistic_v10', '')
process.generator = cms.EDFilter("Pythia8PtGunV2",
PGunParameters = cms.PSet(
AddAntiParticle = cms.bool(True),
MaxEta = cms.double(2.4),
MaxMass = cms.double(1.6),
#MaxMass = cms.double(2.6),
MaxPhi = cms.double(3.14159265359),
MaxPt = cms.double(100.0),
MinEta = cms.double(1.6),
MinMass = cms.double(0.0),
MinPhi = cms.double(-3.14159265359),
MinPt = cms.double(10.0),
ParticleID = cms.vint32(111,111,111,111)
),
PythiaParameters = cms.PSet(
parameterSets = cms.vstring('processParameters'),
processParameters = cms.vstring('111:onMode = off',
'111:onIfMatch = 22 22')
),
Verbosity = cms.untracked.int32(0),
firstRun = cms.untracked.uint32(1),
maxEventsToPrint = cms.untracked.int32(1),
psethack = cms.string('double pi0 pt'),
pythiaHepMCVerbosity = cms.untracked.bool(True),
pythiaPylistVerbosity = cms.untracked.int32(1)
)
process.selectedPhotons = cms.EDFilter("EtaPtMinCandViewSelector",
src = cms.InputTag("gedPhotons"),
generator = cms.EDFilter("Pythia6GeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(1),
# put here the efficiency of your filter (1. if no filter)
filterEfficiency = cms.untracked.double(1.0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
# put here the cross section of your process (in pb)
crossSection = cms.untracked.double(1.0),
maxEventsToPrint = cms.untracked.int32(1),
comEnergy = cms.double(8000.0),
ExternalDecays = cms.PSet(
Tauola = cms.untracked.PSet(
TauolaPolar,
TauolaDefaultInputCards
),
parameterSets = cms.vstring('Tauola')
),
PythiaParameters = cms.PSet(
pythiaUESettingsBlock,
processParameters = cms.vstring('PMAS(25,1)=121.0 !mass of Higgs',
'MSEL=0 ! user selection for process',
'MSUB(102)=0 !ggH',
'MSUB(123)=0 !ZZ fusion to H',
'MSUB(124)=0 !WW fusion to H',
'MSUB(24)=1 !ZH production',
'MSUB(26)=1 !WH production',
'MSUB(121)=0 !gg to ttH',
'MSUB(122)=0 !qq to ttH',
'MDME(210,1)=0 !Higgs decay into dd',
'MDME(211,1)=0 !Higgs decay into uu',
'MDME(212,1)=0 !Higgs decay into ss',
'MDME(213,1)=0 !Higgs decay into cc',
'MDME(214,1)=0 !Higgs decay into bb',
'MDME(215,1)=0 !Higgs decay into tt',
'MDME(216,1)=0 !Higgs decay into',
'MDME(217,1)=0 !Higgs decay into Higgs decay',
'MDME(218,1)=0 !Higgs decay into e nu e',
'MDME(219,1)=0 !Higgs decay into mu nu mu',
'MDME(220,1)=0 !Higgs decay into tau nu tau',
'MDME(221,1)=0 !Higgs decay into Higgs decay',
'MDME(222,1)=0 !Higgs decay into g g',
'MDME(223,1)=1 !Higgs decay into gam gam',
'MDME(224,1)=0 !Higgs decay into gam Z',
'MDME(225,1)=0 !Higgs decay into Z Z',
'MDME(226,1)=0 !Higgs decay into W W'),
# This is a vector of ParameterSet names to be read, in this order
parameterSets = cms.vstring('pythiaUESettings',
'processParameters')
)
)
generator = cms.EDFilter("Pythia8HadronizerFilter",
pythiaPylistVerbosity = cms.untracked.int32(1),
filterEfficiency = cms.untracked.double(1),
pythiaHepMCVerbosity = cms.untracked.bool(False),
comEnergy = cms.double(13000.),
crossSection = cms.untracked.double(CROSS_SECTION),
maxEventsToPrint = cms.untracked.int32(10),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock, # Old PYTHIA tune
# pythia8CP5SettingsBlock,
processParameters = cms.vstring(
'SLHA:useDecayTable = off', # use pythia8 decay mode instead of decays defined in LH accord
'LesHouches:setLifetime = 2',
"1023:new = Zd Zdbar 3 0 0",
"1023:m0 = %s" % MASS_Zd,
# "10231:mWidth = %s" % getGammaEpsilon2(MASS_Zd)*EPSILON_Zd*EPSILON_Zd,
"1023:mWidth = %s" % WIDTH_Zd,
"1023:tau0 = %s" % LIFETIME_Zd_inMM,
"1023:isResonance = on",
"1023:mayDecay = on",
"1023:oneChannel = 1 %s 100 -2 2" %BR_UUBAR,
"1023:addChannel = 1 %s 100 -4 4" %BR_CCBAR,
"1023:addChannel = 1 %s 100 -13 13" %BR_MUMU,
"1023:addChannel = 1 %s 100 -11 11" %BR_EE,
"1023:addChannel = 1 %s 100 -15 15" %BR_TAUTAU,
"1023:addChannel = 1 %s 100 -3 3" %BR_SSBAR,
"1023:addChannel = 1 %s 100 1 -1" %BR_DDBAR,
"1023:addChannel = 1 %s 100 -5 5" %BR_BBBAR,
"1023:addChannel = 1 %s 100 16 -12" %BR_TAUNUTAUNUBAR,
"1023:addChannel = 1 %s 100 14 -12" %BR_MUNUMUNUBAR,
"1023:addChannel = 1 %s 100 12 -12" %BR_ENUENUBAR
),
parameterSets = cms.vstring(
'pythia8CommonSettings',
'pythia8CUEP8M1Settings', #Old Pythia Tune
# 'pythia8CP5Settings',
'processParameters'
)
)
)
generator = cms.EDFilter(
"Pythia6GeneratorFilter",
pythiaHepMCVerbosity=cms.untracked.bool(False),
maxEventsToPrint=cms.untracked.int32(0),
pythiaPylistVerbosity=cms.untracked.int32(1),
filterEfficiency=cms.untracked.double(0.779),
crossSection=cms.untracked.double(7899.),
comEnergy=cms.double(7000.0),
ExternalDecays=cms.PSet(Tauola=cms.untracked.PSet(TauolaPolar,
TauolaDefaultInputCards),
parameterSets=cms.vstring('Tauola')),
UseExternalGenerators=cms.untracked.bool(True),
PythiaParameters=cms.PSet(
pythiaUESettingsBlock,
processParameters=cms.vstring(
'MSEL = 0 !User defined processes',
'MSUB(2) = 1 !W production',
'MDME(190,1) = 0 !W decay into dbar u',
'MDME(191,1) = 0 !W decay into dbar c',
'MDME(192,1) = 0 !W decay into dbar t',
'MDME(194,1) = 0 !W decay into sbar u',
'MDME(195,1) = 0 !W decay into sbar c',
'MDME(196,1) = 0 !W decay into sbar t',
'MDME(198,1) = 0 !W decay into bbar u',
'MDME(199,1) = 0 !W decay into bbar c',
'MDME(200,1) = 0 !W decay into bbar t',
'MDME(205,1) = 0 !W decay into bbar tp',
'MDME(206,1) = 1 !W decay into e+ nu_e',
'MDME(207,1) = 1 !W decay into mu+ nu_mu',
'MDME(208,1) = 1 !W decay into tau+ nu_tau'),
# This is a vector of ParameterSet names to be read, in this order
parameterSets=cms.vstring('pythiaUESettings', 'processParameters')))
from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
generator = cms.EDFilter(
"Pythia8GeneratorFilter",
comEnergy=cms.double(13000.0),
filterEfficiency=cms.untracked.double(1.),
maxEventsToPrint=cms.untracked.int32(1),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(1),
SLHAFileForPythia8=cms.string('%s' % SLHA_TABLE),
ExternalDecays=cms.PSet(Tauola=cms.untracked.PSet(TauolaPolar,
TauolaDefaultInputCards),
parameterSets=cms.vstring('Tauola')),
PythiaParameters=cms.PSet(pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters=cms.vstring(
'SUSY:all off',
'SUSY:gg2squarkantisquark = on',
'SUSY:qqbar2squarkantisquark = on',
'SUSY:idA = 1000006',
'SUSY:idB = 1000006',
),
parameterSets=cms.vstring(
'pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters')))
ProductionFilterSequence = cms.Sequence(generator)
'Configuration/StandardSequences/MagneticField_AutoFromDBCurrent_cff')
process.load('Configuration/StandardSequences/GeometryIdeal_cff')
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.GlobalTag.globaltag = 'GR10_P_V8::All'
process.load("Configuration/StandardSequences/RawToDigi_Data_cff")
process.load("Configuration/StandardSequences/ReconstructionCosmics_cff")
process.load('Configuration/EventContent/EventContentCosmics_cff')
######################################Super Pointing#################################################
process.cosmicMuonsBarrelOnlyFilter = cms.EDFilter(
"HLTMuonPointingFilter",
SALabel=cms.string("cosmicMuonsBarrelOnly"),
PropagatorName=cms.string("SteppingHelixPropagatorAny"),
radius=cms.double(10.0),
maxZ=cms.double(50.0))
process.cosmicMuonsFilter = cms.EDFilter(
"HLTMuonPointingFilter",
SALabel=cms.string("cosmicMuons"),
PropagatorName=cms.string("SteppingHelixPropagatorAny"),
radius=cms.double(10.0),
maxZ=cms.double(50.0))
process.cosmicMuons1LegFilter = cms.EDFilter(
"HLTMuonPointingFilter",
SALabel=cms.string("cosmicMuons1Leg"),
PropagatorName=cms.string("SteppingHelixPropagatorAny"),
radius=cms.double(10.0),
process.GlobalTag.globaltag = '80X_dataRun2_Prompt_v15'
# =============== Import Sequences =====================
#Trigger Selection
### Comment out for the timing being assuming running on secondary dataset with trigger bit selected already
import HLTrigger.HLTfilters.hltHighLevel_cfi
process.hltHM = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
#process.hltHM.HLTPaths = ['HLT_PAFullTracks_Multiplicity280_v*']
process.hltHM.HLTPaths = ['HLT_PAFullTracks_Multiplicity280_v*']
process.hltHM.andOr = cms.bool(True)
process.hltHM.throw = cms.bool(False)
process.PAprimaryVertexFilter = cms.EDFilter(
"VertexSelector",
src=cms.InputTag("offlinePrimaryVertices"),
cut=cms.string(
"!isFake && abs(z) <= 25 && position.Rho <= 2 && tracksSize >= 2"),
# cut = cms.string("!isFake && abs(z) <= 1 && position.Rho <= 2 && tracksSize >= 5"),
filter=cms.bool(True), # otherwise it won't filter the events
)
#Reject beam scraping events standard pp configuration
process.NoScraping = cms.EDFilter("FilterOutScraping",
applyfilter=cms.untracked.bool(True),
debugOn=cms.untracked.bool(False),
numtrack=cms.untracked.uint32(10),
thresh=cms.untracked.double(0.25))
process.PAcollisionEventSelection = cms.Sequence(
process.hfCoincFilter * process.PAprimaryVertexFilter * process.NoScraping)
process.eventFilter_HM = cms.Sequence(
process = cms.Process("VlsbInfo")
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))
process.source = cms.Source(
"HcalTBSource",
fileNames=cms.untracked.vstring("file:/tmp/chenyi/HTB_108538.root"),
streams=cms.untracked.vstring('Chunk699', 'HCAL_Trigger', 'HCAL_SlowData',
'HCAL_QADCTDC', 'HCAL_DCC021'))
process.tbunpack = cms.EDFilter(
"HcalTBObjectUnpacker",
#IncludeUnmatchedHits = cms.untracked.bool(False),
HcalTriggerFED=cms.untracked.int32(1),
HcalVLSBFED=cms.untracked.int32(699),
HcalTDCFED=cms.untracked.int32(8),
HcalQADCFED=cms.untracked.int32(8),
HcalSlowDataFED=cms.untracked.int32(3),
ConfigurationFile=cms.untracked.string('configQADCTDC_TB2009.txt'))
process.vlsbinfo = cms.EDProducer(
"VLSBInformationProducer",
minSample=cms.untracked.uint32(0),
maxSample=cms.untracked.uint32(31),
baselineSamples=cms.untracked.uint32(2),
mip=cms.untracked.string("PercentageCorrectedGeV.txt"),
useMotherBoard0=cms.untracked.bool(True),
useMotherBoard1=cms.untracked.bool(True),
useMotherBoard2=cms.untracked.bool(False),
useMotherBoard3=cms.untracked.bool(True),
fileName=cms.string(outname),
closeFileFast=cms.untracked.bool(True))
'''
################# Sequence ####################
'''
process.l1bits = cms.EDProducer(
"L1TriggerResultsConverter",
src=cms.InputTag("gtStage2Digis"),
# src_ext=cms.InputTag("gtStage2Digis"),
# storeUnprefireableBit=cms.bool(True),
legacyL1=cms.bool(False),
)
process.TnP = cms.EDFilter("TagAndProbeProducer_MC",
muonIDScaleFactors=cms.int32(1),
requireTag=cms.int32(1),
verbose=cms.int32(0))
process.p = cms.Path(process.l1bits + process.TnP)
# DEBUG -- dump the event content
# process.output = cms.OutputModule(
# "PoolOutputModule",
# fileName = cms.untracked.string('edm_output.root'),
# )
# process.output_step = cms.EndPath(process.output)
#
# process.schedule = cms.Schedule(
# process.p,
# process.output_step)
'''
import FWCore.ParameterSet.Config as cms
isoTrackPrefilter = cms.EDFilter("IsoTrackPrefilter",
trackSrc = cms.InputTag("trackSel"),
muonSrc = cms.InputTag("standAloneMuons"),
generatorSrc = cms.InputTag("genParticles"),
minPt = cms.double(25),
coneSize = cms.double(0.3),
isolationThreshold = cms.double(99999.0),
signalPDGId = cms.int32(36))
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
generator = cms.EDFilter("Pythia8GeneratorFilter",
comEnergy=cms.double(13000.0),
crossSection=cms.untracked.double(0.00000782),
filterEfficiency=cms.untracked.double(1),
maxEventsToPrint=cms.untracked.int32(0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(1),
PythiaParameters=cms.PSet(
pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters=cms.vstring(
'ExtraDimensionsG*:ffbar2G* = on',
'ExtraDimensionsG*:kappaMG = 0.0143953282',
'5100039:m0 = 9000', '5100039:onMode = off',
'5100039:onIfAny = 1 2 3 4 5'),
parameterSets=cms.vstring(
'pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters',
)))
ProductionFilterSequence = cms.Sequence(generator)
)
#process.Tracks2Vertex = AssociationMaps.clone(
#
# #Choose which map should be created
# #"VertexToTracks", "TracksToVertex" or "Both"
# AssociationType = cms.InputTag('TracksToVertex'),
#
#)
process.trig = cms.EDFilter(
"HLTHighLevel",
TriggerResultsTag=cms.InputTag("TriggerResults", "", "HLT"),
HLTPaths=cms.vstring(
"HLT_PFJet*"
# HLTPaths = cms.vstring("HLT_ZeroBias_v*"
# HLTPaths = cms.vstring("HLT_*"
),
eventSetupPathsKey=cms.string(''),
andOr=cms.bool(True),
throw=cms.bool(False))
#process.analyzer = cms.EDAnalyzer('TrkEffAnalyzer',
# pionPixel=cms.bool(False),
# kaonPixel=cms.bool(False),
# runMode=cms.int32(1)
#)
process.analyzer = cms.EDAnalyzer(
'LambdaAnalyzer',
doGen=cms.bool(False),
process = cms.Process("ProfileCleanedMIP")
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))
process.source = cms.Source(
"HcalTBSource",
fileNames=cms.untracked.vstring("file:/tmp/chenyi/HTB_108533.root"),
streams=cms.untracked.vstring('Chunk699', 'HCAL_Trigger', 'HCAL_SlowData',
'HCAL_QADCTDC', 'HCAL_DCC021'))
process.tbunpack = cms.EDFilter(
"HcalTBObjectUnpacker",
#IncludeUnmatchedHits = cms.untracked.bool(False),
HcalTriggerFED=cms.untracked.int32(1),
HcalVLSBFED=cms.untracked.int32(699),
HcalTDCFED=cms.untracked.int32(8),
HcalQADCFED=cms.untracked.int32(8),
HcalSlowDataFED=cms.untracked.int32(3),
ConfigurationFile=cms.untracked.string('configQADCTDC_TB2009.txt'))
process.vlsbinfo = cms.EDProducer(
"VLSBInformationProducer",
minSample=cms.untracked.uint32(0),
maxSample=cms.untracked.uint32(31),
baselineSamples=cms.untracked.uint32(2),
mip=cms.untracked.string("MIP_EarlyRejection_Mean.txt"),
# mip = cms.untracked.string("MIP_EarlyRejection_Median.txt"),
useMotherBoard0=cms.untracked.bool(True),
useMotherBoard1=cms.untracked.bool(True),
useMotherBoard2=cms.untracked.bool(False),
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
generator = cms.EDFilter("Pythia8GeneratorFilter",
maxEventsToPrint=cms.untracked.int32(1),
pythiaPylistVerbosity=cms.untracked.int32(1),
filterEfficiency=cms.untracked.double(1.0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
comEnergy=cms.double(13000.),
PythiaParameters=cms.PSet(
pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters=cms.vstring(
'ExtraDimensionsLED:LambdaT = 5000',
'ExtraDimensionsLED:ffbar2llbar = on',
'ExtraDimensionsLED:gg2llbar = on',
'ExtraDimensionsLED:NegInt = 0',
'ExtraDimensionsLED:CutOffmode = 0',
'PhaseSpace:mHatMin = 1300',
'PhaseSpace:mHatMax = 2000',
),
parameterSets=cms.vstring(
'pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters',
)))
import FWCore.ParameterSet.Config as cms
hltHcalSimpleRecHitFilter = cms.EDFilter(
'HLTHcalSimpleRecHitFilter',
saveTags=cms.bool(True),
HFRecHitCollection=cms.InputTag('hltHfreco'),
threshold=cms.double(3),
minNHitsNeg=cms.int32(1),
minNHitsPos=cms.int32(1),
doCoincidence=cms.bool(True),
maskedChannels=cms.vuint32())
