Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:30:39

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 process = cms.Process('HISIGNAL')
0004 
0005 # import of standard configurations
0006 process.load('Configuration.StandardSequences.Services_cff')
0007 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0008 process.load('FWCore.MessageService.MessageLogger_cfi')
0009 process.load('Configuration.StandardSequences.HiEventMixing_cff')
0010 process.load('Configuration.StandardSequences.GeometryExtended_cff')
0011 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0012 process.load('Configuration.StandardSequences.Generator_cff')
0013 process.load('Configuration.StandardSequences.VtxSmeared')
0014 process.load('SimGeneral.MixingModule.himixGEN_cff')
0015 process.load('Configuration.StandardSequences.Sim_cff')
0016 process.load('SimGeneral.MixingModule.himixSIMExtended_cff')
0017 process.load('Configuration.StandardSequences.Digi_cff')
0018 process.load('SimGeneral.MixingModule.himixDIGI_cff')
0019 process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0020 process.load('Configuration.StandardSequences.DigiToRaw_cff')
0021 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0022 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0023 process.load('Configuration.EventContent.EventContentHeavyIons_cff')
0024 
0025 process.configurationMetadata = cms.untracked.PSet(
0026         version = cms.untracked.string('$Revision: 1.168 $'),
0027             annotation = cms.untracked.string('Pyquen_ZeemumuJets_pt10_4TeV_cfi.py nevts:1'),
0028             name = cms.untracked.string('PyReleaseValidation')
0029         )
0030 process.maxEvents = cms.untracked.PSet(
0031         input = cms.untracked.int32(1)
0032         )
0033 process.options = cms.untracked.PSet(
0034 
0035     )
0036 # Input source
0037 process.source = cms.Source("PoolSource",
0038                                 fileNames = cms.untracked.vstring(
0039     'file:/d101/yetkin/FCD18B2A-5A7B-E011-AC06-002618943916.root',
0040     #    '/store/relval/CMSSW_4_2_3/RelValHydjetQ_MinBias_2760GeV/GEN-SIM-DIGI-RAW-HLTDEBUG/STARTHI42_V8-v2/0062/FCD18B2A-5A7B-E011-AC06-002618943916.root'
0041     ),
0042                                 inputCommands = cms.untracked.vstring('drop *',
0043                                                                               'keep *_generator_*_*',
0044                                                                               'keep *_g4SimHits_*_*'),
0045                                 dropDescendantsOfDroppedBranches = cms.untracked.bool(False)
0046                             )
0047 
0048 # Output definition
0049 process.output = cms.OutputModule("PoolOutputModule",
0050                                       splitLevel = cms.untracked.int32(0),
0051                                       outputCommands = process.RAWDEBUGEventContent.outputCommands,
0052                                       fileName = cms.untracked.string('Pyquen_ZeemumuJets_pt10_4TeV_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW.root'),
0053                                       dataset = cms.untracked.PSet(
0054             dataTier = cms.untracked.string('GEN-SIM-RAW'),
0055                     filterName = cms.untracked.string('')
0056                 ),
0057                                       SelectEvents = cms.untracked.PSet(
0058             SelectEvents = cms.vstring('generation_step')
0059                 )
0060                                   )
0061 
0062 # Additional output definition
0063 
0064 # Other statements
0065 process.GlobalTag.globaltag = 'START43_V3::All'
0066 process.hiSignal = cms.EDFilter("PyquenGeneratorFilter",
0067                                     aBeamTarget = cms.double(208.0),
0068                                     comEnergy = cms.double(4000.0),
0069                                     qgpInitialTemperature = cms.double(1.0),
0070                                     doCollisionalEnLoss = cms.bool(False),
0071                                     qgpNumQuarkFlavor = cms.int32(0),
0072                                     qgpProperTimeFormation = cms.double(0.1),
0073                                     numQuarkFlavor = cms.int32(0),
0074                                     hadronFreezoutTemperature = cms.double(0.14),
0075                                     doRadiativeEnLoss = cms.bool(True),
0076                                     backgroundLabel = cms.InputTag("generator","unsmeared"),
0077                                     embeddingMode = cms.bool(True),
0078                                     angularSpectrumSelector = cms.int32(0),
0079                                     doIsospin = cms.bool(True),
0080                                     doQuench = cms.bool(True),
0081                                     cFlag = cms.int32(0),
0082                                     bFixed = cms.double(0.0),
0083                                     bMin = cms.double(0.0),
0084                                     bMax = cms.double(0.0),
0085                                     PythiaParameters = cms.PSet(
0086             pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ',
0087                                                            'BRAT(1035) = 1 ',
0088                                                            'BRAT(1036) = 0 ',
0089                                                            'BRAT(1037) = 0 ',
0090                                                            'BRAT(1038) = 0 ',
0091                                                            'BRAT(1039) = 0 ',
0092                                                            'BRAT(1040) = 0 ',
0093                                                            'BRAT(1041) = 0 ',
0094                                                            'BRAT(1042) = 0 ',
0095                                                            'MDME(1034,1) = 0 ',
0096                                                            'MDME(1035,1) = 1 ',
0097                                                            'MDME(1036,1) = 0 ',
0098                                                            'MDME(1037,1) = 0 ',
0099                                                            'MDME(1038,1) = 0 ',
0100                                                            'MDME(1039,1) = 0 ',
0101                                                            'MDME(1040,1) = 0 ',
0102                                                            'MDME(1041,1) = 0 ',
0103                                                            'MDME(1042,1) = 0 '),
0104                     myParameters = cms.vstring(),
0105                     ppDefault = cms.vstring('MSEL=1   ! QCD hight pT processes',
0106                                                         'CKIN(3)=7.',
0107                                                         'MSTP(81)=0'),
0108                     pythiaZtoElectrons = cms.vstring('MDME(174,1)=0',
0109                                                                  'MDME(175,1)=0',
0110                                                                  'MDME(176,1)=0',
0111                                                                  'MDME(177,1)=0',
0112                                                                  'MDME(178,1)=0',
0113                                                                  'MDME(179,1)=0',
0114                                                                  'MDME(182,1)=1',
0115                                                                  'MDME(183,1)=0',
0116            'MDME(184,1)=0',
0117                                                                  'MDME(185,1)=0',
0118                                                                  'MDME(186,1)=0',
0119                                                                  'MDME(187,1)=0'),
0120                     pythiaZjets = cms.vstring('MSUB(15)=1',
0121                                                           'MSUB(30)=1'),
0122                     pythiaUESettings = cms.vstring('MSTJ(11)=3     ! Choice of the fragmentation function',
0123                                                                'MSTJ(22)=2     ! Decay those unstable particles',
0124                                                                'PARJ(71)=10 .  ! for which ctau  10 mm',
0125                                                                'MSTP(2)=1      ! which order running alphaS',
0126                                                                'MSTP(33)=0     ! no K factors in hard cross sections',
0127                                                                'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
0128                                                                'MSTP(52)=2     ! work with LHAPDF',
0129                                                                'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default',
0130                                                                'MSTP(82)=4     ! Defines the multi-parton model',
0131                                                                'MSTU(21)=1     ! Check on possible errors during program execution',
0132                                                                'PARP(82)=1.8387   ! pt cutoff for multiparton interactions',
0133                                                                'PARP(89)=1960. ! sqrts for which PARP82 is set',
0134                                                                'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter',
0135                                                                'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter',
0136                                                                'PARP(90)=0.16  ! Multiple interactions: rescaling power',
0137                                                                'PARP(67)=2.5    ! amount of initial-state radiation',
0138                                                                'PARP(85)=1.0  ! gluon prod. mechanism in MI',
0139                                                                'PARP(86)=1.0  ! gluon prod. mechanism in MI',
0140                                                                'PARP(62)=1.25   ! ',
0141                                                                'PARP(64)=0.2    ! ',
0142                                                                'MSTP(91)=1      !',
0143                                                                'PARP(91)=2.1   ! kt distribution',
0144                                                                'PARP(93)=15.0  ! '),
0145                     pythiaZtoMuonsAndElectrons = cms.vstring('MDME(174,1)=0',
0146                                                                          'MDME(175,1)=0',
0147                                                                          'MDME(176,1)=0',
0148                                                                          'MDME(177,1)=0',
0149                                                                          'MDME(178,1)=0',
0150                                                                          'MDME(179,1)=0',
0151                                                                          'MDME(182,1)=1',
0152                                                                          'MDME(183,1)=0',
0153                                                                          'MDME(184,1)=1',
0154                                                                          'MDME(185,1)=0',
0155                                                                          'MDME(186,1)=0',
0156                                                                          'MDME(187,1)=0'),
0157                     pythiaPromptPhotons = cms.vstring('MSUB(14)=1',
0158                                                                   'MSUB(18)=1',
0159                                                                   'MSUB(29)=1',
0160                                                                   'MSUB(114)=1',
0161                                                                   'MSUB(115)=1'),
0162                     pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
0163                                                                     'MSUB(422) = 1',
0164                                                                     'MSUB(423) = 1',
0165                                                                     'MSUB(424) = 1',
0166                                                                     'MSUB(425) = 1',
0167                                                                     'MSUB(426) = 1',
0168                                                                     'MSUB(427) = 1',
0169                                                                     'MSUB(428) = 1',
0170                                                                     'MSUB(429) = 1',
0171                                                                     'MSUB(430) = 1',
0172                                                                     'MSUB(431) = 1',
0173                                                                     'MSUB(432) = 1',
0174                                                                     'MSUB(433) = 1',
0175                                                                     'MSUB(434) = 1',
0176                                                                     'MSUB(435) = 1',
0177                                                                     'MSUB(436) = 1',
0178                                                                     'MSUB(437) = 1',
0179                                                                     'MSUB(438) = 1',
0180                                                                     'MSUB(439) = 1'),
0181                     pythiaMuonCandidates = cms.vstring('CKIN(3)=20',
0182                                                                    'MSTJ(22)=2',
0183                                                                    'PARJ(71)=40.'),
0184                     pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16',
0185                                                                       'PARP(142)=0.0119',
0186                                                                       'PARP(143)=0.01',
0187                                                                       'PARP(144)=0.01',
0188                                                                       'PARP(145)=0.05',
0189                                                                       'PARP(146)=9.28',
0190                                                                       'PARP(147)=0.15',
0191                                                                       'PARP(148)=0.02',
0192                                                                       'PARP(149)=0.02',
0193                                                                       'PARP(150)=0.085',
0194                                                                       'PARJ(13)=0.60',
0195                                                                       'PARJ(14)=0.162',
0196                                                                       'PARJ(15)=0.018',
0197                                                                       'PARJ(16)=0.054',
0198                                                                       'MSTP(145)=0',
0199                                                                       'MSTP(146)=0',
0200                                                                       'MSTP(147)=0',
0201                                                                       'MSTP(148)=1',
0202                                                                       'MSTP(149)=1',
0203                                                                       'BRAT(861)=0.202',
0204                                                                       'BRAT(862)=0.798',
0205                                                                       'BRAT(1501)=0.013',
0206                                                                       'BRAT(1502)=0.987',
0207                                                                       'BRAT(1555)=0.356',
0208                                                                       'BRAT(1556)=0.644'),
0209 
0210         pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
0211                                                          'MSUB(462) = 1',
0212                                                          'MSUB(463) = 1',
0213                                                          'MSUB(464) = 1',
0214                                                          'MSUB(465) = 1',
0215                                                          'MSUB(466) = 1',
0216                                                          'MSUB(467) = 1',
0217                                                          'MSUB(468) = 1',
0218                                                          'MSUB(469) = 1',
0219                                                          'MSUB(470) = 1',
0220                                                          'MSUB(471) = 1',
0221                                                          'MSUB(472) = 1',
0222                                                          'MSUB(473) = 1',
0223                                                          'MSUB(474) = 1',
0224                                                          'MSUB(475) = 1',
0225                                                          'MSUB(476) = 1',
0226                                                          'MSUB(477) = 1',
0227                                                          'MSUB(478) = 1',
0228                                                          'MSUB(479) = 1'),
0229                     pythiaWeakBosons = cms.vstring('MSUB(1)=1',
0230                                                                'MSUB(2)=1'),
0231                     pythiaJets = cms.vstring('MSUB(11)=1',
0232                                                          'MSUB(12)=1',
0233                                                          'MSUB(13)=1',
0234                                                          'MSUB(28)=1',
0235                                                          'MSUB(53)=1',
0236                                                          'MSUB(68)=1'),
0237                     pythiaZtoMuons = cms.vstring('MDME(174,1)=0',
0238                                                              'MDME(175,1)=0',
0239                                                              'MDME(176,1)=0',
0240                                                              'MDME(177,1)=0',
0241                                                              'MDME(178,1)=0',
0242                                                              'MDME(179,1)=0',
0243                                                              'MDME(182,1)=0',
0244                                                              'MDME(183,1)=0',
0245                                                              'MDME(184,1)=1',
0246                                                              'MDME(185,1)=0',
0247                                                              'MDME(186,1)=0',
0248                                                              'MDME(187,1)=0'),
0249                     ppJets = cms.vstring('MSEL=1   ! QCD hight pT processes'),
0250                     pythiaHirootDefault = cms.vstring('MSEL=0',
0251                                                                   'MSTU(21)=1',
0252                                                                   'PARU(14)=1.',
0253                                                                   'MSTP(81)=0',
0254                                                                   'PMAS(5,1)=4.8',
0255                                                                   'PMAS(6,1)=175.0',
0256                                                                   'CKIN(3)=7.',
0257                                                                   'MSTJ(22)=2',
0258                                                                   'PARJ(71)=10.',
0259                                                                   'PARP(67)=1.',
0260                                                                   'PARP(82)=1.9',
0261                                                                   'PARP(85)=0.33',
0262                                                                   'PARP(86)=0.66',
0263                                                                   'PARP(89)=1000.',
0264                                                                   'PARP(91)=1.0',
0265                                                                   'MSTJ(11)=3',
0266                                                                   'MSTJ(22)=2'),
0267                     pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ',
0268                                                                 'BRAT(859) = 1 ',
0269                                                                 'BRAT(860) = 0 ',
0270                                                                 'MDME(858,1) = 0 ',
0271                                                                 'MDME(859,1) = 1 ',
0272                                                                 'MDME(860,1) = 0 '),
0273                     parameterSets = cms.vstring('pythiaUESettings',
0274                                                             'pythiaZjets',
0275                                                             'pythiaZtoMuonsAndElectrons',
0276                                                             'kinematics'),
0277                     kinematics = cms.vstring('CKIN(3)=10',
0278                                                          'CKIN(4)=9999',
0279                                                          'CKIN(7)=-2.',
0280                                                          'CKIN(8)=2.')
0281                 )
0282                                 )
0283 
0284 # Path and EndPath definitions
0285 process.generation_step = cms.Path(process.pgen_himix)
0286 process.simulation_step = cms.Path(process.psim)
0287 process.digitisation_step = cms.Path(process.pdigi)
0288 process.L1simulation_step = cms.Path(process.SimL1Emulator)
0289 process.digi2raw_step = cms.Path(process.DigiToRaw)
0290 process.endjob_step = cms.Path(process.endOfProcess)
0291 process.out_step = cms.EndPath(process.output)
0292 
0293 # Schedule definition
0294 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)
0295 
0296 # special treatment in case of production filter sequence
0297 for path in process.paths:
0298     getattr(process,path)._seq = process.hiSignal*getattr(process,path)._seq
0299 
0300