File indexing completed on 2024-11-06 06:06:57
0001
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 ("geometry", "GeometryExtendedRun4D49Reco", VarParsing.multiplicity.singleton, VarParsing.varType.string)
0010 options.register ("conditions", "TDR_600V", 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 from Configuration.AlCa.GlobalTag import GlobalTag
0019 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic', '')
0020
0021 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1) )
0022 process.source = cms.Source("EmptySource")
0023
0024 from SimCalorimetry.HGCalSimAlgos.hgcSensorOpParams_cfi import hgcSiSensorIleak,hgcSiSensorCCE
0025 from SimCalorimetry.HGCalSimProducers.hgcalDigitizer_cfi import HGCAL_ileakParam_toUse, HGCAL_cceParams_toUse
0026
0027 HGCAL_ileakParam_toUse = cms.PSet(
0028 ileakParam = cms.vdouble( hgcSiSensorIleak(options.conditions) )
0029 )
0030
0031 HGCAL_cceParams_toUse = cms.PSet(
0032 cceParamFine = cms.vdouble(hgcSiSensorCCE(120,options.conditions)),
0033 cceParamThin = cms.vdouble(hgcSiSensorCCE(200,options.conditions)),
0034 cceParamThick = cms.vdouble(hgcSiSensorCCE(300,options.conditions)),
0035 )
0036
0037 process.plotter_eol = cms.EDAnalyzer("HGCSiNoiseMapAnalyzer",
0038 scaleByDoseFactor = cms.double(1.),
0039 doseMap = cms.string( options.doseMap ),
0040 doseMapAlgo = cms.uint32(0),
0041 ileakParam = HGCAL_ileakParam_toUse,
0042 cceParams = HGCAL_cceParams_toUse,
0043 aimMIPtoADC = cms.int32(10),
0044 ignoreGainSettings = cms.bool(False)
0045 )
0046
0047 process.plotter_eol_nogain = process.plotter_eol.clone( ignoreGainSettings = cms.bool(True) )
0048
0049 process.plotter_start = process.plotter_eol.clone( doseMapAlgo=cms.uint32(3) )
0050
0051 process.TFileService = cms.Service("TFileService",
0052 fileName = cms.string("dosemap_output_{}_{}.root".format(options.geometry,options.conditions))
0053 )
0054
0055 process.p = cms.Path(process.plotter_eol
0056 *process.plotter_eol_nogain
0057 *process.plotter_start)