| 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_tb10_cff") |
| 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 |
< |
#print "Enter Runnumber: " |
| 15 |
< |
#rn = sys.stdin.readline() |
| 16 |
< |
RUNNUMBER = 2264 |
| 17 |
< |
print "Running on: " + str(RUNNUMBER) |
| 18 |
< |
|
| 19 |
< |
# Use TB source for HTB data fromat |
| 20 |
< |
#process.source = cms.Source("HcalTBSource", |
| 21 |
< |
# fileNames = cms.untracked.vstring( |
| 22 |
< |
#"file:/data/spool/HTB_00" + str(RUNNUMBER) + ".root" |
| 23 |
< |
#) |
| 24 |
< |
#) |
| 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:/tmp_mnt/data0/spool/EcalHcalCombined2010_0000" + str(RUNNUMBER) + ".0.root", |
| 33 |
> |
"file:moe5/EcalHcalCombined2011_{0:08d}.0.root".format(RUNNUMBER), |
| 34 |
|
) |
| 35 |
|
) |
| 36 |
|
|
| 37 |
< |
process.TFileService = cms.Service("TFileService", |
| 38 |
< |
fileName = cms.string('test/tb2_0000' + str(RUNNUMBER) +'.root') |
| 39 |
< |
) |
| 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") |
| 53 |
|
process.ecalEBunpacker.memUnpacking = cms.bool(False) |
| 54 |
|
process.ecalEBunpacker.srpUnpacking = cms.bool(False) |
| 55 |
|
|
| 56 |
< |
process.load("ecalLocalRecoSequence_cff") |
| 57 |
< |
## #uncalibrated rechits |
| 54 |
< |
## import RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi |
| 55 |
< |
## process.ecalUncalibHit = RecoLocalCalo.EcalRecProducers.ecalFixedAlphaBetaFitUncalibRecHit_cfi.ecalFixedAlphaBetaFitUncalibRecHit.clone() |
| 56 |
< |
|
| 57 |
< |
## #rechits |
| 58 |
< |
## process.load("RecoLocalCalo.EcalRecProducers.ecalRecHit_cfi") |
| 59 |
< |
|
| 60 |
< |
## process.ecalUncalibHit.MinAmplBarrel = 12. |
| 61 |
< |
## process.ecalUncalibHit.MinAmplEndcap = 16. |
| 62 |
< |
## process.ecalUncalibHit.EBdigiCollection = 'ecalEBunpacker:ebDigis' |
| 63 |
< |
## process.ecalUncalibHit.EEdigiCollection = 'ecalEBunpacker:eeDigis' |
| 64 |
< |
|
| 65 |
< |
## process.ecalRecHit.EBuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEB' |
| 66 |
< |
## process.ecalRecHit.EEuncalibRecHitCollection = 'ecalUncalibHit:EcalUncalibRecHitsEE' |
| 67 |
< |
|
| 68 |
< |
## process.ecalDataSequence = cms.Sequence(process.ecalEBunpacker*process.ecalUncalibHit*process.ecalRecHit) |
| 69 |
< |
|
| 70 |
< |
## process.ecalUncalibRecHit = cms.EDProducer( |
| 71 |
< |
## "EcalUncalibRecHitProducer", |
| 72 |
< |
## EBdigiCollection = cms.InputTag("ecalEBunpacker","ebDigis"), |
| 73 |
< |
## EEhitCollection = cms.string('EcalUncalibRecHitsEE'), |
| 74 |
< |
## EEdigiCollection = cms.InputTag("ecalEBunpacker","eeDigis"), |
| 75 |
< |
## EBhitCollection = cms.string('EcalUncalibRecHitsEB'), |
| 76 |
< |
## algo = cms.string("EcalUncalibRecHitWorkerWeights") |
| 77 |
< |
## ) |
| 78 |
< |
|
| 79 |
< |
## process.ecalRecHit = cms.EDProducer( |
| 80 |
< |
## "EcalRecHitProducer", |
| 81 |
< |
## EErechitCollection = cms.string(''), |
| 82 |
< |
## EEuncalibRecHitCollection = cms.InputTag(''), |
| 83 |
< |
## EBuncalibRecHitCollection = cms.InputTag('ecalUncalibRecHit', 'EcalUncalibRecHitsEB'), |
| 84 |
< |
## EBrechitCollection = cms.string('EcalRecHitsEB'), |
| 85 |
< |
## algo = cms.string('EcalRecHitWorkerSimple') |
| 86 |
< |
## ) |
| 56 |
> |
process.load("andersj.HcalPlotter.ecalLocalRecoSequence_cff") |
| 57 |
> |
process.ecalDataSequence = cms.Sequence(process.ecalEBunpacker*process.ecalGlobalUncalibRecHit*process.ecalRecHit) |
| 58 |
|
|
| 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 |
< |
process.load("RecoTBCalo.HcalTBObjectUnpacker.HcalTBObjectUnpacker_Normal_cfi") |
| 63 |
> |
process.load("andersj.HcalPlotter.HcalTBObjectUnpacker_Normal_cfi") |
| 64 |
|
|
| 65 |
|
process.horeco.firstSample=2 |
| 66 |
|
process.horeco.samplesToAdd=6 |
| 70 |
|
|
| 71 |
|
## process.plotanal=cms.EDAnalyzer( |
| 72 |
|
## "HcalHOTBPlotAnal", |
| 73 |
< |
## hbheRHtag = cms.untracked.InputTag("hbhereco"), |
| 73 |
> |
## hbheRHtag = cms.untracked.InputTag("hbheprereco"), |
| 74 |
|
## hoRHtag = cms.untracked.InputTag("horeco"), |
| 75 |
|
## hfRHtag = cms.untracked.InputTag("hfreco"), |
| 76 |
|
## hcalDigiTag = cms.untracked.InputTag("hcalDigis"), |
| 77 |
|
## hcalTrigTag = cms.untracked.InputTag("tbunpack"), |
| 78 |
+ |
## ebRHtag = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"), |
| 79 |
+ |
## doBeamCounters = cms.untracked.bool(True), |
| 80 |
|
## calibFC2GeV = cms.untracked.double(0.01), |
| 81 |
|
## HistoParameters = cms.PSet( |
| 82 |
|
## pedGeVlo = cms.double(-15), |
| 90 |
|
## otherGeVlo = cms.double(-5), |
| 91 |
|
## otherGeVhi = cms.double(250), |
| 92 |
|
## beamGeVlo = cms.double(-20), |
| 93 |
< |
## beamGeVhi = cms.double(50), |
| 93 |
> |
## beamGeVhi = cms.double(500), |
| 94 |
|
## #beamGeVhi = cms.double(80*200), |
| 95 |
|
## timeNSlo = cms.double(50), |
| 96 |
|
## timeNShi = cms.double(250) |
| 99 |
|
|
| 100 |
|
process.dump = cms.EDAnalyzer("EventContentAnalyzer") |
| 101 |
|
|
| 102 |
< |
#process.p=cms.Path(process.ecalEBunpacker+process.hcalDigis+process.hbhereco+process.horeco+process.hfreco+process.tbunpack+process.dump) |
| 130 |
< |
process.p=cms.Path(process.ecalEBunpacker+process.ecalGlobalUncalibRecHit+process.ecalRecHit+process.hcalDigis+process.hbhereco+process.horeco+process.hfreco+process.tbunpack+process.dump) |
| 131 |
< |
#process.p=cms.Path(process.ecalEBunpacker+process.ecalLocalRecoSequence+process.hcalDigis+process.hbhereco+process.horeco+process.hfreco+process.tbunpack+process.dump) |
| 102 |
> |
process.p=cms.Path(process.ecalDataSequence+process.hcalDigis+process.hbheprereco+process.horeco+process.hfreco+process.tbunpack+process.dump) |
