CustomPhysics/test/Sexaquark_RunIIFall18GenSim_cfg.py

0001 # Auto generated configuration file
0002 # using:
0003 # Revision: 1.19
0004 # Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
0005 # with command line options: Configuration/GenProduction/python/BPH-RunIIFall18GS-00369-fragment.py --python_filename BPH-RunIIFall18GS-00369_1_cfg.py --eventcontent RAWSIM --datatier GEN-SIM --fileout file:BPH-RunIIFall18GS-00369.root --conditions 102X_upgrade2018_realistic_v11 --beamspot Realistic25ns13TeVEarly2018Collision --step GEN,SIM --geometry DB:Extended --era Run2_2018 --no_exec --mc -n 1000
0006 import FWCore.ParameterSet.Config as cms
0007 from Configuration.StandardSequences.Eras import eras
0008 from FWCore.ParameterSet.VarParsing import VarParsing
0009
0010 process = cms.Process('SIM',eras.Run3_2023)
0011 options = VarParsing ('analysis')
0012 options.outputFile = 'file:sexaq_sim.root'
0013 #options.inputFiles = 'root://cmsxrootd.hep.wisc.edu//store/user/wvetens/crmc_Sexaq/crmc/Sexaquark_13TeV_trial_4_1p8GeV/0/crmc_Sexaq_1.root'
0014 options.inputFiles = 'file:crmc_Sexaq_1.root'
0015 options.maxEvents= 100
0016 options.parseArguments()
0017
0018 # import of standard configurations
0019 process.load('Configuration.StandardSequences.Services_cff')
0020 #process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0021 process.load('FWCore.MessageService.MessageLogger_cfi')
0022 process.load('Configuration.EventContent.EventContent_cff')
0023 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0024 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0025 process.load('Configuration.StandardSequences.GeometrySimDB_cff')
0026 process.load('Configuration.StandardSequences.MagneticField_cff')
0027 process.load('Configuration.StandardSequences.Generator_cff')
0028 # Vtx Smearing done in hepmc 2 gen step
0029 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic25ns13TeVEarly2018Collision_cfi')
0030 #process.load('GeneratorInterface.Core.genFilterSummary_cff')
0031 process.load('Configuration.StandardSequences.SimIdeal_cff')
0032 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0033 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0034 process.load('RecoVertex.BeamSpotProducer.BeamSpot_cfi')
0035
0036 # Lengthy message logs - uncomment to debug
0037 process.MessageLogger = cms.Service("MessageLogger",
0038   destinations = cms.untracked.vstring('cout'),
0039   cout = cms.untracked.PSet(
0040     threshold = cms.untracked.string('INFO')
0041   )
0042 )
0043
0044 process.maxEvents = cms.untracked.PSet(
0045     input = cms.untracked.int32(options.maxEvents)
0046     #input = cms.untracked.int32(-1)
0047 )
0048
0049 # Input source
0050
0051 #process.source = cms.Source("EmptySource")
0052 process.source = cms.Source("PoolSource",
0053     fileNames = cms.untracked.vstring(options.inputFiles),
0054     skipEvents = cms.untracked.uint32(0),
0055     duplicateCheckMode = cms.untracked.string ("noDuplicateCheck")
0056 )
0057
0058
0059 process.options = cms.untracked.PSet(
0060   wantSummary = cms.untracked.bool(True)
0061 )
0062
0063 # Output definition
0064
0065 process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
0066     SelectEvents = cms.untracked.PSet(
0067         SelectEvents = cms.vstring('generation_step')
0068     ),
0069     compressionAlgorithm = cms.untracked.string('LZMA'),
0070     compressionLevel = cms.untracked.int32(1),
0071     dataset = cms.untracked.PSet(
0072         dataTier = cms.untracked.string('GEN-SIM'),
0073         filterName = cms.untracked.string('')
0074     ),
0075     eventAutoFlushCompressedSize = cms.untracked.int32(20971520),
0076     fileName = cms.untracked.string(options.outputFile),
0077     outputCommands = process.RAWSIMEventContent.outputCommands,
0078     splitLevel = cms.untracked.int32(0)
0079 )
0080
0081 # Additional output definition
0082 process.RAWSIMoutput.outputCommands += ("keep *_genParticlesPlusGEANT_*_*",)
0083
0084 # Other statements
0085 #process.XMLFromDBSource.label = cms.string("Extended")
0086 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0087 from Configuration.AlCa.GlobalTag import GlobalTag
0088 #process.GlobalTag = GlobalTag(process.GlobalTag, '102X_upgrade2018_realistic_v11', '')
0089 #process.GlobalTag = GlobalTag(process.GlobalTag, '140X_mcRun3_2024_realistic_v14', '')
0090 process.GlobalTag = GlobalTag(process.GlobalTag, '140X_mcRun3_2023_realistic_v3', '')
0091 #process.GlobalTag = GlobalTag(process.GlobalTag, '140X_mcRun3_2022_realistic_v3', '')
0092
0093
0094
0095 process.genParticlesPlusGEANT = cms.EDProducer("GenPlusSimParticleProducer",
0096   src           = cms.InputTag("g4SimHits"),
0097   setStatus     = cms.int32(8),                 # set status = 8 for GEANT GPs
0098   particleTypes = cms.vstring(),
0099   filter = cms.vstring(),
0100   genParticles  = cms.InputTag("genParticles") # original genParticle list
0101 )
0102
0103 from SimG4Core.CustomPhysics.CustomPhysics_cfi import customPhysicsSetup
0104 process.g4SimHits.Physics.type = cms.string('SimG4Core/Physics/CustomPhysics')
0105 process.g4SimHits.Physics.RHadronDummyFlip = cms.bool(False)
0106 process.g4SimHits.Physics.Verbosity = 1
0107 process.g4SimHits.Physics = cms.PSet(
0108   process.g4SimHits.Physics, #keep all default value and add others
0109   customPhysicsSetup
0110   )
0111 process.g4SimHits.Physics.particlesDef = cms.FileInPath('SimG4Core/CustomPhysics/data/particles_sexaq_1p8_GeV.txt')
0112
0113
0114 ## Vtx Smearing done in hepmc 2 gen step
0115 process.VtxSmeared.src = cms.InputTag("source", "generator")
0116 process.genParticles.src = cms.InputTag("generatorSmeared")
0117 process.g4SimHits.HepMCProductLabel = cms.InputTag("generatorSmeared")
0118 process.g4SimHits.Generator.HepMCProductLabel = cms.InputTag("generatorSmeared")
0119
0120
0121 # Path and EndPath definitions
0122 process.generation_step = cms.Path(process.pgen)
0123 #moved beamspot and vtx smearing to hepmc2gen step
0124 #process.simulation_step = cms.Path(process.offlineBeamSpot*process.generatorSmeared*process.psim*process.genParticlesPlusGEANT)
0125 process.simulation_step = cms.Path(process.psim*process.genParticlesPlusGEANT)
0126 #process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0127 process.endjob_step = cms.EndPath(process.endOfProcess)
0128 process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
0129
0130 # Schedule definition
0131 #process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.endjob_step,process.RAWSIMoutput_step)
0132 process.schedule = cms.Schedule(process.generation_step,process.simulation_step,process.endjob_step,process.RAWSIMoutput_step)
0133
0134 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0135 associatePatAlgosToolsTask(process)
0136 ## filter all path with the production filter sequence
0137 #for path in process.paths:
0138 #   getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
0139
0140 from Configuration.DataProcessing.Utils import addMonitoring
0141 process = addMonitoring(process)
0142
0143 # Customisation from command line
0144
0145 # Add early deletion of temporary data products to reduce peak memory need
0146 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0147 process = customiseEarlyDelete(process)
0148
0149 # End adding early deletion
0150 # For debug:
0151 #print process.dumpPython()