andersj/HcalPlotter/tests_scp.py

import FWCore.ParameterSet.Config as cms
import sys
import os
import subprocess

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)
)
inFileName = 'EcalHcalCombined2011_{0:08d}.0.root'.format(RUNNUMBER)
if not os.access('moe5/{0}'.format(inFileName), os.R_OK):
    cmd = ['scp','-Cp','daq@cmsmoe5:/data0/spool/{0}'.format(inFileName),
           'moe5/{0}'.format(inFileName)]
    print ' '.join(cmd)
    subprocess.call(cmd)
    

process.source = cms.Source("PoolSource",
    fileNames = cms.untracked.vstring(
    "file:moe5/{0}".format(inFileName),
    )
)

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.ecalRatioUncalibRecHit_cfi
process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalRatioUncalibRecHit_cfi.ecalRatioUncalibRecHit.clone()

## 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.5),
        pedADChi   = cms.double(20.5),
        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.plotanal)
Revision:	1.2
Committed:	Thu Aug 18 19:50:17 2011 UTC (13 years, 8 months ago) by andersj
Content type:	text/x-python
Branch:	MAIN
Changes since 1.1:	+1 -1 lines
Log Message:	new gains peds etc
#	User	Rev	Content
1	andersj	1.1	import FWCore.ParameterSet.Config as cms
2			import sys
3			import os
4			import subprocess
5
6			process = cms.Process("Demo")
7
8			process.load("FWCore.MessageService.MessageLogger_cfi")
9
10			process.MessageLogger.cerr.FwkReport.reportEvery = 500
11			process.MessageLogger.suppressWarning.extend(["ecalRecHit"])
12
13			process.load("andersj.HcalPlotter.hcal_tb11_cff")
14
15			# We are getting the runnumber from the command line now. No need to edit this file.
16			if len(sys.argv) > 2:
17			#print sys.argv
18			arg1 = sys.argv[2]
19			RUNNUMBER = int(arg1)
20			#print arg1
21			else:
22			RUNNUMBER = 0
23			if (RUNNUMBER < 1):
24			print "Enter Runnumber: "
25			rn = sys.stdin.readline()
26			RUNNUMBER = int(rn.strip())
27			print "Running on: {0:08d}".format(RUNNUMBER)
28
29			process.maxEvents = cms.untracked.PSet(
30			input = cms.untracked.int32(-1)
31			)
32			inFileName = 'EcalHcalCombined2011_{0:08d}.0.root'.format(RUNNUMBER)
33			if not os.access('moe5/{0}'.format(inFileName), os.R_OK):
34	andersj	1.2	cmd = ['scp','-Cp','daq@cmsmoe5:/data0/spool/{0}'.format(inFileName),
35	andersj	1.1	'moe5/{0}'.format(inFileName)]
36			print ' '.join(cmd)
37			subprocess.call(cmd)
38
39
40			process.source = cms.Source("PoolSource",
41			fileNames = cms.untracked.vstring(
42			"file:moe5/{0}".format(inFileName),
43			)
44			)
45
46			process.TFileService = cms.Service("TFileService",
47			fileName = cms.string('test/tb2011_{0:08d}.root'.format(RUNNUMBER))
48			)
49
50			process.load("EventFilter.EcalRawToDigi.EcalUnpackerData_cfi")
51			process.load("EventFilter.EcalRawToDigi.EcalUnpackerMapping_cfi")
52			process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
53			process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
54			process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
55			process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
56			process.load("Geometry.EcalMapping.EcalMapping_cfi")
57			process.load("Geometry.EcalMapping.EcalMappingRecord_cfi")
58			process.load("CalibCalorimetry.EcalLaserCorrection.ecalLaserCorrectionService_cfi")
59			process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
60
61			process.ecalEBunpacker.FEDs = cms.vint32(610)
62			process.ecalEBunpacker.memUnpacking = cms.bool(False)
63			process.ecalEBunpacker.srpUnpacking = cms.bool(False)
64			process.ecalEBunpacker.syncCheck = cms.bool(False)
65			process.ecalEBunpacker.orderedFedList = cms.vint32(610)
66			process.ecalEBunpacker.silentMode = cms.untracked.bool(False)
67
68			import RecoLocalCalo.EcalRecProducers.ecalRatioUncalibRecHit_cfi
69			process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalRatioUncalibRecHit_cfi.ecalRatioUncalibRecHit.clone()
70
71			## import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi
72
73			## process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone()
74
75			process.ecalUncalibHit.EBdigiCollection = 'ecalEBunpacker:ebDigis'
76			process.ecalUncalibHit.EEdigiCollection = 'ecalEBunpacker:eeDigis'
77
78			process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi")
79			process.ecalRecHit.ebDetIdToBeRecovered = cms.InputTag("","")
80			process.ecalRecHit.eeDetIdToBeRecovered = cms.InputTag("","")
81			process.ecalRecHit.eeFEToBeRecovered = cms.InputTag("","")
82			process.ecalRecHit.ebFEToBeRecovered = cms.InputTag("","")
83			process.ecalRecHit.recoverEBFE = cms.bool(False)
84			process.ecalRecHit.recoverEEFE = cms.bool(False)
85			process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB'
86			process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE'
87
88			process.ecalDataSequence = cms.Sequence(process.ecalEBunpackerprocess.ecalUncalibHitprocess.ecalRecHit)
89
90			process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
91			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hbhe_cfi")
92			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
93			process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hf_cfi")
94			process.load("andersj.HcalPlotter.HcalTBObjectUnpacker_Normal_cfi")
95
96			process.horeco.firstSample=2
97			process.horeco.samplesToAdd=6
98
99			process.hfreco.firstSample=2
100			process.hfreco.samplesToAdd=9
101
102			process.plotanal=cms.EDAnalyzer(
103			"HcalHOTBPlotAnal",
104			hbheRHtag = cms.untracked.InputTag("hbheprereco"),
105			hoRHtag = cms.untracked.InputTag("horeco"),
106			hfRHtag = cms.untracked.InputTag("hfreco"),
107			hcalDigiTag = cms.untracked.InputTag("hcalDigis"),
108			hcalTrigTag = cms.untracked.InputTag("tbunpack"),
109			ebRHtag = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"),
110			doBeamCounters = cms.untracked.bool(True),
111			calibFC2GeV = cms.untracked.double(0.01),
112			HistoParameters = cms.PSet(
113			pedGeVlo = cms.double(-15),
114			pedGeVhi = cms.double(15),
115			pedADClo = cms.double(-0.5),
116			pedADChi = cms.double(20.5),
117			ledGeVlo = cms.double(-5),
118			ledGeVhi = cms.double(250),
119			laserGeVlo = cms.double(-5),
120			laserGeVhi = cms.double(350),
121			otherGeVlo = cms.double(-5),
122			otherGeVhi = cms.double(250),
123			beamGeVlo = cms.double(-20),
124			beamGeVhi = cms.double(500),
125			#beamGeVhi = cms.double(80*200),
126			timeNSlo = cms.double(50),
127			timeNShi = cms.double(250)
128			)
129			)
130
131			##process.dump = cms.EDAnalyzer("EventContentAnalyzer")
132
133			process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.hfreco+process.tbunpack+process.plotanal)