Changing the proton propagator to the Optical Function method

SimPPS/Configuration/test/gluglu_step1_GEN_SIM_2021.py

@@ -7,7 +7,6 @@
 process = cms.Process('SIM',eras.Run3)
 
 # import of standard configurations
-process.load("CondCore.CondDB.CondDB_cfi")
 process.load('Configuration.StandardSequences.Services_cff')
 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
 process.load('FWCore.MessageService.MessageLogger_cfi')
@@ -19,8 +18,7 @@
 process.load('Configuration.StandardSequences.SimIdeal_cff')
 process.load('Configuration.StandardSequences.EndOfProcess_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
-
-process.load('Configuration.Geometry.GeometryExtended2021_cff')
+process.load('Configuration.StandardSequences.GeometrySimDB_cff')
 
 process.RandomNumberGeneratorService.generator.initialSeed = cms.untracked.uint32(random.randint(0,900000000))
 
@@ -33,38 +31,30 @@
         )
 
 process.source = cms.Source("EmptySource")
-"""
-process.source = cms.Source("EmptySource",
-                firstRun = cms.untracked.uint32(324612),               #2018D
-                firstTime = cms.untracked.uint64(6612348794983940096)  
-)
-"""
+
 from Configuration.AlCa.GlobalTag import GlobalTag
 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2021_realistic', '')
-#process.GlobalTag = GlobalTag(process.GlobalTag, "106X_dataRun2_v26")
 
 # generator
-
-
 process.generator = cms.EDFilter("ExhumeGeneratorFilter",
     ExhumeParameters = cms.PSet(
         AlphaEw = cms.double(0.0072974),
         B = cms.double(4.0),
-        BottomMass = cms.double(4.6),
-        CharmMass = cms.double(1.42),
-        HiggsMass = cms.double(120.0),
+        MuonMass = cms.double(0.1057),
+        BottomMass = cms.double(4.65),
+        CharmMass = cms.double(1.28),
+        StrangeMass = cms.double(0.095),
+        TauMass = cms.double(1.78),
+        TopMass = cms.double(172.8),
+        WMass = cms.double(80.38),
+        ZMass = cms.double(91.187),
+        HiggsMass = cms.double(125.1),
         HiggsVev = cms.double(246.0),
         LambdaQCD = cms.double(80.0),
         MinQt2 = cms.double(0.64),
-        MuonMass = cms.double(0.1057),
         PDF = cms.double(11000),
         Rg = cms.double(1.2),
-        StrangeMass = cms.double(0.19),
-        Survive = cms.double(0.03),
-        TauMass = cms.double(1.77),
-        TopMass = cms.double(175.0),
-        WMass = cms.double(80.33),
-        ZMass = cms.double(91.187)
+        Survive = cms.double(0.03)
     ),
     ExhumeProcess = cms.PSet(
         MassRangeHigh = cms.double(2000.0),
@@ -75,7 +65,7 @@
     PythiaParameters = cms.PSet(
         parameterSets = cms.vstring()
     ),
-    comEnergy = cms.double(13000.0),
+    comEnergy = cms.double(14000.0),
     maxEventsToPrint = cms.untracked.int32(2),
     pythiaHepMCVerbosity = cms.untracked.bool(False),
     pythiaPylistVerbosity = cms.untracked.int32(1)
@@ -88,7 +78,6 @@
     version = cms.untracked.string('$Revision: 1.19 $')
 )
 
-
 process.ProductionFilterSequence = cms.Sequence(process.generator)
 
 ############
@@ -100,7 +89,6 @@
 process.generation_step = cms.Path(process.pgen)
 process.simulation_step = cms.Path(process.psim)
 
-
 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
 process.outpath = cms.EndPath(process.o1)
 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.outpath)

SimPPS/Configuration/test/gluglu_step2_DIGI_DIGI2RAW_2021.py

