Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-11-19 23:20:10

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 #
0004 # All units are cm and radians
0005 #
0006 # UNITS:
0007 #
0008 # TimeOffset in nanoseconds
0009 # spacial displacement in cm
0010 
0011 # common parameters
0012 VtxSmearedCommon = cms.PSet(
0013     src = cms.InputTag("generator", "unsmeared"),
0014     readDB = cms.bool(False)
0015 )
0016 # Gaussian smearing
0017 GaussVtxSmearingParameters = cms.PSet(
0018     MeanX = cms.double(0.0),
0019     MeanY = cms.double(0.0),
0020     MeanZ = cms.double(0.0),
0021     SigmaY = cms.double(0.0015),
0022     SigmaX = cms.double(0.0015),
0023     SigmaZ = cms.double(5.3),
0024     TimeOffset = cms.double(0.0)
0025 )
0026 # Gaussian smearing
0027 GaussVtxSigmaZ4cmSmearingParameters = cms.PSet(
0028     MeanX = cms.double(0.0),
0029     MeanY = cms.double(0.0),
0030     MeanZ = cms.double(0.0),
0031     SigmaY = cms.double(0.0015),
0032     SigmaX = cms.double(0.0015),
0033     SigmaZ = cms.double(4.0),
0034     TimeOffset = cms.double(0.0)
0035 )
0036 # Gaussian smearing
0037 # Flat optics for Run3 - Low SigmaZ
0038 # SigmaZ = 4.2 cm
0039 # SigmaX = 11.8 um
0040 # SigmaY = 5.5 um
0041 # BS positions extracted from 2018B 3.8T data, run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
0042 # X0         =  0.09676  [cm]
0043 # Y0         = -0.06245  [cm]
0044 # Z0         = -0.292    [cm]
0045 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0046 # X = 0.0859918 cm
0047 # Y = -0.104172 cm
0048 # Z = -0.327748 cm
0049 Run3FlatOpticsGaussVtxSigmaZ4p2cmSmearingParameters = cms.PSet(
0050     MeanX = cms.double(0.0107682),
0051     MeanY = cms.double(0.041722),
0052     MeanZ = cms.double(0.035748),
0053     SigmaY = cms.double(0.00055),
0054     SigmaX = cms.double(0.00118),
0055     SigmaZ = cms.double(4.2),
0056     TimeOffset = cms.double(0.0)
0057 )
0058 # Gaussian smearing
0059 # Flat optics for Run3 - High SigmaZ
0060 # SigmaZ = 5.3 cm
0061 # SigmaX = 15 um
0062 # SigmaY = 13 um
0063 # BS positions extracted from 2018B 3.8T data, run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
0064 # X0         =  0.09676  [cm]
0065 # Y0         = -0.06245  [cm]
0066 # Z0         = -0.292    [cm]
0067 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0068 # X = 0.0859918 cm
0069 # Y = -0.104172 cm
0070 # Z = -0.327748 cm
0071 Run3FlatOpticsGaussVtxSigmaZ5p3cmSmearingParameters = cms.PSet(
0072     MeanX = cms.double(0.0107682),
0073     MeanY = cms.double(0.041722),
0074     MeanZ = cms.double(0.035748),
0075     SigmaY = cms.double(0.0013),
0076     SigmaX = cms.double(0.0015),
0077     SigmaZ = cms.double(5.3),
0078     TimeOffset = cms.double(0.0)
0079 )
0080 
0081 # Flat Smearing
0082 # Important note: flat independent distributions in Z and T are not correct for physics production
0083 # In reality, if two flat beams interact the real distribution will not be flat with independent Z and T
0084 # but Z and T will be correlated, as example in GaussEvtVtxGenerator.
0085 # Can restore correlation via MinT += (MinZ - MaxZ)/2 and MaxT += (MaxZ - MinZ)/2
0086 # in [ns] units (recall c_light = 29.98cm/ns)
0087 FlatVtxSmearingParameters = cms.PSet(
0088     MaxZ = cms.double(5.3),
0089     MaxX = cms.double(0.0015),
0090     MaxY = cms.double(0.0015),
0091     MinX = cms.double(-0.0015),
0092     MinY = cms.double(-0.0015),
0093     MinZ = cms.double(-5.3),
0094     MaxT = cms.double(0.177),
0095     MinT = cms.double(-0.177)
0096 )
0097 #############################################
0098 # Beta functions smearing (pp 7+7 TeV)
0099 #
0100 # Values taken from LHC optics simulation V6.5:
0101 # see http://proj-lhc-optics-web.web.cern.ch/proj-lhc-optics-web/V6.500/IR5.html
0102 # alpha = angle of the crossing plane 0 degrees means XZ plane
0103 # phi = half-crossing beam angle
0104 #
0105 # Emittance is the no normalized emittance in cm = normalized emittance/gamma (beta=1)
0106 #
0107 # length variables are in cm
0108 #
0109 
0110 # 900 GeV collisions, transverse beam size = 293 microns
0111 Early900GeVCollisionVtxSmearingParameters = cms.PSet(
0112     Phi = cms.double(0.0),
0113     BetaStar = cms.double(1100.0),
0114     Emittance = cms.double(1.564e-06),
0115     Alpha = cms.double(0.0),
0116     SigmaZ = cms.double(7.4),
0117     TimeOffset = cms.double(0.0),
0118     Y0 = cms.double(0.0),
0119     X0 = cms.double(0.0322),
0120     Z0 = cms.double(0.0)
0121 )
0122 #  2.2 TeV collisions, transverse beam size 188 microns
0123 Early2p2TeVCollisionVtxSmearingParameters = cms.PSet(
0124     Phi = cms.double(0.0),
0125     BetaStar = cms.double(1100.0),
0126     Emittance = cms.double(6.4e-07),
0127     Alpha = cms.double(0.0),
0128     SigmaZ = cms.double(5.5),
0129     TimeOffset = cms.double(0.0),
0130     Y0 = cms.double(0.0),
0131     X0 = cms.double(0.0322),
0132     Z0 = cms.double(0.0)
0133 )
0134 #  7 TeV collisions, transverse beam size with betastar=  11m is 105 microns,
0135 Early7TeVCollisionVtxSmearingParameters = cms.PSet(
0136     Phi = cms.double(0.0),
0137     BetaStar = cms.double(1100.0),
0138     Emittance = cms.double(2.0e-07),
0139     Alpha = cms.double(0.0),
0140     SigmaZ = cms.double(4.2),
0141     TimeOffset = cms.double(0.0),
0142     Y0 = cms.double(0.0),
0143     X0 = cms.double(0.0322),
0144     Z0 = cms.double(0.0)
0145 )
0146 #  7 TeV collisions, transverse beam size with betastar=  2m is  45 microns,
0147 Nominal7TeVCollisionVtxSmearingParameters = cms.PSet(
0148     Phi = cms.double(0.0),
0149     BetaStar = cms.double(200.0),
0150     Emittance = cms.double(2.0e-07),
0151     Alpha = cms.double(0.0),
0152     SigmaZ = cms.double(4.2),
0153     TimeOffset = cms.double(0.0),
0154     Y0 = cms.double(0.0),
0155     X0 = cms.double(0.0322),
0156     Z0 = cms.double(0.0)
0157 )
0158 # 900 GeV realistic 2010 collisions, transverse beam size is 200 microns
0159 Realistic900GeVCollisionVtxSmearingParameters = cms.PSet(
0160     Phi = cms.double(0.0),
0161     BetaStar = cms.double(1000.0),
0162     Emittance = cms.double(8.34e-07),
0163     Alpha = cms.double(0.0),
0164     SigmaZ = cms.double(6.17),
0165     TimeOffset = cms.double(0.0),
0166     X0 = cms.double(0.2452),
0167     Y0 = cms.double(0.3993),
0168     Z0 = cms.double(0.8222)
0169 )
0170 # 7 TeV realistic collisions, beamspot width ~28 microns - appropriate for 2nd half of Commissioning10
0171 Realistic7TeVCollisionComm10VtxSmearingParameters = cms.PSet(
0172     Phi = cms.double(0.0),
0173     BetaStar = cms.double(200.0),
0174     Emittance = cms.double(0.804e-07),
0175     Alpha = cms.double(0.0),
0176     SigmaZ = cms.double(3.50),
0177     TimeOffset = cms.double(0.0),
0178     X0 = cms.double(0.2440),
0179     Y0 = cms.double(0.3929),
0180     Z0 = cms.double(0.4145)
0181 )
0182 # 7 TeV realistic collisions, beamspot width ~43 microns - appropriate for 2010A
0183 Realistic7TeVCollisionVtxSmearingParameters = cms.PSet(
0184     Phi = cms.double(0.0),
0185     BetaStar = cms.double(350.0),
0186     Emittance = cms.double(1.072e-07),
0187     Alpha = cms.double(0.0),
0188     SigmaZ = cms.double(6.26),
0189     TimeOffset = cms.double(0.0),
0190     X0 = cms.double(0.2440),
0191     Y0 = cms.double(0.3929),
0192     Z0 = cms.double(0.4145)
0193 )
0194 # 7 TeV realistic collisions, beamspot width ~38 microns - appropriate for 2010B
0195 Realistic7TeVCollision2010BVtxSmearingParameters = cms.PSet(
0196     Phi = cms.double(0.0),
0197     BetaStar = cms.double(350.0),
0198     Emittance = cms.double(0.804e-07),
0199     Alpha = cms.double(0.0),
0200     SigmaZ = cms.double(5.40),
0201     TimeOffset = cms.double(0.0),
0202     X0 = cms.double(0.2440),
0203     Y0 = cms.double(0.3929),
0204     Z0 = cms.double(0.4145)
0205 )
0206 # 7 TeV realistic collisions, updated for 2011
0207 # normalized emittance 2.5 microns, transverse beam size is 32 microns
0208 Realistic7TeV2011CollisionVtxSmearingParameters = cms.PSet(
0209     Phi = cms.double(0.0),
0210     BetaStar = cms.double(150.0),
0211     Emittance = cms.double(0.67e-07),
0212     Alpha = cms.double(0.0),
0213     SigmaZ = cms.double(5.22),
0214     TimeOffset = cms.double(0.0),
0215     X0 = cms.double(0.2440),
0216     Y0 = cms.double(0.3929),
0217     Z0 = cms.double(0.4145)
0218 )
0219 # HI realistic collisions, updated for 2011
0220 # estimated beamspot width 31-35 microns
0221 RealisticHI2011CollisionVtxSmearingParameters = cms.PSet(
0222     Phi = cms.double(0.0),
0223     BetaStar = cms.double(100.0),
0224     Emittance = cms.double(2.04e-07),
0225     Alpha = cms.double(0.0),
0226     SigmaZ = cms.double(7.06),
0227     TimeOffset = cms.double(0.0),
0228     X0 = cms.double(0.2245),
0229     Y0 = cms.double(0.4182),
0230     Z0 = cms.double(0.0847)
0231 )
0232 # 2.76 TeV estimated collisions, 11m beta*
0233 # normalized emittance 2.5 microns, transverse beam size is 140 microns
0234 Realistic2p76TeV2011CollisionVtxSmearingParameters = cms.PSet(
0235     Phi = cms.double(0.0),
0236     BetaStar = cms.double(1100.0),
0237     Emittance = cms.double(1.70e-07),
0238     Alpha = cms.double(0.0),
0239     SigmaZ = cms.double(5.22),
0240     TimeOffset = cms.double(0.0),
0241     X0 = cms.double(0.2440),
0242     Y0 = cms.double(0.3929),
0243     Z0 = cms.double(0.4145)
0244 )
0245 # 2.76 TeV estimated collisions for 2013, 11m beta*
0246 # sigmaZ set to 8 cm
0247 Realistic2p76TeV2013CollisionVtxSmearingParameters = cms.PSet(
0248     Phi = cms.double(0.0),
0249     BetaStar = cms.double(1100.0),
0250     Emittance = cms.double(1.70e-07),
0251     Alpha = cms.double(0.0),
0252     SigmaZ = cms.double(8.0),
0253     TimeOffset = cms.double(0.0),
0254     X0 = cms.double(0.2440),
0255     Y0 = cms.double(0.3929),
0256     Z0 = cms.double(0.4145)
0257 )
0258 # HI realistic pPb collisions, updated for 2013
0259 #
0260 RealisticHIpPb2013CollisionVtxSmearingParameters = cms.PSet(
0261     Phi = cms.double(0.0),
0262     BetaStar = cms.double(80.0),
0263     Emittance = cms.double(6.25e-07),
0264     Alpha = cms.double(0.0),
0265     SigmaZ = cms.double(8.0),
0266     TimeOffset = cms.double(0.0),
0267     X0 = cms.double(0.2440),
0268     Y0 = cms.double(0.3929),
0269     Z0 = cms.double(0.4145)
0270 )
0271 # 7 TeV centered collisions with parameters for 2011
0272 # normalized emittance 2.5 microns, transverse beam size is 32 microns
0273 Centered7TeV2011CollisionVtxSmearingParameters = cms.PSet(
0274     Phi = cms.double(0.0),
0275     BetaStar = cms.double(150.0),
0276     Emittance = cms.double(0.67e-07),
0277     Alpha = cms.double(0.0),
0278     SigmaZ = cms.double(5.22),
0279     TimeOffset = cms.double(0.0),
0280     X0 = cms.double(0.0),
0281     Y0 = cms.double(0.0),
0282     Z0 = cms.double(0.0)
0283 )
0284 # 8 TeV realistic collisions, transverse beam width size is 20 microns
0285 Realistic8TeVCollisionVtxSmearingParameters = cms.PSet(
0286     Phi = cms.double(0.0),
0287     BetaStar = cms.double(70.0),
0288     Emittance = cms.double(0.586e-07),
0289     Alpha = cms.double(0.0),
0290     SigmaZ = cms.double(6.16),
0291     TimeOffset = cms.double(0.0),
0292     X0 = cms.double(0.2440),
0293     Y0 = cms.double(0.3929),
0294     Z0 = cms.double(0.4145)
0295 )
0296 # 8 TeV realistic collisions, transverse beam width size is 20 microns, updated for observed SigmaZ
0297 Realistic8TeV2012CollisionVtxSmearingParameters = cms.PSet(
0298     Phi = cms.double(0.0),
0299     BetaStar = cms.double(70.0),
0300     Emittance = cms.double(0.586e-07),
0301     Alpha = cms.double(0.0),
0302     SigmaZ = cms.double(4.8),
0303     TimeOffset = cms.double(0.0),
0304     X0 = cms.double(0.2440),
0305     Y0 = cms.double(0.3929),
0306     Z0 = cms.double(0.4145)
0307 )
0308 # 10 TeV collisions, transverse beam size = 46 microns
0309 Early10TeVCollisionVtxSmearingParameters = cms.PSet(
0310     Phi = cms.double(0.0),
0311     BetaStar = cms.double(300.0),
0312     Emittance = cms.double(1.406e-07),
0313     Alpha = cms.double(0.0),
0314     SigmaZ = cms.double(3.8),
0315     TimeOffset = cms.double(0.0),
0316     Y0 = cms.double(0.0),
0317     X0 = cms.double(0.0322),
0318     Z0 = cms.double(0.0)
0319 )
0320 # Test offset
0321 Early10TeVX322Y100VtxSmearingParameters = cms.PSet(
0322     Phi = cms.double(0.0),
0323     BetaStar = cms.double(300.0),
0324     Emittance = cms.double(1.406e-07),
0325     Alpha = cms.double(0.0),
0326     SigmaZ = cms.double(3.8),
0327     TimeOffset = cms.double(0.0),
0328     Y0 = cms.double(0.0100),
0329     X0 = cms.double(0.0322),
0330     Z0 = cms.double(0.0)
0331 )
0332 # Test offset
0333 Early10TeVX322Y250VtxSmearingParameters = cms.PSet(
0334     Phi = cms.double(0.0),
0335     BetaStar = cms.double(300.0),
0336     Emittance = cms.double(1.406e-07),
0337     Alpha = cms.double(0.0),
0338     SigmaZ = cms.double(3.8),
0339     TimeOffset = cms.double(0.0),
0340     Y0 = cms.double(0.0250),
0341     X0 = cms.double(0.0322),
0342     Z0 = cms.double(0.0)
0343 )
0344 # Test offset
0345 Early10TeVX322Y500VtxSmearingParameters = cms.PSet(
0346     Phi = cms.double(0.0),
0347     BetaStar = cms.double(300.0),
0348     Emittance = cms.double(1.406e-07),
0349     Alpha = cms.double(0.0),
0350     SigmaZ = cms.double(3.8),
0351     TimeOffset = cms.double(0.0),
0352     Y0 = cms.double(0.0500),
0353     X0 = cms.double(0.0322),
0354     Z0 = cms.double(0.0)
0355 )
0356 # Test offset
0357 Early10TeVX322Y1000VtxSmearingParameters = cms.PSet(
0358     Phi = cms.double(0.0),
0359     BetaStar = cms.double(300.0),
0360     Emittance = cms.double(1.406e-07),
0361     Alpha = cms.double(0.0),
0362     SigmaZ = cms.double(3.8),
0363     TimeOffset = cms.double(0.0),
0364     Y0 = cms.double(0.1),
0365     X0 = cms.double(0.0322),
0366     Z0 = cms.double(0.0)
0367 )
0368 # Test offset
0369 Early10TeVX322Y5000VtxSmearingParameters = cms.PSet(
0370     Phi = cms.double(0.0),
0371     BetaStar = cms.double(300.0),
0372     Emittance = cms.double(1.406e-07),
0373     Alpha = cms.double(0.0),
0374     SigmaZ = cms.double(3.8),
0375     TimeOffset = cms.double(0.0),
0376     Y0 = cms.double(0.5),
0377     X0 = cms.double(0.0322),
0378     Z0 = cms.double(0.0)
0379 )
0380 # Test offset
0381 Early10TeVX322Y10000VtxSmearingParameters = cms.PSet(
0382     Phi = cms.double(0.0),
0383     BetaStar = cms.double(300.0),
0384     Emittance = cms.double(1.406e-07),
0385     Alpha = cms.double(0.0),
0386     SigmaZ = cms.double(3.8),
0387     TimeOffset = cms.double(0.0),
0388     Y0 = cms.double(1.0),
0389     X0 = cms.double(0.0322),
0390     Z0 = cms.double(0.0)
0391 )
0392 
0393 EarlyCollisionVtxSmearingParameters = cms.PSet(
0394     Phi = cms.double(0.0),
0395     BetaStar = cms.double(200.0),
0396     Emittance = cms.double(1.006e-07),
0397     Alpha = cms.double(0.0),
0398     SigmaZ = cms.double(5.3),
0399     TimeOffset = cms.double(0.0),
0400     Y0 = cms.double(0.0),
0401     X0 = cms.double(0.0322),
0402     Z0 = cms.double(0.0)
0403 )
0404 NominalCollisionVtxSmearingParameters = cms.PSet(
0405     Phi = cms.double(0.000142),
0406     BetaStar = cms.double(55.0),
0407     Emittance = cms.double(1.006e-07),
0408     Alpha = cms.double(0.0),
0409     SigmaZ = cms.double(5.3),
0410     TimeOffset = cms.double(0.0),
0411     Y0 = cms.double(0.0),
0412     X0 = cms.double(0.05),
0413     Z0 = cms.double(0.0)
0414 )
0415 NominalCollision1VtxSmearingParameters = cms.PSet(
0416     Phi = cms.double(0.0),
0417     BetaStar = cms.double(55.0),
0418     Emittance = cms.double(1.006e-07),
0419     Alpha = cms.double(0.0),
0420     SigmaZ = cms.double(5.3),
0421     TimeOffset = cms.double(0.0),
0422     Y0 = cms.double(0.025),
0423     X0 = cms.double(0.05),
0424     Z0 = cms.double(0.0)
0425 )
0426 NominalCollision2VtxSmearingParameters = cms.PSet(
0427     Phi = cms.double(0.000142),
0428     BetaStar = cms.double(55.0),
0429     Emittance = cms.double(1.006e-07),
0430     Alpha = cms.double(0.0),
0431     SigmaZ = cms.double(5.3),
0432     TimeOffset = cms.double(0.0),
0433     Y0 = cms.double(0.025),
0434     X0 = cms.double(0.05),
0435     Z0 = cms.double(0.0)
0436 )
0437 NominalCollision3VtxSmearingParameters = cms.PSet(
0438     Phi = cms.double(0.0),
0439     BetaStar = cms.double(55.0),
0440     Emittance = cms.double(1.006e-07),
0441     Alpha = cms.double(0.0),
0442     SigmaZ = cms.double(5.3),
0443     TimeOffset = cms.double(0.0),
0444     Y0 = cms.double(0.025),
0445     X0 = cms.double(0.1),
0446     Z0 = cms.double(0.0)
0447 )
0448 NominalCollision4VtxSmearingParameters = cms.PSet(
0449     Phi = cms.double(0.0),
0450     BetaStar = cms.double(55.0),
0451     Emittance = cms.double(1.006e-07),
0452     Alpha = cms.double(0.0),
0453     SigmaZ = cms.double(5.3),
0454     TimeOffset = cms.double(0.0),
0455     Y0 = cms.double(0.025),
0456     X0 = cms.double(0.2),
0457     Z0 = cms.double(0.0)
0458 )
0459 NominalCollision2015VtxSmearingParameters = cms.PSet(
0460     Phi = cms.double(0.0),
0461     BetaStar = cms.double(65.0),
0462     Emittance = cms.double(5.411e-08),
0463     Alpha = cms.double(0.0),
0464     SigmaZ = cms.double(5.3),
0465     TimeOffset = cms.double(0.0),
0466     X0 = cms.double(0.0322),
0467     Y0 = cms.double(0.0),
0468     Z0 = cms.double(0.0)
0469 )
0470 ZeroTeslaRun247324CollisionVtxSmearingParameters = cms.PSet(
0471     Phi = cms.double(0.0),
0472     BetaStar = cms.double(80.0),
0473     Emittance = cms.double(1.070e-5),
0474     Alpha = cms.double(0.0),
0475     SigmaZ = cms.double(4.125),
0476     TimeOffset = cms.double(0.0),
0477     X0 = cms.double(0.08621),
0478     Y0 = cms.double(0.1657),
0479     Z0 = cms.double(-1.688)
0480 )
0481 
0482 # From 2015A 0T data
0483 # Centroid absolute positions extracted from fills:
0484 # X = 0.059395  cm
0485 # Y = 0.099686  cm
0486 # Z = -1.722240 cm
0487 #
0488 # BPIX absolute position extracted from first collision alignment:
0489 # X = -0.0259503 cm
0490 # Y = -0.07004   cm
0491 # Z = -0.498917  cm
0492 Realistic50ns13TeVCollisionZeroTeslaVtxSmearingParameters = cms.PSet(
0493     Phi = cms.double(0.0),
0494     BetaStar = cms.double(65.0),
0495     Emittance = cms.double(5.411e-08),
0496     Alpha = cms.double(0.0),
0497     SigmaZ = cms.double(5.3),
0498     TimeOffset = cms.double(0.0),
0499     X0 = cms.double(0.08533),
0500     Y0 = cms.double(0.16973),
0501     Z0 = cms.double(-1.2230)
0502 )
0503 
0504 # From 2015B 3.8T data
0505 # Centroid absolute positions extracted from fill 4008:
0506 # X =  0.07798 cm
0507 # Y =  0.09714 cm
0508 # Z = -1.610   cm
0509 #
0510 # BPIX absolute position extracted from PCL-like alignment run after magnet ramp-up:
0511 # X = -0.026837  cm
0512 # Y = -0.0715252 cm
0513 # Z = -0.511453  cm
0514 Realistic50ns13TeVCollisionVtxSmearingParameters = cms.PSet(
0515     Phi = cms.double(0.0),
0516     BetaStar = cms.double(65.0),
0517     Emittance = cms.double(5.411e-08),
0518     Alpha = cms.double(0.0),
0519     SigmaZ = cms.double(5.3),
0520     TimeOffset = cms.double(0.0),
0521     X0 = cms.double(0.10482),
0522     Y0 = cms.double(0.16867),
0523     Z0 = cms.double(-1.0985)
0524 )
0525 
0526 # From 2015B 3.8T data, beta*=90m (700 bunches fills)
0527 # Centroid absolute positions extracted from 700 bunches fills 4499-4511:
0528 # X = 0.068357 cm
0529 # Y = 0.109159 cm
0530 # Z = 0.131811 cm
0531 #
0532 # BPIX absolute position extracted from Prompt Reco alignment of run 259352
0533 # X = -0.041651 cm
0534 # Y = -0.199279 cm
0535 # Z = -0.565093 cm
0536 #
0537 # Emittance has been calculated to match a BeamWidht of O(10um) with: https://lpc.web.cern.ch/lumi2.html
0538 #
0539 Realistic100ns13TeVCollisionBetaStar90mVtxSmearingParameters = cms.PSet(
0540     Phi = cms.double(0.0),
0541     BetaStar = cms.double(9121.0),
0542     Emittance = cms.double(0.12e-7),
0543     Alpha = cms.double(0.0),
0544     SigmaZ = cms.double(4.9),
0545     TimeOffset = cms.double(0.0),
0546     X0 = cms.double(0.11000),
0547     Y0 = cms.double(0.30844),
0548     Z0 = cms.double(0.69690)
0549 )
0550 
0551 # From 2015B 3.8T data, beta*=90m (42/240 bunches fills)
0552 # Centroid absolute positions extracted from 42/240 bunches fills 4495-4496:
0553 # X = 0.064925 cm
0554 # Y = 0.112761 cm
0555 # Z = 0.170413 cm
0556 #
0557 # BPIX absolute position extracted from Prompt Reco alignment of run 259202
0558 # X = -0.041651 cm
0559 # Y = -0.199279 cm
0560 # Z = -0.565093 cm
0561 #
0562 # Emittance has been calculated to match a BeamWidht of O(10um) with: https://lpc.web.cern.ch/lumi2.html
0563 #
0564 Realistic100ns13TeVCollisionBetaStar90mLowBunchesVtxSmearingParameters = cms.PSet(
0565     Phi = cms.double(0.0),
0566     BetaStar = cms.double(9121.0),
0567     Emittance = cms.double(0.12e-7),
0568     Alpha = cms.double(0.0),
0569     SigmaZ = cms.double(5.24),
0570     TimeOffset = cms.double(0.0),
0571     X0 = cms.double(0.10658),
0572     Y0 = cms.double(0.31204),
0573     Z0 = cms.double(0.735506)
0574 )
0575 
0576 # From 2016B 3.8T data
0577 # BS parameters extracted from fills 4895 - 4935:
0578 # X0         = 0.064870 +/- 4.9575E-07 [cm]
0579 # Y0         = 0.093639 +/- 4.9233E-07 [cm]
0580 # Z0         = 0.420085 +/- 4.1102E-04 [cm]
0581 # sigmaZ0    = 3.645533 +/- 2.9064E-04 [cm]
0582 #
0583 # From LHC calculator, emittance is 4.906e-8 cm
0584 # https://lpc.web.cern.ch/lpc/lumi2.html
0585 #
0586 # BPIX absolute position:
0587 # X = -0.0267572 cm
0588 # Y = -0.0759102 cm
0589 # Z = -0.511428  cm
0590 Realistic25ns13TeV2016CollisionVtxSmearingParameters = cms.PSet(
0591     Phi = cms.double(0.0),
0592     BetaStar = cms.double(40.0),
0593     Emittance = cms.double(4.906e-8),
0594     Alpha = cms.double(0.0),
0595     SigmaZ = cms.double(3.65),
0596     TimeOffset = cms.double(0.0),
0597     X0 = cms.double(0.09163),
0598     Y0 = cms.double(0.16955),
0599     Z0 = cms.double(0.9315 )
0600 )
0601 
0602 # From 2017A 3.8T data
0603 # BS parameters extracted from run 295463 (from offline DQM, i.e. PCL):
0604 # X0         =  0.08497  [cm]
0605 # Y0         = -0.03976  [cm]
0606 # Z0         =  1.6      [cm] ==> 0.5 adjusted after cogging tuning by LHC, see  https://hypernews.cern.ch/HyperNews/CMS/get/beamspot/159/1.html
0607 # sigmaZ0    =  3.5      [cm]
0608 #
0609 # From LHC calculator, emittance is 3.319e-8 cm
0610 # https://lpc.web.cern.ch/lpc/lumi2.html
0611 #
0612 # BPIX absolute position (https://hypernews.cern.ch/HyperNews/CMS/get/tif-alignment/657/1/1.html):
0613 # X = 0.109725 cm
0614 # Y = -0.108993 cm
0615 # Z = -0.32054  cm
0616 Realistic25ns13TeVEarly2017CollisionVtxSmearingParameters = cms.PSet(
0617     Phi = cms.double(0.0),
0618     BetaStar = cms.double(40.0),
0619     Emittance = cms.double(3.319e-8),
0620     Alpha = cms.double(0.0),
0621     SigmaZ = cms.double(3.5),
0622     TimeOffset = cms.double(0.0),
0623     X0 = cms.double(-0.024755),
0624     Y0 = cms.double(0.069233 ),
0625     Z0 = cms.double(0.82054  )
0626 )
0627 
0628 # Beam spot extracted from data for 2017 pp run @ 5 TeV
0629 Realistic5TeVppCollision2017VtxSmearingParameters = cms.PSet(
0630     Phi = cms.double(0.0),
0631     BetaStar = cms.double(311),
0632     Emittance = cms.double(3.8e-8),
0633     Alpha = cms.double(0.0),
0634     SigmaZ = cms.double(5.82),
0635     TimeOffset = cms.double(0.0),
0636     X0 = cms.double(-0.0228),
0637     Y0 = cms.double(0.0795),
0638     Z0 = cms.double(0.619)
0639 )
0640 
0641 # Fixed Emittance (X2) in Beam spot extracted from data for 2017 pp run @ 5 TeV
0642 Fixed_EmitRealistic5TeVppCollision2017VtxSmearingParameters = cms.PSet(
0643     Phi = cms.double(0.0),
0644     BetaStar = cms.double(311),
0645     Emittance = cms.double(7.6e-8),
0646     Alpha = cms.double(0.0),
0647     SigmaZ = cms.double(5.82),
0648     TimeOffset = cms.double(0.0),
0649     X0 = cms.double(-0.0228),
0650     Y0 = cms.double(0.0795),
0651     Z0 = cms.double(0.619)
0652 )
0653 
0654 
0655 # From 2018B 3.8T data
0656 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
0657 # X0         =  0.09676  [cm]
0658 # Y0         = -0.06245  [cm]
0659 # Z0         = -0.292    [cm]
0660 # sigmaZ0    =  3.5      [cm] => mean sigmaZ0 in this run is 3.2676
0661 # BeamWidthX 0.0008050
0662 # BeamWidthY 0.0006238
0663 #
0664 # From LHC calculator, emittance is 1.634e-8 cm
0665 # https://lpc.web.cern.ch/lpc/lumi2.html
0666 #
0667 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0668 # X = 0.0859918 cm
0669 # Y = -0.104172 cm
0670 # Z = -0.327748 cm
0671 Realistic25ns13TeVEarly2018CollisionVtxSmearingParameters = cms.PSet(
0672     Phi = cms.double(0.0),
0673     BetaStar = cms.double(30.0),
0674     Emittance = cms.double(1.634e-8),
0675     Alpha = cms.double(0.0),
0676     SigmaZ = cms.double(3.5),
0677     TimeOffset = cms.double(0.0),
0678     X0 = cms.double(0.0107682),
0679     Y0 = cms.double(0.041722 ),
0680     Z0 = cms.double(0.035748 )
0681 )
0682 
0683 # Run3 possible beam parameters
0684 # Round optics - Low SigmaZ = 3.4 cm
0685 # From 2018B 3.8T data
0686 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
0687 # X0         =  0.09676  [cm]
0688 # Y0         = -0.06245  [cm]
0689 # Z0         = -0.292    [cm]
0690 # sigmaZ0    =  3.2676   [cm]
0691 # BeamWidthX 0.0008050
0692 # BeamWidthY 0.0006238
0693 #
0694 # set SigmaZ0 = 3.4 [cm]
0695 # set BeamWidthX = BeamWidthY = 11.5 [um]
0696 # set beta* = 28 cm
0697 # energy = 13 TeV
0698 # From LHC calculator, emittance is 4.762e-8 cm
0699 # https://lpc.web.cern.ch/lpc/lumi2.html
0700 #
0701 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0702 # X = 0.0859918 cm
0703 # Y = -0.104172 cm
0704 # Z = -0.327748 cm
0705 Run3RoundOptics25ns13TeVLowSigmaZVtxSmearingParameters = cms.PSet(
0706     Phi = cms.double(0.0),
0707     BetaStar = cms.double(28.0),
0708     Emittance = cms.double(4.762e-8),
0709     Alpha = cms.double(0.0),
0710     SigmaZ = cms.double(3.4),
0711     TimeOffset = cms.double(0.0),
0712     X0 = cms.double(0.0107682),
0713     Y0 = cms.double(0.041722 ),
0714     Z0 = cms.double(0.035748 )
0715 )
0716 
0717 # Run3 possible beam parameters
0718 # Round optics - High SigmaZ = 5.7 cm
0719 # From 2018B 3.8T data
0720 # BS parameters extracted from run 316199, fill 6675 (from StreamExpressAlignment, HP BS):
0721 # X0         =  0.09676  [cm]
0722 # Y0         = -0.06245  [cm]
0723 # Z0         = -0.292    [cm]
0724 # sigmaZ0    =  3.2676   [cm]
0725 # BeamWidthX 0.0008050
0726 # BeamWidthY 0.0006238
0727 #
0728 # set SigmaZ0 = 5.7 [cm]
0729 # set BeamWidthX = BeamWidthY = 11.5 [um]
0730 # set beta* = 28 cm
0731 # energy = 13 TeV
0732 # From LHC calculator, emittance is 4.762e-8 cm
0733 # https://lpc.web.cern.ch/lpc/lumi2.html
0734 #
0735 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0736 # X = 0.0859918 cm
0737 # Y = -0.104172 cm
0738 # Z = -0.327748 cm
0739 Run3RoundOptics25ns13TeVHighSigmaZVtxSmearingParameters = cms.PSet(
0740     Phi = cms.double(0.0),
0741     BetaStar = cms.double(28.0),
0742     Emittance = cms.double(4.762e-8),
0743     Alpha = cms.double(0.0),
0744     SigmaZ = cms.double(5.7),
0745     TimeOffset = cms.double(0.0),
0746     X0 = cms.double(0.0107682),
0747     Y0 = cms.double(0.041722 ),
0748     Z0 = cms.double(0.035748 )
0749 )
0750 
0751 # From Run3 PilotBeams 2021 3.8T data
0752 # BS parameters extracted from run 346512, fill 7531 (from ExpressPhysics FEVT, Legacy BS):
0753 # X0         =  0.174282 [cm]
0754 # Y0         = -0.187132 [cm]
0755 # Z0         =  0.167616 [cm]
0756 # sigmaZ0    =  6.80728  [cm]
0757 # BeamWidthX 0.0142174
0758 # BeamWidthY 0.0150789
0759 #
0760 # set SigmaZ0 = 6.8 [cm]
0761 # set BeamWidthX = BeamWidthY = 150.0 [um]
0762 # set beta* = 1100 cm
0763 # energy = 900 GeV
0764 # From LHC calculator, emittance is 4.762e-8 cm
0765 # https://lpc.web.cern.ch/lpc/lumi2.html
0766 #
0767 # BPIX absolute position (https://twiki.cern.ch/twiki/bin/view/CMS/TkAlignmentPixelPosition?rev=40#2021):
0768 # X =  0.06076 cm
0769 # Y = -0.14702 cm
0770 # Z = -0.25616 cm
0771 Realistic25ns900GeV2021CollisionVtxSmearingParameters = cms.PSet(
0772     Phi = cms.double(0.0),
0773     BetaStar = cms.double(1100.0),
0774     Emittance = cms.double(4.169e-7),
0775     Alpha = cms.double(0.0),
0776     SigmaZ = cms.double(6.8),
0777     TimeOffset = cms.double(0.0),
0778     X0 = cms.double(0.113522),
0779     Y0 = cms.double(-0.040112),
0780     Z0 = cms.double(0.423776)
0781 )
0782 
0783 # From first Run 3 data at 13.6 TeV and 3.8T
0784 # BS parameters extracted from run 355100, fill 7920:
0785 # X0         =  0.172394 [cm]
0786 # Y0         = -0.180946 [cm]
0787 # Z0         =  0.94181  [cm]
0788 # sigmaZ0    =  3.81941  [cm]
0789 # BeamWidthX = 0.0008772 [cm]
0790 # BeamWidthY = 0.0010078 [cm]
0791 #
0792 # set SigmaZ0 = 3.8 [cm]
0793 # set BeamWidthX = BeamWidthY = 10.0 [um]
0794 # set beta* = 30 cm
0795 # energy = 13.6 TeV
0796 # From LHC calculator, emittance is 6.621e-8 cm
0797 # https://lpc.web.cern.ch/lumiCalc.html
0798 #
0799 # BPIX absolute position (https://twiki.cern.ch/twiki/bin/view/CMS/TkAlignmentPixelPosition?rev=45#Collisions_at_s_13_6_TeV):
0800 # X =  0.0717651 cm
0801 # Y = -0.165951  cm
0802 # Z = -0.356345  cm
0803 Realistic25ns13p6TeVEarly2022CollisionVtxSmearingParameters = cms.PSet(
0804     Phi = cms.double(0.0),
0805     BetaStar = cms.double(30.0),
0806     Emittance = cms.double(6.621e-8),
0807     Alpha = cms.double(0.0),
0808     SigmaZ = cms.double(3.8),
0809     TimeOffset = cms.double(0.0),
0810     X0 = cms.double(0.100629),
0811     Y0 = cms.double(-0.014995),
0812     Z0 = cms.double(1.298155)
0813 )
0814 
0815 # BS parameters extracted from run 360459, Fill 8274:
0816 # X0         =  0.1742 [cm]
0817 # Y0         = -0.1831 [cm]
0818 # Z0         = -0.2531 [cm]
0819 # sigmaZ0    =  3.4019  [cm]
0820 # BeamWidthX = 0.0007519 [cm]
0821 # BeamWidthY = 0.0008636 [cm]
0822 #
0823 # set SigmaZ0 = 3.4 [cm]
0824 # set BeamWidthX = BeamWidthY = 8.0 [um]
0825 # set beta* = 30 cm
0826 # energy = 13.6 TeV
0827 # From LHC calculator, emittance is 4.276-8 cm
0828 # https://lpc.web.cern.ch/lumiCalc.html
0829 #
0830 # BPIX absolute position:
0831 # X =  0.0714025 cm
0832 # Y = -0.166338  cm
0833 # Z = -0.354856  cm
0834 Realistic25ns13p6TeVEOY2022CollisionVtxSmearingParameters = cms.PSet(
0835     Phi = cms.double(0.0),
0836     BetaStar = cms.double(30.0),
0837     Emittance = cms.double(4.276e-8),
0838     Alpha = cms.double(0.0),
0839     SigmaZ = cms.double(3.4),
0840     TimeOffset = cms.double(0.0),
0841     X0 = cms.double(0.1027975),
0842     Y0 = cms.double(-0.016762),
0843     Z0 = cms.double(0.101756)
0844 )
0845 
0846 # BS parameters extracted averaging Fills 8728-8750 (2023C):
0847 # X0         =  0.117154 [cm]
0848 # Y0         = -0.186556 [cm]
0849 # Z0         = -0.431777 [cm]
0850 # sigmaZ0    =  3.599 cm [cm]
0851 # BeamWidthX = 0.0007333 [cm]
0852 # BeamWidthY = 0.0008046 [cm]
0853 #
0854 # set SigmaZ0 = 3.6 [cm]
0855 # set BeamWidthX = BeamWidthY = 7.7 [um]
0856 # set beta* = 30 cm
0857 # energy = 13.6 TeV
0858 # From LHC calculator, emittance is 3.931e-8 cm
0859 # https://lpc.web.cern.ch/lumiCalc.html
0860 #
0861 # BPIX absolute position (from Runs 367094-367589):
0862 # X =  0.0713008 cm
0863 # Y = -0.169590  cm
0864 # Z = -0.356785  cm
0865 Realistic25ns13p6TeVEarly2023CollisionVtxSmearingParameters = cms.PSet(
0866     Phi = cms.double(0.0),
0867     BetaStar = cms.double(30.0),
0868     Emittance = cms.double(3.931e-8),
0869     Alpha = cms.double(0.0),
0870     SigmaZ = cms.double(3.6),
0871     TimeOffset = cms.double(0.0),
0872     X0 = cms.double(0.0458532),
0873     Y0 = cms.double(-0.016966),
0874     Z0 = cms.double(-0.074992)
0875 )
0876 
0877 # Test HF offset
0878 ShiftedCollision2015VtxSmearingParameters = cms.PSet(
0879     Phi = cms.double(0.0),
0880     BetaStar = cms.double(65.0),
0881     Emittance = cms.double(5.411e-08),
0882     Alpha = cms.double(0.0),
0883     SigmaZ = cms.double(5.3),
0884     TimeOffset = cms.double(0.0),
0885     X0 = cms.double(1.0),
0886     Y0 = cms.double(0.0),
0887     Z0 = cms.double(0.0)
0888 )
0889 Shifted5mmCollision2015VtxSmearingParameters = cms.PSet(
0890     Phi = cms.double(0.0),
0891     BetaStar = cms.double(65.0),
0892     Emittance = cms.double(5.411e-08),
0893     Alpha = cms.double(0.0),
0894     SigmaZ = cms.double(5.3),
0895     TimeOffset = cms.double(0.0),
0896     X0 = cms.double(0.5),
0897     Y0 = cms.double(0.0),
0898     Z0 = cms.double(0.0)
0899 )
0900 Shifted15mmCollision2015VtxSmearingParameters = cms.PSet(
0901     Phi = cms.double(0.0),
0902     BetaStar = cms.double(65.0),
0903     Emittance = cms.double(5.411e-08),
0904     Alpha = cms.double(0.0),
0905     SigmaZ = cms.double(5.3),
0906     TimeOffset = cms.double(0.0),
0907     X0 = cms.double(1.5),
0908     Y0 = cms.double(0.0),
0909     Z0 = cms.double(0.0)
0910 )
0911 
0912 # Estimate for 2015 PbPb collisions, based on feedback from accelerator
0913 # Beamspot centroid shifted to match pp expectation for 2015
0914 NominalHICollision2015VtxSmearingParameters = cms.PSet(
0915     Phi = cms.double(0.0),
0916     BetaStar = cms.double(60.0),
0917     Emittance = cms.double(1.70e-07),
0918     Alpha = cms.double(0.0),
0919     SigmaZ = cms.double(7.06),
0920     TimeOffset = cms.double(0.0),
0921     X0 = cms.double(0.0322),
0922     Y0 = cms.double(0.),
0923     Z0 = cms.double(0.)
0924 )
0925 
0926 # updated numbers based on beamspot fits to 2015 PbPb data
0927 # Later found to be incorrect, see following entry
0928 # RealisticHICollision2015VtxSmearingParameters = cms.PSet(
0929 #    Phi = cms.double(0.0),
0930 #    BetaStar = cms.double(60.0),
0931 #    Emittance = cms.double(1.70e-07),
0932 #    Alpha = cms.double(0.0),
0933 #    SigmaZ = cms.double(5.2278),
0934 #    TimeOffset = cms.double(0.0),
0935 #    X0 = cms.double(0.1025),
0936 #    Y0 = cms.double(0.1654),
0937 #    Z0 = cms.double(3.2528)
0938 #)
0939 # updated numbers for 2015 PbPb data with Z centroid from fixed beamspot fits
0940 # See discussion here https://hypernews.cern.ch/HyperNews/CMS/get/hi-general/3968.html
0941 # See plot of difference here: https://www.dropbox.com/s/tsnkgvvpkdqjtyq/vzDataMCOverlay_c_20170420.pdf?dl=0
0942 #
0943 RealisticHICollisionFixZ2015VtxSmearingParameters = cms.PSet(
0944     Phi = cms.double(0.0),
0945     BetaStar = cms.double(60.0),
0946     Emittance = cms.double(1.70e-07),
0947     Alpha = cms.double(0.0),
0948     SigmaZ = cms.double(5.2278),
0949     TimeOffset = cms.double(0.0),
0950     X0 = cms.double(0.1025),
0951     Y0 = cms.double(0.1654),
0952     Z0 = cms.double(0.771)
0953 )
0954 # Numbers based on beamspot fits to 2017 XeXe data
0955 # Documentation here: https://twiki.cern.ch/twiki/pub/CMS/XeXeRereco/IanLRU_AlCaTkAlBS_20171130_approvedByLucaSara.pdf
0956 RealisticXeXeCollision2017VtxSmearingParameters = cms.PSet(
0957     Phi = cms.double(0.0),
0958     BetaStar = cms.double(30.0),
0959     Emittance = cms.double(4.33e-08),
0960     Alpha = cms.double(0.0),
0961     SigmaZ = cms.double(4.64),
0962     TimeOffset = cms.double(0.0),
0963     X0 = cms.double(-0.026),
0964     Y0 = cms.double(0.081),
0965     Z0 = cms.double(0.645)
0966 )
0967 
0968 # From fit to 5 TeV PbPb data
0969 # From 2018 PbPb  data
0970 # BS parameters extracted from run 327211, Fill 7471 (from StreamExpressAlignment, HP, BS):
0971 # X0         =  0.09443  [cm]
0972 # Y0         = -0.06377  [cm]
0973 # Z0         =  0.58067  [cm]
0974 # sigmaZ0    =  4.969    [cm]
0975 # BeamWidthX 0.0014392
0976 # BeamWidthY 0.0011545
0977 #
0978 # Emittance is 3.36e-8 cm. Calculated by  ((BeamWidthX + BeamWidth)/2)^2/BetaStar
0979 #
0980 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0981 # X =  0.084335 cm
0982 # Y = -0.110381 cm
0983 # Z = -0.321904 cm
0984 RealisticPbPbCollision2018VtxSmearingParameters = cms.PSet(
0985     Phi = cms.double(0.0),
0986     BetaStar = cms.double(50),
0987     Emittance = cms.double(3.36e-08),
0988     Alpha = cms.double(0.0),
0989     SigmaZ = cms.double(4.97),
0990     TimeOffset = cms.double(0.0),
0991     X0 = cms.double(0.010),
0992     Y0 = cms.double(0.047),
0993     Z0 = cms.double(0.903)
0994 )
0995 
0996 # 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
0997 Nominal5TeVpp2015VtxSmearingParameters = cms.PSet(
0998     Phi = cms.double(0.0),
0999     BetaStar = cms.double(400.0),
1000     Emittance = cms.double(1.0e-07),
1001     Alpha = cms.double(0.0),
1002     SigmaZ = cms.double(5.5),
1003     TimeOffset = cms.double(0.0),
1004     X0 = cms.double(0.1044),
1005     Y0 = cms.double(0.1676),
1006     Z0 = cms.double(0.6707)
1007 )
1008 
1009 # From fit to 5 TeV pPb data
1010 Realistic5TeVPACollision2016VtxSmearingParameters = cms.PSet(
1011     Phi = cms.double(0.0),
1012     BetaStar = cms.double(1100.0),
1013     Emittance = cms.double(6.75e-08),
1014     Alpha = cms.double(0.0),
1015     SigmaZ = cms.double(6.4891),
1016     TimeOffset = cms.double(0.0),
1017     X0 = cms.double(0.0889),
1018     Y0 = cms.double(0.1820),
1019     Z0 = cms.double(1.6066)
1020 )
1021 
1022 # From fit to 8 TeV pPb data
1023 Realistic8TeVPACollision2016VtxSmearingParameters = cms.PSet(
1024     Phi = cms.double(0.0),
1025     BetaStar = cms.double(60.0),
1026     Emittance = cms.double(6.75e-08),
1027     Alpha = cms.double(0.0),
1028     SigmaZ = cms.double(4.6914),
1029     TimeOffset = cms.double(0.0),
1030     X0 = cms.double(0.0836),
1031     Y0 = cms.double(0.1837),
1032     Z0 = cms.double(1.3577)
1033 )
1034 
1035 # Guess for 2022 PbPb beam conditions, which takes the 2018 PbPb beam width parameters from RealisticPbPbCollision2018VtxSmearingParameters with the current pp MC beam centroid from Realistic25ns13p6TeVEarly2022Collision
1036 Nominal2022PbPbCollisionVtxSmearingParameters = cms.PSet(
1037     Phi = cms.double(0.0),
1038     BetaStar = cms.double(50),
1039     Emittance = cms.double(3.36e-08),
1040     Alpha = cms.double(0.0),
1041     SigmaZ = cms.double(4.97),
1042     TimeOffset = cms.double(0.0),
1043     X0 = cms.double(0.100629),
1044     Y0 = cms.double(-0.014995),
1045     Z0 = cms.double(1.298155)
1046 )
1047 
1048 # From 2022 PbPb test data 362294
1049 Realistic2022PbPbCollisionVtxSmearingParameters = cms.PSet(
1050     Phi = cms.double(0.0),
1051     BetaStar = cms.double(50),
1052     Emittance = cms.double(3.36e-08),
1053     Alpha = cms.double(0.0),
1054     SigmaZ = cms.double(5.01265),
1055     TimeOffset = cms.double(0.0),
1056     X0 = cms.double(0.1017599),
1057     Y0 = cms.double(-0.015602),
1058     Z0 = cms.double(0.131175)
1059 )
1060 
1061 # From 2023 PbPb data run 374719 
1062 Realistic2023PbPbCollisionVtxSmearingParameters = cms.PSet(
1063     Phi = cms.double(0.0),
1064     BetaStar = cms.double(50),
1065     Emittance = cms.double(5.854e-08),
1066     Alpha = cms.double(0.0),
1067     SigmaZ = cms.double(4.60055),
1068     TimeOffset = cms.double(0.0),
1069     X0 = cms.double(0.0394727),
1070     Y0 = cms.double(-0.0157799),
1071     Z0 = cms.double(0.3064731)
1072 )
1073 
1074 # From 2024 ppRef data run 387590-387721
1075 Realistic2024ppRefCollisionVtxSmearingParameters = cms.PSet(
1076     Phi = cms.double(0.0),
1077     BetaStar = cms.double(310),
1078     Emittance = cms.double(8.636e-08),
1079     Alpha = cms.double(0.0),
1080     SigmaZ = cms.double(5.82376),
1081     TimeOffset = cms.double(0.0),
1082     X0 = cms.double(0.0244619),
1083     Y0 = cms.double(-0.0169181),
1084     Z0 = cms.double(0.3513597)
1085 )
1086 
1087 # From 2024 PbPb data run 387998-388425
1088 Realistic2024PbPbCollisionVtxSmearingParameters = cms.PSet(
1089     Phi = cms.double(0.0),
1090     BetaStar = cms.double(50),
1091     Emittance = cms.double(6.684e-08),
1092     Alpha = cms.double(0.0),
1093     SigmaZ = cms.double(4.9068349),
1094     TimeOffset = cms.double(0.0),
1095     X0 = cms.double(0.0302546),
1096     Y0 = cms.double(-0.0170382),
1097     Z0 = cms.double(0.2290316)
1098 )
1099 
1100 # Parameters for HL-LHC operation at 13TeV
1101 HLLHCVtxSmearingParameters = cms.PSet(
1102     MeanXIncm = cms.double(0.),
1103     MeanYIncm = cms.double(0.),
1104     MeanZIncm = cms.double(0.),
1105     TimeOffsetInns = cms.double(0.0),
1106     EprotonInGeV = cms.double(6500.0),
1107     CrossingAngleInurad = cms.double(510.0),
1108     CrabFrequencyInMHz = cms.double(400.0),
1109     RF800 = cms.bool(False),
1110     BetaCrossingPlaneInm = cms.double(0.20),
1111     BetaSeparationPlaneInm = cms.double(0.20),
1112     HorizontalEmittance = cms.double(2.5e-06),
1113     VerticalEmittance = cms.double(2.05e-06),
1114     BunchLengthInm = cms.double(0.090),
1115     CrabbingAngleCrossingInurad = cms.double(380.0),
1116     CrabbingAngleSeparationInurad = cms.double(0.0)
1117 )