/
githubmirror
/
cmssw
Обзор
Документация
Войти
/
githubmirror
/
cmssw
Код
Запросы
0
Пакеты
0
Релизы
0
Аналитика
Безопасность
master
SimGeneral/MixingModule/test/testhimix_cfg.py
298 строк
21 KB
wmtan
Do not update HepMCProduct in place
15 сен 2015, 19:47
15 сен 2015, 19:47
f719245
Код
Авторство
О чём код?
import FWCore.ParameterSet.Config as cms process = cms.Process('HISIGNAL') # import of standard configurations process.load('Configuration.StandardSequences.Services_cff') process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi') process.load('FWCore.MessageService.MessageLogger_cfi') process.load('Configuration.StandardSequences.HiEventMixing_cff') process.load('Configuration.StandardSequences.GeometryExtended_cff') process.load('Configuration.StandardSequences.MagneticField_38T_cff') process.load('Configuration.StandardSequences.Generator_cff') process.load('Configuration.StandardSequences.VtxSmeared') process.load('SimGeneral.MixingModule.himixGEN_cff') process.load('Configuration.StandardSequences.Sim_cff') process.load('SimGeneral.MixingModule.himixSIMExtended_cff') process.load('Configuration.StandardSequences.Digi_cff') process.load('SimGeneral.MixingModule.himixDIGI_cff') process.load('Configuration.StandardSequences.SimL1Emulator_cff') process.load('Configuration.StandardSequences.DigiToRaw_cff') process.load('Configuration.StandardSequences.EndOfProcess_cff') process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff') process.load('Configuration.EventContent.EventContentHeavyIons_cff') process.configurationMetadata = cms.untracked.PSet( version = cms.untracked.string('$Revision: 1.168 $'), annotation = cms.untracked.string('Pyquen_ZeemumuJets_pt10_4TeV_cfi.py nevts:1'), name = cms.untracked.string('PyReleaseValidation') ) process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1) ) process.options = cms.untracked.PSet( ) # Input source process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring( 'file:/d101/yetkin/FCD18B2A-5A7B-E011-AC06-002618943916.root', # '/store/relval/CMSSW_4_2_3/RelValHydjetQ_MinBias_2760GeV/GEN-SIM-DIGI-RAW-HLTDEBUG/STARTHI42_V8-v2/0062/FCD18B2A-5A7B-E011-AC06-002618943916.root' ), inputCommands = cms.untracked.vstring('drop *', 'keep *_generator_*_*', 'keep *_g4SimHits_*_*'), dropDescendantsOfDroppedBranches = cms.untracked.bool(False) ) # Output definition process.output = cms.OutputModule("PoolOutputModule", splitLevel = cms.untracked.int32(0), outputCommands = process.RAWDEBUGEventContent.outputCommands, fileName = cms.untracked.string('Pyquen_ZeemumuJets_pt10_4TeV_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW.root'), dataset = cms.untracked.PSet( dataTier = cms.untracked.string('GEN-SIM-RAW'), filterName = cms.untracked.string('') ), SelectEvents = cms.untracked.PSet( SelectEvents = cms.vstring('generation_step') ) ) # Additional output definition # Other statements process.GlobalTag.globaltag = 'START43_V3::All' process.hiSignal = cms.EDFilter("PyquenGeneratorFilter", aBeamTarget = cms.double(208.0), comEnergy = cms.double(4000.0), qgpInitialTemperature = cms.double(1.0), doCollisionalEnLoss = cms.bool(False), qgpNumQuarkFlavor = cms.int32(0), qgpProperTimeFormation = cms.double(0.1), numQuarkFlavor = cms.int32(0), hadronFreezoutTemperature = cms.double(0.14), doRadiativeEnLoss = cms.bool(True), backgroundLabel = cms.InputTag("generator","unsmeared"), embeddingMode = cms.bool(True), angularSpectrumSelector = cms.int32(0), doIsospin = cms.bool(True), doQuench = cms.bool(True), cFlag = cms.int32(0), bFixed = cms.double(0.0), bMin = cms.double(0.0), bMax = cms.double(0.0), PythiaParameters = cms.PSet( pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ', 'BRAT(1035) = 1 ', 'BRAT(1036) = 0 ', 'BRAT(1037) = 0 ', 'BRAT(1038) = 0 ', 'BRAT(1039) = 0 ', 'BRAT(1040) = 0 ', 'BRAT(1041) = 0 ', 'BRAT(1042) = 0 ', 'MDME(1034,1) = 0 ', 'MDME(1035,1) = 1 ', 'MDME(1036,1) = 0 ', 'MDME(1037,1) = 0 ', 'MDME(1038,1) = 0 ', 'MDME(1039,1) = 0 ', 'MDME(1040,1) = 0 ', 'MDME(1041,1) = 0 ', 'MDME(1042,1) = 0 '), myParameters = cms.vstring(), ppDefault = cms.vstring('MSEL=1 ! QCD hight pT processes', 'CKIN(3)=7.', 'MSTP(81)=0'), pythiaZtoElectrons = cms.vstring('MDME(174,1)=0', 'MDME(175,1)=0', 'MDME(176,1)=0', 'MDME(177,1)=0', 'MDME(178,1)=0', 'MDME(179,1)=0', 'MDME(182,1)=1', 'MDME(183,1)=0', 'MDME(184,1)=0', 'MDME(185,1)=0', 'MDME(186,1)=0', 'MDME(187,1)=0'), pythiaZjets = cms.vstring('MSUB(15)=1', 'MSUB(30)=1'), pythiaUESettings = cms.vstring('MSTJ(11)=3 ! Choice of the fragmentation function', 'MSTJ(22)=2 ! Decay those unstable particles', 'PARJ(71)=10 . ! for which ctau 10 mm', 'MSTP(2)=1 ! which order running alphaS', 'MSTP(33)=0 ! no K factors in hard cross sections', 'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)', 'MSTP(52)=2 ! work with LHAPDF', 'MSTP(81)=1 ! multiple parton interactions 1 is Pythia default', 'MSTP(82)=4 ! Defines the multi-parton model', 'MSTU(21)=1 ! Check on possible errors during program execution', 'PARP(82)=1.8387 ! pt cutoff for multiparton interactions', 'PARP(89)=1960. ! sqrts for which PARP82 is set', 'PARP(83)=0.5 ! Multiple interactions: matter distrbn parameter', 'PARP(84)=0.4 ! Multiple interactions: matter distribution parameter', 'PARP(90)=0.16 ! Multiple interactions: rescaling power', 'PARP(67)=2.5 ! amount of initial-state radiation', 'PARP(85)=1.0 ! gluon prod. mechanism in MI', 'PARP(86)=1.0 ! gluon prod. mechanism in MI', 'PARP(62)=1.25 ! ', 'PARP(64)=0.2 ! ', 'MSTP(91)=1 !', 'PARP(91)=2.1 ! kt distribution', 'PARP(93)=15.0 ! '), pythiaZtoMuonsAndElectrons = cms.vstring('MDME(174,1)=0', 'MDME(175,1)=0', 'MDME(176,1)=0', 'MDME(177,1)=0', 'MDME(178,1)=0', 'MDME(179,1)=0', 'MDME(182,1)=1', 'MDME(183,1)=0', 'MDME(184,1)=1', 'MDME(185,1)=0', 'MDME(186,1)=0', 'MDME(187,1)=0'), pythiaPromptPhotons = cms.vstring('MSUB(14)=1', 'MSUB(18)=1', 'MSUB(29)=1', 'MSUB(114)=1', 'MSUB(115)=1'), pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1', 'MSUB(422) = 1', 'MSUB(423) = 1', 'MSUB(424) = 1', 'MSUB(425) = 1', 'MSUB(426) = 1', 'MSUB(427) = 1', 'MSUB(428) = 1', 'MSUB(429) = 1', 'MSUB(430) = 1', 'MSUB(431) = 1', 'MSUB(432) = 1', 'MSUB(433) = 1', 'MSUB(434) = 1', 'MSUB(435) = 1', 'MSUB(436) = 1', 'MSUB(437) = 1', 'MSUB(438) = 1', 'MSUB(439) = 1'), pythiaMuonCandidates = cms.vstring('CKIN(3)=20', 'MSTJ(22)=2', 'PARJ(71)=40.'), pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16', 'PARP(142)=0.0119', 'PARP(143)=0.01', 'PARP(144)=0.01', 'PARP(145)=0.05', 'PARP(146)=9.28', 'PARP(147)=0.15', 'PARP(148)=0.02', 'PARP(149)=0.02', 'PARP(150)=0.085', 'PARJ(13)=0.60', 'PARJ(14)=0.162', 'PARJ(15)=0.018', 'PARJ(16)=0.054', 'MSTP(145)=0', 'MSTP(146)=0', 'MSTP(147)=0', 'MSTP(148)=1', 'MSTP(149)=1', 'BRAT(861)=0.202', 'BRAT(862)=0.798', 'BRAT(1501)=0.013', 'BRAT(1502)=0.987', 'BRAT(1555)=0.356', 'BRAT(1556)=0.644'), pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1', 'MSUB(462) = 1', 'MSUB(463) = 1', 'MSUB(464) = 1', 'MSUB(465) = 1', 'MSUB(466) = 1', 'MSUB(467) = 1', 'MSUB(468) = 1', 'MSUB(469) = 1', 'MSUB(470) = 1', 'MSUB(471) = 1', 'MSUB(472) = 1', 'MSUB(473) = 1', 'MSUB(474) = 1', 'MSUB(475) = 1', 'MSUB(476) = 1', 'MSUB(477) = 1', 'MSUB(478) = 1', 'MSUB(479) = 1'), pythiaWeakBosons = cms.vstring('MSUB(1)=1', 'MSUB(2)=1'), pythiaJets = cms.vstring('MSUB(11)=1', 'MSUB(12)=1', 'MSUB(13)=1', 'MSUB(28)=1', 'MSUB(53)=1', 'MSUB(68)=1'), pythiaZtoMuons = cms.vstring('MDME(174,1)=0', 'MDME(175,1)=0', 'MDME(176,1)=0', 'MDME(177,1)=0', 'MDME(178,1)=0', 'MDME(179,1)=0', 'MDME(182,1)=0', 'MDME(183,1)=0', 'MDME(184,1)=1', 'MDME(185,1)=0', 'MDME(186,1)=0', 'MDME(187,1)=0'), ppJets = cms.vstring('MSEL=1 ! QCD hight pT processes'), pythiaHirootDefault = cms.vstring('MSEL=0', 'MSTU(21)=1', 'PARU(14)=1.', 'MSTP(81)=0', 'PMAS(5,1)=4.8', 'PMAS(6,1)=175.0', 'CKIN(3)=7.', 'MSTJ(22)=2', 'PARJ(71)=10.', 'PARP(67)=1.', 'PARP(82)=1.9', 'PARP(85)=0.33', 'PARP(86)=0.66', 'PARP(89)=1000.', 'PARP(91)=1.0', 'MSTJ(11)=3', 'MSTJ(22)=2'), pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ', 'BRAT(859) = 1 ', 'BRAT(860) = 0 ', 'MDME(858,1) = 0 ', 'MDME(859,1) = 1 ', 'MDME(860,1) = 0 '), parameterSets = cms.vstring('pythiaUESettings', 'pythiaZjets', 'pythiaZtoMuonsAndElectrons', 'kinematics'), kinematics = cms.vstring('CKIN(3)=10', 'CKIN(4)=9999', 'CKIN(7)=-2.', 'CKIN(8)=2.') ) ) # Path and EndPath definitions process.generation_step = cms.Path(process.pgen_himix) process.simulation_step = cms.Path(process.psim) process.digitisation_step = cms.Path(process.pdigi) process.L1simulation_step = cms.Path(process.SimL1Emulator) process.digi2raw_step = cms.Path(process.DigiToRaw) process.endjob_step = cms.Path(process.endOfProcess) process.out_step = cms.EndPath(process.output) # Schedule definition process.schedule = cms.Schedule(process.generation_step,process.simulation_step,process.digitisation_step,process.L1simulation_step,process.digi2raw_step,process.endjob_step,process.out_step) # special treatment in case of production filter sequence for path in process.paths: getattr(process,path)._seq = process.hiSignal*getattr(process,path)._seq