CMS 3D CMS Logo

EMStep_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
4 from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
5 from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
6 from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer
7 
8 # CLUSTER FILTERING/MASKING
9 
10 filteredLayerClustersEM = _filteredLayerClustersProducer.clone(
11  clusterFilter = "ClusterFilterByAlgoAndSizeAndLayerRange",
12  min_cluster_size = 3, # inclusive
13  max_layerId = 30, # inclusive
14  algo_number = 8,
15  LayerClustersInputMask = 'ticlTrackstersTrkEM',
16  iteration_label = "EM"
17 )
18 
19 # CA - PATTERN RECOGNITION
20 
21 ticlTrackstersEM = _trackstersProducer.clone(
22  filtered_mask = "filteredLayerClustersEM:EM",
23  original_mask = 'ticlTrackstersTrkEM',
24  seeding_regions = "ticlSeedingGlobal",
25  filter_on_categories = [0, 1],
26  pid_threshold = 0.5,
27  energy_em_over_total_threshold = 0.9,
28  max_longitudinal_sigmaPCA = 10,
29  shower_start_max_layer = 5, #inclusive
30  max_out_in_hops = 1,
31  skip_layers = 2,
32  max_missing_layers_in_trackster = 1,
33  min_layers_per_trackster = 10,
34  min_cos_theta = 0.97, # ~14 degrees
35  min_cos_pointing = 0.9, # ~25 degrees
36  max_delta_time = 3.,
37  itername = "EM",
38  algo_verbosity = 0,
39 )
40 
41 # MULTICLUSTERS
42 
43 ticlMultiClustersFromTrackstersEM = _multiClustersFromTrackstersProducer.clone(
44  Tracksters = "ticlTrackstersEM"
45 )
46 
47 ticlEMStepTask = cms.Task(ticlSeedingGlobal
48  ,filteredLayerClustersEM
49  ,ticlTrackstersEM
50  ,ticlMultiClustersFromTrackstersEM)
51 
52 filteredLayerClustersHFNoseEM = filteredLayerClustersEM.clone(
53  LayerClusters = 'hgcalLayerClustersHFNose',
54  LayerClustersInputMask = "hgcalLayerClustersHFNose:InitialLayerClustersMask",
55  iteration_label = "EMn",
56  algo_number = 9
57 #no tracking mask for EM for now
58 )
59 
60 ticlTrackstersHFNoseEM = ticlTrackstersEM.clone(
61  detector = "HFNose",
62  layer_clusters = "hgcalLayerClustersHFNose",
63  layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
64  original_mask = "hgcalLayerClustersHFNose:InitialLayerClustersMask",
65  filtered_mask = "filteredLayerClustersHFNoseEM:EMn",
66  seeding_regions = "ticlSeedingGlobalHFNose",
67  time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
68  min_layers_per_trackster = 6
69 )
70 
71 ticlHFNoseEMStepTask = cms.Task(ticlSeedingGlobalHFNose
72  ,filteredLayerClustersHFNoseEM
73  ,ticlTrackstersHFNoseEM
74 )