RecoTracker/IterativeTracking/python/ElectronSeeds_cff.py - master - githubmirror/cmssw

RecoTracker/IterativeTracking/python/ElectronSeeds_cff.py
342 строки16 KB
Marco Musich
adapt ElectronSeeds to work with the trackingIters01 processModifier
10 сен 2025, 12:03
10 сен 2025, 12:032aec040
import FWCore.ParameterSet.Config as cms

from RecoLocalTracker.SubCollectionProducers.SeedClusterRemover_cfi import seedClusterRemover
initialStepSeedClusterMask = seedClusterRemover.clone(
    trajectories = 'initialStepSeeds',
    oldClusterRemovalInfo = cms.InputTag("pixelLessStepClusters")
)

from Configuration.ProcessModifiers.trackingIters01_cff import trackingIters01
trackingIters01.toModify(initialStepSeedClusterMask,
                         oldClusterRemovalInfo = cms.InputTag('highPtTripletStepClusters')
                         )

from RecoLocalTracker.SubCollectionProducers.seedClusterRemoverPhase2_cfi import seedClusterRemoverPhase2
from Configuration.Eras.Modifier_trackingPhase2PU140_cff import trackingPhase2PU140
trackingPhase2PU140.toReplaceWith(initialStepSeedClusterMask, seedClusterRemoverPhase2.clone(
    trajectories = 'initialStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('highPtTripletStepClusters')
    )
)

highPtTripletStepSeedClusterMask = seedClusterRemover.clone( # for Phase2PU140
    trajectories = 'highPtTripletStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('initialStepSeedClusterMask')
)
pixelPairStepSeedClusterMask = seedClusterRemover.clone(
    trajectories = 'pixelPairStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('initialStepSeedClusterMask')
)

trackingPhase2PU140.toReplaceWith(highPtTripletStepSeedClusterMask, seedClusterRemoverPhase2.clone(
    trajectories = 'highPtTripletStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('initialStepSeedClusterMask')
    )
)
trackingPhase2PU140.toReplaceWith(pixelPairStepSeedClusterMask, seedClusterRemoverPhase2.clone(
    trajectories = 'detachedQuadStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('highPtTripletStepSeedClusterMask')
    )
)

# This is a pure guess to use detachedTripletStep for phase1 here instead of the pixelPair in Run2 configuration
detachedTripletStepSeedClusterMask = seedClusterRemover.clone(
    trajectories = 'lowPtTripletStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('initialStepSeedClusterMask')
)
mixedTripletStepSeedClusterMask = seedClusterRemover.clone(
    trajectories = 'mixedTripletStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('pixelPairStepSeedClusterMask')
)
from Configuration.Eras.Modifier_trackingPhase1_cff import trackingPhase1
trackingPhase1.toModify(mixedTripletStepSeedClusterMask,
    oldClusterRemovalInfo = 'detachedTripletStepSeedClusterMask'
)
pixelLessStepSeedClusterMask = seedClusterRemover.clone(
    trajectories = 'pixelLessStepSeeds',
    oldClusterRemovalInfo = cms.InputTag('mixedTripletStepSeedClusterMask')
)

import RecoTracker.TkSeedingLayers.seedingLayersEDProducer_cfi as _mod

tripletElectronSeedLayers = _mod.seedingLayersEDProducer.clone(
    layerList = ['BPix1+BPix2+BPix3', 
                 'BPix1+BPix2+FPix1_pos', 'BPix1+BPix2+FPix1_neg', 
                 'BPix1+FPix1_pos+FPix2_pos', 'BPix1+FPix1_neg+FPix2_neg'],
    BPix = dict(
        TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4PixelTriplets'),
        HitProducer = cms.string('siPixelRecHits'),
        skipClusters = cms.InputTag('pixelLessStepSeedClusterMask')
    ),
    FPix = dict(
        TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4PixelTriplets'),
        HitProducer = cms.string('siPixelRecHits'),
        skipClusters = cms.InputTag('pixelLessStepSeedClusterMask')
    )
)
_layerListForPhase1 = [
    'BPix1+BPix2+BPix3', 'BPix2+BPix3+BPix4',
    'BPix1+BPix3+BPix4', 'BPix1+BPix2+BPix4',
    'BPix2+BPix3+FPix1_pos', 'BPix2+BPix3+FPix1_neg',
    'BPix1+BPix2+FPix1_pos', 'BPix1+BPix2+FPix1_neg',
    'BPix1+BPix3+FPix1_pos', 'BPix1+BPix3+FPix1_neg',
    'BPix2+FPix1_pos+FPix2_pos', 'BPix2+FPix1_neg+FPix2_neg',
    'BPix1+FPix1_pos+FPix2_pos', 'BPix1+FPix1_neg+FPix2_neg',
    'BPix1+BPix2+FPix2_pos', 'BPix1+BPix2+FPix2_neg',
    'FPix1_pos+FPix2_pos+FPix3_pos', 'FPix1_neg+FPix2_neg+FPix3_neg',
    'BPix1+FPix2_pos+FPix3_pos', 'BPix1+FPix2_neg+FPix3_neg',
    'BPix1+FPix1_pos+FPix3_pos', 'BPix1+FPix1_neg+FPix3_neg'
]
trackingPhase1.toModify(tripletElectronSeedLayers, layerList = _layerListForPhase1)

