File indexing completed on 2023-03-17 10:53:20
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 import sys
0012 import FWCore.ParameterSet.Config as cms
0013 from Configuration.StandardSequences.Eras import eras
0014
0015 process = cms.Process("TEST", eras.Run3)
0016 process.load("Configuration.StandardSequences.GeometryDB_cff")
0017 process.load("CondCore.CondDB.CondDB_cfi")
0018 process.load("EventFilter.L1GlobalTriggerRawToDigi.l1GtUnpack_cfi")
0019 process.l1GtUnpack.DaqGtInputTag = 'source'
0020
0021 process.load('Configuration.StandardSequences.Services_cff')
0022 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0023 process.load('FWCore.MessageService.MessageLogger_cfi')
0024 process.load('Configuration.EventContent.EventContent_cff')
0025 process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
0026 process.load('Configuration.StandardSequences.Reconstruction_Data_cff')
0027 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0028
0029 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0030 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0031
0032
0033 runnumber = sys.argv[2]
0034 rundir = sys.argv[3]
0035 histodir = sys.argv[4]
0036
0037
0038
0039
0040
0041
0042 process.maxEvents = cms.untracked.PSet(
0043
0044 input = cms.untracked.int32(-1)
0045 )
0046
0047
0048
0049
0050
0051
0052
0053
0054
0055
0056
0057
0058
0059 process.TFileService = cms.Service("TFileService",
0060 fileName = cms.string(histodir+'/LED_'+runnumber+'.root')
0061
0062 )
0063
0064
0065
0066
0067
0068 process.source = cms.Source("HcalTBSource",
0069 skipBadFiles=cms.untracked.bool(True),
0070 firstLuminosityBlockForEachRun = cms.untracked.VLuminosityBlockID([]),
0071 firstRun = cms.untracked.uint32(331370),
0072
0073
0074 fileNames = cms.untracked.vstring(
0075 rundir+'/run'+runnumber+'/USC_'+runnumber+'.root'
0076
0077
0078
0079 ),
0080 secondaryFileNames = cms.untracked.vstring()
0081 )
0082
0083
0084 process.Analyzer = cms.EDAnalyzer("CMTRawAnalyzer",
0085
0086 Verbosity = cms.untracked.int32(0),
0087
0088 MapCreation = cms.untracked.int32(1),
0089
0090 recordNtuples = cms.untracked.bool(False),
0091 maxNeventsInNtuple = cms.int32(1),
0092 recordHistoes = cms.untracked.bool(True),
0093 studyRunDependenceHist = cms.untracked.bool(True),
0094 studyCapIDErrorsHist = cms.untracked.bool(True),
0095 studyRMSshapeHist = cms.untracked.bool(True),
0096 studyRatioShapeHist = cms.untracked.bool(True),
0097 studyADCAmplHist = cms.untracked.bool(True),
0098 studyTSmeanShapeHist = cms.untracked.bool(True),
0099 studyTSmaxShapeHist = cms.untracked.bool(True),
0100 studyCalibCellsHist = cms.untracked.bool(True),
0101 studyDiffAmplHist = cms.untracked.bool(True),
0102 studyPedestalsHist = cms.untracked.bool(True),
0103 studyPedestalCorrelations = cms.untracked.bool(True),
0104
0105 ratioHBMin = cms.double(0.70),
0106 ratioHBMax = cms.double(0.94),
0107 ratioHEMin = cms.double(0.60),
0108 ratioHEMax = cms.double(0.95),
0109 ratioHFMin = cms.double(0.45),
0110 ratioHFMax = cms.double(1.02),
0111 ratioHOMin = cms.double(0.40),
0112 ratioHOMax = cms.double(1.04),
0113 rmsHBMin = cms.double(0.7),
0114 rmsHBMax = cms.double(2.5),
0115 rmsHEMin = cms.double(0.7),
0116 rmsHEMax = cms.double(2.2),
0117 rmsHFMin = cms.double(0.1),
0118 rmsHFMax = cms.double(2.6),
0119 rmsHOMin = cms.double(0.1),
0120 rmsHOMax = cms.double(2.8),
0121 TSmeanHBMin = cms.double(2.5),
0122 TSmeanHBMax = cms.double(5.5),
0123 TSmeanHEMin = cms.double(1.0),
0124 TSmeanHEMax = cms.double(5.2),
0125 TSmeanHFMin = cms.double(1.0),
0126 TSmeanHFMax = cms.double(4.2),
0127 TSmeanHOMin = cms.double(1.0),
0128 TSmeanHOMax = cms.double(4.8),
0129 TSpeakHBMin = cms.double(2.2),
0130 TSpeakHBMax = cms.double(5.5),
0131 TSpeakHEMin = cms.double(1.5),
0132 TSpeakHEMax = cms.double(6.5),
0133 TSpeakHFMin = cms.double(0.5),
0134 TSpeakHFMax = cms.double(4.5),
0135 TSpeakHOMin = cms.double(0.5),
0136 TSpeakHOMax = cms.double(7.5),
0137 ADCAmplHBMin = cms.double(10000.),
0138 ADCAmplHBMax = cms.double(300000.),
0139 ADCAmplHEMin = cms.double(20000.),
0140 ADCAmplHEMax = cms.double(300000.),
0141 ADCAmplHFMin = cms.double(50.),
0142 ADCAmplHFMax = cms.double(9000.),
0143 ADCAmplHOMin = cms.double(50.),
0144 ADCAmplHOMax = cms.double(9000.),
0145 pedestalwHBMax = cms.double(0.1),
0146 pedestalwHEMax = cms.double(0.1),
0147 pedestalwHFMax = cms.double(0.4),
0148 pedestalwHOMax = cms.double(0.1),
0149 pedestalHBMax = cms.double(0.1),
0150 pedestalHEMax = cms.double(0.6),
0151 pedestalHFMax = cms.double(0.8),
0152 pedestalHOMax = cms.double(0.1),
0153
0154 calibrADCHBMin = cms.double(1000.),
0155 calibrADCHBMax = cms.double(100000000.),
0156 calibrADCHEMin = cms.double(1000.),
0157 calibrADCHEMax = cms.double(100000000.),
0158 calibrADCHOMin = cms.double(1000.),
0159 calibrADCHOMax = cms.double(100000000.),
0160 calibrADCHFMin = cms.double(100.),
0161 calibrADCHFMax = cms.double(100000000.),
0162 calibrRatioHBMin = cms.double(0.76),
0163 calibrRatioHBMax = cms.double(0.94),
0164 calibrRatioHEMin = cms.double(0.76),
0165 calibrRatioHEMax = cms.double(0.94),
0166 calibrRatioHOMin = cms.double(0.85),
0167 calibrRatioHOMax = cms.double(0.99),
0168 calibrRatioHFMin = cms.double(0.5),
0169 calibrRatioHFMax = cms.double(0.8),
0170 calibrTSmaxHBMin = cms.double(1.50),
0171 calibrTSmaxHBMax = cms.double(2.50),
0172 calibrTSmaxHEMin = cms.double(1.50),
0173 calibrTSmaxHEMax = cms.double(2.50),
0174 calibrTSmaxHOMin = cms.double(1.50),
0175 calibrTSmaxHOMax = cms.double(2.50),
0176 calibrTSmaxHFMin = cms.double(3.50),
0177 calibrTSmaxHFMax = cms.double(4.50),
0178 calibrTSmeanHBMin = cms.double(2.40),
0179 calibrTSmeanHBMax = cms.double(3.70),
0180 calibrTSmeanHEMin = cms.double(2.40),
0181 calibrTSmeanHEMax = cms.double(3.70),
0182 calibrTSmeanHOMin = cms.double(1.50),
0183 calibrTSmeanHOMax = cms.double(2.70),
0184 calibrTSmeanHFMin = cms.double(3.50),
0185 calibrTSmeanHFMax = cms.double(4.50),
0186 calibrWidthHBMin = cms.double(1.30),
0187 calibrWidthHBMax = cms.double(1.90),
0188 calibrWidthHEMin = cms.double(1.30),
0189 calibrWidthHEMax = cms.double(1.90),
0190 calibrWidthHOMin = cms.double(0.70),
0191 calibrWidthHOMax = cms.double(1.65),
0192 calibrWidthHFMin = cms.double(0.30),
0193 calibrWidthHFMax = cms.double(1.50),
0194
0195
0196
0197
0198
0199 flagtoaskrunsorls = cms.int32(0),
0200
0201
0202
0203
0204 flagtodefinebadchannel = cms.int32(0),
0205
0206 howmanybinsonplots = cms.int32(25),
0207
0208
0209 lsmin = cms.int32(1),
0210
0211 lsmax = cms.int32(2600),
0212
0213 flagabortgaprejected = cms.int32(1),
0214 bcnrejectedlow = cms.int32(3446),
0215 bcnrejectedhigh= cms.int32(3564),
0216
0217
0218
0219 flagcpuoptimization = cms.int32(0),
0220
0221
0222
0223 flagestimatornormalization = cms.int32(1),
0224
0225
0226
0227
0228
0229 lsdep_cut1_peak_HBdepth1 = cms.int32(20),
0230 lsdep_cut1_peak_HBdepth2 = cms.int32(7),
0231 lsdep_cut1_peak_HEdepth1 = cms.int32(16),
0232 lsdep_cut1_peak_HEdepth2 = cms.int32(13),
0233 lsdep_cut1_peak_HEdepth3 = cms.int32(4),
0234 lsdep_cut1_peak_HFdepth1 = cms.int32(10),
0235 lsdep_cut1_peak_HFdepth2 = cms.int32(5),
0236 lsdep_cut1_peak_HOdepth4 = cms.int32(45),
0237
0238 lsdep_cut3_max_HBdepth1 = cms.int32(19),
0239 lsdep_cut3_max_HBdepth2 = cms.int32(6),
0240 lsdep_cut3_max_HEdepth1 = cms.int32(15),
0241 lsdep_cut3_max_HEdepth2 = cms.int32(12),
0242 lsdep_cut3_max_HEdepth3 = cms.int32(3),
0243 lsdep_cut3_max_HFdepth1 = cms.int32(9),
0244 lsdep_cut3_max_HFdepth2 = cms.int32(4),
0245 lsdep_cut3_max_HOdepth4 = cms.int32(40),
0246
0247
0248 lsdep_estimator1_HBdepth1 = cms.double(2500.),
0249 lsdep_estimator1_HBdepth2 = cms.double(2500.),
0250 lsdep_estimator1_HBdepth3 = cms.double(2500.),
0251 lsdep_estimator1_HBdepth4 = cms.double(2500.),
0252 lsdep_estimator1_HEdepth1 = cms.double(2500.),
0253 lsdep_estimator1_HEdepth2 = cms.double(2500.),
0254 lsdep_estimator1_HEdepth3 = cms.double(2500.),
0255 lsdep_estimator1_HEdepth4 = cms.double(2500.),
0256 lsdep_estimator1_HEdepth5 = cms.double(2500.),
0257 lsdep_estimator1_HEdepth6 = cms.double(2500.),
0258 lsdep_estimator1_HEdepth7 = cms.double(2500.),
0259 lsdep_estimator1_HFdepth1 = cms.double(2500.),
0260 lsdep_estimator1_HFdepth2 = cms.double(2500.),
0261 lsdep_estimator1_HFdepth3 = cms.double(2500.),
0262 lsdep_estimator1_HFdepth4 = cms.double(2500.),
0263 lsdep_estimator1_HOdepth4 = cms.double(2500.),
0264
0265 lsdep_estimator2_HBdepth1 = cms.double(7.),
0266 lsdep_estimator2_HBdepth2 = cms.double(7.),
0267 lsdep_estimator2_HEdepth1 = cms.double(7.),
0268 lsdep_estimator2_HEdepth2 = cms.double(7.),
0269 lsdep_estimator2_HEdepth3 = cms.double(7.),
0270 lsdep_estimator2_HFdepth1 = cms.double(7.),
0271 lsdep_estimator2_HFdepth2 = cms.double(7.),
0272 lsdep_estimator2_HOdepth4 = cms.double(7.),
0273
0274 lsdep_estimator3_HBdepth1 = cms.double(7.),
0275 lsdep_estimator3_HBdepth2 = cms.double(7.),
0276 lsdep_estimator3_HEdepth1 = cms.double(7.),
0277 lsdep_estimator3_HEdepth2 = cms.double(7.),
0278 lsdep_estimator3_HEdepth3 = cms.double(7.),
0279 lsdep_estimator3_HFdepth1 = cms.double(7.),
0280 lsdep_estimator3_HFdepth2 = cms.double(7.),
0281 lsdep_estimator3_HOdepth4 = cms.double(7.),
0282
0283 lsdep_estimator4_HBdepth1 = cms.double(5.),
0284 lsdep_estimator4_HBdepth2 = cms.double(5.),
0285 lsdep_estimator4_HEdepth1 = cms.double(5.),
0286 lsdep_estimator4_HEdepth2 = cms.double(5.),
0287 lsdep_estimator4_HEdepth3 = cms.double(5.),
0288 lsdep_estimator4_HFdepth1 = cms.double(5.),
0289 lsdep_estimator4_HFdepth2 = cms.double(5.),
0290 lsdep_estimator4_HOdepth4 = cms.double(5.),
0291
0292 lsdep_estimator5_HBdepth1 = cms.double(1.8),
0293 lsdep_estimator5_HBdepth2 = cms.double(1.8),
0294 lsdep_estimator5_HEdepth1 = cms.double(1.8),
0295 lsdep_estimator5_HEdepth2 = cms.double(1.8),
0296 lsdep_estimator5_HEdepth3 = cms.double(1.8),
0297 lsdep_estimator5_HFdepth1 = cms.double(1.8),
0298 lsdep_estimator5_HFdepth2 = cms.double(1.8),
0299 lsdep_estimator5_HOdepth4 = cms.double(1.8),
0300
0301
0302 useADCmassive = cms.untracked.bool(True),
0303 useADCfC = cms.untracked.bool(False),
0304 useADCcounts = cms.untracked.bool(False),
0305
0306
0307
0308 usePedestalSubtraction = cms.untracked.bool(False),
0309
0310
0311
0312
0313 forallestimators_amplitude_bigger = cms.double(-100.),
0314
0315
0316 usecontinuousnumbering = cms.untracked.bool(True),
0317
0318
0319
0320
0321
0322 flagToUseDigiCollectionsORNot = cms.int32(1),
0323
0324
0325
0326
0327
0328 hcalCalibDigiCollectionTag = cms.InputTag('hcalDigis'),
0329 hbheDigiCollectionTag = cms.InputTag('hcalDigis'),
0330 hoDigiCollectionTag = cms.InputTag('hcalDigis'),
0331 hfDigiCollectionTag = cms.InputTag('hcalDigis'),
0332
0333
0334 hbheQIE11DigiCollectionTag = cms.InputTag('hcalDigis'),
0335 hbheQIE10DigiCollectionTag = cms.InputTag('hcalDigis'),
0336
0337
0338
0339
0340
0341
0342
0343
0344
0345
0346
0347
0348
0349
0350
0351
0352
0353
0354 flagupgradeqie1011 = cms.int32(6),
0355
0356
0357 flaguseshunt = cms.int32(6),
0358
0359 flagsipmcorrection = cms.int32(1),
0360
0361
0362
0363
0364 flagLaserRaddam = cms.int32(0),
0365
0366 flagfitshunt1pedorledlowintensity = cms.int32(0),
0367
0368 splashesUpperLimit = cms.int32(10000),
0369
0370
0371
0372
0373
0374
0375
0376 flagIterativeMethodCalibrationGroupDigi = cms.int32(0),
0377
0378
0379 flagIterativeMethodCalibrationGroupReco = cms.int32(0),
0380
0381 hbheInputSignalTag = cms.InputTag('hbherecoMBNZS'),
0382 hbheInputNoiseTag = cms.InputTag('hbherecoNoise'),
0383 hfInputSignalTag = cms.InputTag('hfrecoMBNZS'),
0384 hfInputNoiseTag = cms.InputTag('hfrecoNoise'),
0385
0386
0387
0388
0389
0390
0391
0392
0393 )
0394
0395 process.hcal_db_producer = cms.ESProducer("HcalDbProducer",
0396 dump = cms.untracked.vstring(''),
0397 file = cms.untracked.string('')
0398 )
0399
0400 process.es_hardcode = cms.ESSource("HcalHardcodeCalibrations",
0401 toGet = cms.untracked.vstring('QIEShape',
0402 'QIEData',
0403 'ChannelQuality',
0404 'HcalQIEData',
0405 'Pedestals',
0406 'PedestalWidths',
0407 'Gains',
0408 'GainWidths',
0409 'ZSThresholds',
0410 'RespCorrs')
0411 )
0412
0413 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0414 from Configuration.AlCa.autoCond import autoCond
0415
0416 process.GlobalTag.globaltag = '104X_dataRun2_v1'
0417
0418
0419
0420
0421 process.load('Configuration.StandardSequences.RawToDigi_Data_cff')
0422 process.hcalDigis.FilterDataQuality = cms.bool(False)
0423 process.hcalDigis.InputLabel = cms.InputTag("source")
0424
0425
0426
0427
0428
0429
0430 process.load('EventFilter.HcalRawToDigi.hcalRawToDigi_cfi')
0431 process.hcalDigis= process.hcalRawToDigi.clone(
0432 FilterDataQuality = False,
0433 InputLabel = "source",
0434
0435 )
0436
0437
0438 process.p = cms.Path(process.hcalDigis*process.Analyzer)
0439
0440
0441 process.MessageLogger = cms.Service("MessageLogger",
0442 categories = cms.untracked.vstring(''),
0443 destinations = cms.untracked.vstring('cout'),
0444 debugModules = cms.untracked.vstring('*'),
0445 cout = cms.untracked.PSet(
0446 threshold = cms.untracked.string('WARNING'),
0447 WARNING = cms.untracked.PSet(limit = cms.untracked.int32(0))
0448 )
0449 )
0450
0451
0452
0453
0454
0455