HydjetInterface/test/testHydjet.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("ANA")
0004
0005 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0006 process.load("Configuration.StandardSequences.Services_cff")
0007 process.load("GeneratorInterface.HydjetInterface.hydjetDefault_cfi")
0008
0009 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0010     generator = cms.PSet(
0011         initialSeed = cms.untracked.uint32(123456789),
0012         engineName = cms.untracked.string('HepJamesRandom')
0013     )
0014 )
0015
0016 process.maxEvents = cms.untracked.PSet(
0017     input = cms.untracked.int32(100)
0018     )
0019
0020 process.ana = cms.EDAnalyzer('HydjetAnalyzer',
0021
0022     doHistos        = cms.untracked.bool(True),
0023     userHistos      = cms.untracked.bool(False),
0024
0025     # Settings for USER histos
0026
0027     uStatus     = cms.untracked.int32(2),   # 1 - it's 1,2,3,4,5 of Pythia status; 2 - 11,12,13,14,15; 3 - All
0028     uPDG_1      = cms.untracked.int32(443),
0029     uPDG_2      = cms.untracked.int32(-443),
0030
0031     ### Eta cut for pT dep. dist.
0032     uPTetaCut   = cms.untracked.double(4.),
0033     dPTetaCut   = cms.untracked.double(2.5),
0034
0035     ### Vectors of bins borders (when 0 -  uniform bins would be used)
0036     PtBins      = cms.untracked.vdouble(0.,1.,2.,3.,4.,5.,6.,8.,12.,16.,20.),
0037     EtaBins         = cms.untracked.vdouble(0.),
0038     PhiBins         = cms.untracked.vdouble(0.),
0039     v2EtaBins       = cms.untracked.vdouble(0.),
0040     v2PtBins        = cms.untracked.vdouble(0.,1.,2.,3.,4.,6.,8.,12.,16.,20.),
0041
0042     ### Settings for uniform bins
0043     nintPt      = cms.untracked.int32(1000),
0044     nintEta     = cms.untracked.int32(100),
0045     nintPhi     = cms.untracked.int32(100),
0046     nintV2pt        = cms.untracked.int32(100),
0047     nintV2eta       = cms.untracked.int32(100),
0048
0049     minPt       = cms.untracked.double(0.),
0050     minEta      = cms.untracked.double(-10.),
0051     minPhi      = cms.untracked.double(-3.14159265358979),
0052     minV2pt     = cms.untracked.double(0.),
0053     minV2eta        = cms.untracked.double(-10.),
0054
0055     maxPt       = cms.untracked.double(100.),
0056     maxEta      = cms.untracked.double(10.),
0057     maxPhi      = cms.untracked.double(3.14159265358979),
0058     maxV2pt     = cms.untracked.double(10.),
0059     maxV2eta        = cms.untracked.double(10.),
0060
0061 )
0062
0063
0064 #process.generator.signalVtx = cms.untracked.vdouble(0.,0.,0.,0.) # Signal event vertex option, to set it by hand (instead of smearing)
0065
0066 process.TFileService = cms.Service('TFileService',
0067     fileName = cms.string('Hydjet1_MB_5020GeV.root')
0068 )
0069
0070 process.p = cms.Path(process.generator*process.ana)
0071
0072