@@ -8,9 +8,8 @@
 process = cms.Process('DIGI2RAW',eras.Run3)
 
 # import of standard configurations
-process.load("CondCore.CondDB.CondDB_cfi")
 process.load('Configuration.StandardSequences.Services_cff')
-process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
+#process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
 process.load('FWCore.MessageService.MessageLogger_cfi')
 process.load('Configuration.EventContent.EventContent_cff')
 process.load('Configuration.StandardSequences.MagneticField_cff')
@@ -20,11 +19,7 @@
 process.load('Configuration.StandardSequences.DigiToRaw_cff')
 process.load('Configuration.StandardSequences.EndOfProcess_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
-
-# Since the PPS simulation geometry is not yet in the database, the line below is needed
 process.load('Configuration.StandardSequences.GeometryDB_cff')
-#process.load("Geometry.VeryForwardGeometry.geometryPPS_CMSxz_fromDD_2021_cfi")
-process.load("Geometry.VeryForwardGeometry.geometryRPFromDB_cfi")
 
 # Input source
 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1))
@@ -70,16 +65,15 @@
         filterName = cms.untracked.string('')
     ),
     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
-    fileName = cms.untracked.string('GluGlu_step2_DIGI_DIGI2RAW2021.root'),
-    outputCommands = process.FEVTDEBUGEventContent.outputCommands + ['keep *_CTPPS*_*_*',"keep *_*RP*_*_*",'keep *_LHCTransport_*_*'],
+    fileName = cms.untracked.string('GluGlu_step2_DIGI_DIGI2RAW_2021.root'),
+    outputCommands = process.FEVTDEBUGEventContent.outputCommands + ['keep *_CTPPS*_*_*',"keep *_*RP*_*_*",'keep *_generatorSmeared_*_*'],
     splitLevel = cms.untracked.int32(0)
 )
 
 # Additional output definition
 
 # Other statements
 from Configuration.AlCa.GlobalTag import GlobalTag
-#process.GlobalTag = GlobalTag(process.GlobalTag, '111X_mcRun3_2021_realistic_Candidate_2020_06_03_17_55_10', '')
 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2021_realistic', '')
 
 

SimPPS/Configuration/test/gluglu_step3_RAW2DIGI_L1Reco_RECO_2021.py

@@ -12,26 +12,19 @@
 process.load("Configuration.EventContent.EventContent_cff")
 process.load('Configuration.StandardSequences.MagneticField_cff')
 process.load('Configuration.StandardSequences.EndOfProcess_cff')