(trackingIters01 & ~trackingPhase2PU140).toModify(tripletElectronSeedLayers,
    BPix = dict(skipClusters = cms.InputTag('highPtTripletStepSeedClusterMask')),
    FPix = dict(skipClusters = cms.InputTag('highPtTripletStepSeedClusterMask'))
)

trackingPhase2PU140.toModify(tripletElectronSeedLayers,
    layerList = _layerListForPhase1,
    BPix = dict(skipClusters = 'pixelPairStepSeedClusterMask'),
    FPix = dict(skipClusters = 'pixelPairStepSeedClusterMask')
)

from RecoTracker.TkTrackingRegions.globalTrackingRegionFromBeamSpot_cfi import globalTrackingRegionFromBeamSpot as _globalTrackingRegionFromBeamSpot
tripletElectronTrackingRegions = _globalTrackingRegionFromBeamSpot.clone(RegionPSet = dict(
    ptMin        = 1.0,
    originRadius = 0.02,
    nSigmaZ      = 4.0
))

from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
from RecoTracker.TkTrackingRegions.globalTrackingRegionWithVertices_cff import globalTrackingRegionWithVertices as _globalTrackingRegionWithVertices
pp_on_AA.toReplaceWith(tripletElectronTrackingRegions,
    _globalTrackingRegionWithVertices.clone(
        RegionPSet = dict(
            fixedError   = 0.5,
            ptMin        = 8.0,
            originRadius = 0.02
        )
))

from RecoTracker.TkHitPairs.hitPairEDProducer_cfi import hitPairEDProducer as _hitPairEDProducer
tripletElectronHitDoublets = _hitPairEDProducer.clone(
    seedingLayers   = 'tripletElectronSeedLayers',
    trackingRegions = 'tripletElectronTrackingRegions',
    maxElement      = 50000000,
    produceIntermediateHitDoublets = True,
)
from RecoTracker.PixelSeeding.pixelTripletHLTEDProducer_cfi import pixelTripletHLTEDProducer as _pixelTripletHLTEDProducer
import RecoTracker.PixelLowPtUtilities.LowPtClusterShapeSeedComparitor_cfi
tripletElectronHitTriplets = _pixelTripletHLTEDProducer.clone(
    doublets   = 'tripletElectronHitDoublets',
    maxElement = 1000000,
    produceSeedingHitSets = True,
)
from RecoTracker.TkSeedGenerator.seedCreatorFromRegionConsecutiveHitsEDProducer_cff import seedCreatorFromRegionConsecutiveHitsEDProducer as _seedCreatorFromRegionConsecutiveHitsEDProducer
tripletElectronSeeds = _seedCreatorFromRegionConsecutiveHitsEDProducer.clone(
    seedingHitSets = 'tripletElectronHitTriplets',
)
trackingPhase2PU140.toModify(tripletElectronHitTriplets,
    maxElement = 0,
)

from RecoLocalTracker.SubCollectionProducers.SeedClusterRemover_cfi import seedClusterRemover
tripletElectronClusterMask = seedClusterRemover.clone(
    trajectories = 'tripletElectronSeeds',
    oldClusterRemovalInfo = cms.InputTag('pixelLessStepSeedClusterMask')
)

(trackingIters01 & ~trackingPhase2PU140).toModify(tripletElectronClusterMask,
    oldClusterRemovalInfo = cms.InputTag('highPtTripletStepSeedClusterMask')
)

trackingPhase2PU140.toReplaceWith(tripletElectronClusterMask, seedClusterRemoverPhase2.clone(
    trajectories = 'tripletElectronSeeds',
    oldClusterRemovalInfo = cms.InputTag('pixelLessStepSeedClusterMask')
    )
)

