EcalRecProducers/test/ecalRawDecodingAndMultifit.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.StandardSequences.Eras import eras
0003
0004 process = cms.Process('RECO', eras.Run2_2018)
0005
0006 # import of standard configurations
0007 process.load('Configuration.StandardSequences.Services_cff')
0008 process.load('FWCore.MessageService.MessageLogger_cfi')
0009 process.load('HeterogeneousCore.CUDACore.ProcessAcceleratorCUDA_cfi')
0010 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0011 process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
0012 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0013
0014 # Other statements
0015 from Configuration.AlCa.GlobalTag import GlobalTag
0016 process.GlobalTag = GlobalTag(process.GlobalTag, '102X_dataRun2_HLT_v2', '')
0017
0018
0019 process.maxEvents = cms.untracked.PSet(
0020     input = cms.untracked.int32(100)
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 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0029 process.load("RecoLocalCalo.Configuration.hcalLocalReco_cff")
0030 process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
0031 process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
0032 process.load("RecoLuminosity.LumiProducer.bunchSpacingProducer_cfi")
0033
0034 # load both cpu and gpu plugins
0035 process.load("RecoLocalCalo.EcalRecProducers.ecalUncalibRecHitProducerGPU_cfi")
0036 process.load("RecoLocalCalo.EcalRecProducers.ecalMultiFitUncalibRecHit_cfi")
0037
0038 # for validation of gpu multifit products
0039 process.load("RecoLocalCalo.EcalRecProducers.ecalCPUUncalibRecHitProducer_cfi")
0040 process.load("EventFilter.EcalRawToDigi.ecalCPUDigisProducer_cfi")
0041
0042 process.load("EventFilter.EcalRawToDigi.ecalRawToDigiGPU_cfi")
0043 process.load("EventFilter.EcalRawToDigi.ecalElectronicsMappingGPUESProducer_cfi")
0044
0045 process.load("RecoLocalCalo.EcalRecProducers.ecalPedestalsGPUESProducer_cfi")
0046 process.load("RecoLocalCalo.EcalRecProducers.ecalGainRatiosGPUESProducer_cfi")
0047 process.load("RecoLocalCalo.EcalRecProducers.ecalPulseShapesGPUESProducer_cfi")
0048 process.load("RecoLocalCalo.EcalRecProducers.ecalPulseCovariancesGPUESProducer_cfi")
0049 process.load("RecoLocalCalo.EcalRecProducers.ecalSamplesCorrelationGPUESProducer_cfi")
0050 process.load("RecoLocalCalo.EcalRecProducers.ecalTimeBiasCorrectionsGPUESProducer_cfi")
0051 process.load("RecoLocalCalo.EcalRecProducers.ecalTimeCalibConstantsGPUESProducer_cfi")
0052
0053 process.load("RecoLocalCalo.EcalRecProducers.ecalRechitADCToGeVConstantGPUESProducer_cfi")
0054 process.load("RecoLocalCalo.EcalRecProducers.ecalIntercalibConstantsGPUESProducer_cfi")
0055 process.load("RecoLocalCalo.EcalRecProducers.ecalRechitChannelStatusGPUESProducer_cfi")
0056
0057 process.load("RecoLocalCalo.EcalRecProducers.ecalLaserAPDPNRatiosGPUESProducer_cfi")
0058 process.load("RecoLocalCalo.EcalRecProducers.ecalLaserAPDPNRatiosRefGPUESProducer_cfi")
0059 process.load("RecoLocalCalo.EcalRecProducers.ecalLaserAlphasGPUESProducer_cfi")
0060 process.load("RecoLocalCalo.EcalRecProducers.ecalLinearCorrectionsGPUESProducer_cfi")
0061
0062 # force HLT configuration for ecalMultiFitUncalibRecHit
0063 process.ecalMultiFitUncalibRecHit.algoPSet = cms.PSet(
0064     EBtimeFitLimits_Upper = cms.double(1.4),
0065     EEtimeFitLimits_Lower = cms.double(0.2),
0066     timealgo = cms.string("None"),
0067     EBtimeNconst = cms.double(28.5),
0068     prefitMaxChiSqEE = cms.double(10.0),
0069     outOfTimeThresholdGain12mEB = cms.double(5.0),
0070     outOfTimeThresholdGain12mEE = cms.double(1000.0),
0071     EEtimeFitParameters = cms.vdouble(-2.390548, 3.553628, -17.62341, 67.67538, -133.213, 140.7432, -75.41106, 16.20277),
0072     prefitMaxChiSqEB = cms.double(25.0),
0073     simplifiedNoiseModelForGainSwitch = cms.bool(True),
0074     EBtimeFitParameters = cms.vdouble(-2.015452, 3.130702, -12.3473, 41.88921, -82.83944, 91.01147, -50.35761, 11.05621),
0075     selectiveBadSampleCriteriaEB = cms.bool(False),
0076     dynamicPedestalsEB = cms.bool(False),
0077     useLumiInfoRunHeader = cms.bool(False),
0078     EBamplitudeFitParameters = cms.vdouble(1.138, 1.652),
0079     doPrefitEE = cms.bool(False),
0080     dynamicPedestalsEE = cms.bool(False),
0081     selectiveBadSampleCriteriaEE = cms.bool(False),
0082     outOfTimeThresholdGain61pEE = cms.double(1000.0),
0083     outOfTimeThresholdGain61pEB = cms.double(5.0),
0084     activeBXs = cms.vint32(-5, -4, -3, -2, -1, 0, 1, 2, 3, 4),
0085     doPrefitEB = cms.bool(False),
0086     addPedestalUncertaintyEE = cms.double(0.0),
0087     addPedestalUncertaintyEB = cms.double(0.0),
0088     gainSwitchUseMaxSampleEB = cms.bool(True),
0089     EEtimeNconst = cms.double(31.8),
0090     EEamplitudeFitParameters = cms.vdouble(1.89, 1.4),
0091     outOfTimeThresholdGain12pEB = cms.double(5.0),
0092     gainSwitchUseMaxSampleEE = cms.bool(False),
0093     mitigateBadSamplesEB = cms.bool(False),
0094     outOfTimeThresholdGain12pEE = cms.double(1000.0),
0095     ampErrorCalculation = cms.bool(False),
0096     mitigateBadSamplesEE = cms.bool(False),
0097     amplitudeThresholdEB = cms.double(10.0),
0098     amplitudeThresholdEE = cms.double(10.0),
0099     EBtimeFitLimits_Lower = cms.double(0.2),
0100     EEtimeFitLimits_Upper = cms.double(1.4),
0101     outOfTimeThresholdGain61mEE = cms.double(1000.0),
0102     EEtimeConstantTerm = cms.double(1.0),
0103     EBtimeConstantTerm = cms.double(0.6),
0104     outOfTimeThresholdGain61mEB = cms.double(5.0)
0105 )
0106
0107 process.ecalDigis = process.ecalEBunpacker.clone()
0108 process.ecalDigis.InputLabel = cms.InputTag('rawDataCollector')
0109
0110 process.out = cms.OutputModule("PoolOutputModule",
0111     fileName = cms.untracked.string("test.root")
0112 )
0113
0114 process.finalize = cms.EndPath(process.out)
0115
0116 process.bunchSpacing = cms.Path(
0117     process.bunchSpacingProducer
0118 )
0119
0120 process.digiPath = cms.Path(
0121     process.ecalDigis +
0122     process.ecalRawToDigiGPU +
0123     process.ecalCPUDigisProducer
0124 )
0125
0126 process.recoPath = cms.Path(
0127     process.ecalMultiFitUncalibRecHit +
0128     process.ecalUncalibRecHitProducerGPU +
0129     process.ecalCPUUncalibRecHitProducer
0130 )
0131
0132 process.schedule = cms.Schedule(
0133     process.bunchSpacing,
0134     process.digiPath,
0135     process.recoPath,
0136     process.finalize
0137 )
0138
0139 process.options = cms.untracked.PSet(
0140     numberOfThreads = cms.untracked.uint32(4),
0141     numberOfStreams = cms.untracked.uint32(4),
0142     TryToContinue = cms.untracked.vstring('ProductNotFound'),
0143     wantSummary = cms.untracked.bool(True)
0144 )