ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/andersj/HcalPlotter/tests.py
Revision: 1.4
Committed: Fri Oct 14 17:51:23 2011 UTC (13 years, 6 months ago) by andersj
Content type: text/x-python
Branch: MAIN
Changes since 1.3: +11 -6 lines
Log Message:
latest code

File Contents

# Content
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 process.MessageLogger.suppressWarning.extend(["ecalRecHit"])
10
11 process.load("andersj.HcalPlotter.hcal_tb11_cff")
12
13 # We are getting the runnumber from the command line now. No need to edit this file.
14 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 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 "file:moe5/EcalHcalCombined2011_{0:08d}.0.root".format(RUNNUMBER),
34 )
35 )
36
37 process.TFileService = cms.Service("TFileService",
38 fileName = cms.string('test/tb2011_{0:08d}.root'.format(RUNNUMBER))
39 )
40
41 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 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.ecalRatioUncalibRecHit_cfi
60 process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalRatioUncalibRecHit_cfi.ecalRatioUncalibRecHit.clone()
61
62 ## import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi
63
64 ## process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone()
65
66 process.ecalUncalibHit.EBdigiCollection = 'ecalEBunpacker:ebDigis'
67 process.ecalUncalibHit.EEdigiCollection = 'ecalEBunpacker:eeDigis'
68
69 process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi")
70 process.ecalRecHit.ebDetIdToBeRecovered = cms.InputTag("","")
71 process.ecalRecHit.eeDetIdToBeRecovered = cms.InputTag("","")
72 process.ecalRecHit.eeFEToBeRecovered = cms.InputTag("","")
73 process.ecalRecHit.ebFEToBeRecovered = cms.InputTag("","")
74 process.ecalRecHit.recoverEBFE = cms.bool(False)
75 process.ecalRecHit.recoverEEFE = cms.bool(False)
76 process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB'
77 process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE'
78
79 process.ecalDataSequence = cms.Sequence(process.ecalEBunpacker*process.ecalUncalibHit*process.ecalRecHit)
80
81 process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
82 process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hbhe_cfi")
83 process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
84 #process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hf_cfi")
85 process.load("andersj.HcalPlotter.HcalTBObjectUnpacker_Normal_cfi")
86
87 ## process.horeco.firstSample=2
88 ## process.horeco.samplesToAdd=6
89 process.horeco.firstSample=3
90 process.horeco.samplesToAdd=4
91 process.horeco.correctForPhaseContainment = False
92 process.horeco.correctForTimeslew = False
93
94 #process.hfreco.firstSample=2
95 #process.hfreco.samplesToAdd=9
96
97 process.plotanal=cms.EDAnalyzer(
98 "HcalHOTBPlotAnal",
99 hbheRHtag = cms.untracked.InputTag("hbheprereco"),
100 hoRHtag = cms.untracked.InputTag("horeco"),
101 hfRHtag = cms.untracked.InputTag("hfreco"),
102 hcalDigiTag = cms.untracked.InputTag("hcalDigis"),
103 hcalTrigTag = cms.untracked.InputTag("tbunpack"),
104 ebRHtag = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"),
105 doBeamCounters = cms.untracked.bool(True),
106 calibFC2GeV = cms.untracked.double(0.01),
107 HistoParameters = cms.PSet(
108 pedGeVlo = cms.double(-15),
109 pedGeVhi = cms.double(15),
110 pedADClo = cms.double(0),
111 pedADChi = cms.double(49),
112 ledGeVlo = cms.double(-5),
113 ledGeVhi = cms.double(250),
114 laserGeVlo = cms.double(-5),
115 laserGeVhi = cms.double(350),
116 otherGeVlo = cms.double(-5),
117 otherGeVhi = cms.double(250),
118 beamGeVlo = cms.double(-20),
119 beamGeVhi = cms.double(500),
120 #beamGeVhi = cms.double(80*200),
121 timeNSlo = cms.double(50),
122 timeNShi = cms.double(250)
123 )
124 )
125
126 ##process.dump = cms.EDAnalyzer("EventContentAnalyzer")
127
128 #process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.hfreco+process.tbunpack+process.plotanal)
129 process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.tbunpack+process.plotanal)