Hydjet2Interface/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.Hydjet2Interface.hydjet2Default_cfi")
0008
0009 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
0010
0011 process.ana = cms.EDAnalyzer('Hydjet2Analyzer',
0012
0013         doHistos = cms.untracked.bool(True),
0014                 userHistos = cms.untracked.bool(False),
0015         doAnalysis = cms.untracked.bool(True),
0016                 doTestEvent = cms.untracked.bool(False), # for debuging event output information
0017
0018         ###Settings for USER histos
0019
0020         #status
0021                 uStatus = cms.untracked.int32(1), #1 - it's 1,2,3,4,5 of Pythia status; 2 - 11,12,13,14,15; 3 - All
0022
0023         #up to 3 abs(PDG) for selection, if less needed just comment not used
0024                 uPDG_1 = cms.untracked.int32(211),
0025                 uPDG_2 = cms.untracked.int32(321),
0026                 uPDG_3 = cms.untracked.int32(2212),
0027
0028                 # |eta| cut for pT dep.dist.
0029                 dPTetaCut = cms.untracked.double(0.), #down
0030                 uPTetaCut = cms.untracked.double(0.8), #up
0031
0032                 #Vectors of bins borders(when 0 - uniform bins would be used)
0033                 PtBins = cms.untracked.vdouble(0.), #, 1., 2., 3., 4., 5., 6., 8., 12., 16., 20.),
0034                 EtaBins = cms.untracked.vdouble(0.),
0035                 PhiBins = cms.untracked.vdouble(0.),
0036                 v2EtaBins = cms.untracked.vdouble(0.),
0037                 v2PtBins = cms.untracked.vdouble(0.), #, 1., 2., 3., 4., 6., 8., 12., 16., 20.),
0038
0039                 #Settings for uniform bins
0040         nintPt      = cms.untracked.int32(100),
0041                 nintEta     = cms.untracked.int32(51),
0042         nintPhi     = cms.untracked.int32(100),
0043         nintV2pt    = cms.untracked.int32(100),
0044         nintV2eta   = cms.untracked.int32(100),
0045
0046         minPt       = cms.untracked.double(0.),
0047         minEta      = cms.untracked.double(-10.),
0048         minPhi      = cms.untracked.double(-3.14159265358979),
0049         minV2pt     = cms.untracked.double(0.),
0050         minV2eta    = cms.untracked.double(-10.),
0051
0052         maxPt       = cms.untracked.double(100.),
0053         maxEta      = cms.untracked.double(10.),
0054         maxPhi      = cms.untracked.double(3.14159265358979),
0055         maxV2pt     = cms.untracked.double(10.),
0056         maxV2eta    = cms.untracked.double(10.),
0057
0058 )
0059
0060 #to separate hydro and jet parts of hydjet2
0061 process.generator.separateHydjetComponents = cms.untracked.bool(False)
0062 Debug = None
0063
0064 if Debug:
0065     process.load("FWCore.MessageLogger.MessageLogger_cfi")
0066
0067     process.MessageLogger = cms.Service("MessageLogger",
0068
0069         destinations     = cms.untracked.vstring('LogDebug_Hydjet2'),
0070         categories       = cms.untracked.vstring(
0071                                                 'Hydjet2',
0072                                                 'Hydjet2_array'
0073                                                 ),
0074         LogDebug_Hydjet2 = cms.untracked.PSet(
0075                 threshold =  cms.untracked.string('DEBUG'), #Priority: DEBUG < INFO < WARNING < ERROR
0076                 DEBUG   = cms.untracked.PSet(limit = cms.untracked.int32(-1)),
0077                 INFO    = cms.untracked.PSet(limit = cms.untracked.int32(0)),
0078                 WARNING = cms.untracked.PSet(limit = cms.untracked.int32(0)),
0079                 ERROR   = cms.untracked.PSet(limit = cms.untracked.int32(0)),
0080             #Categores
0081                 Hydjet2  = cms.untracked.PSet(
0082                             limit = cms.untracked.int32(-1), # number of masseges
0083                             timespan = cms.untracked.int32(0)     #time to resete limit counter in seconds
0084                             ),
0085
0086                 Hydjet2_array  = cms.untracked.PSet(
0087                             limit = cms.untracked.int32(-1), # number of masseges
0088                             timespan = cms.untracked.int32(0)     #time to resete limit counter in seconds
0089                             )
0090
0091             ),
0092             debugModules     = cms.untracked.vstring('*')
0093         )
0094
0095 process.TFileService = cms.Service('TFileService',
0096     fileName = cms.string('Hydjet2_5020GeV.root')
0097 )
0098
0099 process.p = cms.Path(process.generator*process.ana)