andersj/HcalPlotter/dump.py

import FWCore.ParameterSet.Config as cms
import sys

process = cms.Process("Demo")

process.load("FWCore.MessageService.MessageLogger_cfi")

process.MessageLogger.cerr.FwkReport.reportEvery = 500
process.MessageLogger.suppressWarning.extend(["ecalRecHit"])

process.load("andersj.HcalPlotter.hcal_tb11_cff")

# We are getting the runnumber from the command line now. No need to edit this file.
if len(sys.argv) > 2:
    #print sys.argv
    arg1 = sys.argv[2]
    RUNNUMBER = int(arg1)
    #print arg1
else:
    RUNNUMBER = 0
if (RUNNUMBER < 1):
    print "Enter Runnumber: "
    rn = sys.stdin.readline()
    RUNNUMBER = int(rn.strip())
print "Running on: {0:08d}".format(RUNNUMBER)

process.maxEvents = cms.untracked.PSet(
    input = cms.untracked.int32(1)
)

process.source = cms.Source("PoolSource",
    fileNames = cms.untracked.vstring(
    "file:moe5/EcalHcalCombined2011_{0:08d}.0.root".format(RUNNUMBER),
    )
)

# process.TFileService = cms.Service("TFileService",
#     fileName = cms.string('test/tb2011_{0:08d}.root'.format(RUNNUMBER))
# )

process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
process.load("EventFilter.EcalRawToDigi.EcalUnpackerMapping_cfi")
process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
process.load("Geometry.EcalMapping.EcalMapping_cfi")
process.load("Geometry.EcalMapping.EcalMappingRecord_cfi")
process.load("CalibCalorimetry.EcalLaserCorrection.ecalLaserCorrectionService_cfi")
process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")

process.ecalEBunpacker.FEDs = cms.vint32(610)
process.ecalEBunpacker.memUnpacking = cms.bool(False)
process.ecalEBunpacker.srpUnpacking = cms.bool(False)

process.load("andersj.HcalPlotter.ecalLocalRecoSequence_cff")
process.ecalDataSequence = cms.Sequence(process.ecalEBunpacker*process.ecalGlobalUncalibRecHit*process.ecalRecHit)

process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hbhe_cfi")
process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hf_cfi")
process.load("andersj.HcalPlotter.HcalTBObjectUnpacker_Normal_cfi")

process.horeco.firstSample=2
process.horeco.samplesToAdd=6

process.hfreco.firstSample=2
process.hfreco.samplesToAdd=9

## process.plotanal=cms.EDAnalyzer(
##     "HcalHOTBPlotAnal",
##     hbheRHtag = cms.untracked.InputTag("hbheprereco"),
##     hoRHtag   = cms.untracked.InputTag("horeco"),
##     hfRHtag   = cms.untracked.InputTag("hfreco"),
##     hcalDigiTag = cms.untracked.InputTag("hcalDigis"),
##     hcalTrigTag = cms.untracked.InputTag("tbunpack"),
##     ebRHtag = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"),
##     doBeamCounters = cms.untracked.bool(True),
##     calibFC2GeV = cms.untracked.double(0.01),
##     HistoParameters = cms.PSet(
##         pedGeVlo   = cms.double(-15),
##         pedGeVhi   = cms.double(15),
##         pedADClo   = cms.double(0),
##         pedADChi   = cms.double(49),
##         ledGeVlo   = cms.double(-5),
##         ledGeVhi   = cms.double(250),
##         laserGeVlo = cms.double(-5),
##         laserGeVhi = cms.double(350),
##         otherGeVlo = cms.double(-5),
##         otherGeVhi = cms.double(250),
##         beamGeVlo  = cms.double(-20),
##         beamGeVhi  = cms.double(500),
##         #beamGeVhi  = cms.double(80*200),
##         timeNSlo   = cms.double(50),
##         timeNShi   = cms.double(250)
##         )
##     )

process.dump = cms.EDAnalyzer("EventContentAnalyzer")

