Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2022-08-10 23:12:12

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 # Test HF offset
0816 ShiftedCollision2015VtxSmearingParameters = cms.PSet(
0817     Phi = cms.double(0.0),
0818     BetaStar = cms.double(65.0),
0819     Emittance = cms.double(5.411e-08),
0820     Alpha = cms.double(0.0),
0821     SigmaZ = cms.double(5.3),
0822     TimeOffset = cms.double(0.0),
0823     X0 = cms.double(1.0),
0824     Y0 = cms.double(0.0),
0825     Z0 = cms.double(0.0)
0826 )
0827 Shifted5mmCollision2015VtxSmearingParameters = cms.PSet(
0828     Phi = cms.double(0.0),
0829     BetaStar = cms.double(65.0),
0830     Emittance = cms.double(5.411e-08),
0831     Alpha = cms.double(0.0),
0832     SigmaZ = cms.double(5.3),
0833     TimeOffset = cms.double(0.0),
0834     X0 = cms.double(0.5),
0835     Y0 = cms.double(0.0),
0836     Z0 = cms.double(0.0)
0837 )
0838 Shifted15mmCollision2015VtxSmearingParameters = cms.PSet(
0839     Phi = cms.double(0.0),
0840     BetaStar = cms.double(65.0),
0841     Emittance = cms.double(5.411e-08),
0842     Alpha = cms.double(0.0),
0843     SigmaZ = cms.double(5.3),
0844     TimeOffset = cms.double(0.0),
0845     X0 = cms.double(1.5),
0846     Y0 = cms.double(0.0),
0847     Z0 = cms.double(0.0)
0848 )
0849 
0850 # Estimate for 2015 PbPb collisions, based on feedback from accelerator
0851 # Beamspot centroid shifted to match pp expectation for 2015
0852 NominalHICollision2015VtxSmearingParameters = cms.PSet(
0853     Phi = cms.double(0.0),
0854     BetaStar = cms.double(60.0),
0855     Emittance = cms.double(1.70e-07),
0856     Alpha = cms.double(0.0),
0857     SigmaZ = cms.double(7.06),
0858     TimeOffset = cms.double(0.0),
0859     X0 = cms.double(0.0322),
0860     Y0 = cms.double(0.),
0861     Z0 = cms.double(0.)
0862 )
0863 
0864 # updated numbers based on beamspot fits to 2015 PbPb data
0865 # Later found to be incorrect, see following entry
0866 # RealisticHICollision2015VtxSmearingParameters = cms.PSet(
0867 #    Phi = cms.double(0.0),
0868 #    BetaStar = cms.double(60.0),
0869 #    Emittance = cms.double(1.70e-07),
0870 #    Alpha = cms.double(0.0),
0871 #    SigmaZ = cms.double(5.2278),
0872 #    TimeOffset = cms.double(0.0),
0873 #    X0 = cms.double(0.1025),
0874 #    Y0 = cms.double(0.1654),
0875 #    Z0 = cms.double(3.2528)
0876 #)
0877 # updated numbers for 2015 PbPb data with Z centroid from fixed beamspot fits
0878 # See discussion here https://hypernews.cern.ch/HyperNews/CMS/get/hi-general/3968.html
0879 # See plot of difference here: https://www.dropbox.com/s/tsnkgvvpkdqjtyq/vzDataMCOverlay_c_20170420.pdf?dl=0
0880 #
0881 RealisticHICollisionFixZ2015VtxSmearingParameters = cms.PSet(
0882     Phi = cms.double(0.0),
0883     BetaStar = cms.double(60.0),
0884     Emittance = cms.double(1.70e-07),
0885     Alpha = cms.double(0.0),
0886     SigmaZ = cms.double(5.2278),
0887     TimeOffset = cms.double(0.0),
0888     X0 = cms.double(0.1025),
0889     Y0 = cms.double(0.1654),
0890     Z0 = cms.double(0.771)
0891 )
0892 # Numbers based on beamspot fits to 2017 XeXe data
0893 # Documentation here: https://twiki.cern.ch/twiki/pub/CMS/XeXeRereco/IanLRU_AlCaTkAlBS_20171130_approvedByLucaSara.pdf
0894 RealisticXeXeCollision2017VtxSmearingParameters = cms.PSet(
0895     Phi = cms.double(0.0),
0896     BetaStar = cms.double(30.0),
0897     Emittance = cms.double(4.33e-08),
0898     Alpha = cms.double(0.0),
0899     SigmaZ = cms.double(4.64),
0900     TimeOffset = cms.double(0.0),
0901     X0 = cms.double(-0.026),
0902     Y0 = cms.double(0.081),
0903     Z0 = cms.double(0.645)
0904 )
0905 
0906 # From fit to 5 TeV PbPb data
0907 # From 2018 PbPb  data
0908 # BS parameters extracted from run 327211, Fill 7471 (from StreamExpressAlignment, HP, BS):
0909 # X0         =  0.09443  [cm]
0910 # Y0         = -0.06377  [cm]
0911 # Z0         =  0.58067  [cm]
0912 # sigmaZ0    =  4.969    [cm]
0913 # BeamWidthX 0.0014392
0914 # BeamWidthY 0.0011545
0915 #
0916 # Emittance is 3.36e-8 cm. Calculated by  ((BeamWidthX + BeamWidth)/2)^2/BetaStar
0917 #
0918 # BPIX absolute position (from https://cms-conddb.cern.ch/cmsDbBrowser/payload_inspector/Prod):
0919 # X =  0.084335 cm
0920 # Y = -0.110381 cm
0921 # Z = -0.321904 cm
0922 RealisticPbPbCollision2018VtxSmearingParameters = cms.PSet(
0923     Phi = cms.double(0.0),
0924     BetaStar = cms.double(50),
0925     Emittance = cms.double(3.36e-08),
0926     Alpha = cms.double(0.0),
0927     SigmaZ = cms.double(4.97),
0928     TimeOffset = cms.double(0.0),
0929     X0 = cms.double(0.010),
0930     Y0 = cms.double(0.047),
0931     Z0 = cms.double(0.903)
0932 )
0933 
0934 # 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
0935 Nominal5TeVpp2015VtxSmearingParameters = cms.PSet(
0936     Phi = cms.double(0.0),
0937     BetaStar = cms.double(400.0),
0938     Emittance = cms.double(1.0e-07),
0939     Alpha = cms.double(0.0),
0940     SigmaZ = cms.double(5.5),
0941     TimeOffset = cms.double(0.0),
0942     X0 = cms.double(0.1044),
0943     Y0 = cms.double(0.1676),
0944     Z0 = cms.double(0.6707)
0945 )
0946 
0947 # From fit to 5 TeV pPb data
0948 Realistic5TeVPACollision2016VtxSmearingParameters = cms.PSet(
0949     Phi = cms.double(0.0),
0950     BetaStar = cms.double(1100.0),
0951     Emittance = cms.double(6.75e-08),
0952     Alpha = cms.double(0.0),
0953     SigmaZ = cms.double(6.4891),
0954     TimeOffset = cms.double(0.0),
0955     X0 = cms.double(0.0889),
0956     Y0 = cms.double(0.1820),
0957     Z0 = cms.double(1.6066)
0958 )
0959 
0960 # From fit to 8 TeV pPb data
0961 Realistic8TeVPACollision2016VtxSmearingParameters = cms.PSet(
0962     Phi = cms.double(0.0),
0963     BetaStar = cms.double(60.0),
0964     Emittance = cms.double(6.75e-08),
0965     Alpha = cms.double(0.0),
0966     SigmaZ = cms.double(4.6914),
0967     TimeOffset = cms.double(0.0),
0968     X0 = cms.double(0.0836),
0969     Y0 = cms.double(0.1837),
0970     Z0 = cms.double(1.3577)
0971 )
0972 
0973 # 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
0974 Nominal2022PbPbCollisionVtxSmearingParameters = cms.PSet(
0975     Phi = cms.double(0.0),
0976     BetaStar = cms.double(50),
0977     Emittance = cms.double(3.36e-08),
0978     Alpha = cms.double(0.0),
0979     SigmaZ = cms.double(4.97),
0980     TimeOffset = cms.double(0.0),
0981     X0 = cms.double(0.100629),
0982     Y0 = cms.double(-0.014995),
0983     Z0 = cms.double(1.298155)
0984 )
0985 
0986 # Parameters for HL-LHC operation at 13TeV
0987 HLLHCVtxSmearingParameters = cms.PSet(
0988     MeanXIncm = cms.double(0.),
0989     MeanYIncm = cms.double(0.),
0990     MeanZIncm = cms.double(0.),
0991     TimeOffsetInns = cms.double(0.0),
0992     EprotonInGeV = cms.double(6500.0),
0993     CrossingAngleInurad = cms.double(510.0),
0994     CrabFrequencyInMHz = cms.double(400.0),
0995     RF800 = cms.bool(False),
0996     BetaCrossingPlaneInm = cms.double(0.20),
0997     BetaSeparationPlaneInm = cms.double(0.20),
0998     HorizontalEmittance = cms.double(2.5e-06),
0999     VerticalEmittance = cms.double(2.05e-06),
1000     BunchLengthInm = cms.double(0.090),
1001     CrabbingAngleCrossingInurad = cms.double(380.0),
1002     CrabbingAngleSeparationInurad = cms.double(0.0)
1003 )
1004 
1005 # Parameters for HL-LHC Crab-kissing operation 13 TeV
1006 HLLHCCrabKissingVtxSmearingParameters = cms.PSet(
1007     MeanXIncm = cms.double(0.),
1008     MeanYIncm = cms.double(0.),
1009     MeanZIncm = cms.double(0.),
1010     TimeOffsetInns = cms.double(0.0),
1011     EprotonInGeV = cms.double(6500.0),
1012     HalfCrossingAngleInurad = cms.double(200.0),
1013     CrabAngleCrossingPlaneInurad = cms.double(200.0),
1014     CrabFrequencyCrossingPlaneInMHz = cms.double(400.0),
1015     NormalizedEmittanceCrossingPlaneInum = cms.double(2.5),
1016     BetaStarCrossingPlaneInm = cms.double(0.30),
1017     CrabAngleParallelPlaneInurad = cms.double(100.0),
1018     CrabFrequencyParallelPlaneInMHz = cms.double(400.0),
1019     NormalizedEmittanceParallelPlaneInum = cms.double(2.5),
1020     BetaStarParallelPlaneInm = cms.double(0.075),
1021     ZsizeInm = cms.double(0.15),
1022     BeamProfile=cms.string("Flat")
1023 )