Merge pull request #40704 from bsunanda/Phase2-sim139

Phase2-sim139 Change scripts to study new Phase2 scenarios in the context of usage and overlaps
cms-sw · Feb 7, 2023 · 184e419 · 184e419
2 parents 8dac892 + 0a0885f
commit 184e419
Show file tree

Hide file tree

Showing 5 changed files with 171 additions and 79 deletions.
diff --git a/SimG4Core/PrintGeomInfo/test/python/g4OverlapCheck2026DD4hep_cfg.py b/SimG4Core/PrintGeomInfo/test/python/g4OverlapCheck2026DD4hep_cfg.py
@@ -2,7 +2,7 @@
 # Way to use this:
 #   cmsRun g4OverlapCheck2026DD4hep_cfg.py geometry=D92 tol=0.1
 #
-#   Options for geometry D88, D92, D93
+#   Options for geometry D88, D92, D93, D94, D95, D96, D98, D99
 #
 ###############################################################################
 import FWCore.ParameterSet.Config as cms
@@ -16,7 +16,7 @@
                  "D88",
                   VarParsing.VarParsing.multiplicity.singleton,
                   VarParsing.VarParsing.varType.string,
-                  "geometry of operations: D88, D92, D93")
+                  "geometry of operations: D88, D92, D93, D94, D95, D96, D98, D99")
 options.register('tol',
                  0.1,
                  VarParsing.VarParsing.multiplicity.singleton,
@@ -32,25 +32,17 @@
 ####################################################################
 # Use the options
 
+from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
 from Configuration.ProcessModifiers.dd4hep_cff import dd4hep
-if (options.geometry == "D92"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    process = cms.Process('PROD',Phase2C11M9,dd4hep)
-    process.load('Configuration.Geometry.GeometryDD4hepExtended2026D92Reco_cff')
-    baseName = 'cms2026D92DD4hep'
-elif (options.geometry == "D93"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    process = cms.Process('PROD',Phase2C11M9,Run3_dd4hep)
-    process.load('Configuration.Geometry.GeometryDD4hepExtended2026D93Reco_cff')
-    baseName = 'cms2026D93DD4hep'
-else:
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    process = cms.Process('PROD',Phase2C11M9,dd4hep)
-    process.load('Configuration.Geometry.GeometryDD4hepExtended2026D88Reco_cff')
-    baseName = 'cms2026D88DD4hep'
+process = cms.Process('PROD',Phase2C17I13M9,dd4hep)
 
-print("Base file Name: ", baseName)
+geomFile = "Configuration.Geometry.GeometryDD4hepExtended2026" + options.geometry + "Reco_cff"
+baseName = "cms2026" + options.geometry + "DD4hep"
 
+print("Geometry file Name: ", geomFile)
+print("Base file Name:     ", baseName)
+
+process.load(geomFile)
 process.load('FWCore.MessageService.MessageLogger_cfi')
 
 #if hasattr(process,'MessageLogger'):

diff --git a/SimG4Core/PrintGeomInfo/test/python/g4OverlapCheck2026DDD_cfg.py b/SimG4Core/PrintGeomInfo/test/python/g4OverlapCheck2026DDD_cfg.py
@@ -2,7 +2,7 @@
 # Way to use this:
 #   cmsRun g4OverlapCheck2026DDD_cfg.py geometry=D88 tol=0.1
 #
-#   Options for geometry D88, D92, D93
+#   Options for geometry D88, D92, D93, D94, D95, D96, D98, D99
 #
 ###############################################################################
 import FWCore.ParameterSet.Config as cms
@@ -16,7 +16,7 @@
                  "D88",
                   VarParsing.VarParsing.multiplicity.singleton,
                   VarParsing.VarParsing.varType.string,
-                  "geometry of operations: D88, D92, D93")
+                  "geometry of operations: D88, D92, D93, D94, D95, D96, D98, D99")
 options.register('tol',
                  0.1,
                  VarParsing.VarParsing.multiplicity.singleton,
@@ -32,24 +32,16 @@
 ####################################################################
 # Use the options
 
-if (options.geometry == "D92"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    process = cms.Process('PROD',Phase2C11M9)
-    process.load('Configuration.Geometry.GeometryExtended2026D92Reco_cff')
-    baseName = 'cms2026D92DDD'
-elif (options.geometry == "D93"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    process = cms.Process('PROD',Phase2C11M9)
-    process.load('Configuration.Geometry.GeometryExtended2026D93Reco_cff')
-    baseName = 'cms2026D93DDD'
-else:
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    process = cms.Process('PROD',Phase2C11M9)
-    process.load('Configuration.Geometry.GeometryExtended2026D88Reco_cff')
-    baseName = 'cms2026D88DDD'
+from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
+process = cms.Process('PROD',Phase2C17I13M9)
 
-print("Base file Name: ", baseName)
+geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
+baseName = "cms2026" + options.geometry + "DDD"
 
+print("Geometry file Name: ", geomFile)
+print("Base file Name:     ", baseName)
+
+process.load(geomFile)
 process.load('FWCore.MessageService.MessageLogger_cfi')
 
 #if hasattr(process,'MessageLogger'):

diff --git a/SimG4Core/PrintGeomInfo/test/python/runPrintG4Solids_cfg.py b/SimG4Core/PrintGeomInfo/test/python/runPrintG4Solids_cfg.py
@@ -1,8 +1,8 @@
 ###############################################################################
 # Way to use this:
-#   cmsRun grunPrintG4Solids_cfg.py geometry=D86 dd4hep=False
+#   cmsRun grunPrintG4Solids_cfg.py geometry=D88 dd4hep=False
 #
-#   Options for geometry D77, D83, D88, D92, D93
+#   Options for geometry D88, D92, D93, D94, D95, D96, D98, D99
 #
 ###############################################################################
 import FWCore.ParameterSet.Config as cms
@@ -13,10 +13,10 @@
 ### SETUP OPTIONS
 options = VarParsing.VarParsing('standard')
 options.register('geometry',
-                 "D92",
+                 "D88",
                   VarParsing.VarParsing.multiplicity.singleton,
                   VarParsing.VarParsing.varType.string,
-                  "geometry of operations: D77, D83, D88, D92, D93")
+                  "geometry of operations: D88, D92, D93, D94, D95, D96, D98, D99")
 options.register('dd4hep',
                  False,
                  VarParsing.VarParsing.multiplicity.singleton,
@@ -31,47 +31,18 @@
 ####################################################################
 # Use the options
 from Configuration.ProcessModifiers.dd4hep_cff import dd4hep
-if (options.geometry == "D83"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    if (options.dd4hep):
-        process = cms.Process('PrintG4Solids',Phase2C11M9,dd4hep)
-        process.load('Configuration.Geometry.GeometryDD4hepExtended2026D83Reco_cff')
-    else:
-        process = cms.Process('PrintG4Solids',Phase2C11M9)
-        process.load('Configuration.Geometry.GeometryExtended2026D83Reco_cff')
-elif (options.geometry == "D77"):
-    from Configuration.Eras.Era_Phase2C11_cff import Phase2C11
-    if (options.dd4hep):
-        process = cms.Process('PrintG4Solids',Phase2C11,dd4hep)
-        process.load('Configuration.Geometry.GeometryDD4hepExtended2026D77Reco_cff')
-    else:
-        process = cms.Process('PrintG4Solids',Phase2C11)
-        process.load('Configuration.Geometry.GeometryExtended2026D77Reco_cff')
-elif (options.geometry == "D92"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    if (options.dd4hep):
-        process = cms.Process('PrintG4Solids',Phase2C11M9,dd4hep)
-        process.load('Configuration.Geometry.GeometryDD4hepExtended2026D92Reco_cff')
-    else:
-        process = cms.Process('PrintG4Solids',Phase2C11M9)
-        process.load('Configuration.Geometry.GeometryExtended2026D92Reco_cff')
-elif (options.geometry == "D93"):
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    if (options.dd4hep):
-        process = cms.Process('PrintG4Solids',Phase2C11M9,dd4hep)
-        process.load('Configuration.Geometry.GeometryDD4hepExtended2026D93Reco_cff')
-    else:
-        process = cms.Process('PrintG4Solids',Phase2C11M9)
-        process.load('Configuration.Geometry.GeometryExtended2026D93Reco_cff')
+from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
+
+if (options.dd4hep):
+    process = cms.Process('PrintG4Solids',Phase2C17I13M9,dd4hep)
+    geomFile = "Configuration.Geometry.GeometryDD4hepExtended2026" + options.geometry + "Reco_cff"
 else:
-    from Configuration.Eras.Era_Phase2C11M9_cff import Phase2C11M9
-    if (options.dd4hep):
-        process = cms.Process('PrintG4Solids',Phase2C11M9,dd4hep)
-        process.load('Configuration.Geometry.GeometryDD4hepExtended2026D88Reco_cff')
-    else:
-        process = cms.Process('PrintG4Solids',Phase2C11M9)
-        process.load('Configuration.Geometry.GeometryExtended2026D88Reco_cff')
+    process = cms.Process('PrintG4Solids',Phase2C17I13M9)
+    geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
+
+print("Geometry file Name: ", geomFile)
 
+process.load(geomFile)
 process.load('FWCore.MessageService.MessageLogger_cfi')
 
 if hasattr(process,'MessageLogger'):

diff --git a/SimG4Core/PrintGeomInfo/test/python/runSens2026_cfg.py b/SimG4Core/PrintGeomInfo/test/python/runSens2026_cfg.py
@@ -0,0 +1,92 @@
+###############################################################################
+# Way to use this:
+#   cmsRun runSens2026_cfg.py geometry=D92
+#
+#   Options for geometry D88, D92, D93, D94, D95, D96, D98, D99
+#
+###############################################################################
+import FWCore.ParameterSet.Config as cms
+import os, sys, imp, re
+import FWCore.ParameterSet.VarParsing as VarParsing
+
+####################################################################
+### SETUP OPTIONS
+options = VarParsing.VarParsing('standard')
+options.register('geometry',
+                 "D92",
+                  VarParsing.VarParsing.multiplicity.singleton,
+                  VarParsing.VarParsing.varType.string,
+                  "geometry of operations: D88, D92, D93, D94, D95, D96, D98, D99")
+
+### get and parse the command line arguments
+options.parseArguments()
+
+print(options)
+
+#####p###############################################################
+# Use the options
+
+from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
+process = cms.Process('PrintSensitive',Phase2C17I13M9)
+
+geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
+
+print("Geometry file Name: ", geomFile)
+
+process.load(geomFile)
+process.load('FWCore.MessageService.MessageLogger_cfi')
+
+process.MessageLogger.cerr.enable = False
+process.MessageLogger.files.SensDet = dict(extension="txt")
+process.MessageLogger.G4cout=dict()
+
+process.maxEvents = cms.untracked.PSet(
+    input = cms.untracked.int32(1)
+)
+
+process.load('SimGeneral.HepPDTESSource.pdt_cfi')
+process.load('IOMC.EventVertexGenerators.VtxSmearedFlat_cfi')
+process.load('GeneratorInterface.Core.generatorSmeared_cfi')
+
+process.source = cms.Source("EmptySource")
+
+process.generator = cms.EDProducer("FlatRandomPtGunProducer",
+    PGunParameters = cms.PSet(
+        PartID = cms.vint32(13),
+        MinEta = cms.double(-2.5),
+        MaxEta = cms.double(2.5),
+        MinPhi = cms.double(-3.14159265359),
+        MaxPhi = cms.double(3.14159265359),
+        MinPt  = cms.double(9.99),
+        MaxPt  = cms.double(10.01)
+    ),
+    AddAntiParticle = cms.bool(False),
+    Verbosity       = cms.untracked.int32(0),
+    firstRun        = cms.untracked.uint32(1)
+)
+
+process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
+    generator = cms.PSet(
+         initialSeed = cms.untracked.uint32(123456789),
+         engineName = cms.untracked.string('HepJamesRandom')
+    ),
+    VtxSmeared = cms.PSet(
+        engineName = cms.untracked.string('HepJamesRandom'),
+        initialSeed = cms.untracked.uint32(98765432)
+    ),
+    g4SimHits = cms.PSet(
+         initialSeed = cms.untracked.uint32(11),
+         engineName = cms.untracked.string('HepJamesRandom')
+    )
+)
+
+process.load('SimG4Core.Application.g4SimHits_cfi')
+
+process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits)
+
+process.g4SimHits.UseMagneticField        = False
+process.g4SimHits.Physics.DefaultCutValue = 10. 
+process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
+	Name           = cms.untracked.string('*'),
+	type           = cms.string('PrintSensitive')
+))
diff --git a/SimG4Core/PrintGeomInfo/test/python/runSummary2026_cfg.py b/SimG4Core/PrintGeomInfo/test/python/runSummary2026_cfg.py
@@ -0,0 +1,45 @@
+###############################################################################
+# Way to use this:
+#   cmsRun runSummary2026_cfg.py geometry=D92
+#
+#   Options for geometry D88, D92, D93, D94, D95, D96, D98, D99
+#
+###############################################################################
+import FWCore.ParameterSet.Config as cms
+import os, sys, imp, re
+import FWCore.ParameterSet.VarParsing as VarParsing
+
+####################################################################
+### SETUP OPTIONS
+options = VarParsing.VarParsing('standard')
+options.register('geometry',
+                 "Run3",
+                  VarParsing.VarParsing.multiplicity.singleton,
+                  VarParsing.VarParsing.varType.string,
+                  "geometry of operations: D88, D92, D93, D94, D95, D96, D98, D99")
+
+### get and parse the command line arguments
+options.parseArguments()
+
+print(options)
+
+#####p###############################################################
+# Use the options
+
+from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
+process = cms.Process('PrintGeometry',Phase2C17I13M9)
+
+geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
+
+print("Geometry file Name: ", geomFile)
+
+process.load(geomFile)
+process.load('FWCore.MessageService.MessageLogger_cfi')
+
+if hasattr(process,'MessageLogger'):
+    process.MessageLogger.G4cout=dict()
+    process.MessageLogger.G4cerr=dict()
+
+from SimG4Core.PrintGeomInfo.g4PrintGeomSummary_cfi import *
+
+process = printGeomSummary(process)