pixelPairElectronSeedLayers = _mod.seedingLayersEDProducer.clone(
    layerList = ['BPix1+BPix2', 'BPix1+BPix3', 'BPix2+BPix3', 
                 'BPix1+FPix1_pos', 'BPix1+FPix1_neg', 
                 'BPix1+FPix2_pos', 'BPix1+FPix2_neg', 
                 'BPix2+FPix1_pos', 'BPix2+FPix1_neg', 
                 'FPix1_pos+FPix2_pos', 'FPix1_neg+FPix2_neg'],
    BPix = dict(
        TTRHBuilder = cms.string('WithTrackAngle'),
        HitProducer = cms.string('siPixelRecHits'),
        skipClusters = cms.InputTag('tripletElectronClusterMask')
    ),
    FPix = dict(
        TTRHBuilder = cms.string('WithTrackAngle'),
        HitProducer = cms.string('siPixelRecHits'),
        skipClusters = cms.InputTag('tripletElectronClusterMask')
    )
)
_layerListForPhase1 = [
        'BPix1+BPix2', 'BPix1+BPix3', 'BPix1+BPix4',
        'BPix2+BPix3', 'BPix2+BPix4',
        'BPix3+BPix4',
        'BPix1+FPix1_pos', 'BPix1+FPix1_neg', 
        'BPix1+FPix2_pos', 'BPix1+FPix2_neg', 
        'BPix2+FPix1_pos', 'BPix2+FPix1_neg', 
        'FPix1_pos+FPix2_pos', 'FPix1_neg+FPix2_neg',
        'FPix1_pos+FPix3_pos', 'FPix1_neg+FPix3_neg',
        'FPix2_pos+FPix3_pos', 'FPix2_neg+FPix3_neg' 
    ]
trackingPhase1.toModify(pixelPairElectronSeedLayers, layerList = _layerListForPhase1)

from RecoTracker.TkTrackingRegions.globalTrackingRegionWithVertices_cff import globalTrackingRegionWithVertices as _globalTrackingRegionWithVertices
pixelPairElectronTrackingRegions = _globalTrackingRegionWithVertices.clone(RegionPSet = dict(
    ptMin        = 1.0,
    originRadius = 0.015,
    fixedError   = 0.03,
))
pp_on_AA.toModify(pixelPairElectronTrackingRegions, RegionPSet = dict(ptMin = 8.0))
from RecoTracker.TkHitPairs.hitPairEDProducer_cfi import hitPairEDProducer as _hitPairEDProducer
pixelPairElectronHitDoublets = _hitPairEDProducer.clone(
    seedingLayers         = 'pixelPairElectronSeedLayers',
    trackingRegions       = 'pixelPairElectronTrackingRegions',
    maxElement            = 1000000,
    produceSeedingHitSets = True,
    maxElementTotal       = 12000000,
)
from RecoTracker.TkSeedGenerator.seedCreatorFromRegionConsecutiveHitsEDProducer_cff import seedCreatorFromRegionConsecutiveHitsEDProducer as _seedCreatorFromRegionConsecutiveHitsEDProducer
pixelPairElectronSeeds = _seedCreatorFromRegionConsecutiveHitsEDProducer.clone(
    seedingHitSets = 'pixelPairElectronHitDoublets',
)

stripPairElectronSeedLayers = _mod.seedingLayersEDProducer.clone(
    layerList = ['TIB1+TIB2', 'TIB1+TID1_pos', 'TIB1+TID1_neg', 'TID2_pos+TID3_pos', 'TID2_neg+TID3_neg',
                 'TEC1_pos+TEC2_pos','TEC2_pos+TEC3_pos','TEC3_pos+TEC4_pos','TEC3_pos+TEC5_pos',
                 'TEC1_neg+TEC2_neg','TEC2_neg+TEC3_neg','TEC3_neg+TEC4_neg','TEC3_neg+TEC5_neg'],
    TIB = dict(
        TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone')),
        matchedRecHits = cms.InputTag('siStripMatchedRecHits','matchedRecHit'),
        skipClusters = cms.InputTag('tripletElectronClusterMask')
    ),
    TID = dict(
        matchedRecHits = cms.InputTag('siStripMatchedRecHits','matchedRecHit'),
        skipClusters = cms.InputTag('tripletElectronClusterMask'),
        useRingSlector = cms.bool(True),
        TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone')),
        minRing = cms.int32(1),
        maxRing = cms.int32(2)
    ),
    TEC = dict(
        matchedRecHits = cms.InputTag('siStripMatchedRecHits','matchedRecHit'),
        skipClusters = cms.InputTag('tripletElectronClusterMask'),
        useRingSlector = cms.bool(True),
        TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone')),
        minRing = cms.int32(1),
        maxRing = cms.int32(2)
    )
)

