ticl_iterations Namespace Reference

Functions
def	TICL_iterations
	TICL_iterations: to be run with local HGCAL reco only i.e. More...
def	TICL_iterations_withReco
	withReco: requires full reco of the event to run this part i.e. More...

Function Documentation

def ticl_iterations.TICL_iterations

(

process

)

TICL_iterations: to be run with local HGCAL reco only i.e.

collections of generalTracks (track-seeded iteration) NOT available

Definition at line 132 of file ticl_iterations.py.

def TICL_iterations(process):
  process.ticlLayerTileProducer = ticlLayerTileProducer.clone()

  process.ticlSeedingGlobal = ticlSeedingRegionProducer.clone(
    algoId = 2
  )

  process.filteredLayerClustersMIP = filteredLayerClustersProducer.clone(
      clusterFilter = "ClusterFilterBySize",
      algo_number = 8,
      max_cluster_size = 2, # inclusive
      iteration_label = "MIP"
  )

  process.trackstersMIP = trackstersProducer.clone(
      filtered_mask = "filteredLayerClustersMIP:MIP",
      seeding_regions = "ticlSeedingGlobal",
      skip_layers = 3,
      min_layers_per_trackster = 15,
      min_cos_theta = 0.99, # ~10 degrees
  )

  process.filteredLayerClusters = filteredLayerClustersProducer.clone(
      clusterFilter = "ClusterFilterByAlgoAndSize",
      min_cluster_size = 2,
      algo_number = 8,
      iteration_label = "algo8"
  )

  process.tracksters = trackstersProducer.clone(
      original_mask = "trackstersMIP",
      filtered_mask = "filteredLayerClusters:algo8",
      seeding_regions = "ticlSeedingGlobal",
      skip_layers = 2,
      min_layers_per_trackster = 15,
      min_cos_theta = 0.94, # ~20 degrees
      min_cos_pointing = 0.7
  )

  process.HGCalUncalibRecHit = HGCalUncalibRecHit
  process.HGCalRecHit = HGCalRecHit
  process.hgcalLayerClusters = hgcalLayerClusters
  process.hgcalMultiClusters = hgcalMultiClusters
  process.TICL_Task = cms.Task(process.HGCalUncalibRecHit,
      process.HGCalRecHit,
      process.hgcalLayerClusters,
      process.filteredLayerClustersMIP,
      process.ticlLayerTileProducer,
      process.ticlSeedingGlobal,
      process.trackstersMIP,
      process.filteredLayerClusters,
      process.tracksters,
      process.hgcalMultiClusters)
  process.schedule = cms.Schedule(process.raw2digi_step,process.FEVTDEBUGHLToutput_step)
  process.schedule.associate(process.TICL_Task)
  return process

def ticl_iterations.TICL_iterations_withReco

(

process

)

withReco: requires full reco of the event to run this part i.e.

collections of generalTracks can be accessed

Definition at line 23 of file ticl_iterations.py.

def TICL_iterations_withReco(process):
  process.FEVTDEBUGHLTEventContent.outputCommands.extend([
    'keep *_ticlCandidateFromTrackstersProducer*_*_*',
    'keep *_pfTICLProducer*_*_*'])

  process.ticlLayerTileProducer = ticlLayerTileProducer.clone()

  process.ticlSeedingTrk = ticlSeedingRegionProducer.clone(
    algoId = 1
  )

  process.filteredLayerClustersTrk = filteredLayerClustersProducer.clone(
    clusterFilter = "ClusterFilterByAlgo",
    algo_number = 8,
    iteration_label = "Trk"
  )

  process.trackstersTrk = trackstersProducer.clone(
    filtered_mask = "filteredLayerClustersTrk:Trk",
    seeding_regions = "ticlSeedingTrk",
    skip_layers = 3,
    min_layers_per_trackster = 5,
    min_cos_theta = 0.99, # ~10 degrees                                              
    min_cos_pointing = 0.9
  )

  process.ticlSeedingGlobal = ticlSeedingRegionProducer.clone(
    algoId = 2
  )

  process.filteredLayerClustersMIP = filteredLayerClustersProducer.clone(
      clusterFilter = "ClusterFilterBySize",
      algo_number = 8,
      max_cluster_size = 2, # inclusive
      iteration_label = "MIP"
  )

  process.trackstersMIP = trackstersProducer.clone(
      filtered_mask = "filteredLayerClustersMIP:MIP",
      seeding_regions = "ticlSeedingGlobal",
      skip_layers = 3,
      min_layers_per_trackster = 15,
      min_cos_theta = 0.99, # ~10 degrees
      min_cos_pointing = 0.9,
      out_in_dfs = False,
  )

  process.filteredLayerClusters = filteredLayerClustersProducer.clone(
      clusterFilter = "ClusterFilterByAlgoAndSize",
      min_cluster_size = 2,
      algo_number = 8,
      iteration_label = "algo8",
      LayerClustersInputMask = "trackstersMIP"
  )

  process.trackstersEM = trackstersProducer.clone(
      max_out_in_hops = 4,
      original_mask = "trackstersMIP",
      filtered_mask = "filteredLayerClusters:algo8",
      seeding_regions = "ticlSeedingGlobal",
      skip_layers = 1,
      min_layers_per_trackster = 10,
      min_cos_theta = 0.984, # ~10 degrees
      min_cos_pointing = 0.9 # ~26 degrees
  )

  process.trackstersHAD = trackstersProducer.clone(
      filtered_mask = "filteredLayerClusters:algo8",
      seeding_regions = "ticlSeedingGlobal",
      skip_layers = 2,
      min_layers_per_trackster = 10,
      min_cos_theta = 0.8, 
      min_cos_pointing = 0.7
  )

  process.ticlCandidateFromTrackstersProducer = ticlCandidateFromTrackstersProducer.clone()

  process.pfTICLProducer = pfTICLProducer.clone()

  process.hgcalMultiClusters = hgcalMultiClusters
  process.TICL_Task = cms.Task(
      process.ticlLayerTileProducer,
      process.ticlSeedingTrk,
      process.filteredLayerClustersTrk,
      process.trackstersTrk,
      process.ticlSeedingGlobal,
      process.filteredLayerClustersMIP,
      process.trackstersMIP,
      process.filteredLayerClusters,
      process.trackstersEM,
      process.trackstersHAD,
      process.ticlCandidateFromTrackstersProducer,
      process.pfTICLProducer)

  process.schedule.associate(process.TICL_Task)

  process.ticlPFValidation = ticlPFValidation
  process.hgcalValidation.insert(-1, process.ticlPFValidation)
  
  if getattr(process,'hgcalValidator'):
    process.hgcalValidator.label_lcl = "hgcalLayerClusters"
    process.hgcalValidator.label_mcl = ["multiClustersFromTrackstersEM:MultiClustersFromTracksterByCA", "multiClustersFromTrackstersHAD:MultiClustersFromTracksterByCA"]
    process.hgcalValidator.domulticlustersPlots = True
    
  return process