EcalTools/python/ecalCosmicsHists_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("COSMICSANALYSIS")
0004 process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
0005
0006 process.load("Geometry.CaloEventSetup.CaloGeometry_cff")
0007
0008 #process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0009 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_noesprefer_cff")
0010 process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
0011 process.GlobalTag.globaltag = 'CRUZET4_V5P::All'
0012
0013 process.EcalTrivialConditionRetriever.producedEcalWeights = False
0014 process.EcalTrivialConditionRetriever.producedEcalPedestals = False
0015 process.EcalTrivialConditionRetriever.producedEcalIntercalibConstants = False
0016 process.EcalTrivialConditionRetriever.producedEcalIntercalibErrors = False
0017 process.EcalTrivialConditionRetriever.producedEcalGainRatios = False
0018 process.EcalTrivialConditionRetriever.producedEcalADCToGeVConstant = False
0019 #process.EcalTrivialConditionRetriever.producedEcalADCToGeVConstant = True
0020 #process.EcalTrivialConditionRetriever.adcToGeVEBConstant = 0.035
0021 #process.EcalTrivialConditionRetriever.adcToGeVEEConstant = 0.060
0022 process.EcalTrivialConditionRetriever.producedEcalLaserCorrection = False
0023 #process.EcalTrivialConditionRetriever.producedEcalChannelStatus = False
0024 process.EcalTrivialConditionRetriever.producedEcalChannelStatus = cms.untracked.bool(False)
0025 #process.EcalTrivialConditionRetriever.channelStatusFile = 'CaloOnlineTools/EcalTools/data/listCRUZET4.v5.hashed.txt'
0026 #process.es_prefer_EcalTrivialConditionRetriever = cms.ESPrefer("EcalTrivialConditionRetriever")
0027
0028 process.load("CalibCalorimetry.EcalLaserCorrection.ecalLaserCorrectionService_cfi")
0029
0030 import CalibTracker.Configuration.Common.PoolDBESSource_cfi
0031 process.siStripPedestalFrontier = CalibTracker.Configuration.Common.PoolDBESSource_cfi.poolDBESSource.clone()
0032 process.siStripPedestalFrontier.connect = 'frontier://PromptProd/CMS_COND_21X_STRIP'
0033 process.siStripPedestalFrontier.toGet = cms.VPSet(cms.PSet(
0034             record = cms.string('SiStripPedestalsRcd'),
0035                         tag = cms.string('SiStripPedestals_TKCC_21X_v3_hlt')
0036                     ))
0037 process.siStripPedestalFrontier.BlobStreamerName = 'TBufferBlobStreamingService'
0038 process.es_prefer_SiStripFake = cms.ESPrefer("PoolDBESSource","siStripPedestalFrontier")
0039
0040 process.load("EventFilter.EcalRawToDigi.EcalUnpackerMapping_cfi")
0041 process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
0042
0043 import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi
0044 process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone()
0045 process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi")
0046
0047 process.load("RecoEcal.EgammaClusterProducers.geometryForClustering_cff")
0048 process.load("RecoEcal.EgammaClusterProducers.cosmicClusteringSequence_cff")
0049
0050 process.load("L1TriggerConfig.L1ScalesProducers.L1MuTriggerScalesConfig_cff")
0051 process.load("L1TriggerConfig.L1ScalesProducers.L1MuTriggerPtScaleConfig_cff")
0052 process.load("L1TriggerConfig.L1GtConfigProducers.L1GtBoardMapsConfig_cff")
0053 process.load("L1TriggerConfig.L1GtConfigProducers.L1GtConfig_cff")
0054 process.load("L1TriggerConfig.L1GtConfigProducers.Luminosity.startup.L1Menu_startup_v3_Unprescaled_cff")
0055 import FWCore.Modules.printContent_cfi
0056 process.dumpEv = FWCore.Modules.printContent_cfi.printContent.clone()
0057
0058 import EventFilter.L1GlobalTriggerRawToDigi.l1GtUnpack_cfi
0059 process.gtDigis = EventFilter.L1GlobalTriggerRawToDigi.l1GtUnpack_cfi.l1GtUnpack.clone()
0060 process.load("CaloOnlineTools.EcalTools.ecalCosmicsHists_cfi")
0061
0062 process.load("Configuration.StandardSequences.MagneticField_cff")
0063
0064 process.load("Geometry.EcalMapping.EcalMapping_cfi")
0065 process.load("Geometry.EcalMapping.EcalMappingRecord_cfi")
0066 process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
0067 process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
0068 process.load("Geometry.CommonTopologies.globalTrackingGeometry_cfi")
0069 process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi")
0070
0071 process.load("TrackingTools.TrackAssociator.DetIdAssociatorESProducer_cff")
0072 process.load("TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorAny_cfi")
0073 process.load("TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorAlong_cfi")
0074 process.load("TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorOpposite_cfi")
0075 process.load("Configuration.GlobalRuns.ForceZeroTeslaField_cff")
0076
0077 process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
0078 process.load("RecoLocalCalo.Configuration.hcalLocalReco_cff")
0079
0080 process.load("HLTrigger.special.TriggerTypeFilter_cfi")
0081
0082 process.source = cms.Source("PoolSource",
0083     skipEvents = cms.untracked.uint32(0),
0084     fileNames = (cms.untracked.vstring(
0085         '/store/data/Commissioning08/BeamHalo/RAW/StuffAlmostToP5_v1/000/061/642/94D3CADF-A47D-DD11-BF3E-000423D94E1C.root'
0086 ))
0087 )
0088
0089 process.maxEvents = cms.untracked.PSet(
0090     input = cms.untracked.int32(-1)
0091 )
0092
0093 process.MessageLogger = cms.Service("MessageLogger",
0094     suppressWarning = cms.untracked.vstring('ecalEBunpacker','ecalUncalibHit','ecalRecHit',
0095                                             'hcalDigis','hcalLocalRecoSequence','ecalCosmicsHists'),
0096     suppressInfo = cms.untracked.vstring('ecalEBunpacker', 'ecalUncalibHit',
0097                                          'ecalRecHit', 'hcalDigis','hcalLocalRecoSequence'),
0098     cout = cms.untracked.PSet(
0099         threshold = cms.untracked.string('ERROR')
0100     ),
0101     categories = cms.untracked.vstring('EcalCosmicsHists','CosmicClusterAlgo','CosmicClusterProducer'),
0102     destinations = cms.untracked.vstring('cout')
0103 )
0104
0105 process.p = cms.Path(process.triggerTypeFilter*process.ecalEBunpacker*process.ecalUncalibHit*process.ecalRecHit*process.cosmicClusteringSequence*process.gtDigis*process.ecalCosmicsHists)
0106
0107 process.ecalUncalibHit.EBdigiCollection = 'ecalEBunpacker:ebDigis'
0108 process.ecalUncalibHit.EEdigiCollection = 'ecalEBunpacker:eeDigis'
0109 process.ecalRecHit.ChannelStatusToBeExcluded = [1]
0110 process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB'
0111 process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE'
0112 process.cosmicBasicClusters.barrelUnHitProducer = "ecalUncalibHit"
0113 process.cosmicBasicClusters.endcapUnHitProducer = "ecalUncalibHit"
0114
0115 process.gtDigis.DaqGtInputTag = 'source'
0116
0117 process.ecalCosmicsHists.fileName = 'EcalCosmicsHists'
0118 process.ecalCosmicsHists.runInFileName = True
0119 #process.ecalCosmicsHists.MinTimingAmpEB = 0.1,   # for adcToGeV=0.009, gain 200
0120 #process.ecalCosmicsHists.MinRecHitAmpEB = 0.027, # for adcToGeV=0.009, gain 200
0121 process.ecalCosmicsHists.MinTimingAmpEB = 0.35   # for adcToGeV=0.035, gain 50
0122 process.ecalCosmicsHists.MinRecHitAmpEB = 0.070  # for adcToGeV=0.035, gain 50
0123 process.ecalCosmicsHists.MinTimingAmpEE = 0.9    # for adcToGeV=0.06
0124 process.ecalCosmicsHists.MinRecHitAmpEE = 0.180  # for adcToGeV=0.06
0125
0126 process.hbhereco = process.hbheprereco.clone()
0127 process.hbhereco.firstSample = 1
0128 process.hbhereco.samplesToAdd = 8
0129 process.hbhereco.correctForTimeslew = True
0130 process.hbhereco.correctForPhaseContainment = True
0131 process.hbhereco.correctionPhaseNS = 10.0
0132 process.horeco.firstSample = 1
0133 process.horeco.samplesToAdd = 8
0134 process.horeco.correctForTimeslew = True
0135 process.horeco.correctForPhaseContainment = True
0136 process.horeco.correctionPhaseNS = 10.
0137 process.hfreco.firstSample = 1
0138 process.hfreco.samplesToAdd = 8
0139 process.hfreco.correctForTimeslew = True
0140 process.hfreco.correctForPhaseContainment = True
0141 process.hfreco.correctionPhaseNS = 10.
0142
0143 process.triggerTypeFilter.SelectedTriggerType = 1
0144