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