ProducerES_cfi.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

TrackTrigger_params = cms.PSet (

    fromDD4hep = cms.bool(False),

  # Parameter to check if configured Tracker Geometry is supported
  SupportedGeometry = cms.PSet (
    XMLLabel    = cms.string ("geomXMLFiles"                                    ), # label of ESProducer/ESSource
    XMLPath     = cms.string ("Geometry/TrackerCommonData/data/PhaseII/"        ), # compared path
    XMLFile     = cms.string ("tracker.xml"                                     ), # compared filen ame
    XMLVersions = cms.vstring("TiltedTracker613", "TiltedTracker613_MB_2019_04", "OuterTracker616_2020_04", "OuterTracker800_2020_07", "Tracker_DD4hep_compatible_2021_02", "Tracker_DD4hep_compatible_OT801_IT640" )  # list of supported versions
  ),

  # Parameter to check if Process History is consistent with process configuration
  ProcessHistory = cms.PSet (
    GeometryConfiguration = cms.string( "XMLIdealGeometryESSource@"                    ), # label of compared GeometryConfiguration
    TTStubAlgorithm       = cms.string( "TTStubAlgorithm_official_Phase2TrackerDigi_@" )  # label of compared TTStubAlgorithm
  ),

  # Common track finding parameter
  TrackFinding = cms.PSet (
    BeamWindowZ      = cms.double( 15. ), # half lumi region size in cm
    MatchedLayers    = cms.int32 (  4  ), # required number of layers a found track has to have in common with a TP to consider it matched to it
    MatchedLayersPS  = cms.int32 (  0  ), # required number of ps layers a found track has to have in common with a TP to consider it matched to it
    UnMatchedStubs   = cms.int32 (  1  ), # allowed number of stubs a found track may have not in common with its matched TP
    UnMatchedStubsPS = cms.int32 (  0  )  # allowed number of PS stubs a found track may have not in common with its matched TP
  ),

  # TMTT specific parameter
  TMTT = cms.PSet (
    MinPt            = cms.double(  3.   ), # cut on stub and TP pt, also defines region overlap shape in GeV
    MaxEta           = cms.double(  2.4  ), # cut on stub eta
    ChosenRofPhi     = cms.double( 67.24 ), # critical radius defining region overlap shape in cm
    NumLayers        = cms.int32 (  7    ), # number of detector layers a reconstructbale particle may cross
    WidthR           = cms.int32 ( 12    ), # number of bits used for stub r - ChosenRofPhi
    WidthPhi         = cms.int32 ( 14    ), # number of bits used for stub phi w.r.t. phi sector centre
    WidthZ           = cms.int32 ( 14    )  # number of bits used for stub z
  ),

  # Hybrid specific parameter
  Hybrid = cms.PSet (
    MinPt        = cms.double(  2.   ),                        # cut on stub and TP pt, also defines region overlap shape in GeV
    MaxEta       = cms.double(  2.5  ),                        # cut on stub eta
    ChosenRofPhi = cms.double( 55.   ),                        # critical radius defining region overlap shape in cm
    NumLayers    = cms.int32 (  4    ),                        # max number of detector layer connected to one DTC
    NumRingsPS   = cms.vint32 ( 11, 11, 8, 8, 8 ),             # number of outer PS rings for disk 1, 2, 3, 4, 5
    WidthsR      = cms.vint32 (   7,     7,    12,      7   ), # number of bits used for stub r w.r.t layer/disk centre for module types (barrelPS, barrel2S, diskPS, disk2S)
    WidthsZ      = cms.vint32 (  12,     8,     7,      7   ), # number of bits used for stub z w.r.t layer/disk centre for module types (barrelPS, barrel2S, diskPS, disk2S)
    WidthsPhi    = cms.vint32 (  14,    17,    14,     14   ), # number of bits used for stub phi w.r.t. region centre for module types (barrelPS, barrel2S, diskPS, disk2S)
    WidthsAlpha  = cms.vint32 (   0,     0,     0,      4   ), # number of bits used for stub row number for module types (barrelPS, barrel2S, diskPS, disk2S)
    WidthsBend   = cms.vint32 (   3,     4,     3,      4   ), # number of bits used for stub bend number for module types (barrelPS, barrel2S, diskPS, disk2S)
    RangesR      = cms.vdouble(   7.5,   7.5, 120. ,    0.  ), # range in stub r which needs to be covered for module types (barrelPS, barrel2S, diskPS, disk2S)
    RangesZ      = cms.vdouble( 240.,  240.,    7.5,    7.5 ), # range in stub z which needs to be covered for module types (barrelPS, barrel2S, diskPS, disk2S)
    RangesAlpha  = cms.vdouble(   0.,    0.,    0.,  1024.  ), # range in stub row which needs to be covered for module types (barrelPS, barrel2S, diskPS, disk2S)
    LayerRs      = cms.vdouble(  24.9316,  37.1777,  52.2656,  68.7598,  86.0156, 108.3105 ), # mean radius of outer tracker barrel layer
    DiskZs       = cms.vdouble( 131.1914, 154.9805, 185.3320, 221.6016, 265.0195           ), # mean z of outer tracker endcap disks
    Disk2SRsSet  = cms.VPSet(                                                                 # center radius of outer tracker endcap 2S diks strips
      cms.PSet( Disk2SRs = cms.vdouble( 66.4391, 71.4391, 76.2750, 81.2750, 82.9550, 87.9550, 93.8150, 98.8150, 99.8160, 104.8160 ) ), # disk 1
      cms.PSet( Disk2SRs = cms.vdouble( 66.4391, 71.4391, 76.2750, 81.2750, 82.9550, 87.9550, 93.8150, 98.8150, 99.8160, 104.8160 ) ), # disk 2
      cms.PSet( Disk2SRs = cms.vdouble( 63.9903, 68.9903, 74.2750, 79.2750, 81.9562, 86.9562, 92.4920, 97.4920, 99.8160, 104.8160 ) ), # disk 3
      cms.PSet( Disk2SRs = cms.vdouble( 63.9903, 68.9903, 74.2750, 79.2750, 81.9562, 86.9562, 92.4920, 97.4920, 99.8160, 104.8160 ) ), # disk 4
      cms.PSet( Disk2SRs = cms.vdouble( 63.9903, 68.9903, 74.2750, 79.2750, 81.9562, 86.9562, 92.4920, 97.4920, 99.8160, 104.8160 ) )  # disk 5
    )
  ),

  # Parameter specifying TrackingParticle used for Efficiency measurements
  TrackingParticle = cms.PSet (
    MaxEta           = cms.double(  2.4 ), # eta cut
    MaxVertR         = cms.double(  1.  ), # cut on vertex pos r in cm
    MaxVertZ         = cms.double( 30.  ), # cut on vertex pos z in cm
    MaxD0            = cms.double(  5.  ), # cut on impact parameter in cm
    MinLayers        = cms.int32 (  4   ), # required number of associated layers to a TP to consider it reconstruct-able
    MinLayersPS      = cms.int32 (  0   )  # required number of associated ps layers to a TP to consider it reconstruct-able
  ),

  # Fimrware specific Parameter
  Firmware = cms.PSet (
    NumFramesInfra = cms.int32 (   6                ), # needed gap between events of emp-infrastructure firmware
    FreqLHC        = cms.double(  40.               ), # LHC bunch crossing rate in MHz
    FreqBE         = cms.double( 360.               ), # processing Frequency of DTC & TFP in MHz, has to be integer multiple of FreqLHC
    TMP_FE         = cms.int32 (   8                ), # number of events collected in front-end
    TMP_TFP        = cms.int32 (  18                ), # time multiplexed period of track finding processor
    SpeedOfLight   = cms.double(   2.99792458       ), # in e8 m/s
    BField         = cms.double(   3.81120228767395 ), # in T
    BFieldError    = cms.double(   1.e-6            ), # accepted difference to EventSetup in T
    OuterRadius    = cms.double( 112.7              ), # outer radius of outer tracker in cm
    InnerRadius    = cms.double(  21.8              ), # inner radius of outer tracker in cm
    HalfLength     = cms.double( 270.               ), # half length of outer tracker in cm
    MaxPitch       = cms.double(    .01             )  # max strip/pixel pitch of outer tracker sensors in cm
  ),

  # Parmeter specifying front-end
  FrontEnd = cms.PSet (
    WidthBend      = cms.int32 (  6      ), # number of bits used for internal stub bend
    WidthCol       = cms.int32 (  5      ), # number of bits used for internal stub column
    WidthRow       = cms.int32 ( 11      ), # number of bits used for internal stub row
    BaseBend       = cms.double(   .25   ), # precision of internal stub bend in pitch units
    BaseCol        = cms.double(  1.     ), # precision of internal stub column in pitch units
    BaseRow        = cms.double(   .5    ), # precision of internal stub row in pitch units
    BaseWindowSize = cms.double(   .5    ), # precision of window sizes in pitch units
    BendCut        = cms.double(  1.3125 )  # used stub bend uncertainty in pitch units
  ),

  # Parmeter specifying DTC 
  DTC = cms.PSet (
    NumRegions            = cms.int32 (  9    ), # number of phi slices the outer tracker readout is organized in
    NumOverlappingRegions = cms.int32 (  2    ), # number of regions a reconstructable particles may cross
    NumATCASlots          = cms.int32 ( 12    ), # number of Slots in used ATCA crates
    NumDTCsPerRegion      = cms.int32 ( 24    ), # number of DTC boards used to readout a detector region, likely constructed to be an integerer multiple of NumSlots_
    NumModulesPerDTC      = cms.int32 ( 72    ), # max number of sensor modules connected to one DTC board
    NumRoutingBlocks      = cms.int32 (  2    ), # number of systiloic arrays in stub router firmware
    DepthMemory           = cms.int32 ( 64    ), # fifo depth in stub router firmware
    WidthRowLUT           = cms.int32 (  4    ), # number of row bits used in look up table
    WidthQoverPt          = cms.int32 (  9    ), # number of bits used for stub qOverPt. lut addr is col + bend = 11 => 1 BRAM -> 18 bits for min and max val -> 9
    OffsetDetIdDSV        = cms.int32 (  1    ), # tk layout det id minus DetSetVec->detId
    OffsetDetIdTP         = cms.int32 ( -1    ), # tk layout det id minus TrackerTopology lower det id
    OffsetLayerDisks      = cms.int32 ( 10    ), # offset in layer ids between barrel layer and endcap disks
    OffsetLayerId         = cms.int32 (  1    )  # offset between 0 and smallest layer id (barrel layer 1)
  ),

  # Parmeter specifying GeometricProcessor
  GeometricProcessor = cms.PSet (
    NumSectorsPhi = cms.int32 (  2  ), # number of phi sectors used in hough transform
    ChosenRofZ    = cms.double( 50. ), # critical radius defining r-z sector shape in cm
    RangeChiZ     = cms.double( 90. ), # range of stub z residual w.r.t. sector center which needs to be covered
    DepthMemory   = cms.int32 ( 64  ), # fifo depth in stub router firmware
    BoundariesEta = cms.vdouble( -2.40, -2.08, -1.68, -1.26, -0.90, -0.62, -0.41, -0.20, 0.0, 0.20, 0.41, 0.62, 0.90, 1.26, 1.68, 2.08, 2.40 ) # defining r-z sector shape
  ),

  # Parmeter specifying HoughTransform
  HoughTransform = cms.PSet (
    NumBinsQoverPt = cms.int32( 16 ), # number of used qOverPt bins
    NumBinsPhiT    = cms.int32( 32 ), # number of used phiT bins
    MinLayers      = cms.int32(  5 ), # required number of stub layers to form a candidate
    DepthMemory    = cms.int32( 32 )  # internal fifo depth
  ),

  # Parmeter specifying MiniHoughTransform
  MiniHoughTransform = cms.PSet (
    NumBinsQoverPt = cms.int32( 2 ), # number of finer qOverPt bins inside HT bin
    NumBinsPhiT    = cms.int32( 2 ), # number of finer phiT bins inside HT bin
    NumDLB         = cms.int32( 2 ), # number of dynamic load balancing steps
    MinLayers      = cms.int32( 5 )  # required number of stub layers to form a candidate
  ),

  # Parmeter specifying SeedFilter
  SeedFilter = cms.PSet (
    PowerBaseCot     = cms.int32( -6 ), # used cot(Theta) bin width = 2 ** this
    BaseDiffZ        = cms.int32(  4 ), # used zT bin width = baseZ * 2 ** this
    MinLayers        = cms.int32(  4 )  # required number of stub layers to form a candidate
  ),

  # Parmeter specifying KalmanFilter
  KalmanFilter = cms.PSet (
    WidthLutInvPhi   = cms.int32(  10 ), # number of bits for internal reciprocal look up
    WidthLutInvZ     = cms.int32(  10 ), # number of bits for internal reciprocal look up
    NumTracks        = cms.int32(  16 ), # cut on number of input candidates
    MinLayers        = cms.int32(   4 ), # required number of stub layers to form a track
    MaxLayers        = cms.int32(   4 ), # maximum number of  layers added to a track
    MaxStubsPerLayer = cms.int32(   4 ), # cut on number of stub per layer for input candidates
    MaxSkippedLayers = cms.int32(   2 ), # maximum allowed skipped layers from inside to outside to form a track
    BaseShiftr0      = cms.int32(  -3 ),
    BaseShiftr02     = cms.int32(  -5 ),
    BaseShiftv0      = cms.int32(  -2 ),
    BaseShiftS00     = cms.int32(  -1 ),
    BaseShiftS01     = cms.int32(  -7 ),
    BaseShiftK00     = cms.int32(  -9 ),
    BaseShiftK10     = cms.int32( -15 ),
    BaseShiftR00     = cms.int32(  -2 ),
    BaseShiftInvR00  = cms.int32( -19 ),
    BaseShiftChi20   = cms.int32(  -5 ),
    BaseShiftC00     = cms.int32(   5 ),
    BaseShiftC01     = cms.int32(  -3 ),
    BaseShiftC11     = cms.int32(  -7 ),
    BaseShiftr1      = cms.int32(   2 ),
    BaseShiftr12     = cms.int32(   5 ),
    BaseShiftv1      = cms.int32(   3 ),
    BaseShiftS12     = cms.int32(  -3 ),
    BaseShiftS13     = cms.int32(  -3 ),
    BaseShiftK21     = cms.int32( -13 ),
    BaseShiftK31     = cms.int32( -14 ),
    BaseShiftR11     = cms.int32(   3 ),
    BaseShiftInvR11  = cms.int32( -21 ),
    BaseShiftChi21   = cms.int32(  -5 ),
    BaseShiftC22     = cms.int32(  -3 ),
    BaseShiftC23     = cms.int32(  -5 ),
    BaseShiftC33     = cms.int32(  -3 ),
    BaseShiftChi2    = cms.int32(  -5 )
  ),

  # Parmeter specifying DuplicateRemoval
  DuplicateRemoval = cms.PSet (
    DepthMemory  = cms.int32( 16 ), # internal memory depth
    WidthPhi0    = cms.int32( 12 ), # number of bist used for phi0
    WidthQoverPt = cms.int32( 15 ), # number of bist used for qOverPt
    WidthCot     = cms.int32( 16 ), # number of bist used for cot(theta)
    WidthZ0      = cms.int32( 12 )  # number of bist used for z0
  )

)