Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:25:16

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 
0041 ticlEMStepTask = cms.Task(ticlSeedingGlobal
0042     ,filteredLayerClustersEM
0043     ,ticlTrackstersEM)
0044 
0045 # HFNOSE CLUSTER FILTERING/MASKING
0046 
0047 filteredLayerClustersHFNoseEM = filteredLayerClustersEM.clone(
0048     LayerClusters = 'hgcalLayerClustersHFNose',
0049     LayerClustersInputMask = 'ticlTrackstersHFNoseTrkEM',
0050     min_cluster_size = 3, # inclusive
0051     algo_number = [9], # reco::CaloCluster::hfnose
0052     iteration_label = "EMn"
0053 )
0054 
0055 # HFNOSE CA - PATTERN RECOGNITION
0056 
0057 ticlTrackstersHFNoseEM = ticlTrackstersEM.clone(
0058     detector = "HFNose",
0059     layer_clusters = "hgcalLayerClustersHFNose",
0060     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0061     original_mask = "ticlTrackstersHFNoseTrkEM",
0062     filtered_mask = "filteredLayerClustersHFNoseEM:EMn",
0063     seeding_regions = "ticlSeedingGlobalHFNose",
0064     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0065     itername = "EMn",
0066     pluginPatternRecognitionByCA = dict(
0067        filter_on_categories = [0, 1],
0068        min_layers_per_trackster = 5,
0069        pid_threshold = 0.,
0070        min_cos_pointing = 0.9845, # ~10 degrees
0071        shower_start_max_layer = 4 ### inclusive
0072     )
0073 )
0074 
0075 ticlHFNoseEMStepTask = cms.Task(ticlSeedingGlobalHFNose
0076                               ,filteredLayerClustersHFNoseEM
0077                               ,ticlTrackstersHFNoseEM
0078 )