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 
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 # Fixed Emittance (X2) in Beam spot extracted from data for 2017 pp run @ 5 TeV
642 Fixed_EmitRealistic5TeVppCollision2017VtxSmearingParameters = cms.PSet(
643  Phi = cms.double(0.0),
644  BetaStar = cms.double(311),
645  Emittance = cms.double(7.6e-8),
646  Alpha = cms.double(0.0),
647  SigmaZ = cms.double(5.82),
648  TimeOffset = cms.double(0.0),
649  X0 = cms.double(-0.0228),
650  Y0 = cms.double(0.0795),
651  Z0 = cms.double(0.619)
652 )
653 
654 
655 # From 2018B 3.8T data
656 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
657 # X0 = 0.09676 [cm]
658 # Y0 = -0.06245 [cm]
659 # Z0 = -0.292 [cm]
660 # sigmaZ0 = 3.5 [cm] => mean sigmaZ0 in this run is 3.2676
661 # BeamWidthX 0.0008050
662 # BeamWidthY 0.0006238
663 #
664 # From LHC calculator, emittance is 1.634e-8 cm
665 # https://lpc.web.cern.ch/lpc/lumi2.html
666 #
667 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
668 # X = 0.0859918 cm
669 # Y = -0.104172 cm
670 # Z = -0.327748 cm
671 Realistic25ns13TeVEarly2018CollisionVtxSmearingParameters = cms.PSet(
672  Phi = cms.double(0.0),
673  BetaStar = cms.double(30.0),
674  Emittance = cms.double(1.634e-8),
675  Alpha = cms.double(0.0),
676  SigmaZ = cms.double(3.5),
677  TimeOffset = cms.double(0.0),
678  X0 = cms.double(0.0107682),
679  Y0 = cms.double(0.041722 ),
680  Z0 = cms.double(0.035748 )
681 )
682 
683 # Run3 possible beam parameters
684 # Round optics - Low SigmaZ = 3.4 cm
685 # From 2018B 3.8T data
686 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
687 # X0 = 0.09676 [cm]
688 # Y0 = -0.06245 [cm]
689 # Z0 = -0.292 [cm]
690 # sigmaZ0 = 3.2676 [cm]
691 # BeamWidthX 0.0008050
692 # BeamWidthY 0.0006238
693 #
694 # set SigmaZ0 = 3.4 [cm]
695 # set BeamWidthX = BeamWidthY = 11.5 [um]
696 # set beta* = 28 cm
697 # energy = 13 TeV
698 # From LHC calculator, emittance is 4.762e-8 cm
699 # https://lpc.web.cern.ch/lpc/lumi2.html
700 #
701 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
702 # X = 0.0859918 cm
703 # Y = -0.104172 cm
704 # Z = -0.327748 cm
705 Run3RoundOptics25ns13TeVLowSigmaZVtxSmearingParameters = cms.PSet(
706  Phi = cms.double(0.0),
707  BetaStar = cms.double(28.0),
708  Emittance = cms.double(4.762e-8),
709  Alpha = cms.double(0.0),
710  SigmaZ = cms.double(3.4),
711  TimeOffset = cms.double(0.0),
712  X0 = cms.double(0.0107682),
713  Y0 = cms.double(0.041722 ),
714  Z0 = cms.double(0.035748 )
715 )
716 
717 # Run3 possible beam parameters
718 # Round optics - High SigmaZ = 5.7 cm
719 # From 2018B 3.8T data
720 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
721 # X0 = 0.09676 [cm]
722 # Y0 = -0.06245 [cm]
723 # Z0 = -0.292 [cm]
724 # sigmaZ0 = 3.2676 [cm]
725 # BeamWidthX 0.0008050
726 # BeamWidthY 0.0006238
727 #
728 # set SigmaZ0 = 5.7 [cm]
729 # set BeamWidthX = BeamWidthY = 11.5 [um]
730 # set beta* = 28 cm
731 # energy = 13 TeV
732 # From LHC calculator, emittance is 4.762e-8 cm
733 # https://lpc.web.cern.ch/lpc/lumi2.html
734 #
735 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
736 # X = 0.0859918 cm
737 # Y = -0.104172 cm
738 # Z = -0.327748 cm
739 Run3RoundOptics25ns13TeVHighSigmaZVtxSmearingParameters = cms.PSet(
740  Phi = cms.double(0.0),
741  BetaStar = cms.double(28.0),
742  Emittance = cms.double(4.762e-8),
743  Alpha = cms.double(0.0),
744  SigmaZ = cms.double(5.7),
745  TimeOffset = cms.double(0.0),
746  X0 = cms.double(0.0107682),
747  Y0 = cms.double(0.041722 ),
748  Z0 = cms.double(0.035748 )
749 )
750 
751 # From Run3 PilotBeams 2021 3.8T data
752 # BS parameters extracted from run 346512, fill 7531 (from ExpressPhysics FEVT, Legacy BS):
753 # X0 = 0.174282 [cm]
754 # Y0 = -0.187132 [cm]
755 # Z0 = 0.167616 [cm]
756 # sigmaZ0 = 6.80728 [cm]
757 # BeamWidthX 0.0142174
758 # BeamWidthY 0.0150789
759 #
760 # set SigmaZ0 = 6.8 [cm]
761 # set BeamWidthX = BeamWidthY = 150.0 [um]
762 # set beta* = 1100 cm
763 # energy = 900 GeV
764 # From LHC calculator, emittance is 4.762e-8 cm
765 # https://lpc.web.cern.ch/lpc/lumi2.html
766 #
767 # BPIX absolute position (https://twiki.cern.ch/twiki/bin/view/CMS/TkAlignmentPixelPosition?rev=40#2021):
768 # X = 0.06076 cm
769 # Y = -0.14702 cm
770 # Z = -0.25616 cm
771 Realistic25ns900GeV2021CollisionVtxSmearingParameters = cms.PSet(
772  Phi = cms.double(0.0),
773  BetaStar = cms.double(1100.0),
774  Emittance = cms.double(4.169e-7),
775  Alpha = cms.double(0.0),
776  SigmaZ = cms.double(6.8),
777  TimeOffset = cms.double(0.0),
778  X0 = cms.double(0.113522),
779  Y0 = cms.double(-0.040112),
780  Z0 = cms.double(0.423776)
781 )
782 
783 # From first Run 3 data at 13.6 TeV and 3.8T
784 # BS parameters extracted from run 355100, fill 7920:
785 # X0 = 0.172394 [cm]
786 # Y0 = -0.180946 [cm]
787 # Z0 = 0.94181 [cm]
788 # sigmaZ0 = 3.81941 [cm]
789 # BeamWidthX = 0.0008772 [cm]
790 # BeamWidthY = 0.0010078 [cm]
791 #
792 # set SigmaZ0 = 3.8 [cm]
793 # set BeamWidthX = BeamWidthY = 10.0 [um]
794 # set beta* = 30 cm
795 # energy = 13.6 TeV
796 # From LHC calculator, emittance is 6.621e-8 cm
797 # https://lpc.web.cern.ch/lumiCalc.html
798 #
799 # BPIX absolute position (https://twiki.cern.ch/twiki/bin/view/CMS/TkAlignmentPixelPosition?rev=45#Collisions_at_s_13_6_TeV):
800 # X = 0.0717651 cm
801 # Y = -0.165951 cm
802 # Z = -0.356345 cm
803 Realistic25ns13p6TeVEarly2022CollisionVtxSmearingParameters = cms.PSet(
804  Phi = cms.double(0.0),
805  BetaStar = cms.double(30.0),
806  Emittance = cms.double(6.621e-8),
807  Alpha = cms.double(0.0),
808  SigmaZ = cms.double(3.8),
809  TimeOffset = cms.double(0.0),
810  X0 = cms.double(0.100629),
811  Y0 = cms.double(-0.014995),
812  Z0 = cms.double(1.298155)
813 )
814 
815 # Test HF offset
816 ShiftedCollision2015VtxSmearingParameters = cms.PSet(
817  Phi = cms.double(0.0),
818  BetaStar = cms.double(65.0),
819  Emittance = cms.double(5.411e-08),
820  Alpha = cms.double(0.0),
821  SigmaZ = cms.double(5.3),
822  TimeOffset = cms.double(0.0),
823  X0 = cms.double(1.0),
824  Y0 = cms.double(0.0),
825  Z0 = cms.double(0.0)
826 )
827 Shifted5mmCollision2015VtxSmearingParameters = cms.PSet(
828  Phi = cms.double(0.0),
829  BetaStar = cms.double(65.0),
830  Emittance = cms.double(5.411e-08),
831  Alpha = cms.double(0.0),
832  SigmaZ = cms.double(5.3),
833  TimeOffset = cms.double(0.0),
834  X0 = cms.double(0.5),
835  Y0 = cms.double(0.0),
836  Z0 = cms.double(0.0)
837 )
838 Shifted15mmCollision2015VtxSmearingParameters = cms.PSet(
839  Phi = cms.double(0.0),
840  BetaStar = cms.double(65.0),
841  Emittance = cms.double(5.411e-08),
842  Alpha = cms.double(0.0),
843  SigmaZ = cms.double(5.3),
844  TimeOffset = cms.double(0.0),
845  X0 = cms.double(1.5),
846  Y0 = cms.double(0.0),
847  Z0 = cms.double(0.0)
848 )
849 
850 # Estimate for 2015 PbPb collisions, based on feedback from accelerator
851 # Beamspot centroid shifted to match pp expectation for 2015
852 NominalHICollision2015VtxSmearingParameters = cms.PSet(
853  Phi = cms.double(0.0),
854  BetaStar = cms.double(60.0),
855  Emittance = cms.double(1.70e-07),
856  Alpha = cms.double(0.0),
857  SigmaZ = cms.double(7.06),
858  TimeOffset = cms.double(0.0),
859  X0 = cms.double(0.0322),
860  Y0 = cms.double(0.),
861  Z0 = cms.double(0.)
862 )
863 
864 # updated numbers based on beamspot fits to 2015 PbPb data
865 # Later found to be incorrect, see following entry
866 # RealisticHICollision2015VtxSmearingParameters = cms.PSet(
867 # Phi = cms.double(0.0),
868 # BetaStar = cms.double(60.0),
869 # Emittance = cms.double(1.70e-07),
870 # Alpha = cms.double(0.0),
871 # SigmaZ = cms.double(5.2278),
872 # TimeOffset = cms.double(0.0),
873 # X0 = cms.double(0.1025),
874 # Y0 = cms.double(0.1654),
875 # Z0 = cms.double(3.2528)
876 #)
877 # updated numbers for 2015 PbPb data with Z centroid from fixed beamspot fits
878 # See discussion here https://hypernews.cern.ch/HyperNews/CMS/get/hi-general/3968.html
879 # See plot of difference here: https://www.dropbox.com/s/tsnkgvvpkdqjtyq/vzDataMCOverlay_c_20170420.pdf?dl=0
880 #
881 RealisticHICollisionFixZ2015VtxSmearingParameters = cms.PSet(
882  Phi = cms.double(0.0),
883  BetaStar = cms.double(60.0),
884  Emittance = cms.double(1.70e-07),
885  Alpha = cms.double(0.0),
886  SigmaZ = cms.double(5.2278),
887  TimeOffset = cms.double(0.0),
888  X0 = cms.double(0.1025),
889  Y0 = cms.double(0.1654),
890  Z0 = cms.double(0.771)
891 )
892 # Numbers based on beamspot fits to 2017 XeXe data
893 # Documentation here: https://twiki.cern.ch/twiki/pub/CMS/XeXeRereco/IanLRU_AlCaTkAlBS_20171130_approvedByLucaSara.pdf
894 RealisticXeXeCollision2017VtxSmearingParameters = cms.PSet(
895  Phi = cms.double(0.0),
896  BetaStar = cms.double(30.0),
897  Emittance = cms.double(4.33e-08),
898  Alpha = cms.double(0.0),
899  SigmaZ = cms.double(4.64),
900  TimeOffset = cms.double(0.0),
901  X0 = cms.double(-0.026),
902  Y0 = cms.double(0.081),
903  Z0 = cms.double(0.645)
904 )
905 
906 # From fit to 5 TeV PbPb data
907 # From 2018 PbPb data
908 # BS parameters extracted from run 327211, Fill 7471 (from StreamExpressAlignment, HP, BS):
909 # X0 = 0.09443 [cm]
910 # Y0 = -0.06377 [cm]
911 # Z0 = 0.58067 [cm]
912 # sigmaZ0 = 4.969 [cm]
913 # BeamWidthX 0.0014392
914 # BeamWidthY 0.0011545
915 #
916 # Emittance is 3.36e-8 cm. Calculated by ((BeamWidthX + BeamWidth)/2)^2/BetaStar
917 #
918 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
919 # X = 0.084335 cm
920 # Y = -0.110381 cm
921 # Z = -0.321904 cm
922 RealisticPbPbCollision2018VtxSmearingParameters = cms.PSet(
923  Phi = cms.double(0.0),
924  BetaStar = cms.double(50),
925  Emittance = cms.double(3.36e-08),
926  Alpha = cms.double(0.0),
927  SigmaZ = cms.double(4.97),
928  TimeOffset = cms.double(0.0),
929  X0 = cms.double(0.010),
930  Y0 = cms.double(0.047),
931  Z0 = cms.double(0.903)
932 )
933 
934 # 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
935 Nominal5TeVpp2015VtxSmearingParameters = cms.PSet(
936  Phi = cms.double(0.0),
937  BetaStar = cms.double(400.0),
938  Emittance = cms.double(1.0e-07),
939  Alpha = cms.double(0.0),
940  SigmaZ = cms.double(5.5),
941  TimeOffset = cms.double(0.0),
942  X0 = cms.double(0.1044),
943  Y0 = cms.double(0.1676),
944  Z0 = cms.double(0.6707)
945 )
946 
947 # From fit to 5 TeV pPb data
948 Realistic5TeVPACollision2016VtxSmearingParameters = cms.PSet(
949  Phi = cms.double(0.0),
950  BetaStar = cms.double(1100.0),
951  Emittance = cms.double(6.75e-08),
952  Alpha = cms.double(0.0),
953  SigmaZ = cms.double(6.4891),
954  TimeOffset = cms.double(0.0),
955  X0 = cms.double(0.0889),
956  Y0 = cms.double(0.1820),
957  Z0 = cms.double(1.6066)
958 )
959 
960 # From fit to 8 TeV pPb data
961 Realistic8TeVPACollision2016VtxSmearingParameters = cms.PSet(
962  Phi = cms.double(0.0),
963  BetaStar = cms.double(60.0),
964  Emittance = cms.double(6.75e-08),
965  Alpha = cms.double(0.0),
966  SigmaZ = cms.double(4.6914),
967  TimeOffset = cms.double(0.0),
968  X0 = cms.double(0.0836),
969  Y0 = cms.double(0.1837),
970  Z0 = cms.double(1.3577)
971 )
972 
973 # Parameters for HL-LHC operation at 13TeV
974 HLLHCVtxSmearingParameters = cms.PSet(
975  MeanXIncm = cms.double(0.),
976  MeanYIncm = cms.double(0.),
977  MeanZIncm = cms.double(0.),
978  TimeOffsetInns = cms.double(0.0),
979  EprotonInGeV = cms.double(6500.0),
980  CrossingAngleInurad = cms.double(510.0),
981  CrabFrequencyInMHz = cms.double(400.0),
982  RF800 = cms.bool(False),
983  BetaCrossingPlaneInm = cms.double(0.20),
984  BetaSeparationPlaneInm = cms.double(0.20),
985  HorizontalEmittance = cms.double(2.5e-06),
986  VerticalEmittance = cms.double(2.05e-06),
987  BunchLengthInm = cms.double(0.090),
988  CrabbingAngleCrossingInurad = cms.double(380.0),
989  CrabbingAngleSeparationInurad = cms.double(0.0)
990 )
991 
992 # Parameters for HL-LHC Crab-kissing operation 13 TeV
993 HLLHCCrabKissingVtxSmearingParameters = cms.PSet(
994  MeanXIncm = cms.double(0.),
995  MeanYIncm = cms.double(0.),
996  MeanZIncm = cms.double(0.),
997  TimeOffsetInns = cms.double(0.0),
998  EprotonInGeV = cms.double(6500.0),
999  HalfCrossingAngleInurad = cms.double(200.0),
1000  CrabAngleCrossingPlaneInurad = cms.double(200.0),
1001  CrabFrequencyCrossingPlaneInMHz = cms.double(400.0),
1002  NormalizedEmittanceCrossingPlaneInum = cms.double(2.5),
1003  BetaStarCrossingPlaneInm = cms.double(0.30),
1004  CrabAngleParallelPlaneInurad = cms.double(100.0),
1005  CrabFrequencyParallelPlaneInMHz = cms.double(400.0),
1006  NormalizedEmittanceParallelPlaneInum = cms.double(2.5),
1007  BetaStarParallelPlaneInm = cms.double(0.075),
1008  ZsizeInm = cms.double(0.15),
1009  BeamProfile=cms.string("Flat")
1010 )