Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:52:39

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.PythiaUESettings_cfi import *
0003 
0004 generator = cms.EDFilter("Pythia6GeneratorFilter",
0005     pythiaHepMCVerbosity = cms.untracked.bool(False),
0006     maxEventsToPrint = cms.untracked.int32(0),
0007     pythiaPylistVerbosity = cms.untracked.int32(1),
0008     filterEfficiency = cms.untracked.double(0.002305),
0009     comEnergy = cms.double(10000.0),
0010     crossSection = cms.untracked.double(51560000000.),
0011     PythiaParameters = cms.PSet(
0012         pythiaUESettingsBlock,
0013         processParameters = cms.vstring(
0014             'MSEL=1           ! User defined processes',
0015             'MSTJ(22)=4       ! Decay unstable particles in a cylinder',
0016             'PARJ(73)=1500.   ! max. radius for MSTJ(22)=4',
0017             'PARJ(74)=3000.   ! max. Z for MSTJ(22)=4',
0018             'MDCY(C130,1)=1   ! decay k0-longs',
0019             'MDCY(C211,1)=1   ! decay pions',
0020             'MDCY(C321,1)=1   ! decay kaons'),
0021 
0022         # This is a vector of ParameterSet names to be read, in this order
0023         parameterSets = cms.vstring('pythiaUESettings', 
0024             'processParameters')
0025     )
0026 )
0027 
0028 mugenfilter = cms.EDFilter("MCSmartSingleParticleFilter",
0029                            MinPt = cms.untracked.vdouble(2.5,2.5),
0030                            MinEta = cms.untracked.vdouble(-2.5,-2.5),
0031                            MaxEta = cms.untracked.vdouble(2.5,2.5),
0032                            ParticleID = cms.untracked.vint32(13,-13),
0033                            Status = cms.untracked.vint32(1,1),
0034                            # Decay cuts are in mm
0035                            MaxDecayRadius = cms.untracked.vdouble(1500.,1500.),
0036                            MinDecayZ = cms.untracked.vdouble(-3000.,-3000.),
0037                            MaxDecayZ = cms.untracked.vdouble(3000.,3000.)
0038 )
0039 
0040 ProductionFilterSequence = cms.Sequence(generator*mugenfilter)