+process.load('Configuration.StandardSequences.RawToDigi_cff')
+process.load('Configuration.StandardSequences.Reconstruction_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
 process.load('Configuration.StandardSequences.GeometryDB_cff')
-process.load("Geometry.VeryForwardGeometry.geometryRPFromDB_cfi")
-#process.load("Geometry.VeryForwardGeometry.geometryPPS_CMSxz_fromDD_2021_cfi")           # CMS frame
 
 
 process.maxEvents = cms.untracked.PSet(
     input = cms.untracked.int32(-1)
 )
 
-process.source = cms.Source("EmptyIOVSource",
-    timetype = cms.string('runnumber'),
-    firstValue = cms.uint64(1),
-    lastValue = cms.uint64(1),
-    interval = cms.uint64(1)
-)
-
 # Input source
 process.source = cms.Source("PoolSource",
-    fileNames = cms.untracked.vstring('file:GluGlu_step2_DIGI_DIGI2RAW2021.root'),
+    fileNames = cms.untracked.vstring('file:GluGlu_step2_DIGI_DIGI2RAW_2021.root'),
     secondaryFileNames = cms.untracked.vstring()
 )
 
@@ -45,31 +38,18 @@
 # Output definition
 
 process.output = cms.OutputModule("PoolOutputModule",
-    fileName = cms.untracked.string('GluGlu_step3_RAW2DIGI_RECO2021.root'),
-    outputCommands = cms.untracked.vstring("drop *","keep SimVertexs_g4SimHits_*_*","keep PSimHits*_*_*_*","keep CTPPS*_*_*_*","keep *_*RP*_*_*",'keep *_LHCTransport_*_*')
+    fileName = cms.untracked.string('GluGlu_step3_RAW2DIGI_RECO_2021.root'),
+    outputCommands = cms.untracked.vstring("drop *","keep SimVertexs_g4SimHits_*_*","keep PSimHits*_*_*_*","keep CTPPS*_*_*_*","keep *_*RP*_*_*",'keep *_generatorSmeared_*_*')
 )
 
-
 # Additional output definition
 # Other statements
 from Configuration.AlCa.GlobalTag import GlobalTag
 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2021_realistic', '')
 
-# modify CTPPS 2018 raw-to-digi modules ONLY FOR PARTICLE GUN, TO AVOID RUN THIS FOR THE WHOLE CMS
-process.load('Configuration.StandardSequences.RawToDigi_cff')
-
-# do not make testID for simulation - keeping the frame
-from EventFilter.CTPPSRawToDigi.totemRPRawToDigi_cfi import totemRPRawToDigi
-totemRPRawToDigi.RawToDigi.testID = cms.uint32(1)
-
-from RecoPPS.Local.totemRPLocalReconstruction_cff import totemRPLocalReconstruction
-process.load('RecoPPS.Local.totemRPLocalReconstruction_cff')
-from RecoPPS.Local.ctppsPixelLocalReconstruction_cff import ctppsPixelLocalReconstruction
-process.load('RecoPPS.Local.ctppsPixelLocalReconstruction_cff')
-
 # Path and EndPath definitions
-process.raw2digi_step = cms.Path(process.ctppsRawToDigi)
-process.reco_step = cms.Path(process.totemRPLocalReconstruction*process.ctppsPixelLocalReconstruction)
+process.raw2digi_step = cms.Path(process.RawToDigi)
+process.reco_step = cms.Path(process.reconstruction)
 process.endjob_step = cms.EndPath(process.endOfProcess)
 process.output_step = cms.EndPath(process.output)
 
@@ -80,4 +60,3 @@
 for path in process.paths:
     #  getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
     getattr(process,path)._seq = getattr(process,path)._seq
-

SimPPS/Configuration/test/pg_step1_GEN_SIM_2021.py

@@ -7,7 +7,6 @@
 process = cms.Process('SIM',eras.Run3)
 
 # import of standard configurations
-process.load("CondCore.CondDB.CondDB_cfi")
 process.load('Configuration.StandardSequences.Services_cff')
 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
 process.load('FWCore.MessageService.MessageLogger_cfi')
@@ -19,8 +18,7 @@
 process.load('Configuration.StandardSequences.SimIdeal_cff')
 process.load('Configuration.StandardSequences.EndOfProcess_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
-
-process.load('Configuration.Geometry.GeometryExtended2021_cff')
+process.load('Configuration.StandardSequences.GeometrySimDB_cff')
 
 process.RandomNumberGeneratorService.generator.initialSeed = cms.untracked.uint32(random.randint(0,900000000))
 
@@ -33,25 +31,30 @@
         )
 
 process.source = cms.Source("EmptySource")
-"""
-process.source = cms.Source("EmptySource",
-                firstRun = cms.untracked.uint32(324612),               #2018D
-                firstTime = cms.untracked.uint64(6612348794983940096)  
-)
-"""
 from Configuration.AlCa.GlobalTag import GlobalTag
 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2021_realistic', '')
-#process.GlobalTag = GlobalTag(process.GlobalTag, "106X_dataRun2_v26")
-
+#process.GlobalTag = GlobalTag(process.GlobalTag, '113X_mcRun3_2021_realistic_Candidate_2021_04_06_19_59_53', '')
 # generator
+"""
+process.GlobalTag.toGet = cms.VPSet(
+                 cms.PSet(
+                     record = cms.string('LHCInfoRcd'),
+                     tag = cms.string("LHCInfo_2021_mc_v1")
+                 ),
+                 cms.PSet(
+                     record = cms.string('CTPPSOpticsRcd'),
+                     tag = cms.string("PPSOpticalFunctions_2021_mc_v1")
+                 )
+     )
+"""
 
 phi_min = -math.pi
 phi_max = math.pi
 t_min   = 0.
 t_max   = 2.
-xi_min  = 0.03
-xi_max  = 0.15
-ecms = 13000.
+xi_min  = 0.02
+xi_max  = 0.20
+ecms = 14000. # ATTENTION: if using HECTOR propagator, currently the energy is hardcoded as 6500 but the optics file is prepared for 7 TeV
 
 process.generator = cms.EDProducer("RandomtXiGunProducer",
         PGunParameters = cms.PSet(
@@ -90,4 +93,3 @@
 # filter all path with the production filter sequence
 for path in process.paths:
     getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
-

SimPPS/Configuration/test/pg_step2_DIGI_DIGI2RAW_2021.py

@@ -8,9 +8,8 @@
 process = cms.Process('DIGI2RAW',eras.Run3)
 
 # import of standard configurations
-process.load("CondCore.CondDB.CondDB_cfi")
 process.load('Configuration.StandardSequences.Services_cff')
-process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
+#process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
 process.load('FWCore.MessageService.MessageLogger_cfi')
 process.load('Configuration.EventContent.EventContent_cff')
 process.load('Configuration.StandardSequences.MagneticField_cff')
@@ -20,11 +19,10 @@
 process.load('Configuration.StandardSequences.DigiToRaw_cff')
 process.load('Configuration.StandardSequences.EndOfProcess_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
+process.load('Configuration.StandardSequences.GeometryDB_cff')
 
 # Since the PPS simulation geometry is not yet in the database, the line below is needed
-process.load('Configuration.StandardSequences.GeometryDB_cff')
 #process.load("Geometry.VeryForwardGeometry.geometryPPS_CMSxz_fromDD_2021_cfi")
-process.load("Geometry.VeryForwardGeometry.geometryRPFromDB_cfi")
 
 # Input source
 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1))
@@ -71,15 +69,15 @@
     ),
     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
     fileName = cms.untracked.string('step2_DIGI_DIGI2RAW2021.root'),
-    outputCommands = process.FEVTDEBUGEventContent.outputCommands + ['keep *_CTPPS*_*_*',"keep *_*RP*_*_*",'keep *_LHCTransport_*_*'],
+    outputCommands = process.FEVTDEBUGEventContent.outputCommands + ['keep *_CTPPS*_*_*',"keep *_*RP*_*_*",'keep *_generatorSmeared_*_*'],
     splitLevel = cms.untracked.int32(0)
 )
 
 # Additional output definition
 
 # Other statements
 from Configuration.AlCa.GlobalTag import GlobalTag
