Skip to content

Commit 1f8416b

Browse files
author
Sunanda
committed
Add a scenario with 1mm production cutoff in HGCal TB
1 parent aed1e1d commit 1f8416b

File tree

4 files changed

+182
-2
lines changed

4 files changed

+182
-2
lines changed

‎Geometry/HGCalCommonData/data/TB181/Oct181/hgcProdCuts.xml‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
<?xml version="1.0"?>
2+
<DDDefinition>
3+
4+
<SpecParSection label="hgcProdCuts.xml" eval="true">
5+
<SpecPar name="hgcProdCuts">
6+
<PartSelector path="//HGCalEE"/>
7+
<PartSelector path="//HGCalHE"/>
8+
<Parameter name="CMSCutsRegion" value="HGCalRegion" eval="false"/>
9+
<Parameter name="ProdCutsForElectrons" value="1.0*mm"/>
10+
<Parameter name="ProdCutsForProtons" value="1.0*mm"/>
11+
<Parameter name="ProdCutsForPositrons" value="1.0*mm"/>
12+
<Parameter name="ProdCutsForGamma" value="1.0*mm"/>
13+
</SpecPar>
14+
</SpecParSection>
15+
16+
</DDDefinition>

‎Geometry/HGCalCommonData/data/TB181/hgcProdCuts.xml‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,13 +1,13 @@
11
<?xml version="1.0"?>
2-
<DDDefinition xmlns="http://www.cern.ch/cms/DDL" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.cern.ch/cms/DDL ../../../../DetectorDescription/Schema/DDLSchema.xsd">
2+
<DDDefinition>
33

44
<SpecParSection label="hgcProdCuts.xml" eval="true">
55
<SpecPar name="hgcProdCuts">
66
<PartSelector path="//HGCalEE"/>
77
<PartSelector path="//HGCalHE"/>
88
<Parameter name="CMSCutsRegion" value="HGCalRegion" eval="false"/>
99
<Parameter name="ProdCutsForElectrons" value="0.03*mm"/>
10-
<Parameter name="ProdCutsForProtons" value="0.03*mm"/>
10+
<Parameter name="ProdCutsForProtons" value="0.03*mm"/>
1111
<Parameter name="ProdCutsForPositrons" value="0.03*mm"/>
1212
<Parameter name="ProdCutsForGamma" value="0.03*mm"/>
1313
</SpecPar>

‎SimG4CMS/HGCalTestBeam/python/HGCalTB181Oct1PC1mmXML_cfi.py‎

Lines changed: 21 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,21 @@
1+
import FWCore.ParameterSet.Config as cms
2+
3+
XMLIdealGeometryESSource = cms.ESSource("XMLIdealGeometryESSource",
4+
geomXMLFiles = cms.vstring('Geometry/CMSCommonData/data/materials.xml',
5+
'Geometry/CMSCommonData/data/rotations.xml',
6+
'Geometry/HGCalCommonData/data/hgcalMaterial/v1/hgcalMaterial.xml',
7+
'Geometry/HGCalCommonData/data/TB181/Oct181/cms.xml',
8+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcal.xml',
9+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcalEE.xml',
10+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcalHE.xml',
11+
'Geometry/HGCalCommonData/data/TB181/ahcal.xml',
12+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcalBeam.xml',
13+
'Geometry/HGCalCommonData/data/hgcalwafer/v7/hgcalwafer.xml',
14+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcalsense.xml',
15+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcProdCuts.xml',
16+
'Geometry/HGCalCommonData/data/TB181/Oct181/hgcalCons.xml'
17+
),
18+
rootNodeName = cms.string('cms:OCMS')
19+
)
20+
21+

‎SimG4CMS/HGCalTestBeam/test/HGCalTBCERN181Oct1PC1mm_cfg.py‎

