EcalDigis/test/SRValidCRAFT.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("EcalSelectiveReadoutValid")
0004
0005 process.maxEvents = cms.untracked.PSet(
0006     input = cms.untracked.int32(1)
0007 )
0008
0009 process.source = cms.Source("PoolSource",
0010                             fileNames = cms.untracked.vstring('to_be_replaced')
0011 )
0012
0013 import FWCore.Modules.printContent_cfi
0014 process.dumpEv = FWCore.Modules.printContent_cfi.printContent.clone()
0015
0016 # MessageLogger:
0017 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0018 process.MessageLogger.cerr.INFO.limit = 10
0019 process.MessageLogger.suppressWarning = ['ecalSelectiveReadoutValidation']
0020
0021 # DQM services
0022 process.load("DQMServices.Core.DQM_cfg")
0023 process.DQM.collectorHost = ''
0024
0025 # ECAL Unpacker:
0026 process.load("EventFilter.EcalRawToDigi.EcalUnpackerMapping_cfi")
0027 process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
0028 process.ecalEBunpacker.silentMode = cms.untracked.bool(True)
0029
0030
0031 # ECAL Geometry:
0032 process.load("Geometry.EcalCommonData.EcalOnly_cfi")
0033 #  Calo geometry service model
0034 process.load("Geometry.CaloEventSetup.CaloGeometry_cff")
0035 process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
0036 process.load("Geometry.EcalMapping.EcalMapping_cfi")
0037 process.load("Geometry.EcalMapping.EcalMappingRecord_cfi")
0038
0039 #Magnetic field:
0040 #process.load("MagneticField.Engine.volumeBasedMagneticField_cfi")
0041 process.load("Configuration.StandardSequences.MagneticField_cff")
0042
0043 process.load("CalibCalorimetry.Configuration.Ecal_FakeConditions_cff")
0044 #process.load("Configuration.StandardSequences.FakeConditions_cff")
0045
0046 process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
0047
0048 #List of collection
0049 #process.load("pgras.ListCollection.ListCollection_cfi")
0050
0051 # Srp validation (analysis module):
0052 #   Defines Ecal seletive readout validation module, ecalSelectiveReadoutValidation:
0053 process.load("Validation.EcalDigis.ecalSelectiveReadoutValidation_cfi")
0054 process.ecalSelectiveReadoutValidation.EbDigiCollection = 'ecalEBunpacker:ebDigis'
0055 process.ecalSelectiveReadoutValidation.EeDigiCollection = 'ecalEBunpacker:eeDigis'
0056 process.ecalSelectiveReadoutValidation.EbSrFlagCollection = 'ecalEBunpacker:'
0057 process.ecalSelectiveReadoutValidation.EeSrFlagCollection = 'ecalEBunpacker:'
0058 process.ecalSelectiveReadoutValidation.TrigPrimCollection = 'ecalEBunpacker:EcalTriggerPrimitives'
0059
0060 process.ecalSelectiveReadoutValidation.tpInGeV = False
0061 process.ecalSelectiveReadoutValidation.ecalDccZs1stSample = 3
0062 process.ecalSelectiveReadoutValidation.dccWeights = [ -1.1865, 0.0195, 0.2900, 0.3477, 0.3008, 0.2266 ]
0063 process.ecalSelectiveReadoutValidation.histDir = ''
0064 process.ecalSelectiveReadoutValidation.histograms = [
0065     "EBEventSize",           #Barrel data volume
0066     "EBHighInterestPayload", #Barrel high interest crystal data payload
0067     "EBLowInterestPayload",  #ECAL Barrel low interest crystal data payload
0068     "EcalMidInterestTriggerTowerFlagMap", #Mid interest trigger tower flags
0069     "EEEventSize",           #Endcap data volume
0070     "EEHighInterestPayload", #Endcap high interest crystal data payload
0071     "EELowInterestPayload",  #Endcap low interest crystal data payload
0072     "EbZeroSupp1SRFlagMap",  #Trigger tower read-out with ZS threshold 1
0073     "EcalChannelOccupancy",  #ECAL crystal channel occupancy after zero suppression
0074     "EcalDccEventSize",      #ECAL DCC event fragment size
0075     "EcalDccEventSizeComputed", #ECAL DCC event fragment size
0076     "EcalEventSize",         #ECAL data volume
0077     "EcalFullReadoutSRFlagMap", #Full Read-out trigger tower
0078     "EcalHighInterestPayload",  #ECAL high interest crystal data payload
0079     "EcalHighInterestTriggerTowerFlagMap", #High interest trigger tower flags
0080     "EcalLowInterestPayload",#ECAL low interest crystal data payload
0081     "EcalLowInterestTriggerTowerFlagMap", #Low interest trigger tower flags
0082     "EcalMidInterestTriggerTowerFlagMap", #Mid interest trigger tower flags
0083     "EcalReadoutUnitForcedBitMap", #ECAL readout unit with forced bit of SR flag on
0084     "EcalTriggerPrimitiveEt",#Trigger primitive TT E_{T}
0085     "EcalTriggerPrimitiveEtMap", #Trigger primitive
0086     "EcalTriggerTowerFlag",  #Trigger primitive TT flag
0087     "hEbEMean",              #EE <E_hit>
0088     "hEbNoZsRecVsSimE",      #Crystal no-zs simulated vs reconstructed energy
0089     "hEbRecE",               #Crystal reconstructed energy
0090     "hEbRecEHitXtal",        #EB rec energy of hit crystals
0091     "hEeEMean",              #EE <E_hit>
0092     "hEeRecE",               #EE crystal reconstructed energy
0093     "tpVsEtSum",             #Trigger primitive Et (TP) vs #sumE_{T}
0094     "ttfVsEtSum",            #Trigger tower flag vs #sumE_{T}
0095     "ttfVsTp",               #Trigger tower flag vs TP
0096     "zsEbHiFIRemu",          #Emulated ouput of ZS FIR filter for EB high interest crystals
0097     "zsEbLiFIRemu",          #Emulated ouput of ZS FIR filter for EB low interest crystals
0098     "zsEeHiFIRemu",          #Emulated ouput of ZS FIR filter for EE high interest crystals
0099     "zsEeLiFIRemu"           #Emulated ouput of ZS FIR filter for EE low interest crystals
0100     ]
0101
0102
0103 process.p = cms.Path(process.ecalSelectiveReadoutValidation)
0104
0105 process.ecalSelectiveReadoutValidation.outputFile = 'run69912hists_.root'
0106
0107 process.source.fileNames = ['/store/data/Commissioning08/Cosmics/RAW/v1/000/069/912/049C2F4D-10AD-DD11-BFEA-000423D174FE.root' ]