Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-10-25 09:59:32

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
0004 from RecoHGCal.TICL.ticlLayerTileProducer_cfi import ticlLayerTileProducer as _ticlLayerTileProducer
0005 from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
0006 from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
0007 
0008 # CLUSTER FILTERING/MASKING
0009 
0010 filteredLayerClustersHAD = _filteredLayerClustersProducer.clone(
0011     clusterFilter = "ClusterFilterByAlgoAndSize",
0012     min_cluster_size = 3, # inclusive
0013     iteration_label = "HAD",
0014     LayerClustersInputMask = "ticlTrackstersTrk"
0015 )
0016 
0017 # CA - PATTERN RECOGNITION
0018 
0019 ticlTrackstersHAD = _trackstersProducer.clone(
0020     filtered_mask = "filteredLayerClustersHAD:HAD",
0021     original_mask = 'ticlTrackstersTrk',
0022     seeding_regions = "ticlSeedingGlobal",
0023     # For the moment we mask everything w/o requirements since we are last
0024 #    filter_on_categories = [5], # filter neutral hadrons
0025 #    pid_threshold = 0.7,
0026     pluginPatternRecognitionByCA = dict(
0027        skip_layers = 1,
0028        min_layers_per_trackster = 12,
0029        min_cos_theta = 0.866,    # ~30 degrees
0030        min_cos_pointing = 0.819, # ~35 degrees
0031        max_delta_time = -1
0032     ),
0033     itername = "HAD"
0034 )
0035 
0036 ticlHADStepTask = cms.Task(ticlSeedingGlobal
0037     ,filteredLayerClustersHAD
0038     ,ticlTrackstersHAD)
0039 
0040 # HFNOSE CLUSTER FILTERING/MASKING
0041 
0042 filteredLayerClustersHFNoseHAD = filteredLayerClustersHAD.clone(
0043     min_cluster_size = 2, # inclusive
0044     algo_number = [9], # reco::CaloCluster::hfnose
0045     iteration_label = "HADn",
0046     LayerClusters = 'hgcalLayerClustersHFNose',
0047     LayerClustersInputMask = "ticlTrackstersHFNoseTrk"
0048 )
0049 
0050 # HFNOSE CA - PATTERN RECOGNITION
0051 
0052 ticlTrackstersHFNoseHAD = ticlTrackstersHAD.clone(
0053     detector = "HFNose",
0054     layer_clusters = "hgcalLayerClustersHFNose",
0055     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0056     original_mask = "ticlTrackstersHFNoseTrk",
0057     filtered_mask = "filteredLayerClustersHFNoseHAD:HADn",
0058     seeding_regions = "ticlSeedingGlobalHFNose",
0059     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0060     pluginPatternRecognitionByCA = dict(
0061        pid_threshold = 0.,
0062        skip_layers = 1,
0063        min_layers_per_trackster = 5,
0064        min_cos_theta = 0.866,    # ~30 degrees
0065        min_cos_pointing = 0.866, # ~30 degrees
0066        max_delta_time = -1
0067     ),
0068     itername = "HADn"
0069     )
0070 
0071 ticlHFNoseHADStepTask = cms.Task(ticlSeedingGlobalHFNose
0072                                  ,filteredLayerClustersHFNoseHAD
0073                                  ,ticlTrackstersHFNoseHAD)