CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PileupJetIDCutParams_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 ###########################################################
4 ## Working points for the 80X training
5 ###########################################################
6 full_80x_chs_wp = cms.PSet(
7  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
8 
9  #Tight Id
10  Pt010_Tight = cms.vdouble( 0.26, -0.34, -0.24, -0.26),
11  Pt1020_Tight = cms.vdouble( 0.26, -0.34, -0.24, -0.26),
12  Pt2030_Tight = cms.vdouble( 0.26, -0.34, -0.24, -0.26),
13  Pt3050_Tight = cms.vdouble( 0.62, -0.21, -0.07, -0.03),
14 
15  #Medium Id
16  Pt010_Medium = cms.vdouble(-0.49, -0.53, -0.44, -0.42),
17  Pt1020_Medium = cms.vdouble(-0.49, -0.53, -0.44, -0.42),
18  Pt2030_Medium = cms.vdouble(-0.49, -0.53, -0.44, -0.42),
19  Pt3050_Medium = cms.vdouble(-0.06, -0.42, -0.3 , -0.23),
20 
21  #Loose Id
22  Pt010_Loose = cms.vdouble(-0.96, -0.64, -0.56, -0.54),
23  Pt1020_Loose = cms.vdouble(-0.96, -0.64, -0.56, -0.54),
24  Pt2030_Loose = cms.vdouble(-0.96, -0.64, -0.56, -0.54),
25  Pt3050_Loose = cms.vdouble(-0.92, -0.56, -0.44, -0.39)
26 )
27 
28 
29 ###########################################################
30 ## Working points for the 76X training
31 ###########################################################
32 full_76x_chs_wp = cms.PSet(
33  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
34 
35  #Tight Id
36  Pt010_Tight = cms.vdouble(0.09,-0.37,-0.24,-0.21),
37  Pt1020_Tight = cms.vdouble(0.09,-0.37,-0.24,-0.21),
38  Pt2030_Tight = cms.vdouble(0.09,-0.37,-0.24,-0.21),
39  Pt3050_Tight = cms.vdouble(0.52,-0.19,-0.06,-0.03),
40 
41  #Medium Id
42  Pt010_Medium = cms.vdouble(-0.58,-0.52,-0.40,-0.36),
43  Pt1020_Medium = cms.vdouble(-0.58,-0.52,-0.40,-0.36),
44  Pt2030_Medium = cms.vdouble(-0.58,-0.52,-0.40,-0.36),
45  Pt3050_Medium = cms.vdouble(-0.20,-0.39,-0.24,-0.19),
46 
47  #Loose Id
48  Pt010_Loose = cms.vdouble(-0.96,-0.62,-0.53,-0.49),
49  Pt1020_Loose = cms.vdouble(-0.96,-0.62,-0.53,-0.49),
50  Pt2030_Loose = cms.vdouble(-0.96,-0.62,-0.53,-0.49),
51  Pt3050_Loose = cms.vdouble(-0.93,-0.52,-0.39,-0.31)
52 )
53 
54 ###########################################################
55 ## Working points for the 74X training
56 ###########################################################
57 full_74x_chs_wp = cms.PSet(
58  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
59 
60  #Tight Id
61  Pt010_Tight = cms.vdouble(-0.1,-0.83,-0.83,-0.98),
62  Pt1020_Tight = cms.vdouble(-0.1,-0.83,-0.83,-0.98),
63  Pt2030_Tight = cms.vdouble(-0.1,-0.83,-0.83,-0.98),
64  Pt3050_Tight = cms.vdouble(-0.5,-0.77,-0.80,-0.98),
65 
66  #Medium Id
67  Pt010_Medium = cms.vdouble(-0.3,-0.87,-0.87,-0.99),
68  Pt1020_Medium = cms.vdouble(-0.3,-0.87,-0.87,-0.99),
69  Pt2030_Medium = cms.vdouble(-0.3,-0.87,-0.87,-0.99),
70  Pt3050_Medium = cms.vdouble(-0.6,-0.85,-0.85,-0.99),
71 
72  #Loose Id
73  Pt010_Loose = cms.vdouble(-0.8,-0.97,-0.97,-0.99),
74  Pt1020_Loose = cms.vdouble(-0.8,-0.97,-0.97,-0.99),
75  Pt2030_Loose = cms.vdouble(-0.8,-0.97,-0.97,-0.99),
76  Pt3050_Loose = cms.vdouble(-0.8,-0.95,-0.97,-0.99)
77 )
78 
79 ###########################################################
80 ## Working points for the 53X training/New Met Dec 21, 2012
81 ###########################################################
82 full_53x_wp = cms.PSet(
83  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
84 
85  #Tight Id
86  Pt010_Tight = cms.vdouble(-0.83,-0.81,-0.74,-0.81),
87  Pt1020_Tight = cms.vdouble(-0.83,-0.81,-0.74,-0.81),
88  Pt2030_Tight = cms.vdouble( 0.73, 0.05,-0.26,-0.42),
89  Pt3050_Tight = cms.vdouble( 0.73, 0.05,-0.26,-0.42),
90 
91  #Medium Id
92  Pt010_Medium = cms.vdouble(-0.83,-0.92,-0.90,-0.92),
93  Pt1020_Medium = cms.vdouble(-0.83,-0.92,-0.90,-0.92),
94  Pt2030_Medium = cms.vdouble( 0.10,-0.36,-0.54,-0.54),
95  Pt3050_Medium = cms.vdouble( 0.10,-0.36,-0.54,-0.54),
96 
97  #Loose Id
98  Pt010_Loose = cms.vdouble(-0.95,-0.96,-0.94,-0.95),
99  Pt1020_Loose = cms.vdouble(-0.95,-0.96,-0.94,-0.95),
100  Pt2030_Loose = cms.vdouble(-0.63,-0.60,-0.55,-0.45),
101  Pt3050_Loose = cms.vdouble(-0.63,-0.60,-0.55,-0.45),
102 
103  #MET
104  Pt010_MET = cms.vdouble( 0. ,-0.6,-0.4,-0.4),
105  Pt1020_MET = cms.vdouble( 0.3 ,-0.2,-0.4,-0.4),
106  Pt2030_MET = cms.vdouble( 0. , 0. , 0. , 0. ),
107  Pt3050_MET = cms.vdouble( 0. , 0. ,-0.1,-0.2)
108  )
109 
110 full_53x_chs_wp = cms.PSet(
111  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
112 
113  #Tight Id
114  Pt010_Tight = cms.vdouble(-0.83,-0.81,-0.74,-0.81),
115  Pt1020_Tight = cms.vdouble(-0.83,-0.81,-0.74,-0.81),
116  Pt2030_Tight = cms.vdouble( 0.78, 0.50, 0.17, 0.17),
117  Pt3050_Tight = cms.vdouble( 0.78, 0.50, 0.17, 0.17),
118 
119  #Medium Id
120  Pt010_Medium = cms.vdouble(-0.83,-0.92,-0.90,-0.92),
121  Pt1020_Medium = cms.vdouble(-0.83,-0.92,-0.90,-0.92),
122  Pt2030_Medium = cms.vdouble(-0.07,-0.09, 0.00,-0.06),
123  Pt3050_Medium = cms.vdouble(-0.07,-0.09, 0.00,-0.06),
124 
125  #Loose Id
126  Pt010_Loose = cms.vdouble(-0.95,-0.96,-0.94,-0.95),
127  Pt1020_Loose = cms.vdouble(-0.95,-0.96,-0.94,-0.95),
128  Pt2030_Loose = cms.vdouble(-0.15,-0.26,-0.16,-0.16),
129  Pt3050_Loose = cms.vdouble(-0.15,-0.26,-0.16,-0.16),
130  )
131 
132 met_53x_wp = cms.PSet(
133 
134  #Tight Id
135  Pt010_Tight = cms.vdouble(-2, -2, -2, -2, -2),
136  Pt1020_Tight = cms.vdouble(-2, -2, -2, -2, -2),
137  Pt2030_Tight = cms.vdouble(-2, -2, -2, -2, -2),
138  Pt3050_Tight = cms.vdouble(-2, -2, -2, -2, -2),
139 
140  #Medium Id
141  Pt010_Medium = cms.vdouble(-2, -2, -2, -2, -2),
142  Pt1020_Medium = cms.vdouble(-2, -2, -2, -2, -2),
143  Pt2030_Medium = cms.vdouble(-2, -2, -2, -2, -2),
144  Pt3050_Medium = cms.vdouble(-2, -2, -2, -2, -2),
145 
146  #Loose Id
147  Pt010_Loose = cms.vdouble(-2, -2, -2, -2, -2),
148  Pt1020_Loose = cms.vdouble(-2, -2, -2, -2, -2),
149  Pt2030_Loose = cms.vdouble(-2, -2, -2, -2, -2),
150  Pt3050_Loose = cms.vdouble(-2, -2, -2, -2, -2),
151 
152  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
153  #MET
154  Pt010_MET = cms.vdouble(-0.2 ,-0.3,-0.5,-0.5),
155  Pt1020_MET = cms.vdouble(-0.2 ,-0.2,-0.5,-0.3),
156  Pt2030_MET = cms.vdouble(-0.2 ,-0.2,-0.2, 0.1),
157  Pt3050_MET = cms.vdouble(-0.2 ,-0.2, 0. , 0.2)
158  )
159 
160 metfull_53x_wp = cms.PSet(
161  #MET
162  Pt010_MET = cms.vdouble(-0.2 ,-0.3,-0.5,-0.5),
163  Pt1020_MET = cms.vdouble(-0.2 ,-0.2,-0.5,-0.3),
164  Pt2030_MET = cms.vdouble( 0. , 0. , 0. , 0. ),
165  Pt3050_MET = cms.vdouble( 0. , 0. ,-0.1,-0.2)
166  )
167 
168 
169 ###########################################################
170 ## Working points for the 5X training
171 ###########################################################
172 full_5x_wp = cms.PSet(
173  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
174 
175  #Tight Id
176  Pt010_Tight = cms.vdouble(-0.47,-0.92,-0.92,-0.94),
177  Pt1020_Tight = cms.vdouble(-0.47,-0.92,-0.92,-0.94),
178  Pt2030_Tight = cms.vdouble(+0.32,-0.49,-0.61,-0.74),
179  Pt3050_Tight = cms.vdouble(+0.32,-0.49,-0.61,-0.74),
180 
181  #Medium Id
182  Pt010_Medium = cms.vdouble(-0.83,-0.96,-0.95,-0.96),
183  Pt1020_Medium = cms.vdouble(-0.83,-0.96,-0.95,-0.96),
184  Pt2030_Medium = cms.vdouble(-0.40,-0.74,-0.76,-0.81),
185  Pt3050_Medium = cms.vdouble(-0.40,-0.74,-0.76,-0.81),
186 
187  #Loose Id
188  Pt010_Loose = cms.vdouble(-0.95,-0.97,-0.97,-0.97),
189  Pt1020_Loose = cms.vdouble(-0.95,-0.97,-0.97,-0.97),
190  Pt2030_Loose = cms.vdouble(-0.80,-0.85,-0.84,-0.85),
191  Pt3050_Loose = cms.vdouble(-0.80,-0.85,-0.84,-0.85)
192 )
193 
194 simple_5x_wp = cms.PSet(
195  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
196 
197  #Tight Id
198  Pt010_Tight = cms.vdouble(-0.54,-0.93,-0.93,-0.94),
199  Pt1020_Tight = cms.vdouble(-0.54,-0.93,-0.93,-0.94),
200  Pt2030_Tight = cms.vdouble(+0.26,-0.54,-0.63,-0.74),
201  Pt3050_Tight = cms.vdouble(+0.26,-0.54,-0.63,-0.74),
202 
203  #Medium Id
204  Pt010_Medium = cms.vdouble(-0.85,-0.96,-0.95,-0.96),
205  Pt1020_Medium = cms.vdouble(-0.85,-0.96,-0.95,-0.96),
206  Pt2030_Medium = cms.vdouble(-0.40,-0.73,-0.74,-0.80),
207  Pt3050_Medium = cms.vdouble(-0.40,-0.73,-0.74,-0.80),
208 
209  #Loose Id
210  Pt010_Loose = cms.vdouble(-0.95,-0.97,-0.96,-0.97),
211  Pt1020_Loose = cms.vdouble(-0.95,-0.97,-0.96,-0.97),
212  Pt2030_Loose = cms.vdouble(-0.80,-0.86,-0.80,-0.84),
213  Pt3050_Loose = cms.vdouble(-0.80,-0.86,-0.80,-0.84)
214 
215 )
216 
217 ###########################################################
218 ## Working points for the 5X_CHS training
219 ###########################################################
220 full_5x_chs_wp = cms.PSet(
221  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
222 
223  #Tight Id
224  Pt010_Tight = cms.vdouble(-0.59,-0.75,-0.78,-0.80),
225  Pt1020_Tight = cms.vdouble(-0.59,-0.75,-0.78,-0.80),
226  Pt2030_Tight = cms.vdouble(+0.41,-0.10,-0.20,-0.45),
227  Pt3050_Tight = cms.vdouble(+0.41,-0.10,-0.20,-0.45),
228 
229  #Medium Id
230  Pt010_Medium = cms.vdouble(-0.94,-0.91,-0.91,-0.92),
231  Pt1020_Medium = cms.vdouble(-0.94,-0.91,-0.91,-0.92),
232  Pt2030_Medium = cms.vdouble(-0.58,-0.65,-0.57,-0.67),
233  Pt3050_Medium = cms.vdouble(-0.58,-0.65,-0.57,-0.67),
234 
235  #Loose Id
236  Pt010_Loose = cms.vdouble(-0.98,-0.95,-0.94,-0.94),
237  Pt1020_Loose = cms.vdouble(-0.98,-0.95,-0.94,-0.94),
238  Pt2030_Loose = cms.vdouble(-0.89,-0.77,-0.69,-0.75),
239  Pt3050_Loose = cms.vdouble(-0.89,-0.77,-0.69,-0.57)
240 )
241 
242 simple_5x_chs_wp = cms.PSet(
243  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
244 
245  #Tight Id
246  Pt010_Tight = cms.vdouble(-0.60,-0.74,-0.78,-0.81),
247  Pt1020_Tight = cms.vdouble(-0.60,-0.74,-0.78,-0.81),
248  Pt2030_Tight = cms.vdouble(-0.47,-0.06,-0.23,-0.47),
249  Pt3050_Tight = cms.vdouble(-0.47,-0.06,-0.23,-0.47),
250 
251  #Medium Id
252  Pt010_Medium = cms.vdouble(-0.95,-0.94,-0.92,-0.91),
253  Pt1020_Medium = cms.vdouble(-0.95,-0.94,-0.92,-0.91),
254  Pt2030_Medium = cms.vdouble(-0.59,-0.65,-0.56,-0.68),
255  Pt3050_Medium = cms.vdouble(-0.59,-0.65,-0.56,-0.68),
256 
257  #Loose Id
258  Pt010_Loose = cms.vdouble(-0.98,-0.96,-0.94,-0.94),
259  Pt1020_Loose = cms.vdouble(-0.98,-0.96,-0.94,-0.94),
260  Pt2030_Loose = cms.vdouble(-0.89,-0.75,-0.72,-0.75),
261  Pt3050_Loose = cms.vdouble(-0.89,-0.75,-0.72,-0.75)
262 )
263 
264 
265 ###########################################################
266 ## Working points for the 4X training
267 ###########################################################
268 PuJetIdOptMVA_wp = cms.PSet(
269  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
270 
271  #Tight Id
272  Pt010_Tight = cms.vdouble(-0.5,-0.2,-0.83,-0.7),
273  Pt1020_Tight = cms.vdouble(-0.5,-0.2,-0.83,-0.7),
274  Pt2030_Tight = cms.vdouble(-0.2, 0., 0., 0.),
275  Pt3050_Tight = cms.vdouble(-0.2, 0., 0., 0.),
276 
277  #Medium Id
278  Pt010_Medium = cms.vdouble(-0.73,-0.89,-0.89,-0.83),
279  Pt1020_Medium = cms.vdouble(-0.73,-0.89,-0.89,-0.83),
280  Pt2030_Medium = cms.vdouble(0.1, -0.4, -0.4, -0.45),
281  Pt3050_Medium = cms.vdouble(0.1, -0.4, -0.4, -0.45),
282 
283  #Loose Id
284  Pt010_Loose = cms.vdouble(-0.9,-0.9, -0.9,-0.9),
285  Pt1020_Loose = cms.vdouble(-0.9,-0.9, -0.9,-0.9),
286  Pt2030_Loose = cms.vdouble(-0.4,-0.85,-0.7,-0.6),
287  Pt3050_Loose = cms.vdouble(-0.4,-0.85,-0.7,-0.6)
288 )
289 
290 PuJetIdMinMVA_wp = cms.PSet(
291  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
292 
293  #Tight Id
294  Pt010_Tight = cms.vdouble(-0.5,-0.2,-0.83,-0.7),
295  Pt1020_Tight = cms.vdouble(-0.5,-0.2,-0.83,-0.7),
296  Pt2030_Tight = cms.vdouble(-0.2, 0., 0., 0.),
297  Pt3050_Tight = cms.vdouble(-0.2, 0., 0., 0.),
298 
299  #Medium Id
300  Pt010_Medium = cms.vdouble(-0.73,-0.89,-0.89,-0.83),
301  Pt1020_Medium = cms.vdouble(-0.73,-0.89,-0.89,-0.83),
302  Pt2030_Medium = cms.vdouble(0.1, -0.4, -0.5, -0.45),
303  Pt3050_Medium = cms.vdouble(0.1, -0.4, -0.5, -0.45),
304 
305  #Loose Id
306  Pt010_Loose = cms.vdouble(-0.9,-0.9, -0.94,-0.9),
307  Pt1020_Loose = cms.vdouble(-0.9,-0.9, -0.94,-0.9),
308  Pt2030_Loose = cms.vdouble(-0.4,-0.85,-0.7,-0.6),
309  Pt3050_Loose = cms.vdouble(-0.4,-0.85,-0.7,-0.6)
310 )
311 
312 EmptyJetIdParams = cms.PSet(
313  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
314 
315  #Tight Id
316  Pt010_Tight = cms.vdouble(-999.,-999.,-999.,-999.),
317  Pt1020_Tight = cms.vdouble(-999.,-999.,-999.,-999.),
318  Pt2030_Tight = cms.vdouble(-999.,-999.,-999.,-999.),
319  Pt3050_Tight = cms.vdouble(-999.,-999.,-999.,-999.),
320 
321  #Medium Id
322  Pt010_Medium = cms.vdouble(-999.,-999.,-999.,-999.),
323  Pt1020_Medium = cms.vdouble(-999.,-999.,-999.,-999.),
324  Pt2030_Medium = cms.vdouble(-999.,-999.,-999.,-999.),
325  Pt3050_Medium = cms.vdouble(-999.,-999.,-999.,-999.),
326 
327  #Loose Id
328  Pt010_Loose = cms.vdouble(-999.,-999.,-999.,-999.),
329  Pt1020_Loose = cms.vdouble(-999.,-999.,-999.,-999.),
330  Pt2030_Loose = cms.vdouble(-999.,-999.,-999.,-999.),
331  Pt3050_Loose = cms.vdouble(-999.,-999.,-999.,-999.)
332 )
333 
334 
335 PuJetIdCutBased_wp = cms.PSet(
336  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
337  #betaStarClassic/log(nvtx-0.64) Values
338  #Tight Id
339  Pt010_BetaStarTight = cms.vdouble( 0.15, 0.15, 999., 999.),
340  Pt1020_BetaStarTight = cms.vdouble( 0.15, 0.15, 999., 999.),
341  Pt2030_BetaStarTight = cms.vdouble( 0.15, 0.15, 999., 999.),
342  Pt3050_BetaStarTight = cms.vdouble( 0.15, 0.15, 999., 999.),
343 
344  #Medium Id => Daniele
345  Pt010_BetaStarMedium = cms.vdouble( 0.2, 0.3, 999., 999.),
346  Pt1020_BetaStarMedium = cms.vdouble( 0.2, 0.3, 999., 999.),
347  Pt2030_BetaStarMedium = cms.vdouble( 0.2, 0.3, 999., 999.),
348  Pt3050_BetaStarMedium = cms.vdouble( 0.2, 0.3, 999., 999.),
349 
350  #Loose Id
351  Pt010_BetaStarLoose = cms.vdouble( 0.2, 0.3, 999., 999.),
352  Pt1020_BetaStarLoose = cms.vdouble( 0.2, 0.3, 999., 999.),
353  Pt2030_BetaStarLoose = cms.vdouble( 0.2, 0.3, 999., 999.),
354  Pt3050_BetaStarLoose = cms.vdouble( 0.2, 0.3, 999., 999.),
355 
356  #RMS variable
357  #Tight Id
358  Pt010_RMSTight = cms.vdouble( 0.06, 0.07, 0.04, 0.05),
359  Pt1020_RMSTight = cms.vdouble( 0.06, 0.07, 0.04, 0.05),
360  Pt2030_RMSTight = cms.vdouble( 0.05, 0.07, 0.03, 0.045),
361  Pt3050_RMSTight = cms.vdouble( 0.05, 0.06, 0.03, 0.04),
362 
363  #Medium Id => Daniele
364  Pt010_RMSMedium = cms.vdouble( 0.06, 0.03, 0.03, 0.04),
365  Pt1020_RMSMedium = cms.vdouble( 0.06, 0.03, 0.03, 0.04),
366  Pt2030_RMSMedium = cms.vdouble( 0.06, 0.03, 0.03, 0.04),
367  Pt3050_RMSMedium = cms.vdouble( 0.06, 0.03, 0.03, 0.04),
368 
369  #Loose Id
370  Pt010_RMSLoose = cms.vdouble( 0.06, 0.05, 0.05, 0.07),
371  Pt1020_RMSLoose = cms.vdouble( 0.06, 0.05, 0.05, 0.07),
372  Pt2030_RMSLoose = cms.vdouble( 0.06, 0.05, 0.05, 0.055),
373  Pt3050_RMSLoose = cms.vdouble( 0.06, 0.05, 0.05, 0.055)
374  )
375 
376 
377 JetIdParams = cms.PSet(
378  #4 Eta Categories 0-2.5 2.5-2.75 2.75-3.0 3.0-5.0
379 
380  #Tight Id
381  Pt010_Tight = cms.vdouble( 0.5,0.6,0.6,0.9),
382  Pt1020_Tight = cms.vdouble(-0.2,0.2,0.2,0.6),
383  Pt2030_Tight = cms.vdouble( 0.3,0.4,0.7,0.8),
384  Pt3050_Tight = cms.vdouble( 0.5,0.4,0.8,0.9),
385 
386  #Medium Id
387  Pt010_Medium = cms.vdouble( 0.2,0.4,0.2,0.6),
388  Pt1020_Medium = cms.vdouble(-0.3,0. ,0. ,0.5),
389  Pt2030_Medium = cms.vdouble( 0.2,0.2,0.5,0.7),
390  Pt3050_Medium = cms.vdouble( 0.3,0.2,0.7,0.8),
391 
392  #Loose Id
393  Pt010_Loose = cms.vdouble( 0. , 0. , 0. ,0.2),
394  Pt1020_Loose = cms.vdouble(-0.4,-0.4,-0.4,0.4),
395  Pt2030_Loose = cms.vdouble( 0. , 0. , 0.2,0.6),
396  Pt3050_Loose = cms.vdouble( 0. , 0. , 0.6,0.2)
397 )
398