Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2025-06-06 01:33:42

0001 ###############################################################################
0002 # Way to use this:
0003 #   cmsRun g4OverlapCheckCalo_cfg.py geometry=2021 tol=0.1
0004 #
0005 #   Options for geometry 2016, 2017, 2021, Run4D102, Run4D103, Run4D104,
0006 #                        Run4D108, Run4D109, Run4D110, Run4D111, Run4D112, 
0007 #                        Run4D113, Run4D114, Run4D115, Run4D116, Run4D117,
0008 #                        Run4D118, Run4D119, Run4D120, Run4D122, Run4D123
0009 #
0010 ###############################################################################
0011 import FWCore.ParameterSet.Config as cms
0012 import os, sys, importlib, re
0013 import FWCore.ParameterSet.VarParsing as VarParsing
0014 
0015 ####################################################################
0016 ### SETUP OPTIONS
0017 options = VarParsing.VarParsing('standard')
0018 options.register('geometry',
0019                  "2021",
0020                   VarParsing.VarParsing.multiplicity.singleton,
0021                   VarParsing.VarParsing.varType.string,
0022                   "geometry of operations: 2016, 2017, 2021, Run4D102, Run4D103, Run4D104, Run4D108, Run4D109, Run4D110, Run4D111, Run4D112, Run4D113, Run4D114, Run4D115, Run4D116, Run4D117, Run4D118, Run4D119, Run4D120, Run4D122, Run4D123")
0023 options.register('tol',
0024                  0.01,
0025                  VarParsing.VarParsing.multiplicity.singleton,
0026                  VarParsing.VarParsing.varType.float,
0027                  "Tolerance for checking overlaps: 0.0, 0.01, 0.1, 1.0"
0028 )
0029 
0030 ### get and parse the command line arguments
0031 options.parseArguments()
0032 
0033 print(options)
0034 
0035 ####################################################################
0036 # Use the options
0037 
0038 if (options.geometry == "2016"):
0039     from Configuration.Eras.Era_Run2_2016_cff import Run2_2016
0040     process = cms.Process('G4PrintGeometry',Run2_2016)
0041     process.load('Configuration.Geometry.GeometryExtended2016Reco_cff')
0042     baseName = 'Calo2016'
0043 elif (options.geometry == "2017"):
0044     from Configuration.Eras.Era_Run2_2017_cff import Run2_2017
0045     process = cms.Process('G4PrintGeometry',Run2_2017)
0046     process.load('Configuration.Geometry.GeometryExtended2017Reco_cff')
0047     baseName = 'Calo2017'
0048 elif (options.geometry == "2021"):
0049     from Configuration.Eras.Era_Run3_DDD_cff import Run3_DDD
0050     process = cms.Process('G4PrintGeometry',Run3_DDD)
0051     process.load('Configuration.Geometry.GeometryExtended2021Reco_cff')
0052     baseName = 'Calo2021'
0053 else:
0054     geomName = options.geometry
0055     geomFile = "Configuration.Geometry.GeometryExtended" + geomName + "Reco_cff"
0056     baseName = "Calo" + geomName
0057     import Configuration.Geometry.defaultPhase2ConditionsEra_cff as _settings
0058     GLOBAL_TAG, ERA = _settings.get_era_and_conditions(geomName)
0059     print("Geometry Name:   ", geomName)
0060     print("Geom file Name:  ", geomFile)
0061     print("Global Tag Name: ", GLOBAL_TAG)
0062     print("Era Name:        ", ERA)
0063     process = cms.Process('G4PrintGeometry',ERA)
0064     process.load(geomFile)
0065 
0066 print("Base file Name: ", baseName)
0067 
0068 from SimG4Core.PrintGeomInfo.g4TestGeometry_cfi import *
0069 process = checkOverlap(process)
0070 
0071 # enable Geant4 overlap check 
0072 process.g4SimHits.CheckGeometry = True
0073 
0074 # Geant4 geometry check 
0075 process.g4SimHits.G4CheckOverlap.OutputBaseName = cms.string(baseName)
0076 process.g4SimHits.G4CheckOverlap.OverlapFlag = cms.bool(True)
0077 process.g4SimHits.G4CheckOverlap.Tolerance  = cms.double(options.tol)
0078 process.g4SimHits.G4CheckOverlap.Resolution = cms.int32(10000)
0079 process.g4SimHits.G4CheckOverlap.Depth      = cms.int32(-1)
0080 # tells if NodeName is G4Region or G4PhysicalVolume
0081 process.g4SimHits.G4CheckOverlap.RegionFlag = cms.bool(True)
0082 # list of names
0083 process.g4SimHits.G4CheckOverlap.NodeNames  = cms.vstring('HcalRegion')
0084 # enable dump gdml file 
0085 process.g4SimHits.G4CheckOverlap.gdmlFlag   = cms.bool(True)
0086 # if defined a G4PhysicsVolume info is printed
0087 process.g4SimHits.G4CheckOverlap.PVname     = ''
0088 # if defined a list of daughter volumes is printed
0089 process.g4SimHits.G4CheckOverlap.LVname     = 'HCAL'
0090 
0091 # extra output files, created if a name is not empty
0092 process.g4SimHits.FileNameField   = ''
0093 process.g4SimHits.FileNameGDML    = ''
0094 process.g4SimHits.FileNameRegions = ''
0095 #