SpecialSeedGenerators/python/SimpleCosmicBONSeeder_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002 import RecoTracker.SpecialSeedGenerators.CombinatorialSeedGeneratorForCosmics_cfi
0003 def makeSimpleCosmicSeedLayers(*layers):
0004     layerList = cms.vstring()
0005     if 'ALL' in layers:
0006         layers = [ 'TOB', 'TEC', 'TOBTEC', 'TECSKIP' ]
0007     if 'TOB' in layers:
0008         layerList += ['MTOB4+MTOB5+MTOB6',
0009                       'MTOB3+MTOB5+MTOB6',
0010                       'MTOB3+MTOB4+MTOB5',
0011                       'MTOB3+MTOB4+MTOB6',
0012                       'TOB2+MTOB4+MTOB5',
0013                       'TOB2+MTOB3+MTOB5']
0014     if 'TEC' in layers:
0015         TECwheelTriplets = [ (i,i+1,i+2) for i in range(7,0,-1)]
0016         layerList += [ 'TEC%d_pos+TEC%d_pos+TEC%d_pos' % ls for ls in TECwheelTriplets ]
0017         layerList += [ 'TEC%d_neg+TEC%d_neg+TEC%d_neg' % ls for ls in TECwheelTriplets ]
0018     if 'TECSKIP' in layers:
0019         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)  ]
0020         layerList += [ 'TEC%d_pos+TEC%d_pos+TEC%d_pos' % ls for ls in TECwheelTriplets ]
0021         layerList += [ 'TEC%d_neg+TEC%d_neg+TEC%d_neg' % ls for ls in TECwheelTriplets ]
0022     if 'TOBTEC' in layers:
0023         layerList += [ 'MTOB6+TEC1_pos+TEC2_pos',
0024                        'MTOB6+TEC1_neg+TEC2_neg',
0025                        'MTOB6+MTOB5+TEC1_pos',
0026                        'MTOB6+MTOB5+TEC1_neg' ]
0027     #print "SEEDING LAYER LIST = ", layerList
0028     return layerList
0029
0030 layerInfo = RecoTracker.SpecialSeedGenerators.CombinatorialSeedGeneratorForCosmics_cfi.layerInfo.clone(
0031     TEC = dict(useSimpleRphiHitsCleaner = False)
0032 )
0033 layerList = makeSimpleCosmicSeedLayers('ALL'),
0034
0035
0036 from RecoTracker.SpecialSeedGenerators.simpleCosmicBONSeeder_cfi import simpleCosmicBONSeeder
0037 simpleCosmicBONSeeds = simpleCosmicBONSeeder.clone(
0038     TTRHBuilder = 'WithTrackAngle',
0039     ClusterCheckPSet = cms.PSet(
0040         doClusterCheck = cms.bool(True),
0041         MaxNumberOfStripClusters = cms.uint32(1000),
0042         ClusterCollectionLabel = cms.InputTag("siStripClusters"),
0043         DontCountDetsAboveNClusters = cms.uint32(20),  # if N > 0, ignore in total the dets with more than N clusters
0044         MaxNumberOfPixelClusters = cms.uint32(1000),
0045         PixelClusterCollectionLabel = cms.InputTag("siPixelClusters")
0046     ),
0047     maxTriplets = 50000,
0048     maxSeeds    = 20000,
0049     RegionPSet = cms.PSet(
0050         originZPosition  = cms.double(0.0),    # \    These three parameters
0051         originRadius     = cms.double(150.0),  #  |-> probably don't change
0052         originHalfLength = cms.double(90.0),   # /    anything at all.
0053         ptMin = cms.double(0.5),               # pt cut, applied both at the triplet finding and at the seeding level
0054         pMin  = cms.double(1.0),               # p  cut, applied only at the seeding level
0055     ),
0056     TripletsSrc = "simpleCosmicBONSeedingLayers",
0057     TripletsDebugLevel = 0,  # debug triplet finding (0 to 3)
0058     seedOnMiddle    = False, # after finding the triplet, add only two hits to the seed
0059     rescaleError    = 1.0, # we don't need it anymore. At least for runs with BON
0060
0061     ClusterChargeCheck = cms.PSet(
0062         checkCharge                 = cms.bool(False), # Apply cuts on cluster charge
0063         matchedRecHitsUseAnd        = cms.bool(True), # Both clusters in the pair should pass the charge cut
0064         Thresholds  = cms.PSet( # Uncorrected thresholds
0065             TIB = cms.int32(0), #
0066             TID = cms.int32(0), # Currenlty not used
0067             TOB = cms.int32(0), #
0068             TEC = cms.int32(0), #
0069         ),
0070     ),
0071     HitsPerModuleCheck = cms.PSet(
0072         checkHitsPerModule = cms.bool(True), # Apply cuts on the number of hits per module
0073         Thresholds  = cms.PSet( #
0074             TIB = cms.int32(20), #
0075             TID = cms.int32(20), # FIXME: to be optimized
0076             TOB = cms.int32(20), #
0077             TEC = cms.int32(20), #
0078         ),
0079     ),
0080
0081     minimumGoodHitsInSeed = 3,   # NO bad hits in the seed (set to '2' to allow one bad hit in the seed)
0082
0083     writeTriplets   = False, # write the triplets to the Event as OwnVector<TrackingRecHit>
0084     helixDebugLevel = 0, # debug FastHelix (0 to 2)
0085     seedDebugLevel  = 0, # debug seed building (0 to 3)
0086     #***top-bottom
0087     PositiveYOnly = False,
0088     NegativeYOnly = False
0089     #***
0090 )
0091