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     LayerClustersInputMask = 'ticlTrackstersTrkEM',
0014     iteration_label = "EM"
0015 )
0016
0017 # CA - PATTERN RECOGNITION
0018
0019 ticlTrackstersEM = _trackstersProducer.clone(
0020     filtered_mask = "filteredLayerClustersEM:EM",
0021     original_mask = 'ticlTrackstersTrkEM',
0022     seeding_regions = "ticlSeedingGlobal",
0023     pluginPatternRecognitionByCA = dict(
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         skip_layers = 2,
0031         max_missing_layers_in_trackster = 1,
0032         min_layers_per_trackster = 10,
0033         min_cos_theta = 0.97,  # ~14 degrees
0034         min_cos_pointing = 0.9, # ~25 degrees
0035         max_delta_time = 3.,
0036         algo_verbosity = 0
0037     ),
0038     itername = "EM"
0039 )
0040 from Configuration.ProcessModifiers.ticl_v5_cff import ticl_v5
0041 ticl_v5.toModify(ticlTrackstersEM.pluginPatternRecognitionByCA, computeLocalTime = cms.bool(True))
0042
0043 ticlEMStepTask = cms.Task(ticlSeedingGlobal
0044     ,filteredLayerClustersEM
0045     ,ticlTrackstersEM)
0046
0047 # HFNOSE CLUSTER FILTERING/MASKING
0048
0049 filteredLayerClustersHFNoseEM = filteredLayerClustersEM.clone(
0050     LayerClusters = 'hgcalLayerClustersHFNose',
0051     LayerClustersInputMask = 'ticlTrackstersHFNoseTrkEM',
0052     min_cluster_size = 3, # inclusive
0053     algo_number = [9], # reco::CaloCluster::hfnose
0054     iteration_label = "EMn"
0055 )
0056
0057 # HFNOSE CA - PATTERN RECOGNITION
0058
0059 ticlTrackstersHFNoseEM = ticlTrackstersEM.clone(
0060     detector = "HFNose",
0061     layer_clusters = "hgcalLayerClustersHFNose",
0062     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0063     original_mask = "ticlTrackstersHFNoseTrkEM",
0064     filtered_mask = "filteredLayerClustersHFNoseEM:EMn",
0065     seeding_regions = "ticlSeedingGlobalHFNose",
0066     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0067     itername = "EMn",
0068     pluginPatternRecognitionByCA = dict(
0069        filter_on_categories = [0, 1],
0070        min_layers_per_trackster = 5,
0071        pid_threshold = 0.,
0072        min_cos_pointing = 0.9845, # ~10 degrees
0073        shower_start_max_layer = 4 ### inclusive
0074     )
0075 )
0076 ticl_v5.toModify(ticlTrackstersHFNoseEM.pluginPatternRecognitionByCA, computeLocalTime = cms.bool(True))
0077
0078 ticlHFNoseEMStepTask = cms.Task(ticlSeedingGlobalHFNose
0079                               ,filteredLayerClustersHFNoseEM
0080                               ,ticlTrackstersHFNoseEM
0081 )