TICL/python/TrkEMStep_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingTrk, ticlSeedingTrkHFNose
0004 from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
0005 from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
0006
0007 # CLUSTER FILTERING/MASKING
0008
0009 filteredLayerClustersTrkEM = _filteredLayerClustersProducer.clone(
0010     clusterFilter = "ClusterFilterByAlgoAndSizeAndLayerRange",
0011     min_cluster_size = 3, # inclusive
0012     max_layerId = 30, # inclusive
0013     algo_number = 8,
0014     iteration_label = "TrkEM"
0015 )
0016
0017 # CA - PATTERN RECOGNITION
0018
0019 ticlTrackstersTrkEM = _trackstersProducer.clone(
0020     filtered_mask = "filteredLayerClustersTrkEM:TrkEM",
0021     seeding_regions = "ticlSeedingTrk",
0022     pluginPatternRecognitionByCA = dict(
0023         algo_verbosity = 0,
0024         filter_on_categories = [0, 1],
0025         pid_threshold = 0.5,
0026         energy_em_over_total_threshold = 0.9,
0027         max_longitudinal_sigmaPCA = 10,
0028         shower_start_max_layer = 5, #inclusive
0029         max_out_in_hops = 1,
0030         max_missing_layers_in_trackster = 2,
0031         skip_layers = 2,
0032         min_layers_per_trackster = 10,
0033         min_cos_theta = 0.97,  # ~14 degrees
0034         min_cos_pointing = 0.94, # ~20 degrees
0035         root_doublet_max_distance_from_seed_squared = 2.5e-3, # dR=0.05
0036         max_delta_time = 3.
0037     ),
0038     itername = "TrkEM",
0039 )
0040
0041 ticlTrkEMStepTask = cms.Task(ticlSeedingTrk
0042     ,filteredLayerClustersTrkEM
0043     ,ticlTrackstersTrkEM)
0044
0045 # HFNOSE CLUSTER FILTERING/MASKING
0046
0047 filteredLayerClustersHFNoseTrkEM = filteredLayerClustersTrkEM.clone(
0048     LayerClusters = 'hgcalLayerClustersHFNose',
0049     LayerClustersInputMask = "hgcalLayerClustersHFNose:InitialLayerClustersMask",
0050     min_cluster_size = 3, # inclusive
0051     algo_number = 9,
0052     iteration_label = "TrkEMn"
0053 )
0054
0055 # HFNOSE CA - PATTERN RECOGNITION
0056
0057 ticlTrackstersHFNoseTrkEM = ticlTrackstersTrkEM.clone(
0058     detector = "HFNose",
0059     layer_clusters = "hgcalLayerClustersHFNose",
0060     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0061     original_mask = "hgcalLayerClustersHFNose:InitialLayerClustersMask",
0062     filtered_mask = "filteredLayerClustersHFNoseTrkEM:TrkEMn",
0063     seeding_regions = "ticlSeedingTrkHFNose",
0064     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0065     itername = "TrkEMn",
0066     pluginPatternRecognitionByCA = dict(
0067         filter_on_categories = [0, 1],
0068         min_layers_per_trackster = 5,
0069         pid_threshold = 0.,
0070         shower_start_max_layer = 5 #inclusive
0071     )
0072 )
0073
0074 ticlHFNoseTrkEMStepTask = cms.Task(ticlSeedingTrkHFNose
0075     ,filteredLayerClustersHFNoseTrkEM
0076     ,ticlTrackstersHFNoseTrkEM)
0077