TICL/python/HADStep_cff.py

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 from Configuration.ProcessModifiers.ticl_v5_cff import ticl_v5
0037 ticl_v5.toModify(ticlTrackstersHAD.pluginPatternRecognitionByCA, computeLocalTime = cms.bool(True))
0038
0039 ticlHADStepTask = cms.Task(ticlSeedingGlobal
0040     ,filteredLayerClustersHAD
0041     ,ticlTrackstersHAD)
0042
0043 # HFNOSE CLUSTER FILTERING/MASKING
0044
0045 filteredLayerClustersHFNoseHAD = filteredLayerClustersHAD.clone(
0046     min_cluster_size = 2, # inclusive
0047     algo_number = [9], # reco::CaloCluster::hfnose
0048     iteration_label = "HADn",
0049     LayerClusters = 'hgcalLayerClustersHFNose',
0050     LayerClustersInputMask = "ticlTrackstersHFNoseTrk"
0051 )
0052
0053 # HFNOSE CA - PATTERN RECOGNITION
0054
0055 ticlTrackstersHFNoseHAD = ticlTrackstersHAD.clone(
0056     detector = "HFNose",
0057     layer_clusters = "hgcalLayerClustersHFNose",
0058     layer_clusters_hfnose_tiles = "ticlLayerTileHFNose",
0059     original_mask = "ticlTrackstersHFNoseTrk",
0060     filtered_mask = "filteredLayerClustersHFNoseHAD:HADn",
0061     seeding_regions = "ticlSeedingGlobalHFNose",
0062     time_layerclusters = "hgcalLayerClustersHFNose:timeLayerCluster",
0063     pluginPatternRecognitionByCA = dict(
0064        pid_threshold = 0.,
0065        skip_layers = 1,
0066        min_layers_per_trackster = 5,
0067        min_cos_theta = 0.866,    # ~30 degrees
0068        min_cos_pointing = 0.866, # ~30 degrees
0069        max_delta_time = -1
0070     ),
0071     itername = "HADn"
0072     )
0073
0074 ticl_v5.toModify(ticlTrackstersHFNoseHAD.pluginPatternRecognitionByCA, computeLocalTime = cms.bool(True))
0075
0076 ticlHFNoseHADStepTask = cms.Task(ticlSeedingGlobalHFNose
0077                                  ,filteredLayerClustersHFNoseHAD
0078                                  ,ticlTrackstersHFNoseHAD)