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     iteration_label = "TrkEM"
0014 )
0015
0016 # CA - PATTERN RECOGNITION
0017
0018 ticlTrackstersTrkEM = _trackstersProducer.clone(
0019     filtered_mask = "filteredLayerClustersTrkEM:TrkEM",
0020     seeding_regions = "ticlSeedingTrk",
0021     pluginPatternRecognitionByCA = dict(
0022         algo_verbosity = 0,
0023         filter_on_categories = [0, 1],
0024         pid_threshold = 0.5,
0025         energy_em_over_total_threshold = 0.9,
0026         max_longitudinal_sigmaPCA = 10,
0027         shower_start_max_layer = 5, #inclusive
0028         max_out_in_hops = 1,
0029         max_missing_layers_in_trackster = 2,
0030         skip_layers = 2,
0031         min_layers_per_trackster = 10,
0032         min_cos_theta = 0.97,  # ~14 degrees
0033         min_cos_pointing = 0.94, # ~20 degrees
0034         root_doublet_max_distance_from_seed_squared = 2.5e-3, # dR=0.05
0035         max_delta_time = 3.
0036     ),
0037     itername = "TrkEM",
0038 )
0039
0040 from Configuration.ProcessModifiers.ticl_v5_cff import ticl_v5
0041 ticl_v5.toModify(ticlTrackstersTrkEM.pluginPatternRecognitionByCA, computeLocalTime = cms.bool(True))
0042
0043 ticlTrkEMStepTask = cms.Task(ticlSeedingTrk
0044     ,filteredLayerClustersTrkEM
0045     ,ticlTrackstersTrkEM)
0046
0047 # HFNOSE CLUSTER FILTERING/MASKING
0048
0049 filteredLayerClustersHFNoseTrkEM = filteredLayerClustersTrkEM.clone(
0050     LayerClusters = 'hgcalLayerClustersHFNose',
0051     LayerClustersInputMask = "hgcalLayerClustersHFNose:InitialLayerClustersMask",
0052     min_cluster_size = 3, # inclusive
0053     algo_number = [9], # reco::CaloCluster::hfnose
0054     iteration_label = "TrkEMn"
0055 )
0056
0057 # HFNOSE CA - PATTERN RECOGNITION
0058
0059 ticlTrackstersHFNoseTrkEM = ticlTrackstersTrkEM.clone(
0060     detector = "HFNose",
0061     layer_clusters = "hgcalLayerClustersHFNose",
0062     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0063     original_mask = "hgcalLayerClustersHFNose:InitialLayerClustersMask",
0064     filtered_mask = "filteredLayerClustersHFNoseTrkEM:TrkEMn",
0065     seeding_regions = "ticlSeedingTrkHFNose",
0066     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0067     itername = "TrkEMn",
0068     pluginPatternRecognitionByCA = dict(
0069         filter_on_categories = [0, 1],
0070         min_layers_per_trackster = 5,
0071         pid_threshold = 0.,
0072         shower_start_max_layer = 5 #inclusive
0073     )
0074 )
0075
0076 ticl_v5.toModify(ticlTrackstersHFNoseTrkEM.pluginPatternRecognitionByCA, computeLocalTime = cms.bool(True))
0077
0078 ticlHFNoseTrkEMStepTask = cms.Task(ticlSeedingTrkHFNose
0079     ,filteredLayerClustersHFNoseTrkEM
0080     ,ticlTrackstersHFNoseTrkEM)
0081