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 # Test HF offset
598 ShiftedCollision2015VtxSmearingParameters = cms.PSet(
599  Phi = cms.double(0.0),
600  BetaStar = cms.double(65.0),
601  Emittance = cms.double(5.411e-08),
602  Alpha = cms.double(0.0),
603  SigmaZ = cms.double(5.3),
604  TimeOffset = cms.double(0.0),
605  X0 = cms.double(1.0),
606  Y0 = cms.double(0.0),
607  Z0 = cms.double(0.0)
608 )
609 Shifted5mmCollision2015VtxSmearingParameters = cms.PSet(
610  Phi = cms.double(0.0),
611  BetaStar = cms.double(65.0),
612  Emittance = cms.double(5.411e-08),
613  Alpha = cms.double(0.0),
614  SigmaZ = cms.double(5.3),
615  TimeOffset = cms.double(0.0),
616  X0 = cms.double(0.5),
617  Y0 = cms.double(0.0),
618  Z0 = cms.double(0.0)
619 )
620 Shifted15mmCollision2015VtxSmearingParameters = cms.PSet(
621  Phi = cms.double(0.0),
622  BetaStar = cms.double(65.0),
623  Emittance = cms.double(5.411e-08),
624  Alpha = cms.double(0.0),
625  SigmaZ = cms.double(5.3),
626  TimeOffset = cms.double(0.0),
627  X0 = cms.double(1.5),
628  Y0 = cms.double(0.0),
629  Z0 = cms.double(0.0)
630 )
631 
632 # Estimate for 2015 PbPb collisions, based on feedback from accelerator
633 # Beamspot centroid shifted to match pp expectation for 2015
634 NominalHICollision2015VtxSmearingParameters = cms.PSet(
635  Phi = cms.double(0.0),
636  BetaStar = cms.double(60.0),
637  Emittance = cms.double(1.70e-07),
638  Alpha = cms.double(0.0),
639  SigmaZ = cms.double(7.06),
640  TimeOffset = cms.double(0.0),
641  X0 = cms.double(0.0322),
642  Y0 = cms.double(0.),
643  Z0 = cms.double(0.)
644 )
645 
646 # updated numbers based on beamspot fits to 2015 PbPb data
647 # Later found to be incorrect, see following entry
648 # RealisticHICollision2015VtxSmearingParameters = cms.PSet(
649 # Phi = cms.double(0.0),
650 # BetaStar = cms.double(60.0),
651 # Emittance = cms.double(1.70e-07),
652 # Alpha = cms.double(0.0),
653 # SigmaZ = cms.double(5.2278),
654 # TimeOffset = cms.double(0.0),
655 # X0 = cms.double(0.1025),
656 # Y0 = cms.double(0.1654),
657 # Z0 = cms.double(3.2528)
658 #)
659 # updated numbers for 2015 PbPb data with Z centroid from fixed beamspot fits
660 # See discussion here https://hypernews.cern.ch/HyperNews/CMS/get/hi-general/3968.html
661 # See plot of difference here: https://www.dropbox.com/s/tsnkgvvpkdqjtyq/vzDataMCOverlay_c_20170420.pdf?dl=0
662 #
663 RealisticHICollisionFixZ2015VtxSmearingParameters = cms.PSet(
664  Phi = cms.double(0.0),
665  BetaStar = cms.double(60.0),
666  Emittance = cms.double(1.70e-07),
667  Alpha = cms.double(0.0),
668  SigmaZ = cms.double(5.2278),
669  TimeOffset = cms.double(0.0),
670  X0 = cms.double(0.1025),
671  Y0 = cms.double(0.1654),
672  Z0 = cms.double(0.771)
673 )
674 # Numbers based on beamspot fits to 2017 XeXe data
675 # Documentation here: https://twiki.cern.ch/twiki/pub/CMS/XeXeRereco/IanLRU_AlCaTkAlBS_20171130_approvedByLucaSara.pdf
676 RealisticXeXeCollision2017VtxSmearingParameters = cms.PSet(
677  Phi = cms.double(0.0),
678  BetaStar = cms.double(30.0),
679  Emittance = cms.double(4.33e-08),
680  Alpha = cms.double(0.0),
681  SigmaZ = cms.double(4.64),
682  TimeOffset = cms.double(0.0),
683  X0 = cms.double(-0.026),
684  Y0 = cms.double(0.081),
685  Z0 = cms.double(0.645)
686 )
687 
688 # 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
689 Nominal5TeVpp2015VtxSmearingParameters = cms.PSet(
690  Phi = cms.double(0.0),
691  BetaStar = cms.double(400.0),
692  Emittance = cms.double(1.0e-07),
693  Alpha = cms.double(0.0),
694  SigmaZ = cms.double(5.5),
695  TimeOffset = cms.double(0.0),
696  X0 = cms.double(0.1044),
697  Y0 = cms.double(0.1676),
698  Z0 = cms.double(0.6707)
699 )
700 
701 # From fit to 5 TeV pPb data
702 Realistic5TeVPACollision2016VtxSmearingParameters = cms.PSet(
703  Phi = cms.double(0.0),
704  BetaStar = cms.double(1100.0),
705  Emittance = cms.double(6.75e-08),
706  Alpha = cms.double(0.0),
707  SigmaZ = cms.double(6.4891),
708  TimeOffset = cms.double(0.0),
709  X0 = cms.double(0.0889),
710  Y0 = cms.double(0.1820),
711  Z0 = cms.double(1.6066)
712 )
713 
714 # From fit to 8 TeV pPb data
715 Realistic8TeVPACollision2016VtxSmearingParameters = cms.PSet(
716  Phi = cms.double(0.0),
717  BetaStar = cms.double(60.0),
718  Emittance = cms.double(6.75e-08),
719  Alpha = cms.double(0.0),
720  SigmaZ = cms.double(4.6914),
721  TimeOffset = cms.double(0.0),
722  X0 = cms.double(0.0836),
723  Y0 = cms.double(0.1837),
724  Z0 = cms.double(1.3577)
725 )
726 
727 # Parameters for HL-LHC operation at 13TeV
728 HLLHCVtxSmearingParameters = cms.PSet(
729  MeanXIncm = cms.double(0.),
730  MeanYIncm = cms.double(0.),
731  MeanZIncm = cms.double(0.),
732  TimeOffsetInns = cms.double(0.0),
733  EprotonInGeV = cms.double(6500.0),
734  CrossingAngleInurad = cms.double(510.0),
735  CrabFrequencyInMHz = cms.double(400.0),
736  RF800 = cms.bool(False),
737  BetaCrossingPlaneInm = cms.double(0.20),
738  BetaSeparationPlaneInm = cms.double(0.20),
739  HorizontalEmittance = cms.double(2.5e-06),
740  VerticalEmittance = cms.double(2.05e-06),
741  BunchLengthInm = cms.double(0.090),
742  CrabbingAngleCrossingInurad = cms.double(380.0),
743  CrabbingAngleSeparationInurad = cms.double(0.0)
744 )
745 
746 # Parameters for HL-LHC Crab-kissing operation 13 TeV
747 HLLHCCrabKissingVtxSmearingParameters = cms.PSet(
748  MeanXIncm = cms.double(0.),
749  MeanYIncm = cms.double(0.),
750  MeanZIncm = cms.double(0.),
751  TimeOffsetInns = cms.double(0.0),
752  EprotonInGeV = cms.double(6500.0),
753  HalfCrossingAngleInurad = cms.double(200.0),
754  CrabAngleCrossingPlaneInurad = cms.double(200.0),
755  CrabFrequencyCrossingPlaneInMHz = cms.double(400.0),
756  NormalizedEmittanceCrossingPlaneInum = cms.double(2.5),
757  BetaStarCrossingPlaneInm = cms.double(0.30),
758  CrabAngleParallelPlaneInurad = cms.double(100.0),
759  CrabFrequencyParallelPlaneInMHz = cms.double(400.0),
760  NormalizedEmittanceParallelPlaneInum = cms.double(2.5),
761  BetaStarParallelPlaneInm = cms.double(0.075),
762  ZsizeInm = cms.double(0.15),
763  BeamProfile=cms.string("Flat")
764 )