HGCalSimAlgos/test/hgchebacksignalscaler_cfg.py

0001 #run with: cmsRun test/hgchebacksignalscaler_cfg.py doseMap=SimCalorimetry/HGCalSimProducers/data/doseParams_3000fb_fluka-3.7.20.txt geom=GeometryExtended2026D49Reco
0002
0003 import FWCore.ParameterSet.Config as cms
0004 from FWCore.ParameterSet.VarParsing import VarParsing
0005 from Configuration.StandardSequences.Eras import eras
0006
0007 options = VarParsing()
0008 options.register("doseMap",  "SimCalorimetry/HGCalSimProducers/data/doseParams_3000fb_fluka-3.7.20.txt",  VarParsing.multiplicity.singleton, VarParsing.varType.string)
0009 options.register("sipmMap",  "SimCalorimetry/HGCalSimProducers/data/sipmParams_geom-10.txt",  VarParsing.multiplicity.singleton, VarParsing.varType.string)
0010 options.register("geometry", "GeometryExtended2026D49Reco",  VarParsing.multiplicity.singleton, VarParsing.varType.string)
0011 options.parseArguments()
0012
0013 from Configuration.Eras.Era_Phase2C11I13M9_cff import Phase2C11I13M9
0014 process = cms.Process('demo',Phase2C11I13M9)
0015
0016 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0017 process.load('Configuration.Geometry.{}_cff'.format(options.geometry))
0018
0019 from Configuration.AlCa.GlobalTag import GlobalTag
0020 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic', '')
0021
0022 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1) )
0023 process.source = cms.Source("EmptySource")
0024
0025
0026 #startup: custom sipm map, no dose or fluence
0027 process.startup = cms.EDAnalyzer("HGCHEbackSignalScalerAnalyzer",
0028                                  doseMap  = cms.string( options.doseMap ),
0029                                  doseMapAlgo = cms.uint32( 2+8+16 ),
0030                                  sipmMap  = cms.string( options.sipmMap ),
0031                                  referenceIdark = cms.double( 0.125 ) )
0032
0033 #end-of-life (same radiation/sipm maps)
0034 process.eol = process.startup.clone( doseMapAlgo = cms.uint32( 2 ),
0035                                      referenceIdark = cms.double( 0.5 ) )
0036
0037 #add tfile service
0038 process.TFileService = cms.Service("TFileService",
0039                                    fileName = cms.string("sipmontile_dosemap_{}.root".format(options.geometry))
0040 )
0041
0042 process.p = cms.Path(
0043     process.startup
0044     *process.eol
0045 )
0046