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