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())