process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.hfreco+process.tbunpack+process.dump)
Revision:	1.2
Committed:	Wed Jul 20 12:19:05 2011 UTC (13 years, 9 months ago) by andersj
Content type:	text/x-python
Branch:	MAIN
Changes since 1.1:	+28 -57 lines
Log Message:	new dump.py
#	User	Rev	Content
1	andersj	1.1	import FWCore.ParameterSet.Config as cms
2			import sys
3
4			process = cms.Process("Demo")
5
6			process.load("FWCore.MessageService.MessageLogger_cfi")
7
8			process.MessageLogger.cerr.FwkReport.reportEvery = 500
9	andersj	1.2	process.MessageLogger.suppressWarning.extend(["ecalRecHit"])
10	andersj	1.1
11	andersj	1.2	process.load("andersj.HcalPlotter.hcal_tb11_cff")
12	andersj	1.1
13			# We are getting the runnumber from the command line now. No need to edit this file.
14	andersj	1.2	if len(sys.argv) > 2:
15			#print sys.argv
16			arg1 = sys.argv[2]
17			RUNNUMBER = int(arg1)
18			#print arg1
19			else:
20			RUNNUMBER = 0
21			if (RUNNUMBER < 1):
22			print "Enter Runnumber: "
23			rn = sys.stdin.readline()
24			RUNNUMBER = int(rn.strip())
25			print "Running on: {0:08d}".format(RUNNUMBER)
26
27	andersj	1.1	process.maxEvents = cms.untracked.PSet(
28			input = cms.untracked.int32(1)
29			)
30
31			process.source = cms.Source("PoolSource",
32			fileNames = cms.untracked.vstring(
33	andersj	1.2	"file:moe5/EcalHcalCombined2011_{0:08d}.0.root".format(RUNNUMBER),
34	andersj	1.1	)
35			)
36
37	andersj	1.2	# process.TFileService = cms.Service("TFileService",
38			# fileName = cms.string('test/tb2011_{0:08d}.root'.format(RUNNUMBER))
39			# )
40
41	andersj	1.1	process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
42			process.load("EventFilter.EcalRawToDigi.EcalUnpackerMapping_cfi")
43			process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
44			process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
45			process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
46			process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
47			process.load("Geometry.EcalMapping.EcalMapping_cfi")
48			process.load("Geometry.EcalMapping.EcalMappingRecord_cfi")
49			process.load("CalibCalorimetry.EcalLaserCorrection.ecalLaserCorrectionService_cfi")
50			process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
51
52			process.ecalEBunpacker.FEDs = cms.vint32(610)
53			process.ecalEBunpacker.memUnpacking = cms.bool(False)
54			process.ecalEBunpacker.srpUnpacking = cms.bool(False)
55
56	andersj	1.2	process.load("andersj.HcalPlotter.ecalLocalRecoSequence_cff")
57			process.ecalDataSequence = cms.Sequence(process.ecalEBunpackerprocess.ecalGlobalUncalibRecHitprocess.ecalRecHit)
58	andersj	1.1
59			process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
60			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hbhe_cfi")
61			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
62			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hf_cfi")
63	andersj	1.2	process.load("andersj.HcalPlotter.HcalTBObjectUnpacker_Normal_cfi")
64	andersj	1.1
65			process.horeco.firstSample=2
66			process.horeco.samplesToAdd=6
67
68			process.hfreco.firstSample=2
69			process.hfreco.samplesToAdd=9
70
71			## process.plotanal=cms.EDAnalyzer(
72			## "HcalHOTBPlotAnal",
73	andersj	1.2	## hbheRHtag = cms.untracked.InputTag("hbheprereco"),
74	andersj	1.1	## hoRHtag = cms.untracked.InputTag("horeco"),
75			## hfRHtag = cms.untracked.InputTag("hfreco"),
76			## hcalDigiTag = cms.untracked.InputTag("hcalDigis"),
77			## hcalTrigTag = cms.untracked.InputTag("tbunpack"),
78	andersj	1.2	## ebRHtag = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"),
79			## doBeamCounters = cms.untracked.bool(True),
80	andersj	1.1	## calibFC2GeV = cms.untracked.double(0.01),
81			## HistoParameters = cms.PSet(
82			## pedGeVlo = cms.double(-15),
83			## pedGeVhi = cms.double(15),
84			## pedADClo = cms.double(0),
85			## pedADChi = cms.double(49),
86			## ledGeVlo = cms.double(-5),
87			## ledGeVhi = cms.double(250),
88			## laserGeVlo = cms.double(-5),
89			## laserGeVhi = cms.double(350),
90			## otherGeVlo = cms.double(-5),
91			## otherGeVhi = cms.double(250),
92			## beamGeVlo = cms.double(-20),
93	andersj	1.2	## beamGeVhi = cms.double(500),
94	andersj	1.1	## #beamGeVhi = cms.double(80*200),
95			## timeNSlo = cms.double(50),
96			## timeNShi = cms.double(250)
97			## )
98			## )
99
100			process.dump = cms.EDAnalyzer("EventContentAnalyzer")
101
102	andersj	1.2	process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.hfreco+process.tbunpack+process.dump)