process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.GlobalTag.globaltag = cms.string('100X_upgrade2018_realistic_v10')
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring('file:../../Core/data/store_mc_RunIISpring18MiniAOD_BulkGravTohhTohbbhbb_narrow_M-2000_13TeV-madgraph_MINIAODSIM_100X_upgrade2018_realistic_v10-v1_30000_24A0230C-B530-E811-ADE3-14187741120B.root')
)
if len(options.inputFiles)>0: process.source.fileNames = options.inputFiles
################### EDProducer ##############################
process.jetImageProducer = cms.EDProducer(prodname,
JetTag = cms.InputTag('slimmedJetsAK8'),
topN = cms.uint32(5),
# imageList = cms.string("../../Core/data/imagenet_classes.txt"),
imageList = cms.string("../../Core/data/top_classes.txt"),
Client = Client,
)
# Let it run
process.p = cms.Path(
process.jetImageProducer
)
process.MessageLogger.cerr.FwkReport.reportEvery = 1
keep_msgs = ['JetImageProducer','TFClientRemote','TFClientLocal']
for msg in keep_msgs:
process.MessageLogger.categories.append(msg)
setattr(process.MessageLogger.cerr,msg,
cms.untracked.PSet(
from RecoJets.JetProducers.JetIDParams_cfi import *
ak5PFmatch = patJetGenJetMatch.clone(
src=cms.InputTag("ak5PFJets"), matched=cms.InputTag("ak5HiGenJetsCleaned"))
ak5PFparton = patJetPartonMatch.clone(src=cms.InputTag("ak5PFJets"))
ak5PFcorr = patJetCorrFactors.clone(
useNPV=False,
# primaryVertices = cms.InputTag("hiSelectedVertex"),
levels=cms.vstring('L2Relative', 'L3Absolute'),
src=cms.InputTag("ak5PFJets"),
payload="AK5PF_hiIterativeTracks")
ak5PFJetID = cms.EDProducer('JetIDProducer',
JetIDParams,
src=cms.InputTag('ak5CaloJets'))
ak5PFclean = heavyIonCleanedGenJets.clone(
src=cms.InputTag('ak5HiGenJetsCleaned'))
ak5PFbTagger = bTaggers("ak5PF")
#create objects locally since they dont load properly otherwise
ak5PFmatch = ak5PFbTagger.match
ak5PFparton = ak5PFbTagger.parton
ak5PFPatJetFlavourAssociation = ak5PFbTagger.PatJetFlavourAssociation
ak5PFJetTracksAssociatorAtVertex = ak5PFbTagger.JetTracksAssociatorAtVertex
ak5PFSimpleSecondaryVertexHighEffBJetTags = ak5PFbTagger.SimpleSecondaryVertexHighEffBJetTags
ak5PFSimpleSecondaryVertexHighPurBJetTags = ak5PFbTagger.SimpleSecondaryVertexHighPurBJetTags
ak5PFCombinedSecondaryVertexBJetTags = ak5PFbTagger.CombinedSecondaryVertexBJetTags
splitLevel=cms.untracked.int32(0))
# Additional output definition
# Other statements
process.genstepfilter.triggerConditions = cms.vstring("generation_step")
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic', '')
process.generator = cms.EDProducer("FlatRandomPtGunProducer",
AddAntiParticle=cms.bool(True),
PGunParameters=cms.PSet(
MaxEta=cms.double(3.0),
MaxPhi=cms.double(3.14159265359),
MaxPt=cms.double(25.01),
MinEta=cms.double(1.3),
MinPhi=cms.double(-3.14159265359),
MinPt=cms.double(24.99),
PartID=cms.vint32(211)),
Verbosity=cms.untracked.int32(0),
firstRun=cms.untracked.uint32(1),
psethack=cms.string('single pi pt 25'))
# Path and EndPath definitions
process.generation_step = cms.Path(process.pgen)
process.simulation_step = cms.Path(process.psim)
process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
process.endjob_step = cms.EndPath(process.endOfProcess)
process.FEVTDEBUGoutput_step = cms.EndPath(process.FEVTDEBUGoutput)
# Schedule definition
import FWCore.ParameterSet.Config as cms
csvWeights = cms.EDProducer(
"CSVWeightProducer",
src=cms.InputTag("catJets"),
minPt=cms.double(20),
maxEta=cms.double(2.41),
csvSFHF=cms.string("csv_rwt_fit_hf_v2_final_2017_1_10test.root"),
csvSFLF=cms.string("csv_rwt_fit_lf_v2_final_2017_1_10test.root"),
)
import FWCore.ParameterSet.Config as cms
externalLHEProducer = cms.EDProducer(
'ExternalLHEProducer',
scriptName=cms.FileInPath(
"GeneratorInterface/LHEInterface/data/run_madgraph_gridpack.sh"),
outputFile=cms.string("W2JetsToLNu_matchingup_8TeV-madgraph_final.lhe"),
numberOfParameters=cms.uint32(10),
args=cms.vstring(
'slc5_ia32_gcc434/madgraph/V5_1.3.30/8TeV_Summer12/W2JetsToLNu_matchingup_8TeV-madgraph/v1',
'W2JetsToLNu_matchingup_8TeV-madgraph', 'false', 'true', 'wjets', '5',
'30', 'true', '2', '99'),
nEvents=cms.uint32(100000))
Multiplex a cut on a PFTauDiscriminator using another PFTauDiscriminator as the
index.
Used by the anti-electron MVA, which needs to choose what cut to apply on the
MVA output depending on what the category is.
'''
import FWCore.ParameterSet.Config as cms
recoTauDiscriminantCutMultiplexer = cms.EDProducer(
"RecoTauDiscriminantCutMultiplexer",
PFTauProducer=cms.InputTag("fixme"),
toMultiplex=cms.InputTag("fixme"),
Prediscriminants=cms.PSet(BooleanOperator=cms.string("and"),
decayMode=cms.PSet(
Producer=cms.InputTag("fixme"),
cut=cms.double(0.))),
key=cms.InputTag("fixme"), # a discriminator
mapping=cms.VPSet(
cms.PSet(
category=cms.uint32(0),
cut=cms.double(0.5),
),
cms.PSet(
category=cms.uint32(1),
cut=cms.double(0.2),
),
))
import RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi
iterativeDetachedTripletTracks = RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi.ctfWithMaterialTracks.clone(
)
iterativeDetachedTripletTracks.src = 'iterativeDetachedTripletTrackCandidates'
iterativeDetachedTripletTracks.TTRHBuilder = 'WithoutRefit'
iterativeDetachedTripletTracks.Fitter = 'KFFittingSmootherSecond'
iterativeDetachedTripletTracks.Propagator = 'PropagatorWithMaterial'
# track merger
#from FastSimulation.Tracking.IterativeSecondTrackMerger_cfi import *
detachedTripletStepTracks = cms.EDProducer(
"FastTrackMerger",
TrackProducers=cms.VInputTag(
cms.InputTag("iterativeDetachedTripletTrackCandidates"),
cms.InputTag("iterativeDetachedTripletTracks")),
RemoveTrackProducers=cms.untracked.VInputTag(
cms.InputTag("initialStepTracks"),
cms.InputTag("lowPtTripletStepTracks"),
cms.InputTag("pixelPairStepTracks")),
trackAlgo=cms.untracked.uint32(7), # iter3
MinNumberOfTrajHits=cms.untracked.uint32(3),
MaxLostTrajHits=cms.untracked.uint32(1))
# TRACK SELECTION AND QUALITY FLAG SETTING.
import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
detachedTripletStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
src='detachedTripletStepTracks',
trackSelectors=cms.VPSet(
RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
name='detachedTripletStepVtxLoose',
chi2n_par=1.6,
res_par=(0.003, 0.001),
def runMETs(process,era):
dBFile = os.path.expandvars(era+".db")
print dBFile
if usePrivateSQlite:
process.jec = cms.ESSource("PoolDBESSource",
CondDBSetup,
connect = cms.string("sqlite_file:"+dBFile),
toGet = cms.VPSet(
cms.PSet(
record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_"+era+"_AK4PF"),
label= cms.untracked.string("AK4PF")
),
cms.PSet(
record = cms.string("JetCorrectionsRecord"),
tag = cms.string("JetCorrectorParametersCollection_"+era+"_AK4PFchs"),
label= cms.untracked.string("AK4PFchs")
),
)
)
process.es_prefer_jec = cms.ESPrefer("PoolDBESSource",'jec')
#===========================================================================================================================#
isMC = False
if 'MC' in dBFile:
isMC = True
##only corrects JEC/JER uncertainties for the baseline collection; means only muon corrected and EG/Muon corrected MET collections are guaranteed to have the correct JEC uncertainties.
runMetCorAndUncFromMiniAOD(process, metType="PF",
recoMetFromPFCs=False,
postfix="",
isData=(not isMC),
)
if isMC or os.environ["CMSSW_VERSION"].count("CMSSW_9"):
process.flashggMets = cms.EDProducer('FlashggMetProducer',
verbose = cms.untracked.bool(False),
metTag = cms.InputTag('slimmedMETs'),
)
process.flashggMetSequence = cms.Sequence(process.flashggMets)
if not isMC and os.environ["CMSSW_VERSION"].count("CMSSW_8_0_28"):
corMETFromMuonAndEG(process,
pfCandCollection="", #not needed
electronCollection="slimmedElectronsBeforeGSFix",
photonCollection="slimmedPhotonsBeforeGSFix",
corElectronCollection="slimmedElectrons",
corPhotonCollection="slimmedPhotons",
allMETEGCorrected=True,
muCorrection=False,
eGCorrection=True,
runOnMiniAOD=True,
postfix="FullMETClean"
)
process.slimmedMETsFullMETClean = process.slimmedMETs.clone()
process.slimmedMETsFullMETClean.src = cms.InputTag("patPFMetT1FullMETClean")
process.slimmedMETsFullMETClean.rawVariation = cms.InputTag("patPFMetRawFullMETClean")
process.slimmedMETsFullMETClean.t1Uncertainties = cms.InputTag("patPFMetT1%sFullMETClean")
del process.slimmedMETsFullMETClean.caloMET
process.egcorrMET = cms.Sequence(
process.cleanedPhotonsFullMETClean+process.cleanedCorPhotonsFullMETClean+
process.matchedPhotonsFullMETClean + process.matchedElectronsFullMETClean +
process.corMETPhotonFullMETClean+process.corMETElectronFullMETClean+
process.patPFMetT1FullMETClean+process.patPFMetRawFullMETClean+
process.patPFMetT1SmearFullMETClean+process.patPFMetT1TxyFullMETClean+
process.patPFMetTxyFullMETClean+process.patPFMetT1JetEnUpFullMETClean+
process.patPFMetT1JetResUpFullMETClean+process.patPFMetT1SmearJetResUpFullMETClean+
process.patPFMetT1ElectronEnUpFullMETClean+process.patPFMetT1PhotonEnUpFullMETClean+
process.patPFMetT1MuonEnUpFullMETClean+process.patPFMetT1TauEnUpFullMETClean+
process.patPFMetT1UnclusteredEnUpFullMETClean+process.patPFMetT1JetEnDownFullMETClean+
process.patPFMetT1JetResDownFullMETClean+process.patPFMetT1SmearJetResDownFullMETClean+
process.patPFMetT1ElectronEnDownFullMETClean+process.patPFMetT1PhotonEnDownFullMETClean+
process.patPFMetT1MuonEnDownFullMETClean+process.patPFMetT1TauEnDownFullMETClean+
process.patPFMetT1UnclusteredEnDownFullMETClean+process.slimmedMETsFullMETClean)
process.flashggMets = cms.EDProducer('FlashggMetProducer',
verbose = cms.untracked.bool(False),
metTag = cms.InputTag('slimmedMETs'),
)
# process.flashggMetsMuons = cms.EDProducer('FlashggMetProducer',
# verbose = cms.untracked.bool(False),
# metTag = cms.InputTag('slimmedMETs'),
# )
# process.flashggMetsEG = cms.EDProducer('FlashggMetProducer',
# verbose = cms.untracked.bool(False),
# metTag = cms.InputTag('slimmedMETsEGClean'),
# )
# process.flashggMets = cms.EDProducer('FlashggMetProducer',
# verbose = cms.untracked.bool(False),
# metTag = cms.InputTag('slimmedMETsFullMETClean'),
# )
# process.flashggMetsEGmuon = cms.EDProducer('FlashggMetProducer',
# verbose = cms.untracked.bool(False),
# metTag = cms.InputTag('slimmedMETsMuEGClean'),
# )
# process.flashggMetsUncorr = cms.EDProducer('FlashggMetProducer',
# verbose = cms.untracked.bool(False),
# metTag = cms.InputTag('slimmedMETsUncorrected'),
# )
# process.flashggMetSequence = cms.Sequence(process.flashggMetsMuons *process.flashggMetsEGmuon*process.flashggMets*process.flashggMetsUncorr)
process.flashggMetSequence = cms.Sequence(process.flashggMets)
PythiaParameters=cms.PSet(
processParameters=cms.vstring(
'Main:timesAllowErrors = 10000',
'ParticleDecays:limitTau0 = on',
'ParticleDecays:tauMax = 10',
'HardQCD:all = on',
'PhaseSpace:pTHatMin = 10 ',
'PhaseSpace:pTHatMax = 20 ',
'Tune:pp 5',
'Tune:ee 3',
),
parameterSets=cms.vstring('processParameters')))
genParticlesForFilter = cms.EDProducer(
"GenParticleProducer",
saveBarCodes=cms.untracked.bool(True),
src=cms.InputTag("generator"),
abortOnUnknownPDGCode=cms.untracked.bool(False))
bctoefilter = cms.EDFilter(
"BCToEFilter",
filterAlgoPSet=cms.PSet(
eTThreshold=cms.double(10),
genParSource=cms.InputTag("genParticlesForFilter")))
emenrichingfilter = cms.EDFilter(
"EMEnrichingFilter",
filterAlgoPSet=cms.PSet(
isoGenParETMin=cms.double(5.),
isoGenParConeSize=cms.double(0.1),
clusterThreshold=cms.double(5.),
],
boolEventSelFlagProducer="BoolEventSelFlagProducer",
pyModuleName=__name__)
produceEventSelFlagsZtoMuTauOS = zToMuTauEventSelConfiguratorOS.configure()
zToMuTauEventSelConfiguratorSS = eventSelFlagProdConfigurator(
[cfgDiTauCandidateForMuTauNonZeroChargeCut],
boolEventSelFlagProducer="BoolEventSelFlagProducer",
pyModuleName=__name__)
produceEventSelFlagsZtoMuTauSS = zToMuTauEventSelConfiguratorSS.configure()
produceEventSelFlagsZtoMuTau = cms.Sequence(produceEventSelFlagsZtoMuTauOS *
produceEventSelFlagsZtoMuTauSS)
isRecZtoMuTau = cms.EDProducer(
"BoolEventSelFlagProducer",
pluginName=cms.string('isRecZtoMuTau'),
pluginType=cms.string('MultiBoolEventSelFlagSelector'),
flags=cms.VInputTag(
cms.InputTag('Trigger'), cms.InputTag('muonTrkIPcut', 'cumulative'),
cms.InputTag('tauElectronVeto', 'cumulative'),
cms.InputTag('diTauCandidateForMuTauZeroChargeCut', 'cumulative'),
cms.InputTag('primaryEventVertexPositionForMuTau'),
cms.InputTag('diMuPairZmumuHypothesisVetoByLooseIsolation'),
cms.InputTag('diMuPairDYmumuHypothesisVeto')))
selectZtoMuTauEvents = cms.Sequence(produceEventSelFlagsZtoMuTau *
isRecZtoMuTau)
)
process.Timing = cms.Service("Timing")
process.SimpleMemoryCheck = cms.Service("SimpleMemoryCheck",
oncePerEventMode = cms.untracked.bool(True),
showMallocInfo = cms.untracked.bool(True),
dump = cms.untracked.bool(True),
ignoreTotal = cms.untracked.int32(1)
)
process.load("IOMC.RandomEngine.IOMC_cff")
process.RandomNumberGeneratorService.generator.initialSeed = 456789
process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
process.rndmStore = cms.EDProducer("RandomEngineStateProducer")
process.TFileService = cms.Service("TFileService",
fileName = cms.string('minbias_FTFP_BERT_EMM.root')
)
# Event output
process.output = cms.OutputModule("PoolOutputModule",
process.FEVTSIMEventContent,
fileName = cms.untracked.string('simevent_minbias_FTFP_BERT_EMM.root')
)
process.generation_step = cms.Path(process.pgen)
process.simulation_step = cms.Path(process.psim)
process.analysis_step = cms.Path(process.CaloSimHitStudy)
process.out_step = cms.EndPath(process.output)
import FWCore.ParameterSet.Config as cms
# link to card:
# https://github.com/cms-sw/genproductions/tree/master/bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/VBF_HH
externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
nEvents = cms.untracked.uint32(5000),
outputFile = cms.string('cmsgrid_final.lhe'),
scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh'),
numberOfParameters = cms.uint32(1),
args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/slc7_amd64_gcc700/14TeV/powheg/V2/ggHH_EWChL_NLO_cHHH_1/ggHH_EWChL_slc7_amd64_gcc700_CMSSW_10_6_0_my_ggHH_EWChL_NLO_cHHH_1.tgz'),
)
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
from Configuration.Generator.Pythia8PowhegEmissionVetoSettings_cfi import *
generator = cms.EDFilter("Pythia8HadronizerFilter",
maxEventsToPrint = cms.untracked.int32(1),
pythiaPylistVerbosity = cms.untracked.int32(1),
filterEfficiency = cms.untracked.double(1.0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
comEnergy = cms.double(14000.),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
pythia8PowhegEmissionVetoSettingsBlock,
processParameters = cms.vstring(
'POWHEG:nFinal = 2', ## Number of final state particles
'23:mMin = 0.05', # the lower limit of the allowed mass range generated by the Breit-Wigner (in GeV)
import FWCore.ParameterSet.Config as cms
hgcalPhotonMVA = cms.EDProducer("HGCalPhotonMVAProducer",
photons = cms.InputTag("photonsFromMultiCl"),
rhoSrc = cms.InputTag("fixedGridRhoFastjetAll"),
electronsFromMultiCl = cms.InputTag("ecalDrivenGsfElectronsFromMultiCl"),
usePAT = cms.bool(False),
gedGsfElectrons = cms.InputTag("gedGsfElectrons"),
conversions = cms.InputTag("conversions"),
beamspot = cms.InputTag("offlineBeamSpot"),
externalVariables = cms.PSet(
sigmaUU = cms.InputTag("hgcPhotonID:sigmaUU"),
sigmaVV = cms.InputTag("hgcPhotonID:sigmaVV"),
e4oEtot = cms.InputTag("hgcPhotonID:e4oEtot"),
layerEfrac10 = cms.InputTag("hgcPhotonID:layerEfrac10"),
layerEfrac90 = cms.InputTag("hgcPhotonID:layerEfrac90"),
seedEnergyFH = cms.InputTag("hgcPhotonID:seedEnergyFH"),
seedEnergyEE = cms.InputTag("hgcPhotonID:seedEnergyEE"),
measuredDepth = cms.InputTag("hgcPhotonID:measuredDepth"),
depthCompatibility = cms.InputTag("hgcPhotonID:depthCompatibility"),
caloIsoRing0 = cms.InputTag("hgcPhotonID:caloIsoRing0"),
caloIsoRing1 = cms.InputTag("hgcPhotonID:caloIsoRing1"),
caloIsoRing2 = cms.InputTag("hgcPhotonID:caloIsoRing2"),
caloIsoRing3 = cms.InputTag("hgcPhotonID:caloIsoRing3"),
caloIsoRing4 = cms.InputTag("hgcPhotonID:caloIsoRing4"),
scRawEnergy = cms.InputTag("hgcPhotonID:scRawEnergy"),
),
barrelTrainingFile = cms.FileInPath("RecoEgamma/Phase2InterimID/data/barrelV4_BDT.weights.xml"),
endcapTrainingFile = cms.FileInPath("RecoEgamma/Phase2InterimID/data/endcapV4_BDT.weights.xml"),
)
import FWCore.ParameterSet.Config as cms
process = cms.Process("TESTMODULEDELETE")
process.maxEvents.input = 1
process.options.deleteNonConsumedUnscheduledModules = False
process.source = cms.Source("EmptySource")
process.intEventProducer = cms.EDProducer("IntProducer", ivalue = cms.int32(1))
process.intEventConsumer = cms.EDAnalyzer("IntTestAnalyzer",
moduleLabel = cms.untracked.InputTag("intEventProducer"),
valueMustMatch = cms.untracked.int32(1)
)
process.intProducerNotConsumed = cms.EDProducer("edmtest::MustRunIntProducer",
ivalue = cms.int32(1),
mustRunEvent = cms.bool(False)
)
process.t = cms.Task(
process.intEventProducer,
process.intProducerNotConsumed
)
process.p = cms.Path(
process.intEventConsumer,
process.t
)
#--------------------------------------------------
# PFElectron Specific
#--------------------------------------------------
process.pfAllElectrons = cms.EDFilter("PdgIdPFCandidateSelector",
pdgId = cms.vint32(11, -11),
#src = cms.InputTag("pfNoPileUp")
src = cms.InputTag("particleFlow")
)
process.gensource = cms.EDProducer("GenParticlePruner",
src = cms.InputTag("genParticles"),
select = cms.vstring('drop *',
# for matching
#'keep+ pdgId = 23',
'keep pdgId = 11',
'keep pdgId = -11'
## for fake rate
#'keep pdgId = 211', # pi+
#'keep pdgId = -211' # pi-
)
)
process.pfPileUp = cms.EDProducer("PFPileUp",
Enable = cms.bool(True),
PFCandidates = cms.InputTag("particleFlow"),
verbose = cms.untracked.bool(False),
#Vertices = cms.InputTag("offlinePrimaryVertices")
Vertices = cms.InputTag("offlinePrimaryVerticesWithBS")
)
#process.pfNoPileUp = cms.EDProducer("TPPileUpPFCandidatesOnPFCandidates",
#--------------------------------------------------------------------------------
# produce combinations of muon + tau-jet pairs
#--------------------------------------------------------------------------------
allMuTauPairs = cms.EDProducer(
"PATMuTauPairProducer",
useLeadingTausOnly=cms.bool(False),
srcLeg1=cms.InputTag('selectedPatMuonsTrkIPcumulative'),
srcLeg2=cms.InputTag('selectedPatTausForMuTauElectronVetoCumulative'),
dRmin12=cms.double(0.3),
srcMET=cms.InputTag('patMETs'),
srcPrimaryVertex=cms.InputTag("offlinePrimaryVerticesWithBS"),
srcBeamSpot=cms.InputTag("offlineBeamSpot"),
srcGenParticles=cms.InputTag('genParticles'),
recoMode=cms.string(""),
doSVreco=cms.bool(True),
nSVfit=cms.PSet(),
scaleFuncImprovedCollinearApprox=cms.string('1'),
doPFMEtSign=cms.bool(True),
pfMEtSign=cms.PSet(srcPFJets=cms.InputTag('ak5PFJets'),
srcPFCandidates=cms.InputTag('particleFlow'),
resolution=METSignificance_params,
dRoverlapPFJet=cms.double(0.3),
dRoverlapPFCandidate=cms.double(0.1)),
doMtautauMin=cms.bool(True),
verbosity=cms.untracked.int32(0))
#--------------------------------------------------------------------------------
# configure (new) SVfit algorithm
# (using combination of PS + MET likelihoods + logM regularization term
# to reconstruct mass of tau lepton pair, as described in CMS AN-11-165)
process.ak8PFJetsCHS = ak8PFJetsCHS.clone(src='puppi', jetPtMin=100.0)
process.ak8PFJetsCHSPruned = ak8PFJetsCHSPruned.clone(src='puppi',
jetPtMin=100.0)
process.ak8PFJetsCHSPrunedMass = ak8PFJetsCHSPrunedMass.clone()
process.ak8PFJetsCHSSoftDrop = ak8PFJetsCHSSoftDrop.clone(src='puppi',
jetPtMin=100.0)
process.ak8PFJetsCHSSoftDropMass = ak8PFJetsCHSSoftDropMass.clone()
process.NjettinessAK8 = cms.EDProducer(
"NjettinessAdder",
src=cms.InputTag("ak8PFJetsCHS"),
Njets=cms.vuint32(1, 2, 3, 4),
# variables for measure definition :
measureDefinition=cms.uint32(0), # CMS default is normalized measure
beta=cms.double(1.0), # CMS default is 1
R0=cms.double(0.8), # CMS default is jet cone size
Rcutoff=cms.double(999.0), # not used by default
# variables for axes definition :
axesDefinition=cms.uint32(6), # CMS default is 1-pass KT axes
nPass=cms.int32(0), # not used by default
akAxesR0=cms.double(-999.0) # not used by default
)
process.substructureSequence = cms.Sequence()
process.substructureSequence += process.puppi
process.substructureSequence += process.ak8PFJetsCHS
process.substructureSequence += process.NjettinessAK8
process.substructureSequence += process.ak8PFJetsCHSPruned
process.substructureSequence += process.ak8PFJetsCHSPrunedMass
process.substructureSequence += process.ak8PFJetsCHSSoftDrop
process.substructureSequence += process.ak8PFJetsCHSSoftDropMass
from PhysicsTools.NanoAOD.jets_cff import *
from PhysicsTools.NanoAOD.muons_cff import *
from PhysicsTools.NanoAOD.taus_cff import *
from PhysicsTools.NanoAOD.electrons_cff import *
from PhysicsTools.NanoAOD.photons_cff import *
from PhysicsTools.NanoAOD.globals_cff import *
from PhysicsTools.NanoAOD.ttbarCategorization_cff import *
from PhysicsTools.NanoAOD.genparticles_cff import *
from PhysicsTools.NanoAOD.particlelevel_cff import *
from PhysicsTools.NanoAOD.vertices_cff import *
from PhysicsTools.NanoAOD.met_cff import *
from PhysicsTools.NanoAOD.triggerObjects_cff import *
from PhysicsTools.NanoAOD.isotracks_cff import *
from PhysicsTools.NanoAOD.NanoAODEDMEventContent_cff import *
nanoMetadata = cms.EDProducer("UniqueStringProducer",
strings=cms.PSet(tag=cms.string("untagged"), ))
linkedObjects = cms.EDProducer(
"PATObjectCrossLinker",
jets=cms.InputTag("finalJets"),
muons=cms.InputTag("finalMuons"),
electrons=cms.InputTag("finalElectrons"),
taus=cms.InputTag("finalTaus"),
photons=cms.InputTag("finalPhotons"),
)
simpleCleanerTable = cms.EDProducer(
"NanoAODSimpleCrossCleaner",
name=cms.string("cleanmask"),
doc=cms.string("simple cleaning mask with priority to leptons"),
jets=cms.InputTag("linkedObjects", "jets"),
# BasicCluster producer
cosmicBasicClusters = cms.EDProducer(
"CosmicClusterProducer",
barrelHits=cms.InputTag('ecalRecHit', 'EcalRecHitsEB'),
endcapHits=cms.InputTag('ecalRecHit', 'EcalRecHitsEE'),
barrelUncalibHits=cms.InputTag('ecalFixedAlphaBetaFitUncalibRecHit',
'EcalUncalibRecHitsEB'),
endcapUncalibHits=cms.InputTag('ecalFixedAlphaBetaFitUncalibRecHit',
'EcalUncalibRecHitsEE'),
barrelClusterCollection=cms.string('CosmicBarrelBasicClusters'),
EndcapSecondThr=cms.double(9.99),
VerbosityLevel=cms.string('ERROR'),
BarrelSingleThr=cms.double(14.99),
BarrelSupThr=cms.double(2.0),
EndcapSupThr=cms.double(3.0),
barrelShapeAssociation=cms.string('CosmicBarrelShapeAssoc'),
clustershapecollectionEE=cms.string('CosmicEndcapShape'),
clustershapecollectionEB=cms.string('CosmicBarrelShape'),
EndcapSingleThr=cms.double(25.99),
endcapClusterCollection=cms.string('CosmicEndcapBasicClusters'),
BarrelSecondThr=cms.double(4.99),
EndcapSeedThr=cms.double(9.99),
BarrelSeedThr=cms.double(4.99),
endcapShapeAssociation=cms.string('CosmicEndcapShapeAssoc'),
posCalcParameters=cms.PSet(T0_barl=cms.double(7.4),
T0_endc=cms.double(3.1),
T0_endcPresh=cms.double(1.2),
LogWeighted=cms.bool(True),
W0=cms.double(4.2),
X0=cms.double(0.89)))
flashggPhotons = cms.EDProducer(
'FlashggPhotonProducer',
photonTag=cms.InputTag('slimmedPhotons'),
pfCandidatesTag=cms.InputTag("packedPFCandidates"),
reducedBarrelRecHitCollection=cms.InputTag('reducedEgamma',
'reducedEBRecHits'),
reducedEndcapRecHitCollection=cms.InputTag('reducedEgamma',
'reducedEERecHits'),
reducedPreshowerRecHitCollection=cms.InputTag('reducedEgamma',
'reducedESRecHits'),
vertexTag=cms.InputTag("offlineSlimmedPrimaryVertices"),
vertexCandidateMapTag=cms.InputTag("flashggVertexMapNonUnique"),
rhoFixedGridCollection=cms.InputTag('fixedGridRhoAll'),
photonIdMVAweightfile_EB=cms.FileInPath(
"flashgg/MicroAOD/data/MVAweights_80X_barrel_ICHEPvtx.xml"),
photonIdMVAweightfile_EE=cms.FileInPath(
"flashgg/MicroAOD/data/MVAweights_80X_endcap_ICHEPvtx.xml"),
useNonZsLazyTools=cms.bool(True),
recomputeNonZsClusterShapes=cms.bool(False),
addRechitFlags=cms.bool(True),
doOverlapRemovalForIsolation=cms.bool(True),
useVtx0ForNeutralIso=cms.bool(True),
extraCaloIsolations=cms.VPSet(),
extraIsolations=cms.VPSet(),
genPhotonTag=cms.InputTag("flashggGenPhotonsExtra"),
maxGenDeltaR=cms.double(0.1),
copyExtraGenInfo=cms.bool(True),
convTag=cms.InputTag('reducedEgamma', 'reducedConversions'),
beamSpotTag=cms.InputTag('offlineBeamSpot'),
elecTag=cms.InputTag("slimmedElectrons"))
mix = cms.EDProducer("MixingModule",
digitizers = cms.PSet(theDigitizers),
LabelPlayback = cms.string(''),
maxBunch = cms.int32(3),
minBunch = cms.int32(-12), ## in terms of 25 nsec
bunchspace = cms.int32(25), ##ns
mixProdStep1 = cms.bool(False),
mixProdStep2 = cms.bool(False),
playback = cms.untracked.bool(False),
useCurrentProcessOnly = cms.bool(False),
input = cms.SecSource("PoolSource",
type = cms.string('probFunction'),
nbPileupEvents = cms.PSet(
probFunctionVariable = cms.vint32(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50),
probValue = cms.vdouble(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025
),
histoFileName = cms.untracked.string('histProbFunction.root'),
),
sequential = cms.untracked.bool(False),
manage_OOT = cms.untracked.bool(True), ## manage out-of-time pileup
## setting this to True means that the out-of-time pileup
## will have a different distribution than in-time, given
## by what is described on the next line:
OOT_type = cms.untracked.string('Poisson'), ## generate OOT with a Poisson matching the number chosen for in-time
#OOT_type = cms.untracked.string('fixed'), ## generate OOT with a fixed distribution
#intFixed_OOT = cms.untracked.int32(2),
fileNames = FileNames
),
mixObjects = cms.PSet(theMixObjects)
)
trajectorySeedProducer = cms.EDProducer("TrajectorySeedProducer",
simTrackSelection = cms.PSet(
# The smallest pT (in GeV) to create a track candidate
pTMin = cms.double(-1),
# skip SimTracks processed in previous iterations
skipSimTrackIds = cms.VInputTag(),
maxZ0 = cms.double(-1),
maxD0 = cms.double(-1),
),
# minimum number of layer crossed (with hits on them) by the simtrack
minLayersCrossed = cms.uint32(0),
#if empty, BS compatibility is skipped
beamSpot = cms.InputTag("offlineBeamSpot"),
#if empty, PV compatibility is skipped
primaryVertex = cms.InputTag(""),
nSigmaZ = cms.double(-1),
originHalfLength= cms.double(-1),
originRadius = cms.double(-1),
ptMin = cms.double(-1),
# Inputs: tracker rechits, beam spot position.
recHits = cms.InputTag("siTrackerGaussianSmearingRecHits","TrackerGSMatchedRecHits"),
layerList = cms.vstring(),
)
#from: ~psilva/public/HZZSkim/;
#FIXME: cf. Configuration.Skimming.PDWG_HZZSkim_cff
import FWCore.ParameterSet.Config as cms
####### MUON SELECTION
MUON_CUT_HZZSKIM = ("abs(eta)<2.5 && (isGlobalMuon || isTrackerMuon)")
DIMUON_MASSCUT_HZZSKIM = ("mass > 40")
DIMUON_KINCUT_HZZSKIM = (
"(max(daughter(0).pt(),daughter(1).pt())>20 && min(daughter(0).pt,daughter(1).pt())>10)"
)
goodHzzMuons = cms.EDFilter("PATMuonRefSelector",
src=cms.InputTag("patMuonsWithTrigger"),
cut=cms.string(MUON_CUT_HZZSKIM))
hzzKinDiMuons = cms.EDProducer("CandViewShallowCloneCombiner",
decay=cms.string("goodHzzMuons goodHzzMuons"),
checkCharge=cms.bool(False),
cut=cms.string(DIMUON_KINCUT_HZZSKIM))
hzzKinDiMuonsFilter = cms.EDFilter("CandViewCountFilter",
src=cms.InputTag("hzzKinDiMuons"),
minNumber=cms.uint32(1))
hzzMassDiMuons = hzzKinDiMuons.clone(cut=cms.string(DIMUON_MASSCUT_HZZSKIM))
hzzMassDiMuonsFilter = hzzKinDiMuonsFilter.clone(
src=cms.InputTag("hzzMassDiMuons"))
####### ELECTRON SELECTION
ELECTRON_CUT_HZZSKIM = ("abs(eta)<2.5")
DIELECTRON_MASSCUT_HZZSKIM = ("mass > 40")
DIELECTRON_KINCUT_HZZSKIM = (
"(max(daughter(0).pt(),daughter(1).pt())>20 && min(daughter(0).pt,daughter(1).pt())>10)"
)
goodHzzElectrons = cms.EDFilter("PATElectronRefSelector",
process.externalLHEAsciiDumper.lheFileName = cms.string('output.lhe')
process.RandomNumberGeneratorService.externalLHEProducer.initialSeed = 1111111
process.source = cms.Source(
"EmptySource",
firstLuminosityBlock=cms.untracked.uint32(23456),
numberEventsInLuminosityBlock=cms.untracked.uint32(10))
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(20000))
process.externalLHEProducer = cms.EDProducer(
'ExternalLHEProducer',
scriptName=cms.FileInPath(
"GeneratorInterface/LHEInterface/data/run_madgraph_tarball.sh"),
outputFile=cms.string(
"TprimeTprimeToBWTZ_M-600_8TeV-madgraph_unweighted_events_final.lhe"),
args=cms.vstring(
'slc5_ia32_gcc434/madgraph/V5_1.3.30/8TeV_Summer12/TprimeTprimetoInclusive/v1',
'TprimeTprimeToBWTZ_M-600_8TeV-madgraph', 'false', 'false', 'ttbar',
'5', '20', 'true', '0', '1'),
nEvents=cms.uint32(process.maxEvents.input.value()))
process.out = cms.OutputModule(
"PoolOutputModule",
splitLevel=cms.untracked.int32(0),
eventAutoFlushCompressedSize=cms.untracked.int32(5242880),
outputCommands=process.LHEEventContent.outputCommands,
fileName=cms.untracked.string('myOutputFile.root'),
dataset=cms.untracked.PSet(filterName=cms.untracked.string(''),
dataTier=cms.untracked.string('LHE')))
process.MessageLogger = cms.Service("MessageLogger",
destinations = cms.untracked.vstring('output'),
threshold = cms.untracked.string('INFO'),
noLineBreaks = cms.untracked.bool(True)
)
process.MessageLogger = cms.Service("MessageLogger",
destinations = cms.untracked.vstring('output'),
threshold = cms.untracked.string('ERROR'),
noLineBreaks = cms.untracked.bool(True)
)
process.MessageLogger = cms.Service("MessageLogger",
destinations = cms.untracked.vstring('output'),
threshold = cms.untracked.string('WARNING'),
noLineBreaks = cms.untracked.bool(True)
)
process.filler = cms.EDProducer("SiPixelFedFillerWordEventNumber",
InputLabel = cms.untracked.string('source'),
InputInstance = cms.untracked.string(''),
SaveFillerWords = cms.bool(False)
)
process.out = cms.EDProducer("PoolOutputModule",
fileName = cms.untracked.string("NewFEDFWs.root")
)
process.producer = cms.Path(process.filler)
process.producer = cms.EndPath(process.out)
#No edmtest::IntProduct's added to this output file
import FWCore.ParameterSet.Config as cms
process = cms.Process("PROD")
process.source = cms.Source("EmptySource",
firstLuminosityBlock=cms.untracked.uint32(2),
firstEvent=cms.untracked.uint32(20))
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(10))
process.c = cms.EDProducer("ThingProducer")
process.d = cms.EDProducer("OtherThingProducer", thingTag=cms.InputTag("c"))
process.o = cms.OutputModule(
"PoolOutputModule",
outputCommands=cms.untracked.vstring("drop *", "keep edmtestThings_*_*_*",
"keep edmtestOtherThings_*_*_*"),
fileName=cms.untracked.string("ref_merge_prod2.root"))
process.p = cms.Path(process.c * process.d)
process.tester = cms.EDAnalyzer("OtherThingAnalyzer",
other=cms.untracked.InputTag(
"d", "testUserTag"))
process.out = cms.EndPath(process.o + process.tester)
maxDeltaR = 0.6
)
akVsSoftDrop6PFparton = patJetPartonMatch.clone(src = cms.InputTag("akVsSoftDrop6PFJets")
)
akVsSoftDrop6PFcorr = patJetCorrFactors.clone(
useNPV = cms.bool(False),
useRho = cms.bool(False),
# primaryVertices = cms.InputTag("hiSelectedVertex"),
levels = cms.vstring('L2Relative','L3Absolute'),
src = cms.InputTag("akVsSoftDrop6PFJets"),
payload = "AK6PF_offline"
)
akVsSoftDrop6PFJetID= cms.EDProducer('JetIDProducer', JetIDParams, src = cms.InputTag('akVsSoftDrop6CaloJets'))
#akVsSoftDrop6PFclean = heavyIonCleanedGenJets.clone(src = cms.InputTag('ak6GenJets'))
akVsSoftDrop6PFbTagger = bTaggers("akVsSoftDrop6PF",0.6)
#create objects locally since they dont load properly otherwise
#akVsSoftDrop6PFmatch = akVsSoftDrop6PFbTagger.match
akVsSoftDrop6PFparton = patJetPartonMatch.clone(src = cms.InputTag("akVsSoftDrop6PFJets"), matched = cms.InputTag("genParticles"))
akVsSoftDrop6PFPatJetFlavourAssociationLegacy = akVsSoftDrop6PFbTagger.PatJetFlavourAssociationLegacy
akVsSoftDrop6PFPatJetPartons = akVsSoftDrop6PFbTagger.PatJetPartons
akVsSoftDrop6PFJetTracksAssociatorAtVertex = akVsSoftDrop6PFbTagger.JetTracksAssociatorAtVertex
akVsSoftDrop6PFJetTracksAssociatorAtVertex.tracks = cms.InputTag("highPurityTracks")
akVsSoftDrop6PFSimpleSecondaryVertexHighEffBJetTags = akVsSoftDrop6PFbTagger.SimpleSecondaryVertexHighEffBJetTags
akVsSoftDrop6PFSimpleSecondaryVertexHighPurBJetTags = akVsSoftDrop6PFbTagger.SimpleSecondaryVertexHighPurBJetTags
akVsSoftDrop6PFCombinedSecondaryVertexBJetTags = akVsSoftDrop6PFbTagger.CombinedSecondaryVertexBJetTags
import FWCore.ParameterSet.Config as cms
from MuonAnalysis.MuonAssociators.muonL1Match_cfi import *
from math import pi
muonL1MatchExtended = cms.EDProducer("L1MatcherExtended",
muons = cms.InputTag("muons"),
l1extra = cms.InputTag("l1extraParticles"),
segmentArbitration = cms.string("SegmentAndTrackArbitration"),
csctfDigis = cms.InputTag("csctfDigis"),
csctfLcts = cms.InputTag("csctfDigis"),
matcherGeom = cms.PSet(
preselection = cms.string("gmtMuonCand.quality > 1"), # FIXME: maybe exclude CSC-only region?
useTrack = cms.string("tracker"),
useState = cms.string("atVertex"),
maxDeltaR = cms.double(1.5), ## FIXME: to be tuned
maxDeltaEta = cms.double(0.3), ## FIXME: to be tuned
l1PhiOffset = cms.double(1.25 * pi/180.),
useSimpleGeometry = cms.bool(True),
fallbackToME1 = cms.bool(True),
)
)
def addUserData(patMuonProducer, matcherLabel='muonL1MatchExtended', addExtraInfo=False):
patMuonProducer.userData.userInts.src += [ cms.InputTag(matcherLabel) ]
if addExtraInfo:
for L in ("cscMode", "canPropagate", "l1q"):
patMuonProducer.userData.userInts.src += [ cms.InputTag(matcherLabel,L) ]
for L in ("deltaR", "deltaEta", "deltaPhi", "l1pt"):
patMuonProducer.userData.userFloats.src += [ cms.InputTag(matcherLabel,L) ]
'4900113:addChannel = 1 0.151 91 15 -15', '4900113:tau0 = 500',
'4900111:onMode = 0'),
parameterSets=cms.vstring('pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'pythia8PSweightsSettings',
'pythia8PowhegEmissionVetoSettings',
'processParameters')))
process.genfilter = cms.EDFilter(
"GenParticleSelector",
src=cms.InputTag("genParticlesForFilter"),
cut=cms.string('(pdgId==25) && pt>140. && status==62'))
process.genParticlesForFilter = cms.EDProducer(
"GenParticleProducer",
saveBarCodes=cms.untracked.bool(True),
src=cms.InputTag("generator"),
abortOnUnknownPDGCode=cms.untracked.bool(False))
process.externalLHEProducer = cms.EDProducer(
"ExternalLHEProducer",
nEvents=cms.untracked.uint32(10),
outputFile=cms.string('cmsgrid_final.lhe'),
scriptName=cms.FileInPath(
'GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh'),
numberOfParameters=cms.uint32(1),
args=cms.vstring(
'/cvmfs/cms.cern.ch/phys_generator/gridpacks/slc6_amd64_gcc481/13TeV/powheg/V2/gg_H_quark-mass-effects_NNPDF30_13TeV_M125/v2/gg_H_quark-mass-effects_NNPDF30_13TeV_M125_tarball.tar.gz'
))
process.ProductionFilterSequence = cms.Sequence(process.generator +
resolveByMatchQuality=cms.bool(True),
maxDeltaR=0.2)
akFlowPuCsSoftDropZ05B152PFparton = patJetPartonMatch.clone(
src=cms.InputTag("akFlowPuCsSoftDropZ05B152PFJets"),
matched=cms.InputTag("hiSignalGenParticles"))
akFlowPuCsSoftDropZ05B152PFcorr = patJetCorrFactors.clone(
useNPV=cms.bool(False),
useRho=cms.bool(False),
levels=cms.vstring('L2Relative'),
src=cms.InputTag("akFlowPuCsSoftDropZ05B152PFJets"),
payload="AK2PF")
akFlowPuCsSoftDropZ05B152PFJetID = cms.EDProducer(
'JetIDProducer',
JetIDParams,
src=cms.InputTag('akFlowPuCsSoftDropZ05B152CaloJets'))
# akFlowPuCsSoftDropZ05B152PFclean = heavyIonCleanedGenJets.clone(
# src = cms.InputTag('ak2HiSignalGenJets'))
akFlowPuCsSoftDropZ05B152PFbTagger = bTaggers("akFlowPuCsSoftDropZ05B152PF",
0.2)
# create objects locally since they dont load properly otherwise
akFlowPuCsSoftDropZ05B152PFPatJetFlavourAssociationLegacy = akFlowPuCsSoftDropZ05B152PFbTagger.PatJetFlavourAssociationLegacy
akFlowPuCsSoftDropZ05B152PFPatJetPartons = akFlowPuCsSoftDropZ05B152PFbTagger.PatJetPartons
akFlowPuCsSoftDropZ05B152PFJetTracksAssociatorAtVertex = akFlowPuCsSoftDropZ05B152PFbTagger.JetTracksAssociatorAtVertex
akFlowPuCsSoftDropZ05B152PFJetTracksAssociatorAtVertex.tracks = cms.InputTag(
"highPurityTracks")
akFlowPuCsSoftDropZ05B152PFSimpleSecondaryVertexHighEffBJetTags = akFlowPuCsSoftDropZ05B152PFbTagger.SimpleSecondaryVertexHighEffBJetTags
