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