Pythia8Interface/test/Py8Had_mgmatching_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("Gen")
0004
0005 process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
0006
0007 process.source = cms.Source("LHESource",
0008     fileNames = cms.untracked.vstring('file:../../Pythia6Interface/test/ttbar_5flavours_xqcut20_10TeV.lhe')
0009     # fileNames = cms.untracked.vstring('file:/storage/local/data1/condor/mrenna/lhe/7TeV_Zbb_run45040_unweighted_events_qcut13_mgPostv2.lhe')
0010     # fileNames = cms.untracked.vstring('file:/uscmst1b_scratch/lpc1/3DayLifetime/recovery_for_julia/storage/local/data1/condor/mrenna/lhe/7TeV_ttbarjets_run621_unweighted_events_qcut40_mgPost.lhe')
0011     # fileNames = cms.untracked.vstring('file:/storage/local/data1/condor/mrenna/lhe/7TeV_avjets_run50000_unweighted_events_qcut15_mgPost.lhe')
0012     # fileNames = cms.untracked.vstring('file:/storage/local/data1/condor/mrenna/lhe/7TeV_zvv_200_HT_inf_run114000_unweighted_events_qcut20_mgPostv2.lhe')
0013     # fileNames = cms.untracked.vstring('/store/user/mrenna/7TeV_ZbbToLL_M_50_run1001to1018_3processes_unweighted_events.lhe')
0014 )
0015
0016
0017 # process.load("Configuration.Generator.Hadronizer_MgmMatchTune4C_7TeV_madgraph_pythia8_cff")
0018
0019 process.generator = cms.EDFilter("Pythia8HadronizerFilter",
0020     maxEventsToPrint = cms.untracked.int32(1),
0021     pythiaPylistVerbosity = cms.untracked.int32(1),
0022     filterEfficiency = cms.untracked.double(1.0),
0023     pythiaHepMCVerbosity = cms.untracked.bool(False),
0024     comEnergy = cms.double(7000.),
0025     jetMatching = cms.untracked.PSet(
0026        scheme = cms.string("Madgraph"),
0027        mode = cms.string("auto"),   # soup, or "inclusive"/"exclusive"
0028        #
0029        # ATTENTION PLEASE !
0030        # One can set some parameters to -1 to make the tool pock it up from LHE file.
0031        # However, -1 is ONLY possible if a givcen parameter is present in LHE file
0032        # - otherwise the code will throw.
0033        # So the user should make sure what it is and what she/he wants to do.
0034        #
0035        MEMAIN_etaclmax = cms.double(5.),
0036        MEMAIN_qcut = cms.double(30.),
0037        MEMAIN_minjets = cms.int32(-1),
0038        MEMAIN_maxjets = cms.int32(-1),
0039        MEMAIN_showerkt = cms.double(0),    # use 1=yes only for pt-ordered showers !
0040        MEMAIN_nqmatch = cms.int32(5),      # PID of the flavor until which the QCD radiation are kept in the matching procedure.
0041                                            # If nqmatch=4, then all showered partons from b's are NOT taken into account.
0042                            # In many cases the D=5
0043        MEMAIN_excres = cms.string(""),
0044        outTree_flag = cms.int32(0)         # 1=yes, write out the tree for future sanity check
0045     ),
0046     PythiaParameters = cms.PSet(
0047         pythia8_mg = cms.vstring(''), # this pset is for very initial testing
0048         # this pset below is actually used in large-scale (production-type) tests
0049     processParameters = cms.vstring(
0050             'Main:timesAllowErrors    = 10000',
0051         'ParticleDecays:limitTau0 = on',
0052             'ParticleDecays:tauMax = 10',
0053         # '15:onMode = off', # tmp turn off tau decays, to process av sample (crash in Tau::decay in Py8)
0054         'Tune:ee 3',
0055         'Tune:pp 5',
0056     'ParticleDecays:sophisticatedTau = 0' ),
0057         parameterSets = cms.vstring('processParameters')
0058     )
0059 )
0060
0061 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0062 process.MessageLogger = cms.Service("MessageLogger",
0063     cerr = cms.untracked.PSet(
0064         enable = cms.untracked.bool(False)
0065     ),
0066     cout = cms.untracked.PSet(
0067         default = cms.untracked.PSet(
0068             limit = cms.untracked.int32(2)
0069         ),
0070         enable = cms.untracked.bool(True)
0071     )
0072 )
0073
0074 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0075     generator = cms.PSet(
0076         initialSeed = cms.untracked.uint32(123456789),
0077     )
0078 )
0079
0080 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000) )
0081
0082 process.GEN = cms.OutputModule("PoolOutputModule",
0083     fileName = cms.untracked.string('Py8Had_mgmatching.root')
0084 )
0085
0086 process.p = cms.Path(process.generator)
0087 process.outpath = cms.EndPath(process.GEN)
0088
0089 #process.schedule = cms.Schedule(process.p, process.outpath)
0090 process.schedule = cms.Schedule(process.p)
0091