-#process.GlobalTag = GlobalTag(process.GlobalTag, '111X_mcRun3_2021_realistic_Candidate_2020_06_03_17_55_10', '')
+#process.GlobalTag = GlobalTag(process.GlobalTag, '113X_mcRun3_2021_realistic_Candidate_2021_04_06_19_59_53', '')
 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2021_realistic', '')
 
 

SimPPS/Configuration/test/pg_step3_MCDB_RAW2DIGI_RECO_2021.py

@@ -13,22 +13,17 @@
 process.load('Configuration.StandardSequences.MagneticField_cff')
 process.load('Configuration.StandardSequences.EndOfProcess_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
+process.load('Configuration.StandardSequences.RawToDigi_cff')
+process.load('Configuration.StandardSequences.Reconstruction_cff')
 process.load('Configuration.StandardSequences.GeometryDB_cff')
-process.load("Geometry.VeryForwardGeometry.geometryRPFromDB_cfi")
-#process.load("Geometry.VeryForwardGeometry.geometryPPS_CMSxz_fromDD_2021_cfi")           # CMS frame
-
+#process.load("Geometry.VeryForwardGeometry.geometryRPFromDB_cfi")
+#process.load("Geometry.VeryForwardGeometry.geometryPPS_CMSxz_fromDD_2018_cfi") # CMS frame, 2021 = 2018
+#process.load('CalibPPS.ESProducers.ppsPixelTopologyESSourceRun2_cfi')          # temporary solution, the 2021 geometry is the same as run2, force the usage of run2 topology 
 
 process.maxEvents = cms.untracked.PSet(
     input = cms.untracked.int32(-1)
 )
 
-process.source = cms.Source("EmptyIOVSource",
-    timetype = cms.string('runnumber'),
-    firstValue = cms.uint64(1),
-    lastValue = cms.uint64(1),
-    interval = cms.uint64(1)
-)
-
 # Input source
 process.source = cms.Source("PoolSource",
     fileNames = cms.untracked.vstring('file:step2_DIGI_DIGI2RAW2021.root'),
@@ -46,30 +41,21 @@
 
 process.output = cms.OutputModule("PoolOutputModule",
     fileName = cms.untracked.string('step3_RAW2DIGI_RECO2021.root'),
-    outputCommands = cms.untracked.vstring("drop *","keep SimVertexs_g4SimHits_*_*","keep PSimHits*_*_*_*","keep CTPPS*_*_*_*","keep *_*RP*_*_*",'keep *_LHCTransport_*_*')
+    outputCommands = cms.untracked.vstring("drop *","keep SimVertexs_g4SimHits_*_*","keep PSimHits*_*_*_*","keep CTPPS*_*_*_*","keep *_*RP*_*_*",'keep *_generatorSmeared_*_*')
 )
 
 
 # Additional output definition
 # Other statements
 from Configuration.AlCa.GlobalTag import GlobalTag
 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2021_realistic', '')
-
-# modify CTPPS 2018 raw-to-digi modules ONLY FOR PARTICLE GUN, TO AVOID RUN THIS FOR THE WHOLE CMS
-process.load('Configuration.StandardSequences.RawToDigi_cff')
-
-# do not make testID for simulation - keeping the frame
-from EventFilter.CTPPSRawToDigi.totemRPRawToDigi_cfi import totemRPRawToDigi
-totemRPRawToDigi.RawToDigi.testID = cms.uint32(1)
-
-from RecoPPS.Local.totemRPLocalReconstruction_cff import totemRPLocalReconstruction
-process.load('RecoPPS.Local.totemRPLocalReconstruction_cff')
-from RecoPPS.Local.ctppsPixelLocalReconstruction_cff import ctppsPixelLocalReconstruction
-process.load('RecoPPS.Local.ctppsPixelLocalReconstruction_cff')
+#process.GlobalTag = GlobalTag(process.GlobalTag, '113X_mcRun3_2021_realistic_Candidate_2021_04_06_19_59_53', '')
 
 # Path and EndPath definitions
-process.raw2digi_step = cms.Path(process.ctppsRawToDigi)
-process.reco_step = cms.Path(process.totemRPLocalReconstruction*process.ctppsPixelLocalReconstruction)
+#process.raw2digi_step = cms.Path(process.RawToDigi)
+#process.reco_step = cms.Path(process.reconstruction)
+process.raw2digi_step = cms.Path(process.ctppsRawToDigi) # to avoid full CMS digitization
+process.reco_step = cms.Path(process.totemRPLocalReconstruction*process.ctppsPixelLocalReconstruction) # to avoid full CMS reconstruction
 process.endjob_step = cms.EndPath(process.endOfProcess)
 process.output_step = cms.EndPath(process.output)
 
@@ -80,4 +66,3 @@
 for path in process.paths:
     #  getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
     getattr(process,path)._seq = getattr(process,path)._seq
-