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.ecalEBunpacker.syncCheck = cms.bool(False)
process.ecalEBunpacker.orderedFedList = cms.vint32(610)
process.ecalEBunpacker.silentMode = cms.untracked.bool(False)

import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi

process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone()

process.ecalUncalibHit.EBdigiCollection = 'ecalEBunpacker:ebDigis'
process.ecalUncalibHit.EEdigiCollection = 'ecalEBunpacker:eeDigis'

process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi")
process.ecalRecHit.ebDetIdToBeRecovered = cms.InputTag("","")
process.ecalRecHit.eeDetIdToBeRecovered = cms.InputTag("","")
process.ecalRecHit.eeFEToBeRecovered = cms.InputTag("","")
process.ecalRecHit.ebFEToBeRecovered = cms.InputTag("","")
process.ecalRecHit.recoverEBFE = cms.bool(False)
process.ecalRecHit.recoverEEFE = cms.bool(False)
process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB'
process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE'

process.ecalDataSequence = cms.Sequence(process.ecalEBunpacker*process.ecalUncalibHit*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.3
Committed:	Thu Jul 21 12:09:22 2011 UTC (13 years, 9 months ago) by andersj
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	HEAD
Changes since 1.2:	+21 -2 lines
Log Message:	Ecal working
#	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	andersj	1.3	process.ecalEBunpacker.syncCheck = cms.bool(False)
56			process.ecalEBunpacker.orderedFedList = cms.vint32(610)
57			process.ecalEBunpacker.silentMode = cms.untracked.bool(False)
58
59			import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi
60
61			process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone()
62
63			process.ecalUncalibHit.EBdigiCollection = 'ecalEBunpacker:ebDigis'
64			process.ecalUncalibHit.EEdigiCollection = 'ecalEBunpacker:eeDigis'
65
66			process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi")
67			process.ecalRecHit.ebDetIdToBeRecovered = cms.InputTag("","")
68			process.ecalRecHit.eeDetIdToBeRecovered = cms.InputTag("","")
69			process.ecalRecHit.eeFEToBeRecovered = cms.InputTag("","")
70			process.ecalRecHit.ebFEToBeRecovered = cms.InputTag("","")
71			process.ecalRecHit.recoverEBFE = cms.bool(False)
72			process.ecalRecHit.recoverEEFE = cms.bool(False)
73			process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB'
74			process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE'
75	andersj	1.1
76	andersj	1.3	process.ecalDataSequence = cms.Sequence(process.ecalEBunpackerprocess.ecalUncalibHitprocess.ecalRecHit)
77	andersj	1.1
78			process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
79			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hbhe_cfi")
80			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
81			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hf_cfi")
82	andersj	1.2	process.load("andersj.HcalPlotter.HcalTBObjectUnpacker_Normal_cfi")
83	andersj	1.1
84			process.horeco.firstSample=2
85			process.horeco.samplesToAdd=6
86
87			process.hfreco.firstSample=2
88			process.hfreco.samplesToAdd=9
89
90			## process.plotanal=cms.EDAnalyzer(
91			## "HcalHOTBPlotAnal",
92	andersj	1.2	## hbheRHtag = cms.untracked.InputTag("hbheprereco"),
93	andersj	1.1	## hoRHtag = cms.untracked.InputTag("horeco"),
94			## hfRHtag = cms.untracked.InputTag("hfreco"),
95			## hcalDigiTag = cms.untracked.InputTag("hcalDigis"),
96			## hcalTrigTag = cms.untracked.InputTag("tbunpack"),
97	andersj	1.2	## ebRHtag = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"),
98			## doBeamCounters = cms.untracked.bool(True),
99	andersj	1.1	## calibFC2GeV = cms.untracked.double(0.01),
100			## HistoParameters = cms.PSet(
101			## pedGeVlo = cms.double(-15),
102			## pedGeVhi = cms.double(15),
103			## pedADClo = cms.double(0),
104			## pedADChi = cms.double(49),
105			## ledGeVlo = cms.double(-5),
106			## ledGeVhi = cms.double(250),
107			## laserGeVlo = cms.double(-5),
108			## laserGeVhi = cms.double(350),
109			## otherGeVlo = cms.double(-5),
110			## otherGeVhi = cms.double(250),
111			## beamGeVlo = cms.double(-20),
112	andersj	1.2	## beamGeVhi = cms.double(500),
113	andersj	1.1	## #beamGeVhi = cms.double(80*200),
114			## timeNSlo = cms.double(50),
115			## timeNShi = cms.double(250)
116			## )
117			## )
118
119			process.dump = cms.EDAnalyzer("EventContentAnalyzer")
120
121	andersj	1.2	process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.hfreco+process.tbunpack+process.dump)