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