Merge pull request #14738 from MilanoBicocca-pix/CMSSW_7_1_X_vtx_smea…

…ring_2016B Vertex smearing parameters for 2016
cms-sw · Jun 7, 2016 · 8c00002 · 8c00002
2 parents 9b0d157 + 08a4f54
commit 8c00002
Show file tree

Hide file tree

Showing 3 changed files with 34 additions and 1 deletion.
diff --git a/Configuration/StandardSequences/python/VtxSmeared.py b/Configuration/StandardSequences/python/VtxSmeared.py
@@ -38,7 +38,7 @@
     'Realistic50ns13TeVCollisionZeroTesla': 'IOMC.EventVertexGenerators.VtxSmearedRealistic50ns13TeVCollisionZeroTesla_cfi',
     'Realistic50ns13TeVCollision': 'IOMC.EventVertexGenerators.VtxSmearedRealistic50ns13TeVCollision_cfi',
     'Nominal5TeVpp2015Collision':    'IOMC.EventVertexGenerators.VtxSmearedNominal5TeVpp2015Collision_cfi',
-
+    'Realistic25ns13TeV2016Collision':    'IOMC.EventVertexGenerators.VtxSmearedRealistic25ns13TeV2016Collision_cfi',
 }
 VtxSmearedDefaultKey='NominalCollision2015'
 VtxSmearedHIDefaultKey='NominalHICollision2015'
diff --git a/IOMC/EventVertexGenerators/python/VtxSmearedParameters_cfi.py b/IOMC/EventVertexGenerators/python/VtxSmearedParameters_cfi.py
@@ -462,6 +462,32 @@
     Z0 = cms.double(-1.0985)
 )
 
+# From 2016B 3.8T data
+# BS parameters extracted from fills 4895 - 4935:
+# X0         = 0.064870 +/- 4.9575E-07 [cm]
+# Y0         = 0.093639 +/- 4.9233E-07 [cm]
+# Z0         = 0.420085 +/- 4.1102E-04 [cm]
+# sigmaZ0    = 3.645533 +/- 2.9064E-04 [cm]
+#
+# From LHC calculator, emittance is 4.906e-8 cm 
+# https://lpc.web.cern.ch/lpc/lumi2.html
+# 
+# BPIX absolute position:
+# X = -0.0267572 cm
+# Y = -0.0759102 cm
+# Z = -0.511428  cm
+Realistic25ns13TeV2016CollisionVtxSmearingParameters = cms.PSet(
+    Phi = cms.double(0.0),
+    BetaStar = cms.double(40.0),
+    Emittance = cms.double(4.906e-8),
+    Alpha = cms.double(0.0),
+    SigmaZ = cms.double(3.65),
+    TimeOffset = cms.double(0.0),
+    X0 = cms.double(0.09163),
+    Y0 = cms.double(0.16955),
+    Z0 = cms.double(0.9315 )
+)
+
 # Test HF offset
 ShiftedCollision2015VtxSmearingParameters = cms.PSet(
     Phi = cms.double(0.0),

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