EcalDigis/test/SRValidCRAFT_SR_TP_reprocess.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
0047 process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
0048
0049 #List of collection
0050 #process.load("pgras.ListCollection.ListCollection_cfi")
0051
0052
0053 # SR emulation:
0054 process.load("SimCalorimetry.EcalSelectiveReadoutProducers.ecalDigis_craft_cfi")
0055 process.simEcalDigis.trigPrimProducer = "ecalTriggerPrimitiveDigis"
0056 process.simEcalDigis.digiProducer = cms.string('ecalEBunpacker')
0057 process.simEcalDigis.EBdigiCollection = cms.string('ebDigis')
0058 process.simEcalDigis.EEdigiCollection = cms.string('eeDigis')
0059
0060 # Srp validation (analysis module):
0061 #   Defines Ecal seletive readout validation module, ecalSelectiveReadoutValidation:
0062 process.load("Validation.EcalDigis.ecalSelectiveReadoutValidation_cfi")
0063 process.ecalSelectiveReadoutValidation.EbDigiCollection = 'simEcalDigis:ebDigis'
0064 process.ecalSelectiveReadoutValidation.EeDigiCollection = 'simEcalDigis:eeDigis'
0065 process.ecalSelectiveReadoutValidation.EbSrFlagCollection = 'simEcalDigis:ebSrFlags'
0066 process.ecalSelectiveReadoutValidation.EeSrFlagCollection = 'simEcalDigis:eeSrFlags'
0067 process.ecalSelectiveReadoutValidation.TrigPrimCollection = 'simEcalTriggerPrimitiveDigis'
0068
0069 #process.ecalSelectiveReadoutValidation.EbDigiCollection = 'ecalEBunpacker:ebDigis'
0070 #process.ecalSelectiveReadoutValidation.EeDigiCollection = 'ecalEBunpacker:eeDigis'
0071 #process.ecalSelectiveReadoutValidation.EbSrFlagCollection = 'ecalEBunpacker:'
0072 #process.ecalSelectiveReadoutValidation.EeSrFlagCollection = 'ecalEBunpacker:'
0073 #process.ecalSelectiveReadoutValidation.TrigPrimCollection = 'ecalEBunpacker:EcalTriggerPrimitives'
0074
0075 process.ecalSelectiveReadoutValidation.tpInGeV = False
0076 process.ecalSelectiveReadoutValidation.ecalDccZs1stSample = 3
0077 process.ecalSelectiveReadoutValidation.dccWeights = [ -1.1865, 0.0195, 0.2900, 0.3477, 0.3008, 0.2266 ]
0078 process.ecalSelectiveReadoutValidation.histDir = ''
0079 process.ecalSelectiveReadoutValidation.histograms = [
0080     "EBEventSize",           #Barrel data volume
0081     "EBHighInterestPayload", #Barrel high interest crystal data payload
0082     "EBLowInterestPayload",  #ECAL Barrel low interest crystal data payload
0083     "EcalMidInterestTriggerTowerFlagMap", #Mid interest trigger tower flags
0084     "EEEventSize",           #Endcap data volume
0085     "EEHighInterestPayload", #Endcap high interest crystal data payload
0086     "EELowInterestPayload",  #Endcap low interest crystal data payload
0087     "EbZeroSupp1SRFlagMap",  #Trigger tower read-out with ZS threshold 1
0088     "EcalChannelOccupancy",  #ECAL crystal channel occupancy after zero suppression
0089     "EcalDccEventSize",      #ECAL DCC event fragment size
0090     "EcalDccEventSizeComputed", #ECAL DCC event fragment size
0091     "EcalEventSize",         #ECAL data volume
0092     "EcalFullReadoutSRFlagMap", #Full Read-out trigger tower
0093     "EcalHighInterestPayload",  #ECAL high interest crystal data payload
0094     "EcalHighInterestTriggerTowerFlagMap", #High interest trigger tower flags
0095     "EcalLowInterestPayload",#ECAL low interest crystal data payload
0096     "EcalLowInterestTriggerTowerFlagMap", #Low interest trigger tower flags
0097     "EcalMidInterestTriggerTowerFlagMap", #Mid interest trigger tower flags
0098     "EcalReadoutUnitForcedBitMap", #ECAL readout unit with forced bit of SR flag on
0099     "EcalTriggerPrimitiveEt",#Trigger primitive TT E_{T}
0100     "EcalTriggerPrimitiveEtMap", #Trigger primitive
0101     "EcalTriggerTowerFlag",  #Trigger primitive TT flag
0102     "hEbEMean",              #EE <E_hit>
0103     "hEbNoZsRecVsSimE",      #Crystal no-zs simulated vs reconstructed energy
0104     "hEbRecE",               #Crystal reconstructed energy
0105     "hEbRecEHitXtal",        #EB rec energy of hit crystals
0106     "hEeEMean",              #EE <E_hit>
0107     "hEeRecE",               #EE crystal reconstructed energy
0108     "tpVsEtSum",             #Trigger primitive Et (TP) vs #sumE_{T}
0109     "ttfVsEtSum",            #Trigger tower flag vs #sumE_{T}
0110     "ttfVsTp",               #Trigger tower flag vs TP
0111     "zsEbHiFIRemu",          #Emulated ouput of ZS FIR filter for EB high interest crystals
0112     "zsEbLiFIRemu",          #Emulated ouput of ZS FIR filter for EB low interest crystals
0113     "zsEeHiFIRemu",          #Emulated ouput of ZS FIR filter for EE high interest crystals
0114     "zsEeLiFIRemu"           #Emulated ouput of ZS FIR filter for EE low interest crystals
0115     ]
0116 process.tpparams = cms.ESSource("EmptyESSource",
0117     recordName = cms.string('EcalTPGLinearizationConstRcd'),
0118     iovIsRunNotTime = cms.bool(True),
0119     firstValid = cms.vuint32(1)
0120 )
0121
0122 process.tpparams2 = cms.ESSource("EmptyESSource",
0123     recordName = cms.string('EcalTPGPedestalsRcd'),
0124     iovIsRunNotTime = cms.bool(True),
0125     firstValid = cms.vuint32(1)
0126 )
0127
0128 process.tpparams3 = cms.ESSource("EmptyESSource",
0129     recordName = cms.string('EcalTPGSlidingWindowRcd'),
0130     iovIsRunNotTime = cms.bool(True),
0131     firstValid = cms.vuint32(1)
0132 )
0133
0134 process.tpparams4 = cms.ESSource("EmptyESSource",
0135     recordName = cms.string('EcalTPGWeightIdMapRcd'),
0136     iovIsRunNotTime = cms.bool(True),
0137     firstValid = cms.vuint32(1)
0138 )
0139
0140 process.tpparams5 = cms.ESSource("EmptyESSource",
0141     recordName = cms.string('EcalTPGWeightGroupRcd'),
0142     iovIsRunNotTime = cms.bool(True),
0143     firstValid = cms.vuint32(1)
0144 )
0145
0146 process.tpparams6 = cms.ESSource("EmptyESSource",
0147     recordName = cms.string('EcalTPGLutGroupRcd'),
0148     iovIsRunNotTime = cms.bool(True),
0149     firstValid = cms.vuint32(1)
0150 )
0151
0152 process.tpparams7 = cms.ESSource("EmptyESSource",
0153     recordName = cms.string('EcalTPGLutIdMapRcd'),
0154     iovIsRunNotTime = cms.bool(True),
0155     firstValid = cms.vuint32(1)
0156 )
0157
0158 process.tpparams8 = cms.ESSource("EmptyESSource",
0159     recordName = cms.string('EcalTPGFineGrainEBIdMapRcd'),
0160     iovIsRunNotTime = cms.bool(True),
0161     firstValid = cms.vuint32(1)
0162 )
0163
0164 process.tpparams9 = cms.ESSource("EmptyESSource",
0165     recordName = cms.string('EcalTPGFineGrainEBGroupRcd'),
0166     iovIsRunNotTime = cms.bool(True),
0167     firstValid = cms.vuint32(1)
0168 )
0169
0170 process.tpparams10 = cms.ESSource("EmptyESSource",
0171     recordName = cms.string('EcalTPGFineGrainStripEERcd'),
0172     iovIsRunNotTime = cms.bool(True),
0173     firstValid = cms.vuint32(1)
0174 )
0175
0176 process.tpparams11 = cms.ESSource("EmptyESSource",
0177     recordName = cms.string('EcalTPGFineGrainTowerEERcd'),
0178     iovIsRunNotTime = cms.bool(True),
0179     firstValid = cms.vuint32(1)
0180 )
0181
0182 process.tpparams12 = cms.ESSource("EmptyESSource",
0183     recordName = cms.string('EcalTPGPhysicsConstRcd'),
0184     iovIsRunNotTime = cms.bool(True),
0185     firstValid = cms.vuint32(1)
0186 )
0187
0188 process.EcalTrigPrimESProducer = cms.ESProducer("EcalTrigPrimESProducer",
0189     DatabaseFile = cms.untracked.string('TPG_startup.txt.gz')
0190 )
0191
0192 process.ecalTriggerPrimitiveDigis = cms.EDProducer("EcalTrigPrimProducer",
0193     InstanceEB = cms.string('ebDigis'),
0194     InstanceEE = cms.string('eeDigis'),
0195     Label = cms.string('ecalEBunpacker'),
0196
0197     BarrelOnly = cms.bool(False),
0198     Famos = cms.bool(False),
0199     TcpOutput = cms.bool(False),
0200     Debug = cms.bool(False),
0201
0202     binOfMaximum = cms.int32(6), ## optional from release 200 on, from 1-10
0203
0204     TTFHighEnergyEB = cms.double(1.0),
0205     TTFHighEnergyEE = cms.double(1.0),
0206     TTFLowEnergyEB = cms.double(1.0), ## this + the following is added from 140_pre4 on
0207     TTFLowEnergyEE = cms.double(1.0)
0208 )
0209
0210
0211 #process.tpAnalyzer = cms.EDAnalyzer("EcalTPGAnalyzer",
0212 #
0213 #    Label = cms.string('ecalEBunpacker'),
0214 #    Producer = cms.string('EcalTriggerPrimitives'),
0215 #    DigiLabel = cms.string('ecalEBunpacker'),
0216 #    DigiProducerEE = cms.string(''),
0217 #    DigiProducerEB = cms.string('ebDigis'),
0218 #    EmulLabel = cms.string('ecalTriggerPrimitiveDigis'),
0219 #    EmulProducer = cms.string(''),
0220 #
0221 #    Print = cms.bool(True),
0222 #    CompareWithPreviousEvent = cms.bool(True),
0223 #    ReadTriggerPrimitives = cms.bool(True),
0224 #    UseEndCap = cms.bool(False),
0225 #    KeepOnlyL1Ecal = cms.bool(True) ,
0226 #
0227 #    ADCCut = cms.int32(3), ## crystal energy contributes to tower energy if above the threshold
0228 #    shapeCut = cms.int32(3), ## pulse profile filled with tower energy if above the threshold
0229 #    occupancyCut = cms.int32(3), ## occupancy plots filled if tpg above the threshold
0230 #
0231 #    TPGEmulatorIndexRef = cms.int32(2) ## must be in [0,4]
0232 #)
0233 #
0234
0235 process.p = cms.Path(process.ecalEBunpacker*process.ecalTriggerPrimitiveDigis*process.simEcalDigis*process.ecalSelectiveReadoutValidation)
0236
0237 process.ecalSelectiveReadoutValidation.outputFile = 'run69912hists_.root'
0238
0239 process.source.fileNames = ['/store/data/Commissioning08/Cosmics/RAW/v1/000/069/912/049C2F4D-10AD-DD11-BFEA-000423D174FE.root' ]