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