CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
fftjetcommon_cfi.py
Go to the documentation of this file.
1 import math
2 import FWCore.ParameterSet.Config as cms
4 
5 # Common definitions for FFTJet interface modules
6 
7 # Some useful constants
8 fftjet_large_int = pow(2,31) - 1
9 
10 # Global phi-to-eta bandwidth ratio
11 fftjet_phi_to_eta_bw_ratio = 1.0
12 
13 # Are we including the complete event at the lowest scale
14 # of the clustering tree?
15 fftjet_insert_complete_event = False
16 fftjet_complete_event_scale = 0.05
17 
18 # The standard eta range for various algos
19 fftjet_standard_eta_range = 5.2
20 
21 # Provide several feasible energy discretization grid configurations.
22 # The width of eta bins is 0.087. We want to choose the binning
23 # so that FFT is efficient. At the same time, we have to prevent
24 # wrap-around energy leakage during convolutions from one eta side
25 # to another. Note that the CMS calorimeter extends up to eta of 5.191.
26 fftjet_grid_256_72 = cms.PSet(
27  nEtaBins = cms.uint32(256),
28  etaMin = cms.double(-11.136),
29  etaMax = cms.double(11.136),
30  nPhiBins = cms.uint32(72),
31  phiBin0Edge = cms.double(0.0),
32  title = cms.untracked.string("256 x 72")
33 )
34 
35 fftjet_grid_192_72 = cms.PSet(
36  nEtaBins = cms.uint32(192),
37  etaMin = cms.double(-8.352),
38  etaMax = cms.double(8.352),
39  nPhiBins = cms.uint32(72),
40  phiBin0Edge = cms.double(0.0),
41  title = cms.untracked.string("192 x 72")
42 )
43 
44 fftjet_grid_144_72 = cms.PSet(
45  nEtaBins = cms.uint32(144),
46  etaMin = cms.double(-6.264),
47  etaMax = cms.double(6.264),
48  nPhiBins = cms.uint32(72),
49  phiBin0Edge = cms.double(0.0),
50  title = cms.untracked.string("144 x 72")
51 )
52 
53 fftjet_grid_128_72 = cms.PSet(
54  nEtaBins = cms.uint32(128),
55  etaMin = cms.double(-5.568),
56  etaMax = cms.double(5.568),
57  nPhiBins = cms.uint32(72),
58  phiBin0Edge = cms.double(0.0),
59  title = cms.untracked.string("128 x 72")
60 )
61 
62 #
63 # Definitions for anomalous towers
64 #
65 fftjet_anomalous_tower_default = anom.AnomalousCellParameters
66 
67 fftjet_anomalous_tower_allpass = cms.PSet(
68  maxBadEcalCells = cms.uint32(fftjet_large_int),
69  maxRecoveredEcalCells = cms.uint32(fftjet_large_int),
70  maxProblematicEcalCells = cms.uint32(fftjet_large_int),
71  maxBadHcalCells = cms.uint32(fftjet_large_int),
72  maxRecoveredHcalCells = cms.uint32(fftjet_large_int),
73  maxProblematicHcalCells = cms.uint32(fftjet_large_int)
74 )
75 
76 #
77 # Peak selectors
78 #
79 fftjet_peak_selector_allpass = cms.PSet(
80  Class = cms.string("AllPeaksPass")
81 )
82 
83 #
84 # 50 scales (49 intervals) from 0.087 to 0.6 in log space correspond
85 # to the pattern recognition kernel bandwidth increase of 4.0% per scale
86 #
87 fftjet_patreco_scales_50 = cms.PSet(
88  Class = cms.string("EquidistantInLogSpace"),
89  minScale = cms.double(0.087),
90  maxScale = cms.double(0.6),
91  nScales = cms.uint32(50)
92 )
93 
94 #
95 # Here, the distance calculator for the tree is a simple eta-phi
96 # distance with fixed bandwidth values in eta and phi. However,
97 # if the "etaDependentScaleFactors" are given in the module
98 # configuration, it makes a lot of sense to use eta-dependent
99 # eta-to-phi bandwidth ratio.
100 #
101 fftjet_fixed_bandwidth_distance = cms.PSet(
102  Class = cms.string("PeakEtaPhiDistance"),
103  etaToPhiBandwidthRatio = cms.double(1.0/fftjet_phi_to_eta_bw_ratio)
104 )
105 
106 #
107 # A placeholder for the tree distance calculator with eta-dependent
108 # eta-to-phi bandwidth ratio (must be modified for meaningful use).
109 # Inside the interpolator, bandwidth ratio points are placed at the
110 # cell centers.
111 #
112 fftjet_variable_bandwidth_distance = cms.PSet(
113  Class = cms.string("PeakEtaDependentDistance"),
114  Interpolator = cms.PSet(
115  xmin = cms.double(-5.2),
116  xmax = cms.double(5.2),
117  flow = cms.double(1.0),
118  fhigh = cms.double(1.0),
119  data = cms.vdouble(1.0, 1.0)
120  )
121 )
122 
123 #
124 # Various jet membership functions
125 #
126 fftjet_jet_membership_cone = cms.PSet(
127  Class = cms.string("Linear2d"),
128  sx = cms.double(math.sqrt(1.0/fftjet_phi_to_eta_bw_ratio)),
129  sy = cms.double(math.sqrt(fftjet_phi_to_eta_bw_ratio)),
130  scalePower = cms.int32(1),
131  kernelParameters = cms.vdouble()
132 )
133 
134 #
135 # Background/noise membership functions
136 #
137 fftjet_noise_membership_smallconst = cms.PSet(
138  Class = cms.string("GaussianNoiseMembershipFcn"),
139  minWeight = cms.double(1.0e-8),
140  prior = cms.double(0.0)
141 )
142 
143 #
144 # Distance between jets for convergence determination
145 #
146 fftjet_convergence_jet_distance = cms.PSet(
147  Class = cms.string("JetConvergenceDistance"),
148  etaToPhiBandwidthRatio = cms.double(1.0/fftjet_phi_to_eta_bw_ratio),
149  relativePtBandwidth = cms.double(10.0)
150 )
151 
152 #
153 # Various peak functors
154 #
155 fftjet_peakfunctor_const_zero = cms.PSet(
156  Class = cms.string("ConstDouble"),
157  value = cms.double(0.0)
158 )
159 
160 fftjet_peakfunctor_const_one = cms.PSet(
161  Class = cms.string("ConstDouble"),
162  value = cms.double(1.0)
163 )
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40