from RecoTracker.TkTrackingRegions.globalTrackingRegionFromBeamSpotFixedZ_cfi import globalTrackingRegionFromBeamSpotFixedZ as _globalTrackingRegionFromBeamSpotFixedZ
stripPairElectronTrackingRegions = _globalTrackingRegionFromBeamSpotFixedZ.clone(RegionPSet = dict(
    ptMin            = 1.0,
    originHalfLength = 12.0,
    originRadius     = 0.4,
))
pp_on_AA.toReplaceWith(stripPairElectronTrackingRegions,
    _globalTrackingRegionWithVertices.clone(
        RegionPSet = dict(
            fixedError   = 0.5,
            ptMin        = 8.0,
            originRadius = 0.4
        )
))
from RecoTracker.TkHitPairs.hitPairEDProducer_cfi import hitPairEDProducer as _hitPairEDProducer
stripPairElectronHitDoublets = _hitPairEDProducer.clone(
    seedingLayers         = 'stripPairElectronSeedLayers',
    trackingRegions       = 'stripPairElectronTrackingRegions',
    maxElement            = 1000000,
    produceSeedingHitSets = True,
    maxElementTotal       = 12000000,
)
from RecoTracker.TkSeedGenerator.seedCreatorFromRegionConsecutiveHitsEDProducer_cff import seedCreatorFromRegionConsecutiveHitsEDProducer as _seedCreatorFromRegionConsecutiveHitsEDProducer
stripPairElectronSeeds = _seedCreatorFromRegionConsecutiveHitsEDProducer.clone(
    seedingHitSets = 'stripPairElectronHitDoublets',
)


###This seed collection is produced for electron reconstruction
import RecoTracker.TkSeedGenerator.GlobalCombinedSeeds_cfi
newCombinedSeeds = RecoTracker.TkSeedGenerator.GlobalCombinedSeeds_cfi.globalCombinedSeeds.clone(
    seedCollections = ['initialStepSeeds',
                       'pixelPairStepSeeds',
                       'mixedTripletStepSeeds',
                       'pixelLessStepSeeds',
                       'tripletElectronSeeds',
                       'pixelPairElectronSeeds',
                       'stripPairElectronSeeds']
)
_seedCollections_Phase1 = [
    'initialStepSeeds',
    'highPtTripletStepSeeds',
    'mixedTripletStepSeeds',
    'pixelLessStepSeeds',
    'tripletElectronSeeds',
    'pixelPairElectronSeeds',
    'stripPairElectronSeeds',
    'lowPtTripletStepSeeds',
    'lowPtQuadStepSeeds',
    'detachedTripletStepSeeds',
    'detachedQuadStepSeeds',
    'pixelPairStepSeeds'
]
trackingPhase1.toModify(newCombinedSeeds, seedCollections = _seedCollections_Phase1)
(trackingPhase2PU140 | trackingIters01).toModify(newCombinedSeeds,
                                                 seedCollections = ['initialStepSeeds',
	                                                            'highPtTripletStepSeeds',
	                                                            'tripletElectronSeeds'] )

from Configuration.Eras.Modifier_fastSim_cff import fastSim
from FastSimulation.Tracking.ElectronSeeds_cff import _newCombinedSeeds
fastSim.toReplaceWith(newCombinedSeeds,_newCombinedSeeds.clone())

electronSeedsSeqTask = cms.Task( initialStepSeedClusterMask,
                                 pixelPairStepSeedClusterMask,
                                 mixedTripletStepSeedClusterMask,
                                 pixelLessStepSeedClusterMask,
                                 tripletElectronSeedLayers,
                                 tripletElectronTrackingRegions,
                                 tripletElectronHitDoublets,
                                 tripletElectronHitTriplets,
                                 tripletElectronSeeds,
                                 tripletElectronClusterMask,
                                 pixelPairElectronSeedLayers,
                                 pixelPairElectronTrackingRegions,
                                 pixelPairElectronHitDoublets,
                                 pixelPairElectronSeeds,
                                 stripPairElectronSeedLayers,
                                 stripPairElectronTrackingRegions,
                                 stripPairElectronHitDoublets,
                                 stripPairElectronSeeds,
                                 newCombinedSeeds)
electronSeedsSeq = cms.Sequence(electronSeedsSeqTask)
_electronSeedsSeqTask_Phase1 = electronSeedsSeqTask.copy()
_electronSeedsSeqTask_Phase1.replace(pixelPairStepSeedClusterMask, detachedTripletStepSeedClusterMask)
trackingPhase1.toReplaceWith(electronSeedsSeqTask, _electronSeedsSeqTask_Phase1 )
trackingIters01.toReplaceWith(electronSeedsSeqTask, cms.Task(
    initialStepSeedClusterMask,
    highPtTripletStepSeedClusterMask,
    tripletElectronSeedLayers,
    tripletElectronTrackingRegions,
    tripletElectronHitDoublets,
    tripletElectronHitTriplets,
    tripletElectronSeeds,
    newCombinedSeeds
))
trackingPhase2PU140.toReplaceWith(electronSeedsSeqTask, cms.Task(
    initialStepSeedClusterMask,
    highPtTripletStepSeedClusterMask,
    pixelPairStepSeedClusterMask,
    tripletElectronSeedLayers,
    tripletElectronTrackingRegions,
    tripletElectronHitDoublets,
    tripletElectronHitTriplets,
    tripletElectronSeeds,
    newCombinedSeeds
))