Hydjet2Interface/test/gensimTest.py

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