GeneratorInterface/Hydjet2Interface/test/testHydjet.py - master - githubmirror/cmssw

GeneratorInterface/Hydjet2Interface/test/testHydjet.py
99 строк4 KB
Andrey
Hydjet2 bugfixes for status codes and units
20 июн 2022, 16:16
20 июн 2022, 16:16b310a57
import FWCore.ParameterSet.Config as cms

process = cms.Process("ANA")

process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.load("Configuration.StandardSequences.Services_cff")
process.load("GeneratorInterface.Hydjet2Interface.hydjet2Default_cfi")

process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))

process.ana = cms.EDAnalyzer('Hydjet2Analyzer',

		doHistos = cms.untracked.bool(True),
                userHistos = cms.untracked.bool(False),
		doAnalysis = cms.untracked.bool(True),
                doTestEvent = cms.untracked.bool(False), # for debuging event output information

		###Settings for USER histos

		#status
                uStatus = cms.untracked.int32(1), #1 - it's 1,2,3,4,5 of Pythia status; 2 - 11,12,13,14,15; 3 - All

		#up to 3 abs(PDG) for selection, if less needed just comment not used
                uPDG_1 = cms.untracked.int32(211),
                uPDG_2 = cms.untracked.int32(321),
                uPDG_3 = cms.untracked.int32(2212),

                # |eta| cut for pT dep.dist.
                dPTetaCut = cms.untracked.double(0.), #down
                uPTetaCut = cms.untracked.double(0.8), #up

                #Vectors of bins borders(when 0 - uniform bins would be used)
                PtBins = cms.untracked.vdouble(0.), #, 1., 2., 3., 4., 5., 6., 8., 12., 16., 20.),
                EtaBins = cms.untracked.vdouble(0.),
                PhiBins = cms.untracked.vdouble(0.),
                v2EtaBins = cms.untracked.vdouble(0.),
                v2PtBins = cms.untracked.vdouble(0.), #, 1., 2., 3., 4., 6., 8., 12., 16., 20.),

                #Settings for uniform bins 
		nintPt 		= cms.untracked.int32(100),
                nintEta 	= cms.untracked.int32(51),
		nintPhi		= cms.untracked.int32(100),
		nintV2pt	= cms.untracked.int32(100),
		nintV2eta	= cms.untracked.int32(100),

		minPt		= cms.untracked.double(0.),
		minEta		= cms.untracked.double(-10.),
		minPhi		= cms.untracked.double(-3.14159265358979),
		minV2pt		= cms.untracked.double(0.),
		minV2eta	= cms.untracked.double(-10.),
	
		maxPt		= cms.untracked.double(100.),
		maxEta		= cms.untracked.double(10.),
		maxPhi		= cms.untracked.double(3.14159265358979),
		maxV2pt		= cms.untracked.double(10.),
		maxV2eta	= cms.untracked.double(10.),

)

#to separate hydro and jet parts of hydjet2	
process.generator.separateHydjetComponents = cms.untracked.bool(False)
Debug = None

if Debug:
	process.load("FWCore.MessageLogger.MessageLogger_cfi")

	process.MessageLogger = cms.Service("MessageLogger",

		destinations     = cms.untracked.vstring('LogDebug_Hydjet2'),
		categories       = cms.untracked.vstring(
                	                        	'Hydjet2',
                        	                	'Hydjet2_array'
                                		        ),
		LogDebug_Hydjet2 = cms.untracked.PSet(
        		threshold =  cms.untracked.string('DEBUG'), #Priority: DEBUG < INFO < WARNING < ERROR
        		DEBUG   = cms.untracked.PSet(limit = cms.untracked.int32(-1)),
        		INFO    = cms.untracked.PSet(limit = cms.untracked.int32(0)),
        		WARNING = cms.untracked.PSet(limit = cms.untracked.int32(0)),
        		ERROR   = cms.untracked.PSet(limit = cms.untracked.int32(0)),
			#Categores
        		Hydjet2  = cms.untracked.PSet(
                        	limit = cms.untracked.int32(-1), # number of masseges 
                        	timespan = cms.untracked.int32(0)     #time to resete limit counter in seconds
                        	),

        		Hydjet2_array  = cms.untracked.PSet(
                        	limit = cms.untracked.int32(-1), # number of masseges 
                        	timespan = cms.untracked.int32(0)     #time to resete limit counter in seconds
                        	)

			),
			debugModules     = cms.untracked.vstring('*')
		)

process.TFileService = cms.Service('TFileService',
	fileName = cms.string('Hydjet2_5020GeV.root')
)

process.p = cms.Path(process.generator*process.ana)