FastTimingCommon/python/mtdDigitizer_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 _barrel_MTDDigitizer = cms.PSet(
0004     digitizerName     = cms.string("BTLDigitizer"),
0005     inputSimHits      = cms.InputTag("g4SimHits:FastTimerHitsBarrel"),
0006     digiCollectionTag = cms.string("FTLBarrel"),
0007     maxSimHitsAccTime = cms.uint32(100),
0008     premixStage1      = cms.bool(False),
0009     premixStage1MinCharge = cms.double(1e-4),
0010     premixStage1MaxCharge = cms.double(1e6),
0011     DeviceSimulation = cms.PSet(
0012         bxTime                   = cms.double(25),      # [ns]
0013         LightYield               = cms.double(40000.),  # [photons/MeV]
0014         LightCollectionEff       = cms.double(0.25),
0015         LightCollectionSlopeR    = cms.double(0.075),   # [ns/cm]
0016         LightCollectionSlopeL    = cms.double(0.075),   # [ns/cm]
0017         PhotonDetectionEff       = cms.double(0.20),
0018         ),
0019     ElectronicsSimulation = cms.PSet(
0020         bxTime                     = cms.double(25),    # [ns]
0021         TestBeamMIPTimeRes         = cms.double(4.293), # This is given by 0.048[ns]*sqrt(8000.), in order to
0022                                                         # rescale the time resolution of 1 MIP = 8000 p.e.
0023         ScintillatorRiseTime       = cms.double(1.1),   # [ns]
0024         ScintillatorDecayTime      = cms.double(40.),   # [ns]
0025         ChannelTimeOffset          = cms.double(0.),    # [ns]
0026         smearChannelTimeOffset     = cms.double(0.),    # [ns]
0027         EnergyThreshold            = cms.double(4.),    # [photo-electrons]
0028         TimeThreshold1             = cms.double(20.),   # [photo-electrons]
0029         TimeThreshold2             = cms.double(50.),   # [photo-electrons]
0030         ReferencePulseNpe          = cms.double(100.),  # [photo-electrons]
0031         DarkCountRate              = cms.double(10.),   # [GHz]
0032         SinglePhotonTimeResolution = cms.double(0.060), # [ns]
0033         SigmaElectronicNoise       = cms.double(1.),    # [p.e.]
0034         SigmaClock                 = cms.double(0.015), # [ns]
0035         CorrelationCoefficient     = cms.double(1.),
0036         SmearTimeForOOTtails       = cms.bool(True),
0037         Npe_to_pC                  = cms.double(0.016), # [pC]
0038         Npe_to_V                   = cms.double(0.0064),# [V]
0039
0040         # n bits for the ADC
0041         adcNbits          = cms.uint32(10),
0042         # n bits for the TDC
0043         tdcNbits          = cms.uint32(10),
0044         # ADC saturation
0045         adcSaturation_MIP = cms.double(600.),           # [pC]
0046         # for different thickness
0047         adcThreshold_MIP   = cms.double(0.064),         # [pC]
0048         # LSB for time of arrival estimate from TDC
0049         toaLSB_ns         = cms.double(0.020),          # [ns]
0050         )
0051
0052
0053 )
0054
0055 _endcap_MTDDigitizer = cms.PSet(
0056     digitizerName     = cms.string("ETLDigitizer"),
0057     inputSimHits      = cms.InputTag("g4SimHits:FastTimerHitsEndcap"),
0058     digiCollectionTag = cms.string("FTLEndcap"),
0059     maxSimHitsAccTime = cms.uint32(100),
0060     premixStage1      = cms.bool(False),
0061     premixStage1MinCharge = cms.double(1e-4),
0062     premixStage1MaxCharge = cms.double(1e6),
0063     DeviceSimulation  = cms.PSet(
0064         bxTime            = cms.double(25),
0065         tofDelay          = cms.double(1),
0066         meVPerMIP         = cms.double(0.085), # from HGCal
0067         ),
0068     ElectronicsSimulation = cms.PSet(
0069         bxTime               = cms.double(25),
0070         IntegratedLuminosity = cms.double(1000.),      # [1/fb]
0071         FluenceVsRadius      = cms.string("1.937*TMath::Power(x,-1.706)"),
0072         LGADGainVsFluence    = cms.string("TMath::Min(15.,30.-x)"),
0073         TimeResolution2      = cms.string("0.0225/x"), # [ns^2]
0074         # n bits for the ADC
0075         adcNbits             = cms.uint32(8),
0076         # n bits for the TDC
0077         tdcNbits             = cms.uint32(11),
0078         # ADC saturation
0079         adcSaturation_MIP  = cms.double(25),
0080         # for different thickness
0081         adcThreshold_MIP   = cms.double(0.025),
0082         # LSB for time of arrival estimate from TDC in ns
0083         toaLSB_ns          = cms.double(0.013),
0084         )
0085 )
0086
0087 from Configuration.Eras.Modifier_phase2_etlV4_cff import phase2_etlV4
0088 phase2_etlV4.toModify(_endcap_MTDDigitizer.DeviceSimulation, meVPerMIP = 0.015 )
0089
0090 from Configuration.ProcessModifiers.premix_stage1_cff import premix_stage1
0091 for _m in [_barrel_MTDDigitizer, _endcap_MTDDigitizer]:
0092     premix_stage1.toModify(_m, premixStage1 = True)
0093
0094 # Fast Timing
0095 mtdDigitizer = cms.PSet(
0096     accumulatorType   = cms.string("MTDDigiProducer"),
0097     makeDigiSimLinks  = cms.bool(False),
0098     verbosity         = cms.untracked.uint32(0),
0099
0100     barrelDigitizer = _barrel_MTDDigitizer,
0101     endcapDigitizer = _endcap_MTDDigitizer
0102 )