EcalRecProducers/test/testEcalRechitProducer_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('FWCore.MessageService.MessageLogger_cfi')
0010 process.load('HeterogeneousCore.CUDACore.ProcessAcceleratorCUDA_cfi')
0011 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0012 process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
0013 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0014
0015 # Other statements
0016 from Configuration.AlCa.GlobalTag import GlobalTag
0017 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2022_realistic', '')
0018
0019 process.maxEvents = cms.untracked.PSet(
0020     input = cms.untracked.int32(1000)
0021 )
0022
0023 # load data using the DAQ source
0024 import sys, os, inspect
0025 sys.path.append(os.path.dirname(os.path.abspath(inspect.getfile(inspect.currentframe()))))
0026 process.load('sourceFromRawCmggpu_cff')
0027
0028 #-----------------------------------------
0029 # CMSSW/Hcal non-DQM Related Module import
0030 #-----------------------------------------
0031 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0032 process.load("RecoLocalCalo.Configuration.hcalLocalReco_cff")
0033 #process.load("RecoLocalCalo.Configuration.ecalLocalRecoSequence_cff")
0034 process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
0035 process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
0036 process.load("RecoLuminosity.LumiProducer.bunchSpacingProducer_cfi")
0037
0038 # load both cpu and gpu plugins
0039 #
0040 # ../cfipython/slc7_amd64_gcc700/RecoLocalCalo/EcalRecProducers/ecalUncalibRecHitProducerGPU_cfi.py
0041 #
0042 process.load("RecoLocalCalo.EcalRecProducers.ecalUncalibRecHitProducerGPU_cfi")
0043 process.load("RecoLocalCalo.EcalRecProducers.ecalMultiFitUncalibRecHit_cfi")
0044
0045 # for validation of gpu multifit products
0046 process.load("RecoLocalCalo.EcalRecProducers.ecalCPUUncalibRecHitProducer_cfi")
0047 process.load("EventFilter.EcalRawToDigi.ecalCPUDigisProducer_cfi")
0048
0049 process.load("EventFilter.EcalRawToDigi.ecalRawToDigiGPU_cfi")
0050 process.load("EventFilter.EcalRawToDigi.ecalElectronicsMappingGPUESProducer_cfi")
0051
0052 #process.ecalUncalibRecHitProducerGPU.kernelsVersion = 0
0053 #process.ecalUncalibRecHitProducerGPU.kernelMinimizeThreads = cms.vuint32(16, 1, 1)
0054
0055 process.load("RecoLocalCalo.EcalRecProducers.ecalPedestalsGPUESProducer_cfi")
0056 process.load("RecoLocalCalo.EcalRecProducers.ecalGainRatiosGPUESProducer_cfi")
0057 process.load("RecoLocalCalo.EcalRecProducers.ecalPulseShapesGPUESProducer_cfi")
0058 process.load("RecoLocalCalo.EcalRecProducers.ecalPulseCovariancesGPUESProducer_cfi")
0059 process.load("RecoLocalCalo.EcalRecProducers.ecalSamplesCorrelationGPUESProducer_cfi")
0060 process.load("RecoLocalCalo.EcalRecProducers.ecalTimeBiasCorrectionsGPUESProducer_cfi")
0061 process.load("RecoLocalCalo.EcalRecProducers.ecalTimeCalibConstantsGPUESProducer_cfi")
0062
0063 #process.ecalMultiFitUncalibRecHitgpu.algoPSet.threads = cms.vint32(256, 1, 1)
0064
0065 #from RecoLocalCalo.EcalRecProducers.ecalMultifitParametersGPUESProducer_cfi import ecalMultifitParametersGPUESProducer
0066 process.load("RecoLocalCalo.EcalRecProducers.ecalMultifitParametersGPUESProducer_cfi")
0067
0068 #
0069 #
0070 #   No "JobConfigurationGPURecord" record found in the EventSetup.n
0071 #    #--->
0072 #
0073 process.load("RecoLocalCalo.EcalRecProducers.ecalRecHitParametersGPUESProducer_cfi")
0074 #ecalRecHitParametersGPUESProducer_cfi.py
0075
0076
0077 ##
0078 ## force HLT configuration for ecalMultiFitUncalibRecHit
0079 ##
0080
0081 process.ecalMultiFitUncalibRecHit.algoPSet = cms.PSet(
0082       EBtimeFitLimits_Upper = cms.double( 1.4 ),
0083       EEtimeFitLimits_Lower = cms.double( 0.2 ),
0084       timealgo = cms.string( "None" ),
0085       EBtimeNconst = cms.double( 28.5 ),
0086       prefitMaxChiSqEE = cms.double( 10.0 ),
0087       outOfTimeThresholdGain12mEB = cms.double( 5.0 ),
0088       outOfTimeThresholdGain12mEE = cms.double( 1000.0 ),
0089       EEtimeFitParameters = cms.vdouble( -2.390548, 3.553628, -17.62341, 67.67538, -133.213, 140.7432, -75.41106, 16.20277 ),
0090       prefitMaxChiSqEB = cms.double( 25.0 ),
0091       simplifiedNoiseModelForGainSwitch = cms.bool( True ),
0092       EBtimeFitParameters = cms.vdouble( -2.015452, 3.130702, -12.3473, 41.88921, -82.83944, 91.01147, -50.35761, 11.05621 ),
0093       selectiveBadSampleCriteriaEB = cms.bool( False ),
0094       dynamicPedestalsEB = cms.bool( False ),
0095       useLumiInfoRunHeader = cms.bool( False ),
0096       EBamplitudeFitParameters = cms.vdouble( 1.138, 1.652 ),
0097       doPrefitEE = cms.bool( False ),
0098       dynamicPedestalsEE = cms.bool( False ),
0099       selectiveBadSampleCriteriaEE = cms.bool( False ),
0100       outOfTimeThresholdGain61pEE = cms.double( 1000.0 ),
0101       outOfTimeThresholdGain61pEB = cms.double( 5.0 ),
0102       activeBXs = cms.vint32( -5, -4, -3, -2, -1, 0, 1, 2, 3, 4 ),
0103       doPrefitEB = cms.bool( False ),
0104       addPedestalUncertaintyEE = cms.double( 0.0 ),
0105       addPedestalUncertaintyEB = cms.double( 0.0 ),
0106       gainSwitchUseMaxSampleEB = cms.bool( True ),
0107       EEtimeNconst = cms.double( 31.8 ),
0108       EEamplitudeFitParameters = cms.vdouble( 1.89, 1.4 ),
0109       outOfTimeThresholdGain12pEB = cms.double( 5.0 ),
0110       gainSwitchUseMaxSampleEE = cms.bool( False ),
0111       mitigateBadSamplesEB = cms.bool( False ),
0112       outOfTimeThresholdGain12pEE = cms.double( 1000.0 ),
0113       ampErrorCalculation = cms.bool( False ),
0114       mitigateBadSamplesEE = cms.bool( False ),
0115       amplitudeThresholdEB = cms.double( 10.0 ),
0116       amplitudeThresholdEE = cms.double( 10.0 ),
0117       EBtimeFitLimits_Lower = cms.double( 0.2 ),
0118       EEtimeFitLimits_Upper = cms.double( 1.4 ),
0119       outOfTimeThresholdGain61mEE = cms.double( 1000.0 ),
0120       EEtimeConstantTerm = cms.double( 1.0 ),
0121       EBtimeConstantTerm = cms.double( 0.6 ),
0122       outOfTimeThresholdGain61mEB = cms.double( 5.0 )
0123 )
0124
0125 ##
0126
0127
0128
0129 process.load('Configuration.StandardSequences.Reconstruction_cff')
0130 #process.ecalRecHit
0131
0132
0133 process.load("RecoLocalCalo.EcalRecProducers.ecalRechitADCToGeVConstantGPUESProducer_cfi")
0134 process.load("RecoLocalCalo.EcalRecProducers.ecalRechitChannelStatusGPUESProducer_cfi")
0135 #process.load("RecoLocalCalo.EcalRecProducers.ecalADCToGeVConstantGPUESProducer_cfi")
0136 #process.load("RecoLocalCalo.EcalRecProducers.ecalChannelStatusGPUESProducer_cfi")
0137 process.load("RecoLocalCalo.EcalRecProducers.ecalIntercalibConstantsGPUESProducer_cfi")
0138
0139 process.load("RecoLocalCalo.EcalRecProducers.ecalLaserAPDPNRatiosGPUESProducer_cfi")
0140 process.load("RecoLocalCalo.EcalRecProducers.ecalLaserAPDPNRatiosRefGPUESProducer_cfi")
0141 process.load("RecoLocalCalo.EcalRecProducers.ecalLaserAlphasGPUESProducer_cfi")
0142 process.load("RecoLocalCalo.EcalRecProducers.ecalLinearCorrectionsGPUESProducer_cfi")
0143
0144 process.load("RecoLocalCalo.EcalRecProducers.ecalRecHitGPU_cfi")
0145 process.ecalRecHitProducerGPU = process.ecalRecHitGPU.clone()
0146
0147
0148 process.load("RecoLocalCalo.EcalRecProducers.ecalCPURecHitProducer_cfi")
0149
0150
0151 #
0152 # AM : TEST to see if the number of rechits matches
0153 #
0154 #process.ecalRecHit.killDeadChannels = cms.bool(False)
0155 #
0156 #process.ecalRecHit.recoverEBFE = cms.bool(False)
0157 #process.ecalRecHit.recoverEBIsolatedChannels = cms.bool(False)
0158 #process.ecalRecHit.recoverEBVFE = cms.bool(False)
0159 ##
0160 #process.ecalRecHit.recoverEEFE = cms.bool(False)
0161 #process.ecalRecHit.recoverEEIsolatedChannels = cms.bool(False)
0162 #process.ecalRecHit.recoverEEVFE = cms.bool(False)
0163 #
0164 #process.ecalRecHit.skipTimeCalib = cms.bool(True)
0165 #
0166 #process.ecalRecHitProducerGPU.killDeadChannels = cms.bool(False)
0167 #
0168 #
0169 #process.ecalRecHitProducerGPU.recoverEBFE = cms.bool(False)
0170 #process.ecalRecHitProducerGPU.recoverEBIsolatedChannels = cms.bool(False)
0171 #process.ecalRecHitProducerGPU.recoverEBVFE = cms.bool(False)
0172 #process.ecalRecHitProducerGPU.recoverEEFE = cms.bool(False)
0173 #process.ecalRecHitProducerGPU.recoverEEIsolatedChannels = cms.bool(False)
0174 #process.ecalRecHitProducerGPU.recoverEEVFE = cms.bool(False)
0175 #
0176 #
0177 #
0178 # TEST
0179 #
0180 #process.ecalRecHit.ChannelStatusToBeExcluded = cms.vstring(
0181                                                           #'kDAC',
0182                                                           #'kNoisy',
0183                                                           #'kNNoisy',
0184                                                           #'kFixedG6',
0185                                                           #'kFixedG1',
0186                                                           #'kFixedG0',
0187                                                           #'kNonRespondingIsolated',
0188                                                           #'kDeadVFE',
0189                                                           #'kDeadFE',
0190                                                           #'kNoDataNoTP'
0191                                                           #)
0192 #process.ecalRecHitProducerGPU.ChannelStatusToBeExcluded = cms.vstring(
0193                                                           #'kDAC',
0194                                                           #'kNoisy',
0195                                                           #'kNNoisy',
0196                                                           #'kFixedG6',
0197                                                           #'kFixedG1',
0198                                                           #'kFixedG0',
0199                                                           #'kNonRespondingIsolated',
0200                                                           #'kDeadVFE',
0201                                                           #'kDeadFE',
0202                                                           #'kNoDataNoTP'
0203                                                           #)
0204 #
0205 #
0206
0207     #ChannelStatusToBeExcluded = cms.vstring(
0208         #'kDAC',
0209         #'kNoisy',
0210         #'kNNoisy',
0211         #'kFixedG6',
0212         #'kFixedG1',
0213         #'kFixedG0',
0214         #'kNonRespondingIsolated',
0215         #'kDeadVFE',
0216         #'kDeadFE',
0217         #'kNoDataNoTP'
0218     #),
0219
0220
0221
0222 #process.hcalDigis.silent = cms.untracked.bool(False)
0223 #process.hcalDigis.InputLabel = rawTag
0224 process.ecalDigis = process.ecalEBunpacker.clone()
0225 process.ecalDigis.InputLabel = cms.InputTag('rawDataCollector')
0226 #process.hbheprerecogpu.processQIE11 = cms.bool(True)
0227
0228 process.out = cms.OutputModule(
0229     "PoolOutputModule",
0230     fileName = cms.untracked.string("testRechit.root")
0231 )
0232
0233 #process.out = cms.OutputModule("AsciiOutputModule",
0234 #    outputCommands = cms.untracked.vstring(
0235 #        'keep *_ecalMultiFitUncalibRecHit_*_*',
0236 #    ),
0237 #    verbosity = cms.untracked.uint32(0)
0238 #)
0239 process.finalize = cms.EndPath(process.out)
0240
0241 process.bunchSpacing = cms.Path(
0242     process.bunchSpacingProducer
0243 )
0244
0245 process.digiPath = cms.Path(
0246     #process.hcalDigis
0247     process.ecalDigis
0248     *process.ecalRawToDigiGPU
0249     *process.ecalCPUDigisProducer
0250 )
0251
0252 process.recoPath = cms.Path(
0253     (process.ecalMultiFitUncalibRecHit+process.ecalDetIdToBeRecovered)
0254     #process.ecalMultiFitUncalibRecHit
0255     *process.ecalRecHit
0256 #   gpu
0257     *process.ecalUncalibRecHitProducerGPU
0258     *process.ecalCPUUncalibRecHitProducer
0259     *process.ecalRecHitProducerGPU
0260     *process.ecalCPURecHitProducer
0261 )
0262
0263 process.schedule = cms.Schedule(
0264     process.bunchSpacing,
0265     process.digiPath,
0266     process.recoPath,
0267 #    process.ecalecalLocalRecoSequence
0268     process.finalize
0269 )
0270
0271 process.options = cms.untracked.PSet(
0272     numberOfThreads = cms.untracked.uint32(4),
0273     numberOfStreams = cms.untracked.uint32(4),
0274     TryToContinue = cms.untracked.vstring('ProductNotFound'),
0275     wantSummary = cms.untracked.bool(True)
0276 )
0277
0278 #
0279 #process.DependencyGraph = cms.Service("DependencyGraph")
0280
0281