CMSSW Reference Manual

00001 # Particle gun whose pt and rapidity distributions can be specified
00002 # by TGraphs, and later handled by pythia for decays. The example here
00003 # is jpsi, with kinematic distributions from hep-ph/0310274v1
00004 # and chosen to decay into two muons
00005 
00006 
00007 import FWCore.ParameterSet.Config as cms
00008 
00009 from GeneratorInterface.Pythia6Interface.pythiaDefault_cff import *
00010 
00011 source = cms.Source("PythiaSource",
00012                     pythiaHepMCVerbosity = cms.untracked.bool(False),
00013                     pythiaPylistVerbosity = cms.untracked.int32(0),
00014                     
00015                     ParticleID = cms.untracked.int32(443),
00016                     DoubleParticle = cms.untracked.bool(False),
00017                     kinematicsFile = cms.untracked.string('HeavyIonsAnalysis/Configuration/data/jpsipbpb.root'),
00018                     
00019                     ptBinning = cms.untracked.int32(100000),
00020                     Ptmin = cms.untracked.double(0.0),
00021                     Ptmax = cms.untracked.double(100.0),
00022                     yBinning = cms.untracked.int32(500),
00023                     ymin = cms.untracked.double(-10.0),
00024                     ymax = cms.untracked.double(10.0),
00025                     
00026                     PythiaParameters = cms.PSet(pythiaDefaultBlock,
00027                                                 parameterSets = cms.vstring('pythiaDefault','jpsiDecay'),
00028                                                 jpsiDecay = cms.vstring('BRAT(858) = 0 ! switch off',
00029                                                                         'BRAT(859) = 1 ! switch on',
00030                                                                         'BRAT(860) = 0 ! switch off',
00031                                                                         'MDME(858,1) = 0 ! switch off',
00032                                                                         'MDME(859,1) = 1 ! switch on',
00033                                                                         'MDME(860,1) = 0 ! switch off'),
00034                                                 upsilonDecay = cms.vstring('BRAT(1034) = 0 ! switch off',
00035                                                                            'BRAT(1035) = 1 ! switch on',
00036                                                                            'BRAT(1036) = 0 ! switch off',
00037                                                                            'BRAT(1037) = 0 ! switch off',
00038                                                                            'BRAT(1038) = 0 ! switch off',
00039                                                                            'BRAT(1039) = 0 ! switch off',
00040                                                                            'BRAT(1040) = 0 ! switch off',
00041                                                                            'BRAT(1041) = 0 ! switch off',
00042                                                                            'BRAT(1042) = 0 ! switch off',
00043                                                                            'MDME(1034,1) = 0 ! switch off',
00044                                                                            'MDME(1035,1) = 1 ! switch on',
00045                                                                            'MDME(1036,1) = 0 ! switch off',
00046                                                                            'MDME(1037,1) = 0 ! switch off',
00047                                                                            'MDME(1038,1) = 0 ! switch off',
00048                                                                            'MDME(1039,1) = 0 ! switch off',
00049                                                                            'MDME(1040,1) = 0 ! switch off',
00050                                                                            'MDME(1041,1) = 0 ! switch off',
00051                                                                            'MDME(1042,1) = 0 ! switch off')
00052                                                 )
00053                     )
00054 
00055 # For upsilon generation, add in your configuration the uncommented lines :
00056 #
00057 #process.PythiaSource.kinematicsFile = cms.untracked.string('HeavyIonsAnalysis/Configuration/data/upsipbpb.root')
00058 #process.PythiaSource.ParticleID = cms.untracked.int32(553)
00059
DecayGun_cfi.py
