Merge pull request #43288 from francescobrivio/HI_2023_VtxSmearing_13…

…_2_X [13_2_X] Vertex smearing for PbPb MC based on 2023 PbPb data
cms-sw · Nov 16, 2023 · 5772b7a · 5772b7a
2 parents 74c17f1 + e5c8666
commit 5772b7a
Show file tree

Hide file tree

Showing 3 changed files with 21 additions and 0 deletions.
diff --git a/Configuration/StandardSequences/python/VtxSmeared.py b/Configuration/StandardSequences/python/VtxSmeared.py
@@ -71,6 +71,7 @@
     'Realistic25ns13p6TeVEarly2023Collision' : 'IOMC.EventVertexGenerators.VtxSmearedRealistic25ns13p6TeVEarly2023Collision_cfi',
     'Nominal2022PbPbCollision' : 'IOMC.EventVertexGenerators.VtxSmearedNominal2022PbPbCollision_cfi',
     'Realistic2022PbPbCollision' : 'IOMC.EventVertexGenerators.VtxSmearedRealistic2022PbPbCollision_cfi',
+    'Realistic2023PbPbCollision' : 'IOMC.EventVertexGenerators.VtxSmearedRealistic2023PbPbCollision_cfi',
 }
 VtxSmearedDefaultKey='Realistic50ns13TeVCollision'
 VtxSmearedHIDefaultKey='RealisticPbPbCollision2018'
diff --git a/IOMC/EventVertexGenerators/python/VtxSmearedParameters_cfi.py b/IOMC/EventVertexGenerators/python/VtxSmearedParameters_cfi.py
@@ -1058,6 +1058,19 @@
     Z0 = cms.double(0.131175)
 )
 
+# From 2023 PbPb data run 374719 
+Realistic2023PbPbCollisionVtxSmearingParameters = cms.PSet(
+    Phi = cms.double(0.0),
+    BetaStar = cms.double(50),
+    Emittance = cms.double(5.854e-08),
+    Alpha = cms.double(0.0),
+    SigmaZ = cms.double(4.60055),
+    TimeOffset = cms.double(0.0),
+    X0 = cms.double(0.0394727),
+    Y0 = cms.double(-0.0157799),
+    Z0 = cms.double(0.3064731)
+)
+
 # Parameters for HL-LHC operation at 13TeV
 HLLHCVtxSmearingParameters = cms.PSet(
     MeanXIncm = cms.double(0.),

diff --git a/IOMC/EventVertexGenerators/python/VtxSmearedRealistic2023PbPbCollision_cfi.py b/IOMC/EventVertexGenerators/python/VtxSmearedRealistic2023PbPbCollision_cfi.py
@@ -0,0 +1,7 @@
+import FWCore.ParameterSet.Config as cms
+
+from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import Realistic2023PbPbCollisionVtxSmearingParameters,VtxSmearedCommon
+VtxSmeared = cms.EDProducer("BetafuncEvtVtxGenerator",
+    Realistic2023PbPbCollisionVtxSmearingParameters,
+    VtxSmearedCommon
+)