HcalPlugins/python/HBHEDarkening_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 # energy = center of mass energy (beam energy x2), GeV

0004 # [intlumi] = fb-1

0005 # [lumirate] = fb-1/hr

0006 # dose maps only produced for 8 and 14 TeV - used also for 7 and 13 TeV, respectively

0007
0008 # total for phase0/1 ~ 300 fb-1

0009 _years_LHC = cms.VPSet(
0010     cms.PSet(year = cms.string("2011"), intlumi = cms.double(5.6), lumirate = cms.double(0.005), energy = cms.int32(8)),
0011     cms.PSet(year = cms.string("2012"), intlumi = cms.double(23.3), lumirate = cms.double(0.013), energy = cms.int32(8)),
0012     cms.PSet(year = cms.string("2015"), intlumi = cms.double(4.1), lumirate = cms.double(0.009), energy = cms.int32(14)),
0013     cms.PSet(year = cms.string("2016"), intlumi = cms.double(41.0), lumirate = cms.double(0.026), energy = cms.int32(14)),
0014     cms.PSet(year = cms.string("2017"), intlumi = cms.double(45.0), lumirate = cms.double(0.043), energy = cms.int32(14)),
0015     cms.PSet(year = cms.string("2018"), intlumi = cms.double(45.0), lumirate = cms.double(0.043), energy = cms.int32(14)),
0016     cms.PSet(year = cms.string("2021"), intlumi = cms.double(45.0), lumirate = cms.double(0.05), energy = cms.int32(14)),
0017     cms.PSet(year = cms.string("2022"), intlumi = cms.double(45.0), lumirate = cms.double(0.05), energy = cms.int32(14)),
0018     cms.PSet(year = cms.string("2023"), intlumi = cms.double(50.0), lumirate = cms.double(0.05), energy = cms.int32(14)),
0019 )
0020
0021 # total for phase2 nominal = 3000 fb-1 (including phase0/1) @ 5.0E34/cm^2/s

0022 _years_HLLHC_nominal = cms.VPSet(
0023     cms.PSet(year = cms.string("2038"), intlumi = cms.double(2700), lumirate = cms.double(0.15), energy = cms.int32(14)),
0024 )
0025
0026 # total for phase2 ultimate = 4500 fb-1 (including phase0/1) @ 7.5E34/cm^2/s

0027 _years_HLLHC_ultimate = cms.VPSet(
0028     cms.PSet(year = cms.string("2029"), intlumi = cms.double(700), lumirate = cms.double(0.15), energy = cms.int32(14)),
0029     cms.PSet(year = cms.string("2039"), intlumi = cms.double(3500), lumirate = cms.double(0.225), energy = cms.int32(14)),
0030 )
0031
0032 HEDarkeningEP = cms.ESSource("HBHEDarkeningEP",
0033     appendToDataLabel = cms.string("HE"),
0034     ieta_shift = cms.int32(16),
0035     # parameters taken from https://indico.cern.ch/event/721565/contributions/2966213/attachments/1632432/2603148/2018.04.12_Updated_DvsR_fits_and_predictions.pdf, slide 3

0036     drdA = cms.double(2.7383),
0037     drdB = cms.double(0.37471),
0038     dosemaps = cms.VPSet(
0039         cms.PSet(energy = cms.int32(8), file = cms.FileInPath("CalibCalorimetry/HcalPlugins/data/dosemapHE_4TeV.txt")),
0040         cms.PSet(energy = cms.int32(14), file = cms.FileInPath("CalibCalorimetry/HcalPlugins/data/dosemapHE_7TeV.txt")),
0041     ),
0042     years = _years_LHC + _years_HLLHC_nominal,
0043 )
0044
0045 HBDarkeningEP = HEDarkeningEP.clone(
0046     appendToDataLabel = cms.string("HB"),
0047     ieta_shift = cms.int32(1),
0048     dosemaps = cms.VPSet(
0049         cms.PSet(energy = cms.int32(8), file = cms.FileInPath("CalibCalorimetry/HcalPlugins/data/dosemapHB_4TeV.txt")),
0050         cms.PSet(energy = cms.int32(14), file = cms.FileInPath("CalibCalorimetry/HcalPlugins/data/dosemapHB_7TeV.txt")),
0051     ),
0052 )