Lines changed: 143 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,143 @@
1+
import FWCore.ParameterSet.Config as cms
2+
3+
process = cms.Process('SIM')
4+
5+
# import of standard configurations
6+
process.load("FWCore.MessageService.MessageLogger_cfi")
7+
process.load('Configuration.StandardSequences.Services_cff')
8+
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
9+
process.load('Configuration.EventContent.EventContent_cff')
10+
process.load('SimGeneral.MixingModule.mixNoPU_cfi')
11+
process.load('SimG4CMS.HGCalTestBeam.HGCalTB181Oct1PC1mmXML_cfi')
12+
process.load('Geometry.HGCalCommonData.hgcalNumberingInitialization_cfi')
13+
process.load('Geometry.HGCalCommonData.hgcalParametersInitialization_cfi')
14+
process.load('Geometry.HcalTestBeamData.hcalTB06Parameters_cff')
15+
process.load('Configuration.StandardSequences.MagneticField_0T_cff')
16+
process.load('Configuration.StandardSequences.Generator_cff')
17+
process.load('IOMC.EventVertexGenerators.VtxSmearedFlat_cfi')
18+
process.load('GeneratorInterface.Core.genFilterSummary_cff')
19+
process.load('Configuration.StandardSequences.SimIdeal_cff')
20+
process.load('Configuration.StandardSequences.EndOfProcess_cff')
21+
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
22+
process.load('SimG4CMS.HGCalTestBeam.HGCalTBAnalyzer_cfi')
23+
process.load('SimG4CMS.HGCalTestBeam.HGCalTBCheckGunPosition_cfi')
24+
25+
process.maxEvents = cms.untracked.PSet(
26+
input = cms.untracked.int32(10)
27+
)
28+
29+
#if 'MessageLogger' in process.__dict__:
30+
# process.MessageLogger.HGCalGeom=dict()
31+
# process.MessageLogger.SimG4CoreGeometry=dict()
32+
# process.MessageLogger.HGCSim=dict()
33+
# process.MessageLogger.HcalSim=dict()
34+
# process.MessageLogger.HcalTB06BeamSD=dict()
35+
36+
# Input source
37+
process.source = cms.Source("EmptySource")
38+
39+
process.options = cms.untracked.PSet(
40+
)
41+
42+
# Production Info
43+
process.configurationMetadata = cms.untracked.PSet(
44+
annotation = cms.untracked.string('SingleMuonE200_cfi nevts:10'),
45+
name = cms.untracked.string('Applications'),
46+
version = cms.untracked.string('$Revision: 1.19 $')
47+
)
48+
49+
# Output definition
50+
51+
process.FEVTDEBUGoutput = cms.OutputModule("PoolOutputModule",
52+
SelectEvents = cms.untracked.PSet(
53+
SelectEvents = cms.vstring('generation_step')
54+
),
55+
dataset = cms.untracked.PSet(
56+
dataTier = cms.untracked.string('GEN-SIM'),
57+
filterName = cms.untracked.string('')
58+
),
59+
fileName = cms.untracked.string('file:TBGenSim181Oct1.root'),
60+
outputCommands = process.FEVTDEBUGEventContent.outputCommands,
61+
splitLevel = cms.untracked.int32(0)
62+
)
63+
64+
# Additional output definition
65+
process.TFileService = cms.Service("TFileService",
66+
fileName = cms.string('TBGenSim.root')
67+
)
68+
69+
# Other statements
70+
process.genstepfilter.triggerConditions=cms.vstring("generation_step")
71+
from Configuration.AlCa.GlobalTag import GlobalTag
72+
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
73+
74+
process.generator = cms.EDProducer("FlatRandomEThetaGunProducer",
75+
AddAntiParticle = cms.bool(False),
76+
PGunParameters = cms.PSet(
77+
MinE = cms.double(99.99),
78+
MaxE = cms.double(100.01),
79+
MinTheta = cms.double(0.0),
80+
MaxTheta = cms.double(0.0),
81+
MinPhi = cms.double(-3.14159265359),
82+
MaxPhi = cms.double(3.14159265359),
83+
PartID = cms.vint32(11)
84+
),
85+
Verbosity = cms.untracked.int32(0),
86+
firstRun = cms.untracked.uint32(1),
87+
psethack = cms.string('single muon E 100')
88+
)
89+
process.VtxSmeared.MinZ = -800.0
90+
process.VtxSmeared.MaxZ = -800.0
91+
process.VtxSmeared.MinX = -7.5
92+
process.VtxSmeared.MaxX = 7.5
93+
process.VtxSmeared.MinY = -7.5
94+
process.VtxSmeared.MaxY = 7.5
95+
process.g4SimHits.HGCSD.RejectMouseBite = True
96+
process.g4SimHits.HGCSD.RotatedWafer = True
97+
process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
98+
HGCPassive = cms.PSet(
99+
LVNames = cms.vstring('HGCalEE','HGCalHE','HGCalAH', 'HGCalBeam', 'CMSE'),
100+
MotherName = cms.string('OCMS'),
101+
),
102+
type = cms.string('HGCPassive'),
103+
)
104+
)
105+
process.HGCalTBAnalyzer.doDigis = False
106+
process.HGCalTBAnalyzer.doRecHits = False
107+
process.HGCalTBAnalyzer.useFH = True
108+
process.HGCalTBAnalyzer.useBH = True
109+
process.HGCalTBAnalyzer.useBeam = True
110+
process.HGCalTBAnalyzer.zFrontEE = 1110.0
111+
process.HGCalTBAnalyzer.zFrontFH = 1176.5
112+
process.HGCalTBAnalyzer.zFrontFH = 1307.5
113+
process.HGCalTBAnalyzer.maxDepth = 39
114+
process.HGCalTBAnalyzer.deltaZ = 26.2
115+
process.HGCalTBAnalyzer.zFirst = 22.8
116+
process.HGCalTBAnalyzer.doPassive = True
117+
process.HGCalTBAnalyzer.doPassiveEE = True
118+
process.HGCalTBAnalyzer.doPassiveHE = True
119+
process.HGCalTBAnalyzer.doPassiveBH = True
120+
121+
# Path and EndPath definitions
122+
process.generation_step = cms.Path(process.pgen)
123+
process.gunfilter_step = cms.Path(process.HGCalTBCheckGunPostion)
124+
process.simulation_step = cms.Path(process.psim)
125+
process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
126+
process.analysis_step = cms.Path(process.HGCalTBAnalyzer)
127+
process.endjob_step = cms.EndPath(process.endOfProcess)
128+
process.FEVTDEBUGoutput_step = cms.EndPath(process.FEVTDEBUGoutput)
129+
130+
# Schedule definition
131+
process.schedule = cms.Schedule(process.generation_step,
132+
process.genfiltersummary_step,
133+
process.simulation_step,
134+
process.gunfilter_step,
135+
process.analysis_step,
136+
process.endjob_step,
137+
process.FEVTDEBUGoutput_step
138+
)
139+
# filter all path with the production filter sequence
140+
for path in process.paths:
141+
getattr(process,path)._seq = process.generator * getattr(process,path)._seq
142+
143+

0 commit comments

Comments
 (0)