RecoTracker/SpecialSeedGenerators/python/SimpleCosmicBONSeeder_cfi.py - master - githubmirror/cmssw

RecoTracker/SpecialSeedGenerators/python/SimpleCosmicBONSeeder_cfi.py
90 строк4 KB
jlidrych
increase threshold of max number of strips clusters for cosmics
13 мар 2025, 19:45
13 мар 2025, 19:45d77e23c
import FWCore.ParameterSet.Config as cms
import RecoTracker.SpecialSeedGenerators.CombinatorialSeedGeneratorForCosmics_cfi
def makeSimpleCosmicSeedLayers(*layers):
    layerList = cms.vstring()
    if 'ALL' in layers: 
        layers = [ 'TOB', 'TEC', 'TOBTEC', 'TECSKIP' ]
    if 'TOB' in layers:
        layerList += ['MTOB4+MTOB5+MTOB6',
                      'MTOB3+MTOB5+MTOB6',
                      'MTOB3+MTOB4+MTOB5',
                      'MTOB3+MTOB4+MTOB6',
                      'TOB2+MTOB4+MTOB5',
                      'TOB2+MTOB3+MTOB5']
    if 'TEC' in layers:
        TECwheelTriplets = [ (i,i+1,i+2) for i in range(7,0,-1)]
        layerList += [ 'TEC%d_pos+TEC%d_pos+TEC%d_pos' % ls for ls in TECwheelTriplets ]
        layerList += [ 'TEC%d_neg+TEC%d_neg+TEC%d_neg' % ls for ls in TECwheelTriplets ]
    if 'TECSKIP' in layers:
        TECwheelTriplets = [ (i-1,i+1,i+2) for i in range(7,1,-1)] + [ (i-1,i,i+2) for i in range(7,1,-1)  ]
        layerList += [ 'TEC%d_pos+TEC%d_pos+TEC%d_pos' % ls for ls in TECwheelTriplets ]
        layerList += [ 'TEC%d_neg+TEC%d_neg+TEC%d_neg' % ls for ls in TECwheelTriplets ]
    if 'TOBTEC' in layers:
        layerList += [ 'MTOB6+TEC1_pos+TEC2_pos',
                       'MTOB6+TEC1_neg+TEC2_neg',
                       'MTOB6+MTOB5+TEC1_pos',
                       'MTOB6+MTOB5+TEC1_neg' ]
    #print "SEEDING LAYER LIST = ", layerList
    return layerList

layerInfo = RecoTracker.SpecialSeedGenerators.CombinatorialSeedGeneratorForCosmics_cfi.layerInfo.clone(
    TEC = dict(useSimpleRphiHitsCleaner = False)
)
layerList = makeSimpleCosmicSeedLayers('ALL'),


from RecoTracker.SpecialSeedGenerators.simpleCosmicBONSeeder_cfi import simpleCosmicBONSeeder 
simpleCosmicBONSeeds = simpleCosmicBONSeeder.clone(
    TTRHBuilder = 'WithTrackAngle',
    ClusterCheckPSet = cms.PSet(
        doClusterCheck = cms.bool(True),
        MaxNumberOfStripClusters = cms.uint32(1000),
        ClusterCollectionLabel = cms.InputTag("siStripClusters"),
        DontCountDetsAboveNClusters = cms.uint32(20),  # if N > 0, ignore in total the dets with more than N clusters
        MaxNumberOfPixelClusters = cms.uint32(1000),
        PixelClusterCollectionLabel = cms.InputTag("siPixelClusters")
    ),
    maxTriplets = 50000,
    maxSeeds    = 20000,
    RegionPSet = cms.PSet(
        originZPosition  = cms.double(0.0),    # \    These three parameters
        originRadius     = cms.double(150.0),  #  |-> probably don't change
        originHalfLength = cms.double(90.0),   # /    anything at all.
        ptMin = cms.double(0.5),               # pt cut, applied both at the triplet finding and at the seeding level
        pMin  = cms.double(1.0),               # p  cut, applied only at the seeding level
    ),
    TripletsSrc = "simpleCosmicBONSeedingLayers",
    TripletsDebugLevel = 0,  # debug triplet finding (0 to 3)
    seedOnMiddle    = False, # after finding the triplet, add only two hits to the seed
    rescaleError    = 1.0, # we don't need it anymore. At least for runs with BON

    ClusterChargeCheck = cms.PSet(
        checkCharge                 = cms.bool(False), # Apply cuts on cluster charge
        matchedRecHitsUseAnd        = cms.bool(True), # Both clusters in the pair should pass the charge cut
        Thresholds  = cms.PSet( # Uncorrected thresholds
            TIB = cms.int32(0), #
            TID = cms.int32(0), # Currenlty not used
            TOB = cms.int32(0), # 
            TEC = cms.int32(0), #
        ),
    ),
    HitsPerModuleCheck = cms.PSet(
        checkHitsPerModule = cms.bool(True), # Apply cuts on the number of hits per module 
        Thresholds  = cms.PSet( # 
            TIB = cms.int32(20), #
            TID = cms.int32(20), # FIXME: to be optimized
            TOB = cms.int32(20), # 
            TEC = cms.int32(20), #
        ),
    ),
    
    minimumGoodHitsInSeed = 3,   # NO bad hits in the seed (set to '2' to allow one bad hit in the seed)
                                      
    writeTriplets   = False, # write the triplets to the Event as OwnVector<TrackingRecHit>
    helixDebugLevel = 0, # debug FastHelix (0 to 2)
    seedDebugLevel  = 0, # debug seed building (0 to 3)
    #***top-bottom
    PositiveYOnly = False,
    NegativeYOnly = False
    #***
)