-
Notifications
You must be signed in to change notification settings - Fork 1
/
run_only_CTPPS_cfg_DIG_100.py
121 lines (94 loc) · 4.72 KB
/
run_only_CTPPS_cfg_DIG_100.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
process = cms.Process('CTPPS',eras.Run2_2017)
# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('SimGeneral.MixingModule.mixNoPU_cfi')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic50ns13TeVCollision_cfi')
process.load('GeneratorInterface.Core.genFilterSummary_cff')
process.load('Configuration.StandardSequences.SimIdeal_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
############### using only CTPPS geometry
from Configuration.Geometry.geometry_CTPPS_cfi import XMLIdealGeometryESSource_CTPPS
process.XMLIdealGeometryESSource = XMLIdealGeometryESSource_CTPPS.clone()
#process.load("Configuration.Geometry.geometry_CTPPS_cfi")
##########SimTransport###########
process.load('SimTransport.HectorProducer.HectorTransportCTPPS_cfi')
process.load("SimG4Core.Application.g4SimHits_cfi")
process.g4SimHits.Generator.ApplyPCuts = False
process.g4SimHits.Generator.ApplyPhiCuts = False
process.g4SimHits.Generator.ApplyEtaCuts = False
process.g4SimHits.Generator.HepMCProductLabel = 'LHCTransport'
process.g4SimHits.Generator.MinEtaCut = -13.0
process.g4SimHits.Generator.MaxEtaCut = 13.0
process.g4SimHits.Generator.Verbosity = 0
process.g4SimHits.G4TrackingManagerVerbosity = cms.untracked.int32(3)
process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(2000.0)
process.g4SimHits.StackingAction.MaxTrackTime = cms.double(2000.0)
process.load("IOMC.RandomEngine.IOMC_cff")
process.RandomNumberGeneratorService.generator.initialSeed = 456789
process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
process.RandomNumberGeneratorService.RPixDetDigitizer = cms.PSet(initialSeed =cms.untracked.uint32(137137))
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(100)
)
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2017_realistic', '')
process.load('PhysicsTools.HepMCCandAlgos.genParticles_cfi')
import math
phi_min = -math.pi
phi_max = math.pi
t_min = 0.
t_max = 2.
xi_min = 0.01
xi_max = 0.3
ecms = 13000.
process.generator = cms.EDProducer("RandomtXiGunProducer",
PGunParameters = cms.PSet(
PartID = cms.vint32(2212),
MinPhi = cms.double(phi_min),
MaxPhi = cms.double(phi_max),
ECMS = cms.double(ecms),
Mint = cms.double(t_min),
Maxt = cms.double(t_max),
MinXi = cms.double(xi_min),
MaxXi = cms.double(xi_max)
),
Verbosity = cms.untracked.int32(0),
psethack = cms.string('single protons'),
FireBackward = cms.bool(True),
FireForward = cms.bool(True)
)
process.source = cms.Source("EmptySource", firstRun=cms.untracked.uint32(294700))
#process.ProductionFilterSequence = cms.Sequence(process.generator*process.VtxSmeared)
process.ProductionFilterSequence = cms.Sequence(process.generator)
############
process.o1 = cms.OutputModule("PoolOutputModule",
# process.FEVTSIMEventContent,
outputCommands = cms.untracked.vstring('keep *'),
fileName = cms.untracked.string('simevent_CTPPS_DIG_100.root')
)
process.common_maximum_timex = cms.PSet( # need to be localy redefined
MaxTrackTime = cms.double(2000.0), # need to be localy redefined
MaxTimeNames = cms.vstring('ZDCRegion'), # need to be localy redefined
MaxTrackTimes = cms.vdouble(10000.0), # need to be localy redefined
DeadRegions = cms.vstring()
)
process.generation_step = cms.Path(process.pgen)
process.simulation_step = cms.Path(process.psim)
process.g4Simhits_step = cms.Path(process.g4SimHits)
# Transport
process.transport_step = cms.Path(process.LHCTransport)
process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
process.outpath = cms.EndPath(process.o1)
process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.transport_step,process.g4Simhits_step,process.outpath)
# filter all path with the production filter sequence
for path in process.paths:
getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq