mtdDigitizer_cfi.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

_barrel_tile_MTDDigitizer = cms.PSet(
    digitizerName     = cms.string("BTLTileDigitizer"),
    inputSimHits      = cms.InputTag("g4SimHits:FastTimerHitsBarrel"),
    digiCollectionTag = cms.string("FTLBarrel"),
    maxSimHitsAccTime = cms.uint32(100),
    premixStage1      = cms.bool(False),
    premixStage1MinCharge = cms.double(1e-4), # These give relative precision of 1.1 %,
    premixStage1MaxCharge = cms.double(1e6),  # can the range be squeezed?
    DeviceSimulation = cms.PSet(
        bxTime                   = cms.double(25),      # [ns] 
        LightYield               = cms.double(40000.),  # [photons/MeV]
        LightCollectionEff       = cms.double(0.25),
        LightCollectionTime      = cms.double(0.2),     # [ns]
        smearLightCollectionTime = cms.double(0.),      # [ns]
        PhotonDetectionEff       = cms.double(0.20),
        ),
    ElectronicsSimulation = cms.PSet(
        bxTime                     = cms.double(25),    # [ns]
        TestBeamMIPTimeRes         = cms.double(4.293), # This is given by 0.048[ns]*sqrt(8000.), in order to
                                                        # rescale the time resolution of 1 MIP = 8000 p.e.
        ScintillatorRiseTime       = cms.double(1.1),   # [ns]
        ScintillatorDecayTime      = cms.double(40.),   # [ns]
        ChannelTimeOffset          = cms.double(0.),    # [ns]
        smearChannelTimeOffset     = cms.double(0.),    # [ns]
        EnergyThreshold            = cms.double(4.),    # [photo-electrons]
        TimeThreshold1             = cms.double(20.),   # [photo-electrons]
        TimeThreshold2             = cms.double(50.),   # [photo-electrons]
        ReferencePulseNpe          = cms.double(100.),  # [photo-electrons]
        DarkCountRate              = cms.double(10.),   # [GHz]
        SinglePhotonTimeResolution = cms.double(0.060), # [ns]
        SigmaElectronicNoise       = cms.double(1.),    # [p.e.]
        SigmaClock                 = cms.double(0.015), # [ns]
        CorrelationCoefficient     = cms.double(1.),
        Npe_to_pC                  = cms.double(0.016), # [pC] 
        Npe_to_V                   = cms.double(0.0064),# [V] 

        # n bits for the ADC 
        adcNbits          = cms.uint32(10),
        # n bits for the TDC
        tdcNbits          = cms.uint32(10),
        # ADC saturation
        adcSaturation_MIP = cms.double(600.),           # [pC]
        # for different thickness
        adcThreshold_MIP   = cms.double(0.064),         # [pC]
        # LSB for time of arrival estimate from TDC
        toaLSB_ns         = cms.double(0.020),          # [ns]
        )
    )

_barrel_bar_MTDDigitizer = cms.PSet(
    digitizerName     = cms.string("BTLBarDigitizer"),
    inputSimHits      = cms.InputTag("g4SimHits:FastTimerHitsBarrel"),
    digiCollectionTag = cms.string("FTLBarrel"),
    maxSimHitsAccTime = cms.uint32(100),
    premixStage1      = cms.bool(False),
    premixStage1MinCharge = cms.double(1e-4),
    premixStage1MaxCharge = cms.double(1e6),
    DeviceSimulation = cms.PSet(
        bxTime                   = cms.double(25),      # [ns]
        LightYield               = cms.double(40000.),  # [photons/MeV]
        LightCollectionEff       = cms.double(0.25),
        LightCollectionSlopeR    = cms.double(0.075),   # [ns/cm]
        LightCollectionSlopeL    = cms.double(0.075),   # [ns/cm]
        PhotonDetectionEff       = cms.double(0.20),
        ),
    # The readout electronics is the same for both geometries:
    ElectronicsSimulation = _barrel_tile_MTDDigitizer.ElectronicsSimulation
)

_endcap_MTDDigitizer = cms.PSet(
    digitizerName     = cms.string("ETLDigitizer"),
    inputSimHits      = cms.InputTag("g4SimHits:FastTimerHitsEndcap"),
    digiCollectionTag = cms.string("FTLEndcap"),
    maxSimHitsAccTime = cms.uint32(100),
    premixStage1      = cms.bool(False),
    premixStage1MinCharge = cms.double(1e-4),
    premixStage1MaxCharge = cms.double(1e6),
    DeviceSimulation  = cms.PSet(
        bxTime            = cms.double(25),
        tofDelay          = cms.double(1),
        meVPerMIP         = cms.double(0.085), # from HGCal
        ),
    ElectronicsSimulation = cms.PSet(
        bxTime             = cms.double(25),
        etaResolution      = cms.string("0.03+0.0025*x"), # This is just a dummy dependence on eta.
        # n bits for the ADC 
        adcNbits           = cms.uint32(8),
        # n bits for the TDC
        tdcNbits           = cms.uint32(11),
        # ADC saturation
        adcSaturation_MIP  = cms.double(25),
        # for different thickness
        adcThreshold_MIP   = cms.double(0.025),
        # LSB for time of arrival estimate from TDC in ns
        toaLSB_ns          = cms.double(0.013),
        )
)

from Configuration.ProcessModifiers.premix_stage1_cff import premix_stage1
for _m in [_barrel_tile_MTDDigitizer, _barrel_bar_MTDDigitizer, _endcap_MTDDigitizer]:
    premix_stage1.toModify(_m, premixStage1 = True)

# Fast Timing
mtdDigitizer = cms.PSet( 
    accumulatorType   = cms.string("MTDDigiProducer"),
    makeDigiSimLinks  = cms.bool(False),
    verbosity         = cms.untracked.uint32(0),

    barrelDigitizer = _barrel_tile_MTDDigitizer,
    endcapDigitizer = _endcap_MTDDigitizer
)