CMS 3D CMS Logo

VtxSmearedParameters_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 #
4 # All units are cm and radians
5 #
6 # UNITS:
7 #
8 # TimeOffset in nanoseconds
9 # spacial displacement in cm
10 
11 # common parameters
12 VtxSmearedCommon = cms.PSet(
13  src = cms.InputTag("generator", "unsmeared"),
14  readDB = cms.bool(False)
15 )
16 # Gaussian smearing
17 GaussVtxSmearingParameters = cms.PSet(
18  MeanX = cms.double(0.0),
19  MeanY = cms.double(0.0),
20  MeanZ = cms.double(0.0),
21  SigmaY = cms.double(0.0015),
22  SigmaX = cms.double(0.0015),
23  SigmaZ = cms.double(5.3),
24  TimeOffset = cms.double(0.0)
25 )
26 # Gaussian smearing
27 GaussVtxSigmaZ4cmSmearingParameters = cms.PSet(
28  MeanX = cms.double(0.0),
29  MeanY = cms.double(0.0),
30  MeanZ = cms.double(0.0),
31  SigmaY = cms.double(0.0015),
32  SigmaX = cms.double(0.0015),
33  SigmaZ = cms.double(4.0),
34  TimeOffset = cms.double(0.0)
35 )
36 # Gaussian smearing
37 # Flat optics for Run3 - Low SigmaZ
38 # SigmaZ = 4.2 cm
39 # SigmaX = 11.8 um
40 # SigmaY = 5.5 um
41 # BS positions extracted from 2018B 3.8T data, run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
42 # X0 = 0.09676 [cm]
43 # Y0 = -0.06245 [cm]
44 # Z0 = -0.292 [cm]
45 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
46 # X = 0.0859918 cm
47 # Y = -0.104172 cm
48 # Z = -0.327748 cm
49 Run3FlatOpticsGaussVtxSigmaZ4p2cmSmearingParameters = cms.PSet(
50  MeanX = cms.double(0.0107682),
51  MeanY = cms.double(0.041722),
52  MeanZ = cms.double(0.035748),
53  SigmaY = cms.double(0.00055),
54  SigmaX = cms.double(0.00118),
55  SigmaZ = cms.double(4.2),
56  TimeOffset = cms.double(0.0)
57 )
58 # Gaussian smearing
59 # Flat optics for Run3 - High SigmaZ
60 # SigmaZ = 5.3 cm
61 # SigmaX = 15 um
62 # SigmaY = 13 um
63 # BS positions extracted from 2018B 3.8T data, run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
64 # X0 = 0.09676 [cm]
65 # Y0 = -0.06245 [cm]
66 # Z0 = -0.292 [cm]
67 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
68 # X = 0.0859918 cm
69 # Y = -0.104172 cm
70 # Z = -0.327748 cm
71 Run3FlatOpticsGaussVtxSigmaZ5p3cmSmearingParameters = cms.PSet(
72  MeanX = cms.double(0.0107682),
73  MeanY = cms.double(0.041722),
74  MeanZ = cms.double(0.035748),
75  SigmaY = cms.double(0.0013),
76  SigmaX = cms.double(0.0015),
77  SigmaZ = cms.double(5.3),
78  TimeOffset = cms.double(0.0)
79 )
80 
81 # Flat Smearing
82 # Important note: flat independent distributions in Z and T are not correct for physics production
83 # In reality, if two flat beams interact the real distribution will not be flat with independent Z and T
84 # but Z and T will be correlated, as example in GaussEvtVtxGenerator.
85 # Can restore correlation via MinT += (MinZ - MaxZ)/2 and MaxT += (MaxZ - MinZ)/2
86 # in [ns] units (recall c_light = 29.98cm/ns)
87 FlatVtxSmearingParameters = cms.PSet(
88  MaxZ = cms.double(5.3),
89  MaxX = cms.double(0.0015),
90  MaxY = cms.double(0.0015),
91  MinX = cms.double(-0.0015),
92  MinY = cms.double(-0.0015),
93  MinZ = cms.double(-5.3),
94  MaxT = cms.double(0.177),
95  MinT = cms.double(-0.177)
96 )
97 #############################################
98 # Beta functions smearing (pp 7+7 TeV)
99 #
100 # Values taken from LHC optics simulation V6.5:
101 # see http://proj-lhc-optics-web.web.cern.ch/proj-lhc-optics-web/V6.500/IR5.html
102 # alpha = angle of the crossing plane 0 degrees means XZ plane
103 # phi = half-crossing beam angle
104 #
105 # Emittance is the no normalized emittance in cm = normalized emittance/gamma (beta=1)
106 #
107 # length variables are in cm
108 #
109 
110 # 900 GeV collisions, transverse beam size = 293 microns
111 Early900GeVCollisionVtxSmearingParameters = cms.PSet(
112  Phi = cms.double(0.0),
113  BetaStar = cms.double(1100.0),
114  Emittance = cms.double(1.564e-06),
115  Alpha = cms.double(0.0),
116  SigmaZ = cms.double(7.4),
117  TimeOffset = cms.double(0.0),
118  Y0 = cms.double(0.0),
119  X0 = cms.double(0.0322),
120  Z0 = cms.double(0.0)
121 )
122 # 2.2 TeV collisions, transverse beam size 188 microns
123 Early2p2TeVCollisionVtxSmearingParameters = cms.PSet(
124  Phi = cms.double(0.0),
125  BetaStar = cms.double(1100.0),
126  Emittance = cms.double(6.4e-07),
127  Alpha = cms.double(0.0),
128  SigmaZ = cms.double(5.5),
129  TimeOffset = cms.double(0.0),
130  Y0 = cms.double(0.0),
131  X0 = cms.double(0.0322),
132  Z0 = cms.double(0.0)
133 )
134 # 7 TeV collisions, transverse beam size with betastar= 11m is 105 microns,
135 Early7TeVCollisionVtxSmearingParameters = cms.PSet(
136  Phi = cms.double(0.0),
137  BetaStar = cms.double(1100.0),
138  Emittance = cms.double(2.0e-07),
139  Alpha = cms.double(0.0),
140  SigmaZ = cms.double(4.2),
141  TimeOffset = cms.double(0.0),
142  Y0 = cms.double(0.0),
143  X0 = cms.double(0.0322),
144  Z0 = cms.double(0.0)
145 )
146 # 7 TeV collisions, transverse beam size with betastar= 2m is 45 microns,
147 Nominal7TeVCollisionVtxSmearingParameters = cms.PSet(
148  Phi = cms.double(0.0),
149  BetaStar = cms.double(200.0),
150  Emittance = cms.double(2.0e-07),
151  Alpha = cms.double(0.0),
152  SigmaZ = cms.double(4.2),
153  TimeOffset = cms.double(0.0),
154  Y0 = cms.double(0.0),
155  X0 = cms.double(0.0322),
156  Z0 = cms.double(0.0)
157 )
158 # 900 GeV realistic 2010 collisions, transverse beam size is 200 microns
159 Realistic900GeVCollisionVtxSmearingParameters = cms.PSet(
160  Phi = cms.double(0.0),
161  BetaStar = cms.double(1000.0),
162  Emittance = cms.double(8.34e-07),
163  Alpha = cms.double(0.0),
164  SigmaZ = cms.double(6.17),
165  TimeOffset = cms.double(0.0),
166  X0 = cms.double(0.2452),
167  Y0 = cms.double(0.3993),
168  Z0 = cms.double(0.8222)
169 )
170 # 7 TeV realistic collisions, beamspot width ~28 microns - appropriate for 2nd half of Commissioning10
171 Realistic7TeVCollisionComm10VtxSmearingParameters = cms.PSet(
172  Phi = cms.double(0.0),
173  BetaStar = cms.double(200.0),
174  Emittance = cms.double(0.804e-07),
175  Alpha = cms.double(0.0),
176  SigmaZ = cms.double(3.50),
177  TimeOffset = cms.double(0.0),
178  X0 = cms.double(0.2440),
179  Y0 = cms.double(0.3929),
180  Z0 = cms.double(0.4145)
181 )
182 # 7 TeV realistic collisions, beamspot width ~43 microns - appropriate for 2010A
183 Realistic7TeVCollisionVtxSmearingParameters = cms.PSet(
184  Phi = cms.double(0.0),
185  BetaStar = cms.double(350.0),
186  Emittance = cms.double(1.072e-07),
187  Alpha = cms.double(0.0),
188  SigmaZ = cms.double(6.26),
189  TimeOffset = cms.double(0.0),
190  X0 = cms.double(0.2440),
191  Y0 = cms.double(0.3929),
192  Z0 = cms.double(0.4145)
193 )
194 # 7 TeV realistic collisions, beamspot width ~38 microns - appropriate for 2010B
195 Realistic7TeVCollision2010BVtxSmearingParameters = cms.PSet(
196  Phi = cms.double(0.0),
197  BetaStar = cms.double(350.0),
198  Emittance = cms.double(0.804e-07),
199  Alpha = cms.double(0.0),
200  SigmaZ = cms.double(5.40),
201  TimeOffset = cms.double(0.0),
202  X0 = cms.double(0.2440),
203  Y0 = cms.double(0.3929),
204  Z0 = cms.double(0.4145)
205 )
206 # 7 TeV realistic collisions, updated for 2011
207 # normalized emittance 2.5 microns, transverse beam size is 32 microns
208 Realistic7TeV2011CollisionVtxSmearingParameters = cms.PSet(
209  Phi = cms.double(0.0),
210  BetaStar = cms.double(150.0),
211  Emittance = cms.double(0.67e-07),
212  Alpha = cms.double(0.0),
213  SigmaZ = cms.double(5.22),
214  TimeOffset = cms.double(0.0),
215  X0 = cms.double(0.2440),
216  Y0 = cms.double(0.3929),
217  Z0 = cms.double(0.4145)
218 )
219 # HI realistic collisions, updated for 2011
220 # estimated beamspot width 31-35 microns
221 RealisticHI2011CollisionVtxSmearingParameters = cms.PSet(
222  Phi = cms.double(0.0),
223  BetaStar = cms.double(100.0),
224  Emittance = cms.double(2.04e-07),
225  Alpha = cms.double(0.0),
226  SigmaZ = cms.double(7.06),
227  TimeOffset = cms.double(0.0),
228  X0 = cms.double(0.2245),
229  Y0 = cms.double(0.4182),
230  Z0 = cms.double(0.0847)
231 )
232 # 2.76 TeV estimated collisions, 11m beta*
233 # normalized emittance 2.5 microns, transverse beam size is 140 microns
234 Realistic2p76TeV2011CollisionVtxSmearingParameters = cms.PSet(
235  Phi = cms.double(0.0),
236  BetaStar = cms.double(1100.0),
237  Emittance = cms.double(1.70e-07),
238  Alpha = cms.double(0.0),
239  SigmaZ = cms.double(5.22),
240  TimeOffset = cms.double(0.0),
241  X0 = cms.double(0.2440),
242  Y0 = cms.double(0.3929),
243  Z0 = cms.double(0.4145)
244 )
245 # 2.76 TeV estimated collisions for 2013, 11m beta*
246 # sigmaZ set to 8 cm
247 Realistic2p76TeV2013CollisionVtxSmearingParameters = cms.PSet(
248  Phi = cms.double(0.0),
249  BetaStar = cms.double(1100.0),
250  Emittance = cms.double(1.70e-07),
251  Alpha = cms.double(0.0),
252  SigmaZ = cms.double(8.0),
253  TimeOffset = cms.double(0.0),
254  X0 = cms.double(0.2440),
255  Y0 = cms.double(0.3929),
256  Z0 = cms.double(0.4145)
257 )
258 # HI realistic pPb collisions, updated for 2013
259 #
260 RealisticHIpPb2013CollisionVtxSmearingParameters = cms.PSet(
261  Phi = cms.double(0.0),
262  BetaStar = cms.double(80.0),
263  Emittance = cms.double(6.25e-07),
264  Alpha = cms.double(0.0),
265  SigmaZ = cms.double(8.0),
266  TimeOffset = cms.double(0.0),
267  X0 = cms.double(0.2440),
268  Y0 = cms.double(0.3929),
269  Z0 = cms.double(0.4145)
270 )
271 # 7 TeV centered collisions with parameters for 2011
272 # normalized emittance 2.5 microns, transverse beam size is 32 microns
273 Centered7TeV2011CollisionVtxSmearingParameters = cms.PSet(
274  Phi = cms.double(0.0),
275  BetaStar = cms.double(150.0),
276  Emittance = cms.double(0.67e-07),
277  Alpha = cms.double(0.0),
278  SigmaZ = cms.double(5.22),
279  TimeOffset = cms.double(0.0),
280  X0 = cms.double(0.0),
281  Y0 = cms.double(0.0),
282  Z0 = cms.double(0.0)
283 )
284 # 8 TeV realistic collisions, transverse beam width size is 20 microns
285 Realistic8TeVCollisionVtxSmearingParameters = cms.PSet(
286  Phi = cms.double(0.0),
287  BetaStar = cms.double(70.0),
288  Emittance = cms.double(0.586e-07),
289  Alpha = cms.double(0.0),
290  SigmaZ = cms.double(6.16),
291  TimeOffset = cms.double(0.0),
292  X0 = cms.double(0.2440),
293  Y0 = cms.double(0.3929),
294  Z0 = cms.double(0.4145)
295 )
296 # 8 TeV realistic collisions, transverse beam width size is 20 microns, updated for observed SigmaZ
297 Realistic8TeV2012CollisionVtxSmearingParameters = cms.PSet(
298  Phi = cms.double(0.0),
299  BetaStar = cms.double(70.0),
300  Emittance = cms.double(0.586e-07),
301  Alpha = cms.double(0.0),
302  SigmaZ = cms.double(4.8),
303  TimeOffset = cms.double(0.0),
304  X0 = cms.double(0.2440),
305  Y0 = cms.double(0.3929),
306  Z0 = cms.double(0.4145)
307 )
308 # 10 TeV collisions, transverse beam size = 46 microns
309 Early10TeVCollisionVtxSmearingParameters = cms.PSet(
310  Phi = cms.double(0.0),
311  BetaStar = cms.double(300.0),
312  Emittance = cms.double(1.406e-07),
313  Alpha = cms.double(0.0),
314  SigmaZ = cms.double(3.8),
315  TimeOffset = cms.double(0.0),
316  Y0 = cms.double(0.0),
317  X0 = cms.double(0.0322),
318  Z0 = cms.double(0.0)
319 )
320 # Test offset
321 Early10TeVX322Y100VtxSmearingParameters = cms.PSet(
322  Phi = cms.double(0.0),
323  BetaStar = cms.double(300.0),
324  Emittance = cms.double(1.406e-07),
325  Alpha = cms.double(0.0),
326  SigmaZ = cms.double(3.8),
327  TimeOffset = cms.double(0.0),
328  Y0 = cms.double(0.0100),
329  X0 = cms.double(0.0322),
330  Z0 = cms.double(0.0)
331 )
332 # Test offset
333 Early10TeVX322Y250VtxSmearingParameters = cms.PSet(
334  Phi = cms.double(0.0),
335  BetaStar = cms.double(300.0),
336  Emittance = cms.double(1.406e-07),
337  Alpha = cms.double(0.0),
338  SigmaZ = cms.double(3.8),
339  TimeOffset = cms.double(0.0),
340  Y0 = cms.double(0.0250),
341  X0 = cms.double(0.0322),
342  Z0 = cms.double(0.0)
343 )
344 # Test offset
345 Early10TeVX322Y500VtxSmearingParameters = cms.PSet(
346  Phi = cms.double(0.0),
347  BetaStar = cms.double(300.0),
348  Emittance = cms.double(1.406e-07),
349  Alpha = cms.double(0.0),
350  SigmaZ = cms.double(3.8),
351  TimeOffset = cms.double(0.0),
352  Y0 = cms.double(0.0500),
353  X0 = cms.double(0.0322),
354  Z0 = cms.double(0.0)
355 )
356 # Test offset
357 Early10TeVX322Y1000VtxSmearingParameters = cms.PSet(
358  Phi = cms.double(0.0),
359  BetaStar = cms.double(300.0),
360  Emittance = cms.double(1.406e-07),
361  Alpha = cms.double(0.0),
362  SigmaZ = cms.double(3.8),
363  TimeOffset = cms.double(0.0),
364  Y0 = cms.double(0.1),
365  X0 = cms.double(0.0322),
366  Z0 = cms.double(0.0)
367 )
368 # Test offset
369 Early10TeVX322Y5000VtxSmearingParameters = cms.PSet(
370  Phi = cms.double(0.0),
371  BetaStar = cms.double(300.0),
372  Emittance = cms.double(1.406e-07),
373  Alpha = cms.double(0.0),
374  SigmaZ = cms.double(3.8),
375  TimeOffset = cms.double(0.0),
376  Y0 = cms.double(0.5),
377  X0 = cms.double(0.0322),
378  Z0 = cms.double(0.0)
379 )
380 # Test offset
381 Early10TeVX322Y10000VtxSmearingParameters = cms.PSet(
382  Phi = cms.double(0.0),
383  BetaStar = cms.double(300.0),
384  Emittance = cms.double(1.406e-07),
385  Alpha = cms.double(0.0),
386  SigmaZ = cms.double(3.8),
387  TimeOffset = cms.double(0.0),
388  Y0 = cms.double(1.0),
389  X0 = cms.double(0.0322),
390  Z0 = cms.double(0.0)
391 )
392 
393 EarlyCollisionVtxSmearingParameters = cms.PSet(
394  Phi = cms.double(0.0),
395  BetaStar = cms.double(200.0),
396  Emittance = cms.double(1.006e-07),
397  Alpha = cms.double(0.0),
398  SigmaZ = cms.double(5.3),
399  TimeOffset = cms.double(0.0),
400  Y0 = cms.double(0.0),
401  X0 = cms.double(0.0322),
402  Z0 = cms.double(0.0)
403 )
404 NominalCollisionVtxSmearingParameters = cms.PSet(
405  Phi = cms.double(0.000142),
406  BetaStar = cms.double(55.0),
407  Emittance = cms.double(1.006e-07),
408  Alpha = cms.double(0.0),
409  SigmaZ = cms.double(5.3),
410  TimeOffset = cms.double(0.0),
411  Y0 = cms.double(0.0),
412  X0 = cms.double(0.05),
413  Z0 = cms.double(0.0)
414 )
415 NominalCollision1VtxSmearingParameters = cms.PSet(
416  Phi = cms.double(0.0),
417  BetaStar = cms.double(55.0),
418  Emittance = cms.double(1.006e-07),
419  Alpha = cms.double(0.0),
420  SigmaZ = cms.double(5.3),
421  TimeOffset = cms.double(0.0),
422  Y0 = cms.double(0.025),
423  X0 = cms.double(0.05),
424  Z0 = cms.double(0.0)
425 )
426 NominalCollision2VtxSmearingParameters = cms.PSet(
427  Phi = cms.double(0.000142),
428  BetaStar = cms.double(55.0),
429  Emittance = cms.double(1.006e-07),
430  Alpha = cms.double(0.0),
431  SigmaZ = cms.double(5.3),
432  TimeOffset = cms.double(0.0),
433  Y0 = cms.double(0.025),
434  X0 = cms.double(0.05),
435  Z0 = cms.double(0.0)
436 )
437 NominalCollision3VtxSmearingParameters = cms.PSet(
438  Phi = cms.double(0.0),
439  BetaStar = cms.double(55.0),
440  Emittance = cms.double(1.006e-07),
441  Alpha = cms.double(0.0),
442  SigmaZ = cms.double(5.3),
443  TimeOffset = cms.double(0.0),
444  Y0 = cms.double(0.025),
445  X0 = cms.double(0.1),
446  Z0 = cms.double(0.0)
447 )
448 NominalCollision4VtxSmearingParameters = cms.PSet(
449  Phi = cms.double(0.0),
450  BetaStar = cms.double(55.0),
451  Emittance = cms.double(1.006e-07),
452  Alpha = cms.double(0.0),
453  SigmaZ = cms.double(5.3),
454  TimeOffset = cms.double(0.0),
455  Y0 = cms.double(0.025),
456  X0 = cms.double(0.2),
457  Z0 = cms.double(0.0)
458 )
459 NominalCollision2015VtxSmearingParameters = cms.PSet(
460  Phi = cms.double(0.0),
461  BetaStar = cms.double(65.0),
462  Emittance = cms.double(5.411e-08),
463  Alpha = cms.double(0.0),
464  SigmaZ = cms.double(5.3),
465  TimeOffset = cms.double(0.0),
466  X0 = cms.double(0.0322),
467  Y0 = cms.double(0.0),
468  Z0 = cms.double(0.0)
469 )
470 ZeroTeslaRun247324CollisionVtxSmearingParameters = cms.PSet(
471  Phi = cms.double(0.0),
472  BetaStar = cms.double(80.0),
473  Emittance = cms.double(1.070e-5),
474  Alpha = cms.double(0.0),
475  SigmaZ = cms.double(4.125),
476  TimeOffset = cms.double(0.0),
477  X0 = cms.double(0.08621),
478  Y0 = cms.double(0.1657),
479  Z0 = cms.double(-1.688)
480 )
481 
482 # From 2015A 0T data
483 # Centroid absolute positions extracted from fills:
484 # X = 0.059395 cm
485 # Y = 0.099686 cm
486 # Z = -1.722240 cm
487 #
488 # BPIX absolute position extracted from first collision alignment:
489 # X = -0.0259503 cm
490 # Y = -0.07004 cm
491 # Z = -0.498917 cm
492 Realistic50ns13TeVCollisionZeroTeslaVtxSmearingParameters = cms.PSet(
493  Phi = cms.double(0.0),
494  BetaStar = cms.double(65.0),
495  Emittance = cms.double(5.411e-08),
496  Alpha = cms.double(0.0),
497  SigmaZ = cms.double(5.3),
498  TimeOffset = cms.double(0.0),
499  X0 = cms.double(0.08533),
500  Y0 = cms.double(0.16973),
501  Z0 = cms.double(-1.2230)
502 )
503 
504 # From 2015B 3.8T data
505 # Centroid absolute positions extracted from fill 4008:
506 # X = 0.07798 cm
507 # Y = 0.09714 cm
508 # Z = -1.610 cm
509 #
510 # BPIX absolute position extracted from PCL-like alignment run after magnet ramp-up:
511 # X = -0.026837 cm
512 # Y = -0.0715252 cm
513 # Z = -0.511453 cm
514 Realistic50ns13TeVCollisionVtxSmearingParameters = cms.PSet(
515  Phi = cms.double(0.0),
516  BetaStar = cms.double(65.0),
517  Emittance = cms.double(5.411e-08),
518  Alpha = cms.double(0.0),
519  SigmaZ = cms.double(5.3),
520  TimeOffset = cms.double(0.0),
521  X0 = cms.double(0.10482),
522  Y0 = cms.double(0.16867),
523  Z0 = cms.double(-1.0985)
524 )
525 
526 # From 2015B 3.8T data, beta*=90m (700 bunches fills)
527 # Centroid absolute positions extracted from 700 bunches fills 4499-4511:
528 # X = 0.068357 cm
529 # Y = 0.109159 cm
530 # Z = 0.131811 cm
531 #
532 # BPIX absolute position extracted from Prompt Reco alignment of run 259352
533 # X = -0.041651 cm
534 # Y = -0.199279 cm
535 # Z = -0.565093 cm
536 #
537 # Emittance has been calculated to match a BeamWidht of O(10um) with: https://lpc.web.cern.ch/lumi2.html
538 #
539 Realistic100ns13TeVCollisionBetaStar90mVtxSmearingParameters = cms.PSet(
540  Phi = cms.double(0.0),
541  BetaStar = cms.double(9121.0),
542  Emittance = cms.double(0.12e-7),
543  Alpha = cms.double(0.0),
544  SigmaZ = cms.double(4.9),
545  TimeOffset = cms.double(0.0),
546  X0 = cms.double(0.11000),
547  Y0 = cms.double(0.30844),
548  Z0 = cms.double(0.69690)
549 )
550 
551 # From 2015B 3.8T data, beta*=90m (42/240 bunches fills)
552 # Centroid absolute positions extracted from 42/240 bunches fills 4495-4496:
553 # X = 0.064925 cm
554 # Y = 0.112761 cm
555 # Z = 0.170413 cm
556 #
557 # BPIX absolute position extracted from Prompt Reco alignment of run 259202
558 # X = -0.041651 cm
559 # Y = -0.199279 cm
560 # Z = -0.565093 cm
561 #
562 # Emittance has been calculated to match a BeamWidht of O(10um) with: https://lpc.web.cern.ch/lumi2.html
563 #
564 Realistic100ns13TeVCollisionBetaStar90mLowBunchesVtxSmearingParameters = cms.PSet(
565  Phi = cms.double(0.0),
566  BetaStar = cms.double(9121.0),
567  Emittance = cms.double(0.12e-7),
568  Alpha = cms.double(0.0),
569  SigmaZ = cms.double(5.24),
570  TimeOffset = cms.double(0.0),
571  X0 = cms.double(0.10658),
572  Y0 = cms.double(0.31204),
573  Z0 = cms.double(0.735506)
574 )
575 
576 # From 2016B 3.8T data
577 # BS parameters extracted from fills 4895 - 4935:
578 # X0 = 0.064870 +/- 4.9575E-07 [cm]
579 # Y0 = 0.093639 +/- 4.9233E-07 [cm]
580 # Z0 = 0.420085 +/- 4.1102E-04 [cm]
581 # sigmaZ0 = 3.645533 +/- 2.9064E-04 [cm]
582 #
583 # From LHC calculator, emittance is 4.906e-8 cm
584 # https://lpc.web.cern.ch/lpc/lumi2.html
585 #
586 # BPIX absolute position:
587 # X = -0.0267572 cm
588 # Y = -0.0759102 cm
589 # Z = -0.511428 cm
590 Realistic25ns13TeV2016CollisionVtxSmearingParameters = cms.PSet(
591  Phi = cms.double(0.0),
592  BetaStar = cms.double(40.0),
593  Emittance = cms.double(4.906e-8),
594  Alpha = cms.double(0.0),
595  SigmaZ = cms.double(3.65),
596  TimeOffset = cms.double(0.0),
597  X0 = cms.double(0.09163),
598  Y0 = cms.double(0.16955),
599  Z0 = cms.double(0.9315 )
600 )
601 
602 # From 2017A 3.8T data
603 # BS parameters extracted from run 295463 (from offline DQM, i.e. PCL):
604 # X0 = 0.08497 [cm]
605 # Y0 = -0.03976 [cm]
606 # Z0 = 1.6 [cm] ==> 0.5 adjusted after cogging tuning by LHC, see https://hypernews.cern.ch/HyperNews/CMS/get/beamspot/159/1.html
607 # sigmaZ0 = 3.5 [cm]
608 #
609 # From LHC calculator, emittance is 3.319e-8 cm
610 # https://lpc.web.cern.ch/lpc/lumi2.html
611 #
612 # BPIX absolute position (https://hypernews.cern.ch/HyperNews/CMS/get/tif-alignment/657/1/1.html):
613 # X = 0.109725 cm
614 # Y = -0.108993 cm
615 # Z = -0.32054 cm
616 Realistic25ns13TeVEarly2017CollisionVtxSmearingParameters = cms.PSet(
617  Phi = cms.double(0.0),
618  BetaStar = cms.double(40.0),
619  Emittance = cms.double(3.319e-8),
620  Alpha = cms.double(0.0),
621  SigmaZ = cms.double(3.5),
622  TimeOffset = cms.double(0.0),
623  X0 = cms.double(-0.024755),
624  Y0 = cms.double(0.069233 ),
625  Z0 = cms.double(0.82054 )
626 )
627 
628 # Beam spot extracted from data for 2017 pp run @ 5 TeV
629 Realistic5TeVppCollision2017VtxSmearingParameters = cms.PSet(
630  Phi = cms.double(0.0),
631  BetaStar = cms.double(311),
632  Emittance = cms.double(3.8e-8),
633  Alpha = cms.double(0.0),
634  SigmaZ = cms.double(5.82),
635  TimeOffset = cms.double(0.0),
636  X0 = cms.double(-0.0228),
637  Y0 = cms.double(0.0795),
638  Z0 = cms.double(0.619)
639 )
640 
641 # From 2018B 3.8T data
642 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
643 # X0 = 0.09676 [cm]
644 # Y0 = -0.06245 [cm]
645 # Z0 = -0.292 [cm]
646 # sigmaZ0 = 3.5 [cm] => mean sigmaZ0 in this run is 3.2676
647 # BeamWidthX 0.0008050
648 # BeamWidthY 0.0006238
649 #
650 # From LHC calculator, emittance is 1.634e-8 cm
651 # https://lpc.web.cern.ch/lpc/lumi2.html
652 #
653 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
654 # X = 0.0859918 cm
655 # Y = -0.104172 cm
656 # Z = -0.327748 cm
657 Realistic25ns13TeVEarly2018CollisionVtxSmearingParameters = cms.PSet(
658  Phi = cms.double(0.0),
659  BetaStar = cms.double(30.0),
660  Emittance = cms.double(1.634e-8),
661  Alpha = cms.double(0.0),
662  SigmaZ = cms.double(3.5),
663  TimeOffset = cms.double(0.0),
664  X0 = cms.double(0.0107682),
665  Y0 = cms.double(0.041722 ),
666  Z0 = cms.double(0.035748 )
667 )
668 
669 # Run3 possible beam parameters
670 # Round optics - Low SigmaZ = 3.4 cm
671 # From 2018B 3.8T data
672 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
673 # X0 = 0.09676 [cm]
674 # Y0 = -0.06245 [cm]
675 # Z0 = -0.292 [cm]
676 # sigmaZ0 = 3.2676 [cm]
677 # BeamWidthX 0.0008050
678 # BeamWidthY 0.0006238
679 #
680 # set SigmaZ0 = 3.4 [cm]
681 # set BeamWidthX = BeamWidthY = 11.5 [um]
682 # set beta* = 28 cm
683 # energy = 13 TeV
684 # From LHC calculator, emittance is 4.762e-8 cm
685 # https://lpc.web.cern.ch/lpc/lumi2.html
686 #
687 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
688 # X = 0.0859918 cm
689 # Y = -0.104172 cm
690 # Z = -0.327748 cm
691 Run3RoundOptics25ns13TeVLowSigmaZVtxSmearingParameters = cms.PSet(
692  Phi = cms.double(0.0),
693  BetaStar = cms.double(28.0),
694  Emittance = cms.double(4.762e-8),
695  Alpha = cms.double(0.0),
696  SigmaZ = cms.double(3.4),
697  TimeOffset = cms.double(0.0),
698  X0 = cms.double(0.0107682),
699  Y0 = cms.double(0.041722 ),
700  Z0 = cms.double(0.035748 )
701 )
702 
703 # Run3 possible beam parameters
704 # Round optics - High SigmaZ = 5.7 cm
705 # From 2018B 3.8T data
706 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
707 # X0 = 0.09676 [cm]
708 # Y0 = -0.06245 [cm]
709 # Z0 = -0.292 [cm]
710 # sigmaZ0 = 3.2676 [cm]
711 # BeamWidthX 0.0008050
712 # BeamWidthY 0.0006238
713 #
714 # set SigmaZ0 = 5.7 [cm]
715 # set BeamWidthX = BeamWidthY = 11.5 [um]
716 # set beta* = 28 cm
717 # energy = 13 TeV
718 # From LHC calculator, emittance is 4.762e-8 cm
719 # https://lpc.web.cern.ch/lpc/lumi2.html
720 #
721 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
722 # X = 0.0859918 cm
723 # Y = -0.104172 cm
724 # Z = -0.327748 cm
725 Run3RoundOptics25ns13TeVHighSigmaZVtxSmearingParameters = cms.PSet(
726  Phi = cms.double(0.0),
727  BetaStar = cms.double(28.0),
728  Emittance = cms.double(4.762e-8),
729  Alpha = cms.double(0.0),
730  SigmaZ = cms.double(5.7),
731  TimeOffset = cms.double(0.0),
732  X0 = cms.double(0.0107682),
733  Y0 = cms.double(0.041722 ),
734  Z0 = cms.double(0.035748 )
735 )
736 
737 # Test HF offset
738 ShiftedCollision2015VtxSmearingParameters = cms.PSet(
739  Phi = cms.double(0.0),
740  BetaStar = cms.double(65.0),
741  Emittance = cms.double(5.411e-08),
742  Alpha = cms.double(0.0),
743  SigmaZ = cms.double(5.3),
744  TimeOffset = cms.double(0.0),
745  X0 = cms.double(1.0),
746  Y0 = cms.double(0.0),
747  Z0 = cms.double(0.0)
748 )
749 Shifted5mmCollision2015VtxSmearingParameters = cms.PSet(
750  Phi = cms.double(0.0),
751  BetaStar = cms.double(65.0),
752  Emittance = cms.double(5.411e-08),
753  Alpha = cms.double(0.0),
754  SigmaZ = cms.double(5.3),
755  TimeOffset = cms.double(0.0),
756  X0 = cms.double(0.5),
757  Y0 = cms.double(0.0),
758  Z0 = cms.double(0.0)
759 )
760 Shifted15mmCollision2015VtxSmearingParameters = cms.PSet(
761  Phi = cms.double(0.0),
762  BetaStar = cms.double(65.0),
763  Emittance = cms.double(5.411e-08),
764  Alpha = cms.double(0.0),
765  SigmaZ = cms.double(5.3),
766  TimeOffset = cms.double(0.0),
767  X0 = cms.double(1.5),
768  Y0 = cms.double(0.0),
769  Z0 = cms.double(0.0)
770 )
771 
772 # Estimate for 2015 PbPb collisions, based on feedback from accelerator
773 # Beamspot centroid shifted to match pp expectation for 2015
774 NominalHICollision2015VtxSmearingParameters = cms.PSet(
775  Phi = cms.double(0.0),
776  BetaStar = cms.double(60.0),
777  Emittance = cms.double(1.70e-07),
778  Alpha = cms.double(0.0),
779  SigmaZ = cms.double(7.06),
780  TimeOffset = cms.double(0.0),
781  X0 = cms.double(0.0322),
782  Y0 = cms.double(0.),
783  Z0 = cms.double(0.)
784 )
785 
786 # updated numbers based on beamspot fits to 2015 PbPb data
787 # Later found to be incorrect, see following entry
788 # RealisticHICollision2015VtxSmearingParameters = cms.PSet(
789 # Phi = cms.double(0.0),
790 # BetaStar = cms.double(60.0),
791 # Emittance = cms.double(1.70e-07),
792 # Alpha = cms.double(0.0),
793 # SigmaZ = cms.double(5.2278),
794 # TimeOffset = cms.double(0.0),
795 # X0 = cms.double(0.1025),
796 # Y0 = cms.double(0.1654),
797 # Z0 = cms.double(3.2528)
798 #)
799 # updated numbers for 2015 PbPb data with Z centroid from fixed beamspot fits
800 # See discussion here https://hypernews.cern.ch/HyperNews/CMS/get/hi-general/3968.html
801 # See plot of difference here: https://www.dropbox.com/s/tsnkgvvpkdqjtyq/vzDataMCOverlay_c_20170420.pdf?dl=0
802 #
803 RealisticHICollisionFixZ2015VtxSmearingParameters = cms.PSet(
804  Phi = cms.double(0.0),
805  BetaStar = cms.double(60.0),
806  Emittance = cms.double(1.70e-07),
807  Alpha = cms.double(0.0),
808  SigmaZ = cms.double(5.2278),
809  TimeOffset = cms.double(0.0),
810  X0 = cms.double(0.1025),
811  Y0 = cms.double(0.1654),
812  Z0 = cms.double(0.771)
813 )
814 # Numbers based on beamspot fits to 2017 XeXe data
815 # Documentation here: https://twiki.cern.ch/twiki/pub/CMS/XeXeRereco/IanLRU_AlCaTkAlBS_20171130_approvedByLucaSara.pdf
816 RealisticXeXeCollision2017VtxSmearingParameters = cms.PSet(
817  Phi = cms.double(0.0),
818  BetaStar = cms.double(30.0),
819  Emittance = cms.double(4.33e-08),
820  Alpha = cms.double(0.0),
821  SigmaZ = cms.double(4.64),
822  TimeOffset = cms.double(0.0),
823  X0 = cms.double(-0.026),
824  Y0 = cms.double(0.081),
825  Z0 = cms.double(0.645)
826 )
827 
828 # From fit to 5 TeV PbPb data
829 # From 2018 PbPb data
830 # BS parameters extracted from run 327211, Fill 7471 (from StreamExpressAlignment, HP, BS):
831 # X0 = 0.09443 [cm]
832 # Y0 = -0.06377 [cm]
833 # Z0 = 0.58067 [cm]
834 # sigmaZ0 = 4.969 [cm]
835 # BeamWidthX 0.0014392
836 # BeamWidthY 0.0011545
837 #
838 # Emittance is 3.36e-8 cm. Calculated by ((BeamWidthX + BeamWidth)/2)^2/BetaStar
839 #
840 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
841 # X = 0.084335 cm
842 # Y = -0.110381 cm
843 # Z = -0.321904 cm
844 RealisticPbPbCollision2018VtxSmearingParameters = cms.PSet(
845  Phi = cms.double(0.0),
846  BetaStar = cms.double(50),
847  Emittance = cms.double(3.36e-08),
848  Alpha = cms.double(0.0),
849  SigmaZ = cms.double(4.97),
850  TimeOffset = cms.double(0.0),
851  X0 = cms.double(0.010),
852  Y0 = cms.double(0.047),
853  Z0 = cms.double(0.903)
854 )
855 
856 # Estimate for 2015 pp collisions at 5.02 TeV, based on feedback from accelerator: beta* ~ 400cm, normalized emittance = 2.5 um, SigmaZ similar to RunIIWinter15GS
857 Nominal5TeVpp2015VtxSmearingParameters = cms.PSet(
858  Phi = cms.double(0.0),
859  BetaStar = cms.double(400.0),
860  Emittance = cms.double(1.0e-07),
861  Alpha = cms.double(0.0),
862  SigmaZ = cms.double(5.5),
863  TimeOffset = cms.double(0.0),
864  X0 = cms.double(0.1044),
865  Y0 = cms.double(0.1676),
866  Z0 = cms.double(0.6707)
867 )
868 
869 # From fit to 5 TeV pPb data
870 Realistic5TeVPACollision2016VtxSmearingParameters = cms.PSet(
871  Phi = cms.double(0.0),
872  BetaStar = cms.double(1100.0),
873  Emittance = cms.double(6.75e-08),
874  Alpha = cms.double(0.0),
875  SigmaZ = cms.double(6.4891),
876  TimeOffset = cms.double(0.0),
877  X0 = cms.double(0.0889),
878  Y0 = cms.double(0.1820),
879  Z0 = cms.double(1.6066)
880 )
881 
882 # From fit to 8 TeV pPb data
883 Realistic8TeVPACollision2016VtxSmearingParameters = cms.PSet(
884  Phi = cms.double(0.0),
885  BetaStar = cms.double(60.0),
886  Emittance = cms.double(6.75e-08),
887  Alpha = cms.double(0.0),
888  SigmaZ = cms.double(4.6914),
889  TimeOffset = cms.double(0.0),
890  X0 = cms.double(0.0836),
891  Y0 = cms.double(0.1837),
892  Z0 = cms.double(1.3577)
893 )
894 
895 # Parameters for HL-LHC operation at 13TeV
896 HLLHCVtxSmearingParameters = cms.PSet(
897  MeanXIncm = cms.double(0.),
898  MeanYIncm = cms.double(0.),
899  MeanZIncm = cms.double(0.),
900  TimeOffsetInns = cms.double(0.0),
901  EprotonInGeV = cms.double(6500.0),
902  CrossingAngleInurad = cms.double(510.0),
903  CrabFrequencyInMHz = cms.double(400.0),
904  RF800 = cms.bool(False),
905  BetaCrossingPlaneInm = cms.double(0.20),
906  BetaSeparationPlaneInm = cms.double(0.20),
907  HorizontalEmittance = cms.double(2.5e-06),
908  VerticalEmittance = cms.double(2.05e-06),
909  BunchLengthInm = cms.double(0.090),
910  CrabbingAngleCrossingInurad = cms.double(380.0),
911  CrabbingAngleSeparationInurad = cms.double(0.0)
912 )
913 
914 # Parameters for HL-LHC Crab-kissing operation 13 TeV
915 HLLHCCrabKissingVtxSmearingParameters = cms.PSet(
916  MeanXIncm = cms.double(0.),
917  MeanYIncm = cms.double(0.),
918  MeanZIncm = cms.double(0.),
919  TimeOffsetInns = cms.double(0.0),
920  EprotonInGeV = cms.double(6500.0),
921  HalfCrossingAngleInurad = cms.double(200.0),
922  CrabAngleCrossingPlaneInurad = cms.double(200.0),
923  CrabFrequencyCrossingPlaneInMHz = cms.double(400.0),
924  NormalizedEmittanceCrossingPlaneInum = cms.double(2.5),
925  BetaStarCrossingPlaneInm = cms.double(0.30),
926  CrabAngleParallelPlaneInurad = cms.double(100.0),
927  CrabFrequencyParallelPlaneInMHz = cms.double(400.0),
928  NormalizedEmittanceParallelPlaneInum = cms.double(2.5),
929  BetaStarParallelPlaneInm = cms.double(0.075),
930  ZsizeInm = cms.double(0.15),
931  BeamProfile=cms.string("Flat")
932 )