EcalRecProducers/test/testEcalUncalibRechitProducer_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Configuration.StandardSequences.Eras import eras
0004
0005 process = cms.Process('RECO', eras.Run2_2018)
0006
0007 # import of standard configurations
0008 process.load('Configuration.StandardSequences.Services_cff')
0009 #process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0010 process.load('FWCore.MessageService.MessageLogger_cfi')
0011 process.load('HeterogeneousCore.CUDACore.ProcessAcceleratorCUDA_cfi')
0012 #process.load('Configuration.EventContent.EventContent_cff')
0013 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0014 process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
0015 #process.load('Configuration.StandardSequences.RawToDigi_Data_cff')
0016 #process.load('Configuration.StandardSequences.Reconstruction_Data_cff')
0017 #process.load('DQMOffline.Configuration.DQMOffline_cff')
0018 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0019
0020
0021
0022
0023
0024
0025 # Other statements
0026 from Configuration.AlCa.GlobalTag import GlobalTag
0027 process.GlobalTag = GlobalTag(process.GlobalTag, '102X_dataRun2_HLT_v2', '')
0028
0029
0030 process.maxEvents = cms.untracked.PSet(
0031     #input = cms.untracked.int32(100)
0032     input = cms.untracked.int32(1000)
0033 )
0034
0035 # load data using the DAQ source
0036 import sys, os, inspect
0037 sys.path.append(os.path.dirname(os.path.abspath(inspect.getfile(inspect.currentframe()))))
0038 process.load('sourceFromRawCmggpu_cff')
0039
0040 #-----------------------------------------
0041 # CMSSW/Hcal non-DQM Related Module import
0042 #-----------------------------------------
0043 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0044 process.load("RecoLocalCalo.Configuration.hcalLocalReco_cff")
0045 #process.load("RecoLocalCalo.Configuration.ecalLocalRecoSequence_cff")
0046 process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
0047 process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
0048 process.load("RecoLuminosity.LumiProducer.bunchSpacingProducer_cfi")
0049
0050 # load both cpu and gpu plugins
0051 #
0052 # ../cfipython/slc7_amd64_gcc700/RecoLocalCalo/EcalRecProducers/ecalUncalibRecHitProducerGPU_cfi.py
0053 #
0054 process.load("RecoLocalCalo.EcalRecProducers.ecalUncalibRecHitProducerGPU_cfi")
0055 process.load("RecoLocalCalo.EcalRecProducers.ecalMultiFitUncalibRecHit_cfi")
0056
0057 # for validation of gpu multifit products
0058 process.load("RecoLocalCalo.EcalRecProducers.ecalCPUUncalibRecHitProducer_cfi")
0059 #
0060 # ../cfipython/slc7_amd64_gcc700/RecoLocalCalo/EcalRecProducers/ecalCPUUncalibRecHitProducer_cfi.py
0061 #
0062
0063 process.load("EventFilter.EcalRawToDigi.ecalRawToDigiGPU_cfi")
0064 process.load("EventFilter.EcalRawToDigi.ecalElectronicsMappingGPUESProducer_cfi")
0065
0066 #process.ecalUncalibRecHitProducerGPU.kernelsVersion = 0
0067 #process.ecalUncalibRecHitProducerGPU.kernelMinimizeThreads = cms.vuint32(16, 1, 1)
0068 #
0069 # process.ecalUncalibRecHitProducerGPU.shouldRunTimingComputation = cms.bool(False)
0070 #
0071
0072
0073 process.load("RecoLocalCalo.EcalRecProducers.ecalPedestalsGPUESProducer_cfi")
0074 process.load("RecoLocalCalo.EcalRecProducers.ecalGainRatiosGPUESProducer_cfi")
0075 process.load("RecoLocalCalo.EcalRecProducers.ecalPulseShapesGPUESProducer_cfi")
0076 process.load("RecoLocalCalo.EcalRecProducers.ecalPulseCovariancesGPUESProducer_cfi")
0077 process.load("RecoLocalCalo.EcalRecProducers.ecalSamplesCorrelationGPUESProducer_cfi")
0078 process.load("RecoLocalCalo.EcalRecProducers.ecalTimeBiasCorrectionsGPUESProducer_cfi")
0079 process.load("RecoLocalCalo.EcalRecProducers.ecalTimeCalibConstantsGPUESProducer_cfi")
0080
0081 #process.ecalMultiFitUncalibRecHitgpu.algoPSet.threads = cms.vint32(256, 1, 1)
0082
0083
0084 process.load("RecoLocalCalo.EcalRecProducers.ecalMultifitParametersGPUESProducer_cfi")
0085
0086
0087
0088 ##
0089 ## force HLT configuration for ecalMultiFitUncalibRecHit
0090 ##
0091
0092 process.ecalMultiFitUncalibRecHit.algoPSet = cms.PSet(
0093       EBtimeFitLimits_Upper = cms.double( 1.4 ),
0094       EEtimeFitLimits_Lower = cms.double( 0.2 ),
0095       timealgo = cms.string( "None" ),   # ----> no timing computation for CPU version
0096       EBtimeNconst = cms.double( 28.5 ),
0097       prefitMaxChiSqEE = cms.double( 10.0 ),
0098       outOfTimeThresholdGain12mEB = cms.double( 5.0 ),
0099       outOfTimeThresholdGain12mEE = cms.double( 1000.0 ),
0100       EEtimeFitParameters = cms.vdouble( -2.390548, 3.553628, -17.62341, 67.67538, -133.213, 140.7432, -75.41106, 16.20277 ),
0101       prefitMaxChiSqEB = cms.double( 25.0 ),
0102       simplifiedNoiseModelForGainSwitch = cms.bool( True ),
0103       EBtimeFitParameters = cms.vdouble( -2.015452, 3.130702, -12.3473, 41.88921, -82.83944, 91.01147, -50.35761, 11.05621 ),
0104       selectiveBadSampleCriteriaEB = cms.bool( False ),
0105       dynamicPedestalsEB = cms.bool( False ),
0106       useLumiInfoRunHeader = cms.bool( False ),
0107       EBamplitudeFitParameters = cms.vdouble( 1.138, 1.652 ),
0108       doPrefitEE = cms.bool( False ),
0109       dynamicPedestalsEE = cms.bool( False ),
0110       selectiveBadSampleCriteriaEE = cms.bool( False ),
0111       outOfTimeThresholdGain61pEE = cms.double( 1000.0 ),
0112       outOfTimeThresholdGain61pEB = cms.double( 5.0 ),
0113       activeBXs = cms.vint32( -5, -4, -3, -2, -1, 0, 1, 2, 3, 4 ),
0114       doPrefitEB = cms.bool( False ),
0115       addPedestalUncertaintyEE = cms.double( 0.0 ),
0116       addPedestalUncertaintyEB = cms.double( 0.0 ),
0117       gainSwitchUseMaxSampleEB = cms.bool( True ),
0118       EEtimeNconst = cms.double( 31.8 ),
0119       EEamplitudeFitParameters = cms.vdouble( 1.89, 1.4 ),
0120       outOfTimeThresholdGain12pEB = cms.double( 5.0 ),
0121       gainSwitchUseMaxSampleEE = cms.bool( False ),
0122       mitigateBadSamplesEB = cms.bool( False ),
0123       outOfTimeThresholdGain12pEE = cms.double( 1000.0 ),
0124       ampErrorCalculation = cms.bool( False ),
0125       mitigateBadSamplesEE = cms.bool( False ),
0126       amplitudeThresholdEB = cms.double( 10.0 ),
0127       amplitudeThresholdEE = cms.double( 10.0 ),
0128       EBtimeFitLimits_Lower = cms.double( 0.2 ),
0129       EEtimeFitLimits_Upper = cms.double( 1.4 ),
0130       outOfTimeThresholdGain61mEE = cms.double( 1000.0 ),
0131       EEtimeConstantTerm = cms.double( 1.0 ),
0132       EBtimeConstantTerm = cms.double( 0.6 ),
0133       outOfTimeThresholdGain61mEB = cms.double( 5.0 )
0134 )
0135
0136 ##
0137
0138
0139
0140 #process.load('Configuration.StandardSequences.Reconstruction_cff')
0141 #process.ecalRecHit
0142
0143
0144
0145 #process.load("RecoLocalCalo.EcalRecProducers.ecalRecHitGPU_cfi")
0146 #process.ecalRecHitGPU
0147
0148
0149
0150 #process.hcalDigis.silent = cms.untracked.bool(False)
0151 #process.hcalDigis.InputLabel = rawTag
0152 process.ecalDigis = process.ecalEBunpacker.clone()
0153 process.ecalDigis.InputLabel = cms.InputTag('rawDataCollector')
0154 #process.hbheprerecogpu.processQIE11 = cms.bool(True)
0155
0156 process.out = cms.OutputModule(
0157     "PoolOutputModule",
0158     fileName = cms.untracked.string("test_uncalib.root")
0159 )
0160
0161 #process.out = cms.OutputModule("AsciiOutputModule",
0162 #    outputCommands = cms.untracked.vstring(
0163 #        'keep *_ecalMultiFitUncalibRecHit_*_*',
0164 #    ),
0165 #    verbosity = cms.untracked.uint32(0)
0166 #)
0167 process.finalize = cms.EndPath(process.out)
0168
0169 process.bunchSpacing = cms.Path(
0170     process.bunchSpacingProducer
0171 )
0172
0173 process.digiPath = cms.Path(
0174     #process.hcalDigis
0175     process.ecalDigis
0176     *process.ecalRawToDigiGPU
0177 )
0178
0179 process.recoPath = cms.Path(
0180     #(process.ecalMultiFitUncalibRecHit+process.ecalDetIdToBeRecovered)
0181     process.ecalMultiFitUncalibRecHit
0182     #*process.ecalRecHit
0183 #   gpu
0184     *process.ecalUncalibRecHitProducerGPU
0185     *process.ecalCPUUncalibRecHitProducer
0186     #*process.ecalRecHitGPU
0187 )
0188
0189 process.schedule = cms.Schedule(
0190     process.bunchSpacing,
0191     process.digiPath,
0192     process.recoPath,
0193 #    process.ecalecalLocalRecoSequence
0194     process.finalize
0195 )
0196
0197 process.options = cms.untracked.PSet(
0198     numberOfThreads = cms.untracked.uint32(8),
0199     numberOfStreams = cms.untracked.uint32(8),
0200     TryToContinue = cms.untracked.vstring('ProductNotFound'),
0201     wantSummary = cms.untracked.bool(True)
0202 )