Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-09-23 22:53:26

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 # This object modifies hcalSimParameters for different scenarios
0004 
0005 hcalSimParameters = cms.PSet(
0006     #  In HF, the SimHits energy is actually
0007     # the number of photoelectrons from the shower
0008     # library.  However, we need a lot more
0009     # smearing, because the first level of the PMT
0010     # only gives ~6 electrons per pe, which
0011     # comes out to a 40% smearing of the single pe peak!
0012     #
0013     hf1 = cms.PSet(
0014         readoutFrameSize = cms.int32(4),
0015         binOfMaximum = cms.int32(3),
0016         samplingFactor = cms.double(0.383),
0017         doPhotoStatistics = cms.bool(True),
0018         photoelectronsToAnalog = cms.double(2.79),
0019         simHitToPhotoelectrons = cms.double(6.0),
0020         syncPhase = cms.bool(True),
0021         timePhase = cms.double(14.0),
0022         doSiPMSmearing = cms.bool(False),
0023         sipmTau = cms.double(0.),
0024         threshold_currentTDC = cms.double(-999.),
0025         delayQIE = cms.int32(-999),
0026     ),
0027     hf2 = cms.PSet(
0028         readoutFrameSize = cms.int32(4),
0029         binOfMaximum = cms.int32(3),
0030         samplingFactor = cms.double(0.368),
0031         doPhotoStatistics = cms.bool(True),
0032         photoelectronsToAnalog = cms.double(1.843),
0033         simHitToPhotoelectrons = cms.double(6.0),
0034         syncPhase = cms.bool(True),
0035         timePhase = cms.double(13.0),
0036         doSiPMSmearing = cms.bool(False),
0037         sipmTau = cms.double(0.),
0038         threshold_currentTDC = cms.double(-999.),
0039         delayQIE = cms.int32(-999),
0040     ),
0041     ho = cms.PSet(
0042         readoutFrameSize = cms.int32(10),
0043         firstRing = cms.int32(1),
0044         binOfMaximum = cms.int32(5),
0045         doPhotoStatistics = cms.bool(True),
0046         simHitToPhotoelectrons = cms.double(4000.0), # is not actually used
0047         samplingFactors = cms.vdouble(231.0, 231.0, 231.0, 231.0,
0048             360.0, 360.0, 360.0, 360.0, 360.0, 360.0,
0049             360.0, 360.0, 360.0, 360.0, 360.0),
0050         syncPhase = cms.bool(True),
0051         timePhase = cms.double(5.0),
0052         timeSmearing = cms.bool(False),
0053         # 0 is HPD, 1 is SiPM, 2 fetches HPD/Zecotek/Hamamatsu from DB
0054         siPMCode = cms.int32(2),
0055         doSiPMSmearing = cms.bool(False),
0056         sipmTau = cms.double(5.),
0057         threshold_currentTDC = cms.double(-999.),
0058         delayQIE = cms.int32(-999),
0059     ),
0060     hb = cms.PSet(
0061         readoutFrameSize = cms.int32(10),
0062         firstRing = cms.int32(1),
0063         binOfMaximum = cms.int32(5),
0064         doPhotoStatistics = cms.bool(True),
0065         simHitToPhotoelectrons = cms.double(2000.0),
0066         samplingFactors = cms.vdouble(
0067             125.44, 125.54, 125.32, 125.13, 124.46,
0068             125.01, 125.22, 125.48, 124.45, 125.90,
0069             125.83, 127.01, 126.82, 129.73, 131.83,
0070             143.52),            
0071         syncPhase = cms.bool(True),
0072         timePhase = cms.double(6.0),
0073         timeSmearing = cms.bool(True),
0074         doSiPMSmearing = cms.bool(False),
0075         sipmTau = cms.double(0.),
0076         threshold_currentTDC = cms.double(-999.),
0077         delayQIE = cms.int32(-999),
0078     ),
0079     he = cms.PSet(
0080         readoutFrameSize = cms.int32(10),
0081         firstRing = cms.int32(16),
0082         binOfMaximum = cms.int32(5),
0083         doPhotoStatistics = cms.bool(True),
0084         simHitToPhotoelectrons = cms.double(2000.0),
0085         samplingFactors = cms.vdouble(
0086             210.55, 197.93, 186.12, 189.64, 189.63,
0087             190.28, 189.61, 189.60, 190.12, 191.22,
0088             190.90, 193.06, 188.42, 188.42),
0089         syncPhase = cms.bool(True),
0090         timePhase = cms.double(6.0),
0091         timeSmearing = cms.bool(True),
0092         doSiPMSmearing = cms.bool(False),
0093         sipmTau = cms.double(0.),
0094         threshold_currentTDC = cms.double(-999.),
0095         delayQIE = cms.int32(-999),
0096     ),
0097     zdc = cms.PSet(
0098         readoutFrameSize = cms.int32(10),
0099         binOfMaximum = cms.int32(5),
0100         samplingFactor = cms.double(1.000),
0101         doPhotoStatistics = cms.bool(True),
0102         photoelectronsToAnalog = cms.double(1.843),
0103         simHitToPhotoelectrons = cms.double(6.0),
0104         syncPhase = cms.bool(True),
0105         timePhase = cms.double(-4.0),
0106         doSiPMSmearing = cms.bool(False),
0107         sipmTau = cms.double(0.),
0108         threshold_currentTDC = cms.double(-999.),
0109         delayQIE = cms.int32(-999),
0110     ),
0111 )
0112 
0113 hcalSimParameters.hoZecotek = hcalSimParameters.ho.clone()
0114 
0115 hcalSimParameters.hoHamamatsu = hcalSimParameters.ho.clone()
0116 
0117 # Customises the HCal digitiser for post LS1 running
0118 from Configuration.Eras.Modifier_run2_common_cff import run2_common
0119 run2_common.toModify( hcalSimParameters, 
0120     ho = dict(
0121         siPMCode = cms.int32(1)
0122     ),
0123     hf1 = dict( samplingFactor = cms.double(0.335) ),
0124     hf2 = dict( samplingFactor = cms.double(0.335) )
0125 )
0126 
0127 from Configuration.Eras.Modifier_run2_HE_2017_cff import run2_HE_2017
0128 run2_HE_2017.toModify( hcalSimParameters,
0129     he = dict(
0130         doSiPMSmearing = cms.bool(True),
0131         sipmTau = cms.double(10.),
0132     )
0133 )
0134 
0135 from Configuration.Eras.Modifier_run2_HF_2017_cff import run2_HF_2017
0136 run2_HF_2017.toModify( hcalSimParameters,
0137     hf1 = dict(
0138                readoutFrameSize = cms.int32(3), 
0139                binOfMaximum     = cms.int32(2),
0140                threshold_currentTDC = cms.double(3.),
0141               ),
0142     hf2 = dict(
0143                readoutFrameSize = cms.int32(3), 
0144                binOfMaximum     = cms.int32(2),
0145                threshold_currentTDC = cms.double(3.),
0146               )
0147 )
0148 
0149 from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
0150 run2_HB_2018.toModify( hcalSimParameters,
0151     hb = dict(
0152                readoutFrameSize = cms.int32(8), 
0153                binOfMaximum     = cms.int32(4)
0154               )
0155 )
0156 from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
0157 run2_HE_2018.toModify( hcalSimParameters,
0158     he = dict(
0159                readoutFrameSize = cms.int32(8), 
0160                binOfMaximum     = cms.int32(4),
0161                threshold_currentTDC = cms.double(18.7)
0162               )
0163 )
0164 
0165 
0166 from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
0167 run3_HB.toModify( hcalSimParameters,
0168     hb = dict(
0169         doSiPMSmearing = cms.bool(True),
0170         threshold_currentTDC = cms.double(18.7),
0171         sipmTau = cms.double(10.),
0172     )
0173 )
0174 
0175 from Configuration.Eras.Modifier_run3_common_cff import run3_common
0176 run3_common.toModify( hcalSimParameters, 
0177     hb = dict(
0178                readoutFrameSize = cms.int32(10), 
0179                binOfMaximum     = cms.int32(6)
0180               ),
0181     he = dict(
0182                readoutFrameSize = cms.int32(10), 
0183                binOfMaximum     = cms.int32(6)
0184               ),
0185     hf1 = dict( samplingFactor = 0.37,
0186                 timePhase = 9.0 
0187                ),
0188     hf2 = dict( samplingFactor = 0.37,
0189                 timePhase = 8.0
0190                )
0191 ) 
0192 
0193 
0194 _newFactors = cms.vdouble(
0195     210.55, 197.93, 186.12, 189.64, 189.63,
0196     189.96, 190.03, 190.11, 190.18, 190.25,
0197     190.32, 190.40, 190.47, 190.54, 190.61,
0198     190.69, 190.83, 190.94, 190.94, 190.94,
0199     190.94, 190.94, 190.94, 190.94, 190.94,
0200     190.94, 190.94, 190.94, 190.94, 190.94,
0201     190.94, 190.94, 190.94, 190.94, 190.94,
0202     190.94, 190.94, 190.94, 190.94, 190.94,
0203     190.94, 190.94, 190.94, 190.94, 190.94,
0204     190.94, 190.94, 190.94, 190.94, 190.94,
0205     190.94, 190.94, 190.94, 190.94, 190.94,
0206     190.94, 190.94, 190.94, 190.94, 190.94,
0207     190.94, 190.94, 190.94, 190.94, 190.94,
0208     190.94, 190.94, 190.94, 190.94, 190.94,
0209     190.94, 190.94, 190.94, 190.94, 190.94,
0210     190.94, 190.94, 190.94, 190.94, 190.94,
0211     190.94, 190.94, 190.94, 190.94, 190.94 )
0212 
0213 from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
0214 phase2_hcal.toModify( hcalSimParameters,
0215     he = dict(
0216         samplingFactors = _newFactors,
0217     )
0218 )