Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:24:24

0001 import FWCore.ParameterSet.Config as cms
0002 process = cms.Process("testEcalRecoLocal")
0003 
0004 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0005 
0006 # get timing service up for profiling
0007 process.TimerService = cms.Service("TimerService")
0008 process.options = cms.untracked.PSet(
0009     wantSummary = cms.untracked.bool(True)
0010 )
0011 
0012 
0013 # unpacking
0014 process.load("EventFilter.EcalRawToDigi.EcalUnpackerMapping_cfi")
0015 process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
0016 
0017 # load trivial conditions
0018 process.EcalTrivialConditionRetriever = cms.ESSource("EcalTrivialConditionRetriever",
0019                                                          adcToGeVEBConstant = cms.untracked.double(0.035),
0020                                                          adcToGeVEEConstant = cms.untracked.double(0.06),
0021                                                          pedWeights = cms.untracked.vdouble(0.333, 0.333, 0.333, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0),
0022                                                          amplWeights = cms.untracked.vdouble(-0.333, -0.333, -0.333, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0),
0023                                                          jittWeights = cms.untracked.vdouble(0.041, 0.041, 0.041, 0.0, 1.325, -0.05, -0.504, -0.502, -0.390, 0.0)
0024                                                      )
0025 
0026 
0027 # get uncalibrechits with weights method
0028 import RecoLocalCalo.EcalRecProducers.ecalWeightUncalibRecHit_cfi
0029 process.ecalUncalibHitWeights = RecoLocalCalo.EcalRecProducers.ecalWeightUncalibRecHit_cfi.ecalWeightUncalibRecHit.clone()
0030 process.ecalUncalibHitWeights.EBdigiCollection = 'ecalEBunpacker:ebDigis'
0031 process.ecalUncalibHitWeights.EEdigiCollection = 'ecalEBunpacker:eeDigis'
0032 
0033 # get uncalibrechits with fit method
0034 import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi
0035 process.ecalUncalibHitFixedAlphaBetaFit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone()
0036 process.ecalUncalibHitFixedAlphaBetaFit.EBdigiCollection = 'ecalEBunpacker:ebDigis'
0037 process.ecalUncalibHitFixedAlphaBetaFit.EEdigiCollection = 'ecalEBunpacker:eeDigis'
0038 
0039 # get uncalibrechits with ratio method
0040 import RecoLocalCalo.EcalRecProducers.ecalRatioUncalibRecHit_cfi
0041 process.ecalUncalibHitRatio = RecoLocalCalo.EcalRecProducers.ecalRatioUncalibRecHit_cfi.ecalRatioUncalibRecHit.clone()
0042 process.ecalUncalibHitRatio.EBdigiCollection = 'ecalEBunpacker:ebDigis'
0043 process.ecalUncalibHitRatio.EEdigiCollection = 'ecalEBunpacker:eeDigis'
0044 
0045 
0046 # get uncalibrechits with ratio method
0047 import RecoLocalCalo.EcalRecProducers.ecalGlobalUncalibRecHit_cfi
0048 process.ecalUncalibHitGlobal = RecoLocalCalo.EcalRecProducers.ecalGlobalUncalibRecHit_cfi.ecalGlobalUncalibRecHit.clone()
0049 process.ecalUncalibHitGlobal.EBdigiCollection = 'ecalEBunpacker:ebDigis'
0050 process.ecalUncalibHitGlobal.EEdigiCollection = 'ecalEBunpacker:eeDigis'
0051 
0052 
0053 # get rechits e.g. from the weights
0054 process.load("CalibCalorimetry.EcalLaserCorrection.ecalLaserCorrectionService_cfi")
0055 process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi")
0056 process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB'
0057 process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE'
0058 
0059 
0060 process.maxEvents = cms.untracked.PSet(  input = cms.untracked.int32(10) )
0061 process.source = cms.Source("PoolSource",
0062             fileNames = cms.untracked.vstring('/store/data/Commissioning08/Cosmics/RAW/v1/000/067/838/006945C8-40A5-DD11-BD7E-001617DBD556.root')
0063             #fileNames = cms.untracked.vstring('file:/data/franzoni/data/Commissioning08_RAW_v1_67838-006945C8-40A5-DD11-BD7E-001617DBD556.root')
0064             #fileNames = cms.untracked.vstring('file:/data/franzoni/data/h4b-sm6/h4b.00016428.A.0.0.root')
0065                 )
0066 
0067 
0068 process.outputmodule = cms.OutputModule("PoolOutputModule",
0069                                           outputCommands = cms.untracked.vstring(
0070             'drop *',
0071             'keep *_ecalUncalibHit*_*_*',
0072             'keep *_ecalRecHit_*_*'
0073           ),
0074   fileName = cms.untracked.string('testEcalLocalRecoA.root')
0075 )
0076 
0077 process.dumpEv = cms.EDAnalyzer("EventContentAnalyzer")
0078 
0079 process.ecalTestRecoLocal = cms.Sequence(process.ecalEBunpacker
0080                                          *process.ecalUncalibHitWeights
0081                                          *process.ecalUncalibHitFixedAlphaBetaFit
0082                                          *process.ecalUncalibHitRatio
0083                                          *process.ecalUncalibHitGlobal
0084                                          *process.ecalRecHit
0085                                          *process.outputmodule
0086                                          #*process.dumpEv
0087                                         )
0088 
0089 process.p = cms.Path(process.ecalTestRecoLocal)