| 1 |
yilmaz |
1.1 |
|
| 2 |
|
|
|
| 3 |
|
|
RUN="___RUN___"
|
| 4 |
|
|
|
| 5 |
|
|
import FWCore.ParameterSet.Config as cms
|
| 6 |
|
|
|
| 7 |
|
|
process = cms.Process("HCAL")
|
| 8 |
|
|
process.load("FWCore.MessageLogger.MessageLogger_cfi")
|
| 9 |
|
|
|
| 10 |
|
|
|
| 11 |
|
|
process.load("Configuration.StandardSequences.Geometry_cff")
|
| 12 |
|
|
process.load("RecoLocalCalo.Configuration.RecoLocalCalo_Cosmics_cff")
|
| 13 |
|
|
process.load("CondCore.DBCommon.CondDBSetup_cfi")
|
| 14 |
|
|
process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
|
| 15 |
|
|
process.load("EventFilter.HcalRawToDigi.HcalRawToDigi_cfi")
|
| 16 |
|
|
process.load("Configuration.StandardSequences.RawToDigi_Data_cff")
|
| 17 |
|
|
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
|
| 18 |
|
|
|
| 19 |
|
|
# Event selection
|
| 20 |
|
|
process.load("EventFilter.HcalRawToDigi.HcalCalibTypeFilter_cfi")
|
| 21 |
|
|
process.hcalCalibTypeFilter.CalibTypes = [2]
|
| 22 |
|
|
process.hcalCalibTypeFilter.FilterSummary = True
|
| 23 |
|
|
process.hcalCalibTypeFilter.InputLabel = "hltHcalCalibrationRaw"
|
| 24 |
|
|
|
| 25 |
|
|
# Unpack HF data
|
| 26 |
|
|
import EventFilter.HcalRawToDigi.HcalRawToDigi_cfi
|
| 27 |
|
|
process.hcalDigis = EventFilter.HcalRawToDigi.HcalRawToDigi_cfi.hcalDigis.clone()
|
| 28 |
|
|
#process.load("Configuration.StandardSequences.ReconstructionCosmics_cff")
|
| 29 |
|
|
process.hcalDigis.InputLabel = "hltHcalCalibrationRaw"
|
| 30 |
|
|
|
| 31 |
|
|
process.options = cms.untracked.PSet(
|
| 32 |
|
|
SkipEvent = cms.untracked.vstring('ProductNotFound'),
|
| 33 |
|
|
wantSummary = cms.untracked.bool(True)
|
| 34 |
|
|
)
|
| 35 |
|
|
|
| 36 |
|
|
process.GlobalTag.globaltag = 'GR_P_V37::All'
|
| 37 |
|
|
process.prefer("GlobalTag")
|
| 38 |
|
|
process.Timing = cms.Service("Timing")
|
| 39 |
|
|
|
| 40 |
|
|
process.maxEvents = cms.untracked.PSet(
|
| 41 |
|
|
input = cms.untracked.int32(-1)
|
| 42 |
|
|
)
|
| 43 |
|
|
|
| 44 |
|
|
process.source = cms.Source("PoolSource",
|
| 45 |
|
|
fileNames = cms.untracked.vstring(
|
| 46 |
|
|
"rfio:/castor/cern.ch/cms/store/data/Run2012C/TestEnablesEcalHcalDT/RAW/v1/000/201/136/44A5E464-5BE8-E111-8BAF-001D09F251FE.root"
|
| 47 |
|
|
# "/store/data/Run2012C/TestEnablesEcalHcalDT/RAW/v1/000/201/136/44A5E464-5BE8-E111-8BAF-001D09F251FE.root"
|
| 48 |
|
|
),
|
| 49 |
|
|
)
|
| 50 |
|
|
|
| 51 |
|
|
process.hcal_db_producer = cms.ESProducer("HcalDbProducer",
|
| 52 |
|
|
dump = cms.untracked.vstring(''),
|
| 53 |
|
|
file = cms.untracked.string('')
|
| 54 |
|
|
)
|
| 55 |
|
|
|
| 56 |
|
|
process.hcalLaserReco = cms.EDProducer("HcalLaserReco")
|
| 57 |
|
|
|
| 58 |
|
|
process.myanalysis = cms.EDAnalyzer("HcalRaddamAnalyzer",
|
| 59 |
|
|
hcalLaserLabel = cms.InputTag("hcalLaserReco"),
|
| 60 |
|
|
nevtsample = cms.untracked.int32(100),
|
| 61 |
|
|
STreshold = cms.untracked.double(20.0),
|
| 62 |
|
|
ADCTreshold = cms.untracked.double(20.0),
|
| 63 |
|
|
DHFL_Treshold = cms.untracked.double(3.8),
|
| 64 |
|
|
p4noise = cms.untracked.double(69.0),
|
| 65 |
|
|
DHFS_Treshold = cms.untracked.double(3.8),
|
| 66 |
|
|
EventsToStore = cms.untracked.double(100.0),
|
| 67 |
|
|
rsigma = cms.untracked.double(3.0),
|
| 68 |
|
|
ADCdigiTS = cms.untracked.double(4.0),
|
| 69 |
|
|
HFL_Treshold = cms.untracked.double(0.65),
|
| 70 |
|
|
outputFileName = cms.string('NtupleLaser'+RUN+'.root'),
|
| 71 |
|
|
textoutputFileName = cms.string('newraddam'+RUN+'.txt'),
|
| 72 |
|
|
radfib2 = cms.untracked.int32(4),
|
| 73 |
|
|
radfib1 = cms.untracked.int32(3),
|
| 74 |
|
|
LASERON = cms.bool(True),
|
| 75 |
|
|
hfperchan = cms.bool(True),
|
| 76 |
|
|
HFS_Treshold = cms.untracked.double(0.9),
|
| 77 |
|
|
dsigma = cms.untracked.double(12.0),
|
| 78 |
|
|
)
|
| 79 |
|
|
|
| 80 |
|
|
process.TFileService = cms.Service("TFileService",
|
| 81 |
|
|
fileName = cms.string('HistoLaser'+RUN+'.root')
|
| 82 |
|
|
)
|
| 83 |
|
|
|
| 84 |
|
|
process.RawToDigi = cms.Sequence(process.hcalDigis)
|
| 85 |
|
|
process.allPath = cms.Path(
|
| 86 |
|
|
process.hcalCalibTypeFilter*
|
| 87 |
|
|
process.hcalDigis*process.hcalLaserReco*process.myanalysis)
|
