Pythia8Interface/test/FXFXewk_test.py

0001 # Auto generated configuration file
0002 # using:
0003 # Revision: 1.19
0004 # Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
0005 # with command line options: GeneratorInterface/Pythia8Interface/python/fxfxEWK.py --python_filename cfg_newLHE_newMatching.py --eventcontent NANOAODSIM --customise Configuration/DataProcessing/Utils.addMonitoring --datatier NANOAOD --fileout file:NanoGen_newLHE_newMatching.root --conditions 106X_upgrade2018_realistic_v4 --beamspot Realistic25ns13TeVEarly2018Collision --customise_commands process.source.numberEventsInLuminosityBlock=cms.untracked.uint32(161) --step LHE,GEN,NANOGEN --geometry DB:Extended --era Run2_2018 --no_exec --mc -n 1000
0006 import FWCore.ParameterSet.Config as cms
0007
0008 from Configuration.Eras.Era_Run2_2018_cff import Run2_2018
0009
0010 process = cms.Process('NANOGEN',Run2_2018)
0011
0012 # import of standard configurations
0013 process.load('Configuration.StandardSequences.Services_cff')
0014 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0015 process.load('FWCore.MessageService.MessageLogger_cfi')
0016 process.load('Configuration.EventContent.EventContent_cff')
0017 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0018 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0019 process.load('Configuration.StandardSequences.MagneticField_cff')
0020 process.load('Configuration.StandardSequences.Generator_cff')
0021 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic25ns13TeVEarly2018Collision_cfi')
0022 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0023 process.load('PhysicsTools.NanoAOD.nanogen_cff')
0024 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0025 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0026
0027 process.maxEvents = cms.untracked.PSet(
0028     input = cms.untracked.int32(1000),
0029     output = cms.optional.untracked.allowed(cms.int32,cms.PSet)
0030 )
0031
0032 # Input source
0033 process.source = cms.Source("EmptySource")
0034
0035 process.options = cms.untracked.PSet(
0036     FailPath = cms.untracked.vstring(),
0037     IgnoreCompletely = cms.untracked.vstring(),
0038     Rethrow = cms.untracked.vstring(),
0039     SkipEvent = cms.untracked.vstring(),
0040     accelerators = cms.untracked.vstring('*'),
0041     allowUnscheduled = cms.obsolete.untracked.bool,
0042     canDeleteEarly = cms.untracked.vstring(),
0043     deleteNonConsumedUnscheduledModules = cms.untracked.bool(True),
0044     dumpOptions = cms.untracked.bool(False),
0045     emptyRunLumiMode = cms.obsolete.untracked.string,
0046     eventSetup = cms.untracked.PSet(
0047         forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0048             allowAnyLabel_=cms.required.untracked.uint32
0049         ),
0050         numberOfConcurrentIOVs = cms.untracked.uint32(0)
0051     ),
0052     fileMode = cms.untracked.string('FULLMERGE'),
0053     forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0054     makeTriggerResults = cms.obsolete.untracked.bool,
0055     numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(0),
0056     numberOfConcurrentRuns = cms.untracked.uint32(1),
0057     numberOfStreams = cms.untracked.uint32(0),
0058     numberOfThreads = cms.untracked.uint32(1),
0059     printDependencies = cms.untracked.bool(False),
0060     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0061     throwIfIllegalParameter = cms.untracked.bool(True),
0062     wantSummary = cms.untracked.bool(False)
0063 )
0064
0065 # Production Info
0066 process.configurationMetadata = cms.untracked.PSet(
0067     annotation = cms.untracked.string('GeneratorInterface/Pythia8Interface/python/fxfxEWK.py nevts:1000'),
0068     name = cms.untracked.string('Applications'),
0069     version = cms.untracked.string('$Revision: 1.19 $')
0070 )
0071
0072 # Output definition
0073
0074 process.NANOAODSIMoutput = cms.OutputModule("NanoAODOutputModule",
0075     SelectEvents = cms.untracked.PSet(
0076         SelectEvents = cms.vstring('generation_step')
0077     ),
0078     compressionAlgorithm = cms.untracked.string('LZMA'),
0079     compressionLevel = cms.untracked.int32(9),
0080     dataset = cms.untracked.PSet(
0081         dataTier = cms.untracked.string('NANOAOD'),
0082         filterName = cms.untracked.string('')
0083     ),
0084     fileName = cms.untracked.string('file:NanoGen_newLHE_newMatching.root'),
0085     outputCommands = process.NANOAODSIMEventContent.outputCommands
0086 )
0087
0088 # Additional output definition
0089
0090 # Other statements
0091 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0092 from Configuration.AlCa.GlobalTag import GlobalTag
0093 process.GlobalTag = GlobalTag(process.GlobalTag, '106X_upgrade2018_realistic_v4', '')
0094
0095 process.generator = cms.EDFilter("Pythia8HadronizerFilter",
0096     PythiaParameters = cms.PSet(
0097         parameterSets = cms.vstring(
0098             'pythia8CommonSettings',
0099             'pythia8CP5Settings',
0100             'pythia8aMCatNLOSettings',
0101             'pythia8PSweightsSettings',
0102             'processParameters'
0103         ),
0104         processParameters = cms.vstring(
0105             'JetMatching:setMad = off',
0106             'JetMatching:scheme = 1',
0107             'JetMatching:doVeto = off',
0108             'JetMatching:merge = on',
0109             'JetMatching:jetAlgorithm = 2',
0110             'JetMatching:etaJetMax = 1000',
0111             'JetMatching:coneRadius = 1.',
0112             'JetMatching:slowJetPower = 1',
0113             'JetMatching:qCut = 42.',
0114             'JetMatching:doFxFx = on',
0115             'JetMatching:qCutME = 20.',
0116             'JetMatching:nQmatch = 5',
0117             'JetMatching:nJetMax = 1'
0118         ),
0119         pythia8CP5Settings = cms.vstring(
0120             'Tune:pp 14',
0121             'Tune:ee 7',
0122             'MultipartonInteractions:ecmPow=0.03344',
0123             'MultipartonInteractions:bProfile=2',
0124             'MultipartonInteractions:pT0Ref=1.41',
0125             'MultipartonInteractions:coreRadius=0.7634',
0126             'MultipartonInteractions:coreFraction=0.63',
0127             'ColourReconnection:range=5.176',
0128             'SigmaTotal:zeroAXB=off',
0129             'SpaceShower:alphaSorder=2',
0130             'SpaceShower:alphaSvalue=0.118',
0131             'SigmaProcess:alphaSvalue=0.118',
0132             'SigmaProcess:alphaSorder=2',
0133             'MultipartonInteractions:alphaSvalue=0.118',
0134             'MultipartonInteractions:alphaSorder=2',
0135             'TimeShower:alphaSorder=2',
0136             'TimeShower:alphaSvalue=0.118',
0137             'SigmaTotal:mode = 0',
0138             'SigmaTotal:sigmaEl = 21.89',
0139             'SigmaTotal:sigmaTot = 100.309',
0140             'PDF:pSet=LHAPDF6:NNPDF31_nnlo_as_0118'
0141         ),
0142         pythia8CommonSettings = cms.vstring(
0143             'Tune:preferLHAPDF = 2',
0144             'Main:timesAllowErrors = 10000',
0145             'Check:epTolErr = 0.01',
0146             'Beams:setProductionScalesFromLHEF = off',
0147             'SLHA:minMassSM = 1000.',
0148             'ParticleDecays:limitTau0 = on',
0149             'ParticleDecays:tau0Max = 10',
0150             'ParticleDecays:allowPhotonRadiation = on'
0151         ),
0152         pythia8PSweightsSettings = cms.vstring(
0153             'UncertaintyBands:doVariations = on',
0154             'UncertaintyBands:List = {isrRedHi isr:muRfac=0.707,fsrRedHi fsr:muRfac=0.707,isrRedLo isr:muRfac=1.414,fsrRedLo fsr:muRfac=1.414,isrDefHi isr:muRfac=0.5,fsrDefHi fsr:muRfac=0.5,isrDefLo isr:muRfac=2.0,fsrDefLo fsr:muRfac=2.0,isrConHi isr:muRfac=0.25,fsrConHi fsr:muRfac=0.25,isrConLo isr:muRfac=4.0,fsrConLo fsr:muRfac=4.0,fsr_G2GG_muR_dn fsr:G2GG:muRfac=0.5,fsr_G2GG_muR_up fsr:G2GG:muRfac=2.0,fsr_G2QQ_muR_dn fsr:G2QQ:muRfac=0.5,fsr_G2QQ_muR_up fsr:G2QQ:muRfac=2.0,fsr_Q2QG_muR_dn fsr:Q2QG:muRfac=0.5,fsr_Q2QG_muR_up fsr:Q2QG:muRfac=2.0,fsr_X2XG_muR_dn fsr:X2XG:muRfac=0.5,fsr_X2XG_muR_up fsr:X2XG:muRfac=2.0,fsr_G2GG_cNS_dn fsr:G2GG:cNS=-2.0,fsr_G2GG_cNS_up fsr:G2GG:cNS=2.0,fsr_G2QQ_cNS_dn fsr:G2QQ:cNS=-2.0,fsr_G2QQ_cNS_up fsr:G2QQ:cNS=2.0,fsr_Q2QG_cNS_dn fsr:Q2QG:cNS=-2.0,fsr_Q2QG_cNS_up fsr:Q2QG:cNS=2.0,fsr_X2XG_cNS_dn fsr:X2XG:cNS=-2.0,fsr_X2XG_cNS_up fsr:X2XG:cNS=2.0,isr_G2GG_muR_dn isr:G2GG:muRfac=0.5,isr_G2GG_muR_up isr:G2GG:muRfac=2.0,isr_G2QQ_muR_dn isr:G2QQ:muRfac=0.5,isr_G2QQ_muR_up isr:G2QQ:muRfac=2.0,isr_Q2QG_muR_dn isr:Q2QG:muRfac=0.5,isr_Q2QG_muR_up isr:Q2QG:muRfac=2.0,isr_X2XG_muR_dn isr:X2XG:muRfac=0.5,isr_X2XG_muR_up isr:X2XG:muRfac=2.0,isr_G2GG_cNS_dn isr:G2GG:cNS=-2.0,isr_G2GG_cNS_up isr:G2GG:cNS=2.0,isr_G2QQ_cNS_dn isr:G2QQ:cNS=-2.0,isr_G2QQ_cNS_up isr:G2QQ:cNS=2.0,isr_Q2QG_cNS_dn isr:Q2QG:cNS=-2.0,isr_Q2QG_cNS_up isr:Q2QG:cNS=2.0,isr_X2XG_cNS_dn isr:X2XG:cNS=-2.0,isr_X2XG_cNS_up isr:X2XG:cNS=2.0}',
0155             'UncertaintyBands:nFlavQ = 4',
0156             'UncertaintyBands:MPIshowers = on',
0157             'UncertaintyBands:overSampleFSR = 10.0',
0158             'UncertaintyBands:overSampleISR = 10.0',
0159             'UncertaintyBands:FSRpTmin2Fac = 20',
0160             'UncertaintyBands:ISRpTmin2Fac = 1'
0161         ),
0162         pythia8aMCatNLOSettings = cms.vstring(
0163             'SpaceShower:pTmaxMatch = 1',
0164             'SpaceShower:pTmaxFudge = 1',
0165             'SpaceShower:MEcorrections = off',
0166             'TimeShower:pTmaxMatch = 1',
0167             'TimeShower:pTmaxFudge = 1',
0168             'TimeShower:MEcorrections = off',
0169             'TimeShower:globalRecoil = on',
0170             'TimeShower:limitPTmaxGlobal = on',
0171             'TimeShower:nMaxGlobalRecoil = 1',
0172             'TimeShower:globalRecoilMode = 2',
0173             'TimeShower:nMaxGlobalBranch = 1',
0174             'TimeShower:weightGluonToQuark = 1'
0175         )
0176     ),
0177     UserCustomization = cms.VPSet(cms.PSet(
0178         pluginName = cms.string('JetMatchingEWKFxFx')
0179     )),
0180     comEnergy = cms.double(13000.0),
0181     filterEfficiency = cms.untracked.double(1.0),
0182     maxEventsToPrint = cms.untracked.int32(1),
0183     pythiaHepMCVerbosity = cms.untracked.bool(False),
0184     pythiaPylistVerbosity = cms.untracked.int32(1)
0185 )
0186
0187
0188 process.externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
0189     args = cms.vstring('/afs/cern.ch/user/c/cvicovil/public/TTWJetsToLNu_5f_NLO_FXFX_slc7_amd64_gcc900_CMSSW_12_0_2_tarball.tar.xz'),
0190     nEvents = cms.untracked.uint32(1000),
0191     numberOfParameters = cms.uint32(1),
0192     outputFile = cms.string('cmsgrid_final.lhe'),
0193     scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh')
0194 )
0195
0196
0197 # Path and EndPath definitions
0198 process.lhe_step = cms.Path(process.externalLHEProducer)
0199 process.generation_step = cms.Path(process.pgen)
0200 process.nanoAOD_step = cms.Path(process.nanogenSequence)
0201 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0202 process.endjob_step = cms.EndPath(process.endOfProcess)
0203 process.NANOAODSIMoutput_step = cms.EndPath(process.NANOAODSIMoutput)
0204
0205 # Schedule definition
0206 process.schedule = cms.Schedule(process.lhe_step,process.generation_step,process.genfiltersummary_step,process.nanoAOD_step,process.endjob_step,process.NANOAODSIMoutput_step)
0207 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0208 associatePatAlgosToolsTask(process)
0209
0210 #Setup FWK for multithreaded
0211 process.options.numberOfConcurrentLuminosityBlocks = 1
0212 process.options.eventSetup.numberOfConcurrentIOVs = 1
0213 # filter all path with the production filter sequence
0214 for path in process.paths:
0215     if path in ['lhe_step']: continue
0216     getattr(process,path).insert(0, process.generator)
0217
0218 # customisation of the process.
0219
0220 # Automatic addition of the customisation function from PhysicsTools.NanoAOD.nanogen_cff
0221 from PhysicsTools.NanoAOD.nanogen_cff import customizeNanoGEN
0222
0223 #call to customisation function customizeNanoGEN imported from PhysicsTools.NanoAOD.nanogen_cff
0224 process = customizeNanoGEN(process)
0225
0226 # Automatic addition of the customisation function from Configuration.DataProcessing.Utils
0227 from Configuration.DataProcessing.Utils import addMonitoring
0228
0229 #call to customisation function addMonitoring imported from Configuration.DataProcessing.Utils
0230 process = addMonitoring(process)
0231
0232 # End of customisation functions
0233
0234
0235 # Customisation from command line
0236
0237 process.source.numberEventsInLuminosityBlock=cms.untracked.uint32(161)
0238 # Add early deletion of temporary data products to reduce peak memory need
0239 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0240 process = customiseEarlyDelete(process)
0241 # End adding early deletion