TICL/python/EMStep_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
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 filteredLayerClustersEM = _filteredLayerClustersProducer.clone(
0010     clusterFilter = "ClusterFilterByAlgoAndSizeAndLayerRange",
0011     min_cluster_size = 3, # inclusive
0012     max_layerId = 30, # inclusive
0013     algo_number = 8,
0014     LayerClustersInputMask = 'ticlTrackstersTrkEM',
0015     iteration_label = "EM"
0016 )
0017
0018 # CA - PATTERN RECOGNITION
0019
0020 ticlTrackstersEM = _trackstersProducer.clone(
0021     filtered_mask = "filteredLayerClustersEM:EM",
0022     original_mask = 'ticlTrackstersTrkEM',
0023     seeding_regions = "ticlSeedingGlobal",
0024     pluginPatternRecognitionByCA = dict(
0025         filter_on_categories = [0, 1],
0026         pid_threshold = 0.5,
0027         energy_em_over_total_threshold = 0.9,
0028         max_longitudinal_sigmaPCA = 10,
0029         shower_start_max_layer = 5, #inclusive
0030         max_out_in_hops = 1,
0031         skip_layers = 2,
0032         max_missing_layers_in_trackster = 1,
0033         min_layers_per_trackster = 10,
0034         min_cos_theta = 0.97,  # ~14 degrees
0035         min_cos_pointing = 0.9, # ~25 degrees
0036         max_delta_time = 3.,
0037         algo_verbosity = 0
0038     ),
0039     itername = "EM"
0040 )
0041
0042 ticlEMStepTask = cms.Task(ticlSeedingGlobal
0043     ,filteredLayerClustersEM
0044     ,ticlTrackstersEM)
0045
0046 # HFNOSE CLUSTER FILTERING/MASKING
0047
0048 filteredLayerClustersHFNoseEM = filteredLayerClustersEM.clone(
0049     LayerClusters = 'hgcalLayerClustersHFNose',
0050     LayerClustersInputMask = 'ticlTrackstersHFNoseTrkEM',
0051     min_cluster_size = 3, # inclusive
0052     algo_number = 9,
0053     iteration_label = "EMn"
0054 )
0055
0056 # HFNOSE CA - PATTERN RECOGNITION
0057
0058 ticlTrackstersHFNoseEM = ticlTrackstersEM.clone(
0059     detector = "HFNose",
0060     layer_clusters = "hgcalLayerClustersHFNose",
0061     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0062     original_mask = "ticlTrackstersHFNoseTrkEM",
0063     filtered_mask = "filteredLayerClustersHFNoseEM:EMn",
0064     seeding_regions = "ticlSeedingGlobalHFNose",
0065     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0066     itername = "EMn",
0067     pluginPatternRecognitionByCA = dict(
0068        filter_on_categories = [0, 1],
0069        min_layers_per_trackster = 5,
0070        pid_threshold = 0.,
0071        min_cos_pointing = 0.9845, # ~10 degrees
0072        shower_start_max_layer = 4 ### inclusive
0073     )
0074 )
0075
0076 ticlHFNoseEMStepTask = cms.Task(ticlSeedingGlobalHFNose
0077                               ,filteredLayerClustersHFNoseEM
0078                               ,ticlTrackstersHFNoseEM
0079 )