Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-11-07 23:58:25

0001 #run with: cmsRun hgcsiNoiseMapTester_cfg.py doseMap=SimCalorimetry/HGCalSimProducers/data/doseParams_3000fb_fluka-3.7.20.txt geometry=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 ("geometry",   "GeometryExtended2026D49Reco",  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)