TauolaInterface/test/H130GGgluonfusion_8TeV_Tauola_TauSpinner_Example_cfi.py

0001 # Auto generated configuration file
0002 # using:
0003 # Revision: 1.381.2.28
0004 # Source: /local/reps/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
0005 # with command line options: H130GGgluonfusion_8TeV_tauola_cfi.py --conditions auto:startup -s GEN,VALIDATION:genvalid_all --datatier GEN --relval 1000000,20000 -n 1000 --eventcontent RAWSIM
0006 import FWCore.ParameterSet.Config as cms
0007
0008 process = cms.Process('VALIDATION')
0009
0010 # import of standard configurations
0011 process.load('Configuration.StandardSequences.Services_cff')
0012 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0013 process.load('FWCore.MessageService.MessageLogger_cfi')
0014 process.load('Configuration.EventContent.EventContent_cff')
0015 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0016 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0017 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0018 process.load('Configuration.StandardSequences.Generator_cff')
0019 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
0020 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0021 process.load('Configuration.StandardSequences.Validation_cff')
0022 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0023 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0024 process.load("GeneratorInterface.TauolaInterface.TauSpinner_cfi")
0025 process.load("GeneratorInterface.TauolaInterface.TauSpinnerFilter_cfi")
0026
0027 process.maxEvents = cms.untracked.PSet(
0028     input = cms.untracked.int32(10000)
0029 )
0030
0031 # Input source
0032 process.source = cms.Source("EmptySource")
0033
0034 process.options = cms.untracked.PSet(
0035
0036 )
0037
0038 # Production Info
0039 process.configurationMetadata = cms.untracked.PSet(
0040     version = cms.untracked.string('$Revision: 1.381.2.28 $'),
0041     annotation = cms.untracked.string('H130GGgluonfusion_8TeV_tauola_cfi.py nevts:1000'),
0042     name = cms.untracked.string('PyReleaseValidation')
0043 )
0044
0045 # Output definition
0046
0047 process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
0048     splitLevel = cms.untracked.int32(0),
0049     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0050     outputCommands = process.RAWSIMEventContent.outputCommands,
0051     fileName = cms.untracked.string('file:step1.root'),
0052     dataset = cms.untracked.PSet(
0053         filterName = cms.untracked.string(''),
0054         dataTier = cms.untracked.string('GEN')
0055     ),
0056     SelectEvents = cms.untracked.PSet(
0057         SelectEvents = cms.vstring('generation_step')
0058     )
0059 )
0060
0061 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0062                            generator = cms.PSet(initialSeed = cms.untracked.uint32(12345)),
0063                            TauSpinnerGen  = cms.PSet(initialSeed = cms.untracked.uint32(12345)),
0064                            TauSpinnerZHFilter = cms.PSet(initialSeed = cms.untracked.uint32(429842)),
0065                            VtxSmeared = cms.PSet(initialSeed = cms.untracked.uint32(275744))
0066                            )
0067
0068
0069
0070 # Other statements
0071 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0072 process.mix.playback = True
0073 from Configuration.AlCa.GlobalTag import GlobalTag
0074 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:startup', '')
0075
0076
0077 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0078                  ExternalDecays = cms.PSet(
0079         Tauola = cms.untracked.PSet(
0080     UseTauolaPolarization = cms.bool(False),
0081     InputCards = cms.PSet(
0082     mdtau = cms.int32(0),
0083     pjak2 = cms.int32(4),
0084     pjak1 = cms.int32(4)
0085     )
0086         ),
0087         parameterSets = cms.vstring('Tauola')
0088     ),
0089                  maxEventsToPrint = cms.untracked.int32(3),
0090                  pythiaPylistVerbosity = cms.untracked.int32(1),
0091                  filterEfficiency = cms.untracked.double(1.0),
0092                  pythiaHepMCVerbosity = cms.untracked.bool(False),
0093                  comEnergy = cms.double(13000.0),
0094                  crossSection = cms.untracked.double(0.05),
0095                  UseExternalGenerators = cms.untracked.bool(True),
0096                  PythiaParameters = cms.PSet(
0097         pythiaUESettings = cms.vstring('MSTU(21)=1     ! Check on possible errors during program execution',
0098                        'MSTJ(22)=2     ! Decay those unstable particles',
0099                        'PARJ(71)=10 .  ! for which ctau  10 mm',
0100                        'MSTP(33)=0     ! no K factors in hard cross sections',
0101                        'MSTP(2)=1      ! which order running alphaS',
0102                        'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
0103                        'MSTP(52)=2     ! work with LHAPDF',
0104                        'PARP(82)=1.921 ! pt cutoff for multiparton interactions',
0105                        'PARP(89)=1800. ! sqrts for which PARP82 is set',
0106                        'PARP(90)=0.227 ! Multiple interactions: rescaling power',
0107                        'MSTP(95)=6     ! CR (color reconnection parameters)',
0108                        'PARP(77)=1.016 ! CR',
0109                        'PARP(78)=0.538 ! CR',
0110                        'PARP(80)=0.1   ! Prob. colored parton from BBR',
0111                        'PARP(83)=0.356 ! Multiple interactions: matter distribution parameter',
0112                        'PARP(84)=0.651 ! Multiple interactions: matter distribution parameter',
0113                        'PARP(62)=1.025 ! ISR cutoff',
0114                        'MSTP(91)=1     ! Gaussian primordial kT',
0115                        'PARP(93)=10.0  ! primordial kT-max',
0116                        'MSTP(81)=21    ! multiple parton interactions 1 is Pythia default',
0117                        'MSTP(82)=4     ! Defines the multi-parton model'),
0118         processParameters = cms.vstring('PMAS(25,1)=130.0      !mass of Higgs',
0119                     'MSEL=0                  ! user selection for process',
0120                     'MSUB(102)=1             !ggH',
0121                     'MSUB(123)=0             !ZZ fusion to H',
0122                     'MSUB(124)=0             !WW fusion to H',
0123                     'MSUB(24)=0              !ZH production',
0124                     'MSUB(26)=0              !WH production',
0125                     'MSUB(121)=0             !gg to ttH',
0126                     'MSUB(122)=0             !qq to ttH',
0127                     'MDME(210,1)=0           !Higgs decay into dd',
0128                     'MDME(211,1)=0           !Higgs decay into uu',
0129                     'MDME(212,1)=0           !Higgs decay into ss',
0130                     'MDME(213,1)=0           !Higgs decay into cc',
0131                     'MDME(214,1)=0           !Higgs decay into bb',
0132                     'MDME(215,1)=0           !Higgs decay into tt',
0133                     'MDME(216,1)=0           !Higgs decay into',
0134                     'MDME(217,1)=0           !Higgs decay into Higgs decay',
0135                     'MDME(218,1)=0           !Higgs decay into e nu e',
0136                     'MDME(219,1)=0           !Higgs decay into mu nu mu',
0137                     'MDME(220,1)=1           !Higgs decay into tau nu tau',
0138                     'MDME(221,1)=0           !Higgs decay into Higgs decay',
0139                     'MDME(222,1)=0           !Higgs decay into g g',
0140                     'MDME(223,1)=0           !Higgs decay into gam gam',
0141                     'MDME(224,1)=0           !Higgs decay into gam Z',
0142                     'MDME(225,1)=0           !Higgs decay into Z Z',
0143                     'MDME(226,1)=0           !Higgs decay into W W'),
0144         parameterSets = cms.vstring('pythiaUESettings',
0145                     'processParameters')
0146     )
0147                  )
0148
0149
0150 process.ProductionFilterSequence = cms.Sequence(process.generator)
0151
0152 # Path and EndPath definitions
0153 process.generation_step = cms.Path(process.pgen+process.TauSpinnerGen+process.TauSpinnerZHFilter)
0154 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0155 process.validation_step = cms.EndPath(process.genstepfilter+process.genvalid_all)
0156 process.endjob_step = cms.EndPath(process.endOfProcess)
0157 process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
0158
0159 # Schedule definition
0160 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.validation_step,process.endjob_step,process.RAWSIMoutput_step)
0161 # filter all path with the production filter sequence
0162 for path in process.paths:
0163     getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
0164