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     IgnoreCompletely = cms.untracked.vstring(),
0037     Rethrow = cms.untracked.vstring(),
0038     TryToContinue = cms.untracked.vstring(),
0039     accelerators = cms.untracked.vstring('*'),
0040     allowUnscheduled = cms.obsolete.untracked.bool,
0041     canDeleteEarly = cms.untracked.vstring(),
0042     deleteNonConsumedUnscheduledModules = cms.untracked.bool(True),
0043     dumpOptions = cms.untracked.bool(False),
0044     emptyRunLumiMode = cms.obsolete.untracked.string,
0045     eventSetup = cms.untracked.PSet(
0046         forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0047             allowAnyLabel_=cms.required.untracked.uint32
0048         ),
0049         numberOfConcurrentIOVs = cms.untracked.uint32(0)
0050     ),
0051     fileMode = cms.untracked.string('FULLMERGE'),
0052     forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0053     makeTriggerResults = cms.obsolete.untracked.bool,
0054     numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(0),
0055     numberOfConcurrentRuns = cms.untracked.uint32(1),
0056     numberOfStreams = cms.untracked.uint32(0),
0057     numberOfThreads = cms.untracked.uint32(1),
0058     printDependencies = cms.untracked.bool(False),
0059     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0060     throwIfIllegalParameter = cms.untracked.bool(True),
0061     wantSummary = cms.untracked.bool(False)
0062 )
0063
0064 # Production Info
0065 process.configurationMetadata = cms.untracked.PSet(
0066     annotation = cms.untracked.string('GeneratorInterface/Pythia8Interface/python/fxfxEWK.py nevts:1000'),
0067     name = cms.untracked.string('Applications'),
0068     version = cms.untracked.string('$Revision: 1.19 $')
0069 )
0070
0071 # Output definition
0072
0073 process.NANOAODSIMoutput = cms.OutputModule("NanoAODOutputModule",
0074     SelectEvents = cms.untracked.PSet(
0075         SelectEvents = cms.vstring('generation_step')
0076     ),
0077     compressionAlgorithm = cms.untracked.string('LZMA'),
0078     compressionLevel = cms.untracked.int32(9),
0079     dataset = cms.untracked.PSet(
0080         dataTier = cms.untracked.string('NANOAOD'),
0081         filterName = cms.untracked.string('')
0082     ),
0083     fileName = cms.untracked.string('file:NanoGen_newLHE_newMatching.root'),
0084     outputCommands = process.NANOAODSIMEventContent.outputCommands
0085 )
0086
0087 # Additional output definition
0088
0089 # Other statements
0090 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0091 from Configuration.AlCa.GlobalTag import GlobalTag
0092 process.GlobalTag = GlobalTag(process.GlobalTag, '106X_upgrade2018_realistic_v4', '')
0093
0094 process.generator = cms.EDFilter("Pythia8HadronizerFilter",
0095     PythiaParameters = cms.PSet(
0096         parameterSets = cms.vstring(
0097             'pythia8CommonSettings',
0098             'pythia8CP5Settings',
0099             'pythia8aMCatNLOSettings',
0100             'pythia8PSweightsSettings',
0101             'processParameters'
0102         ),
0103         processParameters = cms.vstring(
0104             'JetMatching:setMad = off',
0105             'JetMatching:scheme = 1',
0106             'JetMatching:doVeto = off',
0107             'JetMatching:merge = on',
0108             'JetMatching:jetAlgorithm = 2',
0109             'JetMatching:etaJetMax = 1000',
0110             'JetMatching:coneRadius = 1.',
0111             'JetMatching:slowJetPower = 1',
0112             'JetMatching:qCut = 42.',
0113             'JetMatching:doFxFx = on',
0114             'JetMatching:qCutME = 20.',
0115             'JetMatching:nQmatch = 5',
0116             'JetMatching:nJetMax = 1'
0117         ),
0118         pythia8CP5Settings = cms.vstring(
0119             'Tune:pp 14',
0120             'Tune:ee 7',
0121             'MultipartonInteractions:ecmPow=0.03344',
0122             'MultipartonInteractions:bProfile=2',
0123             'MultipartonInteractions:pT0Ref=1.41',
0124             'MultipartonInteractions:coreRadius=0.7634',
0125             'MultipartonInteractions:coreFraction=0.63',
0126             'ColourReconnection:range=5.176',
0127             'SigmaTotal:zeroAXB=off',
0128             'SpaceShower:alphaSorder=2',
0129             'SpaceShower:alphaSvalue=0.118',
0130             'SigmaProcess:alphaSvalue=0.118',
0131             'SigmaProcess:alphaSorder=2',
0132             'MultipartonInteractions:alphaSvalue=0.118',
0133             'MultipartonInteractions:alphaSorder=2',
0134             'TimeShower:alphaSorder=2',
0135             'TimeShower:alphaSvalue=0.118',
0136             'SigmaTotal:mode = 0',
0137             'SigmaTotal:sigmaEl = 21.89',
0138             'SigmaTotal:sigmaTot = 100.309',
0139             'PDF:pSet=LHAPDF6:NNPDF31_nnlo_as_0118'
0140         ),
0141         pythia8CommonSettings = cms.vstring(
0142             'Tune:preferLHAPDF = 2',
0143             'Main:timesAllowErrors = 10000',
0144             'Check:epTolErr = 0.01',
0145             'Beams:setProductionScalesFromLHEF = off',
0146             'SLHA:minMassSM = 1000.',
0147             'ParticleDecays:limitTau0 = on',
0148             'ParticleDecays:tau0Max = 10',
0149             'ParticleDecays:allowPhotonRadiation = on'
0150         ),
0151         pythia8PSweightsSettings = cms.vstring(
0152             'UncertaintyBands:doVariations = on',
0153             '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}',
0154             'UncertaintyBands:nFlavQ = 4',
0155             'UncertaintyBands:MPIshowers = on',
0156             'UncertaintyBands:overSampleFSR = 10.0',
0157             'UncertaintyBands:overSampleISR = 10.0',
0158             'UncertaintyBands:FSRpTmin2Fac = 20',
0159             'UncertaintyBands:ISRpTmin2Fac = 1'
0160         ),
0161         pythia8aMCatNLOSettings = cms.vstring(
0162             'SpaceShower:pTmaxMatch = 1',
0163             'SpaceShower:pTmaxFudge = 1',
0164             'SpaceShower:MEcorrections = off',
0165             'TimeShower:pTmaxMatch = 1',
0166             'TimeShower:pTmaxFudge = 1',
0167             'TimeShower:MEcorrections = off',
0168             'TimeShower:globalRecoil = on',
0169             'TimeShower:limitPTmaxGlobal = on',
0170             'TimeShower:nMaxGlobalRecoil = 1',
0171             'TimeShower:globalRecoilMode = 2',
0172             'TimeShower:nMaxGlobalBranch = 1',
0173             'TimeShower:weightGluonToQuark = 1'
0174         )
0175     ),
0176     UserCustomization = cms.VPSet(cms.PSet(
0177         pluginName = cms.string('JetMatchingEWKFxFx')
0178     )),
0179     comEnergy = cms.double(13000.0),
0180     filterEfficiency = cms.untracked.double(1.0),
0181     maxEventsToPrint = cms.untracked.int32(1),
0182     pythiaHepMCVerbosity = cms.untracked.bool(False),
0183     pythiaPylistVerbosity = cms.untracked.int32(1)
0184 )
0185
0186
0187 process.externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
0188     args = cms.vstring('/afs/cern.ch/user/c/cvicovil/public/TTWJetsToLNu_5f_NLO_FXFX_slc7_amd64_gcc900_CMSSW_12_0_2_tarball.tar.xz'),
0189     nEvents = cms.untracked.uint32(1000),
0190     numberOfParameters = cms.uint32(1),
0191     outputFile = cms.string('cmsgrid_final.lhe'),
0192     scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh')
0193 )
0194
0195
0196 # Path and EndPath definitions
0197 process.lhe_step = cms.Path(process.externalLHEProducer)
0198 process.generation_step = cms.Path(process.pgen)
0199 process.nanoAOD_step = cms.Path(process.nanogenSequence)
0200 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0201 process.endjob_step = cms.EndPath(process.endOfProcess)
0202 process.NANOAODSIMoutput_step = cms.EndPath(process.NANOAODSIMoutput)
0203
0204 # Schedule definition
0205 process.schedule = cms.Schedule(process.lhe_step,process.generation_step,process.genfiltersummary_step,process.nanoAOD_step,process.endjob_step,process.NANOAODSIMoutput_step)
0206 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0207 associatePatAlgosToolsTask(process)
0208
0209 #Setup FWK for multithreaded
0210 process.options.numberOfConcurrentLuminosityBlocks = 1
0211 process.options.eventSetup.numberOfConcurrentIOVs = 1
0212 # filter all path with the production filter sequence
0213 for path in process.paths:
0214     if path in ['lhe_step']: continue
0215     getattr(process,path).insert(0, process.generator)
0216
0217 # customisation of the process.
0218
0219 # Automatic addition of the customisation function from PhysicsTools.NanoAOD.nanogen_cff
0220 from PhysicsTools.NanoAOD.nanogen_cff import customizeNanoGEN
0221
0222 #call to customisation function customizeNanoGEN imported from PhysicsTools.NanoAOD.nanogen_cff
0223 process = customizeNanoGEN(process)
0224
0225 # Automatic addition of the customisation function from Configuration.DataProcessing.Utils
0226 from Configuration.DataProcessing.Utils import addMonitoring
0227
0228 #call to customisation function addMonitoring imported from Configuration.DataProcessing.Utils
0229 process = addMonitoring(process)
0230
0231 # End of customisation functions
0232
0233
0234 # Customisation from command line
0235
0236 process.source.numberEventsInLuminosityBlock=cms.untracked.uint32(161)
0237 # Add early deletion of temporary data products to reduce peak memory need
0238 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0239 process = customiseEarlyDelete(process)
0240 # End adding early deletion