File indexing completed on 2023-03-17 11:10:14
0001 from __future__ import print_function
0002
0003
0004
0005
0006
0007
0008
0009 import os,sys
0010 from FWCore.ParameterSet.VarParsing import VarParsing
0011 options = VarParsing ('analysis')
0012 options.parseArguments()
0013
0014 if True:
0015 print('-> You are using a 2 process file to unzip/untar events on the fly')
0016
0017
0018 from subprocess import call, check_output
0019 import threading
0020 import time
0021 def call_exe(cmd):
0022 print("-> Executing cmd: '"+cmd+"'")
0023 st=call(cmd,shell=True)
0024 print("-> End cmd: status=",st)
0025 return
0026 cmd="rm /tmp/"+os.environ['USER']+"/hepmc10K.dat"
0027 call(cmd,shell=True)
0028 cmd="mkfifo /tmp/"+os.environ['USER']+"/hepmc10K.dat"
0029 call(cmd,shell=True)
0030 exe="cat /tmp/"+os.environ['USER']+"/hepmc.dat.tgz | gunzip -c > /tmp/"+os.environ['USER']+"/hepmc10K.dat &"
0031 t = threading.Thread(target=call_exe, args= ( [exe] ) )
0032 t.start()
0033 print("(sleep 1s to allow start of pipes)")
0034 time.sleep(1)
0035
0036
0037 import FWCore.ParameterSet.Config as cms
0038
0039 process = cms.Process("GEN")
0040
0041
0042 process.source = cms.Source("MCFileSource",
0043
0044 fileNames = cms.untracked.vstring('file:/tmp/'+os.environ['USER']+'/hepmc10K.dat'),
0045 )
0046
0047 maxEvents=options.maxEvents
0048 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(maxEvents))
0049
0050
0051 process.load("FWCore.MessageService.MessageLogger_cfi")
0052 process.MessageLogger.cerr.threshold = 'INFO'
0053
0054 process.GEN = cms.OutputModule("PoolOutputModule",
0055 fileName = cms.untracked.string('HepMC_GEN.root')
0056 )
0057
0058
0059 process.load('Configuration.StandardSequences.Services_cff')
0060 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0061 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0062 process.load('Configuration.StandardSequences.Generator_cff')
0063 process.genParticles.src= cms.InputTag("source","generator")
0064
0065
0066
0067
0068 from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import GaussVtxSmearingParameters,VtxSmearedCommon
0069 VtxSmearedCommon.src=cms.InputTag("source","generator")
0070 process.generatorSmeared = cms.EDProducer("GaussEvtVtxGenerator",
0071 GaussVtxSmearingParameters,
0072 VtxSmearedCommon
0073 )
0074 process.load('Configuration.StandardSequences.Services_cff')
0075 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0076 generatorSmeared = cms.PSet( initialSeed = cms.untracked.uint32(1243987),
0077 engineName = cms.untracked.string('TRandom3'),
0078 ),
0079 )
0080
0081
0082
0083 process.p = cms.Path(process.genParticles * process.generatorSmeared)
0084 process.outpath = cms.EndPath(process.GEN)
0085
0086
0087
0088
0089
0090
0091
0092
0093
0094
0095
0096
0097
0098
0099
0100
0101
0102
0103
0104
0105
0106
0107
0108
0109
0110
0111
0112
0113
0114
0115
0116