Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Phase2-hgx359Z Update the scripts in Validation/HGCalValidation for the changes in Configuration/Geometry #45849

Merged
merged 1 commit into from
Sep 6, 2024
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
62 changes: 31 additions & 31 deletions Validation/HGCalValidation/macros/MakePlots.C
Original file line number Diff line number Diff line change
@@ -17,12 +17,12 @@
//
// where
// fanme (std::string) Input ROOT file name
// ("hfnRecHitD31tt.root" for noseRecHitPlots,
// "roots/hgcSimHitD83tt.root" hgcalStudyPlots,
// "roots/hgcGeomCheckD83.root" hgcalGeomCheck,
// "roots/hgcBHValidD83.root" hgcalBHValid,
// "roots/hgcSilValidD86.root" hgcalSiliconAnalysis,
// "missedRecHitD88.root" hgcalMissedHitPlots,
// ("hfnRecHitD115tt.root" for noseRecHitPlots,
// "roots/hgcSimHitD95tt.root" hgcalStudyPlots,
// "roots/hgcGeomCheckD95.root" hgcalGeomCheck,
// "roots/hgcBHValidD95.root" hgcalBHValid,
// "roots/hgcSilValidD100.root" hgcalSiliconAnalysis,
// "missedRecHitD105.root" hgcalMissedHitPlots,
// "ecalHitdd4hep.root" ehcalPlot)
// dirnm (std::string) Directory name
// ("hfnoseRecHitStudy" for noseRecHitPlots
@@ -33,17 +33,17 @@
// type (int) Type: 0 SimHit; 1 Digi; 2 RecHit (0)
// tag (std::string) Name of the tag for the canvas name
// ("HFNose" for noseRecHitPlots,
// "SimHitD83" for hgcalStudyPlots,
// "GeomChkD83" for hgcalGeomCheck,
// "BHValidD83" for hgcalBHValid,
// "SilValidD86" for hgcalSiliconAnalysisPlots,
// "MissedHitD88" for hgcalMissedHitPlots,
// "SimHitD95" for hgcalStudyPlots,
// "GeomChkD95" for hgcalGeomCheck,
// "BHValidD95" for hgcalBHValid,
// "SilValidD100" for hgcalSiliconAnalysisPlots,
// "MissedHitD105" for hgcalMissedHitPlots,
// "DD4Hep" for ehcalPlots)
// dtype (std::string) Data type added for canvas name
// ("ttbar D83" for hgcalStudyPlots,
// "#mu D83" for hgcalGeomCheck, hgcalBHValid',
// "ttbar D86" for hgcalSiliconAnalysisPlots,
// "#mu D88" for hgcalMissedHitPlots)
// ("ttbar D95" for hgcalStudyPlots,
// "#mu D95" for hgcalGeomCheck, hgcalBHValid',
// "ttbar D100" for hgcalSiliconAnalysisPlots,
// "#mu D105" for hgcalMissedHitPlots)
// save (bool) Flag to save the canvas (false)
// debug (bool) Debug flag (false)
//
@@ -76,7 +76,7 @@
#include <iostream>
#include <fstream>

void noseRecHitPlots(std::string fname = "hfnRecHitD31tt.root",
void noseRecHitPlots(std::string fname = "hfnRecHitD115tt.root",
std::string dirnm = "hfnoseRecHitStudy",
std::string tag = "HFNose",
bool save = false) {
@@ -176,10 +176,10 @@ void noseRecHitPlots(std::string fname = "hfnRecHitD31tt.root",
}
}

void hgcalStudyPlots(std::string fname = "roots/hgcSimHitD83tt.root",
void hgcalStudyPlots(std::string fname = "roots/hgcSimHitD95tt.root",
int type = 0,
std::string tag = "SimHitD83",
std::string dtype = "ttbar D83",
std::string tag = "SimHitD95",
std::string dtype = "ttbar D95",
bool save = false,
bool debug = false) {
int ndir[3] = {1, 3, 3};
@@ -290,9 +290,9 @@ void hgcalStudyPlots(std::string fname = "roots/hgcSimHitD83tt.root",
}
}

void hgcalGeomCheckPlots(std::string fname = "roots/hgcGeomCheckD83.root",
std::string tag = "GeomChkD83",
std::string dtype = "#mu D83",
void hgcalGeomCheckPlots(std::string fname = "roots/hgcGeomCheckD95.root",
std::string tag = "GeomChkD95",
std::string dtype = "#mu D95",
bool statbox = true,
bool save = false,
bool debug = false) {
@@ -375,9 +375,9 @@ void hgcalGeomCheckPlots(std::string fname = "roots/hgcGeomCheckD83.root",
}
}

void hgcalBHValidPlots(std::string fname = "roots/hgcBHValidD83.root",
std::string tag = "BHValidD83",
std::string dtype = "#mu D83",
void hgcalBHValidPlots(std::string fname = "roots/hgcBHValidD95.root",
std::string tag = "BHValidD95",
std::string dtype = "#mu D95",
bool save = false,
bool debug = false) {
std::string dirnm = "hgcalBHAnalysis";
@@ -481,9 +481,9 @@ void hgcalBHValidPlots(std::string fname = "roots/hgcBHValidD83.root",
}
}

void hgcalSiliconAnalysisPlots(std::string fname = "roots/hgcSilValidD86.root",
std::string tag = "SilValidD86",
std::string dtype = "ttbar D86",
void hgcalSiliconAnalysisPlots(std::string fname = "roots/hgcSilValidD100.root",
std::string tag = "SilValidD100",
std::string dtype = "ttbar D100",
bool save = false,
bool debug = false) {
std::string dirnm[2] = {"hgcalSiliconAnalysisEE", "hgcalSiliconAnalysisHEF"};
@@ -590,9 +590,9 @@ void hgcalSiliconAnalysisPlots(std::string fname = "roots/hgcSilValidD86.root",
}
}

void hgcalMissedHitPlots(std::string fname = "missedRecHitD88.root",
std::string tag = "MissedHitD88",
std::string dtype = "#mu D88",
void hgcalMissedHitPlots(std::string fname = "missedRecHitD105.root",
std::string tag = "MissedHitD105",
std::string dtype = "#mu D105",
bool save = false,
bool debug = false) {
std::string dirnm = "hgcMissingRecHit";
2 changes: 1 addition & 1 deletion Validation/HGCalValidation/scripts/muonTomo.sh
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
#!/bin/bash
geom=D92
geom=D110
cmsRun testHGCalSingleMuonPt100_cfg.py geometry=$geom
cmsRun testHGCalCellHitSum_cfg.py geometry=$geom
12 changes: 6 additions & 6 deletions Validation/HGCalValidation/scripts/testHGCalCellHitSum_cfg.py
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
###############################################################################
# Way to use this:
# cmsRun testHGCalCellHitSum_cfg.py geometry=D92 layers=1 type=mu
# cmsRun testHGCalCellHitSum_cfg.py geometry=D110 layers=1 type=mu
#
# Options for geometry D88, D92, D93
# Options for geometry D105, D110, D114
# layers '1', '1,2', any combination from 1..47
# type mu, tt
#
@@ -15,10 +15,10 @@
### SETUP OPTIONS
options = VarParsing.VarParsing('standard')
options.register('geometry',
"D92",
"D110",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"geometry of operations: D88, D92, D93")
"geometry of operations: D105, D110, D114")

options.register('layers',
"1",
@@ -48,9 +48,9 @@
inputFile = "file:step1" + options.geometry + options.type + ".root"
outputFile = "file:geantoutput" + options.geometry + options.type + ".root"

if (options.geometry == "D88"):
if (options.geometry == "D105"):
geomFile = 'Validation/HGCalValidation/data/wafer_v16.csv'
elif (options.geometry == "D93"):
elif (options.geometry == "D114"):
geomFile = 'Validation/HGCalValidation/data/wafer_v17.csv'
else:
geomFile = 'Validation/HGCalValidation/data/wafer_v17.csv'
33 changes: 21 additions & 12 deletions Validation/HGCalValidation/scripts/testHGCalDigi_cfg.py
Original file line number Diff line number Diff line change
@@ -1,9 +1,9 @@
###############################################################################
# Way to use this:
# cmsRun testHGCalDigi_cfg.py geometry=D92 type=DDD data=mu noise=none
# cmsRun testHGCalDigi_cfg.py geometry=D111 type=DDD data=mu noise=none
# threshold=none
#
# Options for geometry: D88, D92, D93, V17Shift, V18
# Options for geometry: D105, D111, D114, V17Shift, D104
# type: DDD, DD4hep
# data: mu, tt
# noise: none, ok
@@ -18,10 +18,10 @@
### SETUP OPTIONS
options = VarParsing.VarParsing('standard')
options.register('geometry',
"D92",
"D111",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"geometry of operations: D88, D92, D93, V17Shift, V18")
"geometry of operations: D105, D111, D114, V17Shift, D104")
options.register('type',
"DDD",
VarParsing.VarParsing.multiplicity.singleton,
@@ -50,26 +50,35 @@

####################################################################
# Use the options
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
if (options.type == "DD4hep"):
from Configuration.ProcessModifiers.dd4hep_cff import dd4hep
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"

globalTag = "auto:phase2_realistic_T21"
globalTag = "auto:phase2_realistic_T33"
inFile = "file:step1" + options.type + options.geometry + options.data + ".root"
outFile = "file:step2" + options.type + options.geometry + options.data + ".root"

46 changes: 23 additions & 23 deletions Validation/HGCalValidation/scripts/testHGCalMTRecoStudy_cfg.py
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
###############################################################################
# Way to use this:
# cmsRun runHGCalRecHitStudy_cfg.py geometry=D82
# cmsRun runHGCalRecHitStudy_cfg.py geometry=D110
#
# Options for geometry D88, D92, D93
# Options for geometry D105, D110, D114
#
###############################################################################
import FWCore.ParameterSet.Config as cms
@@ -11,15 +11,15 @@

####################################################################
### Run as
## $cmsRun testHGCalMTRecoStudy_cfg.py geometry=D88 layers=1
## $cmsRun testHGCalMTRecoStudy_cfg.py geometry=D105 layers=1

### SETUP OPTIONS
options = VarParsing.VarParsing('standard')
options.register('geometry',
"D92",
"D110",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"geometry of operations: D88, D92, D93")
"geometry of operations: D105, D110, D114")

options.register('layers',
"1",
@@ -39,23 +39,23 @@

fileInput = "file:step3.root"

if (options.geometry == "D88"):
from Configuration.Eras.Era_Phase2C11I13M9_cff import Phase2C11I13M9
process = cms.Process('HGCalMTReco',Phase2C11I13M9)
process.load('Configuration.Geometry.GeometryExtended2026D88Reco_cff')
outputFile = 'file:recoutputD88.root'
elif (options.geometry == "D93"):
from Configuration.Eras.Era_Phase2C11I13M9_cff import Phase2C11I13M9
process = cms.Process('HGCalMTReco',Phase2C11I13M9)
process.load('Configuration.Geometry.GeometryExtended2026D93Reco_cff')
outputFile = 'file:recoutputD93.root'
elif (options.geometry == "D92"):
from Configuration.Eras.Era_Phase2C11I13M9_cff import Phase2C11I13M9
process = cms.Process('HGCalMTReco',Phase2C11I13M9)
process.load('Configuration.Geometry.GeometryExtended2026D92Reco_cff')
outputFile = 'file:recoutputD92.root'
if (options.geometry == "D105"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('HGCalMTReco',Phase2C17I13M9)
process.load('Configuration.Geometry.GeometryExtended2026D105Reco_cff')
outputFile = 'file:recoutputD105.root'
elif (options.geometry == "D114"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('HGCalMTReco',Phase2C17I13M9)
process.load('Configuration.Geometry.GeometryExtended2026D114Reco_cff')
outputFile = 'file:recoutputD114.root'
elif (options.geometry == "D110"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('HGCalMTReco',Phase2C17I13M9)
process.load('Configuration.Geometry.GeometryExtended2026D110Reco_cff')
outputFile = 'file:recoutputD110.root'
else:
print("Please select a valid geometry version e.g. D88, D92, D93....")
print("Please select a valid geometry version e.g. D105, D110, D114....")

print("Input file: ", fileInput)
print("Output file: ", outputFile)
@@ -67,7 +67,7 @@
process.load('Configuration.EventContent.EventContent_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T21', '')
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T33', '')

process.MessageLogger.cerr.FwkReport.reportEvery = 100

@@ -76,7 +76,7 @@
fileNames = cms.untracked.vstring(fileInput),
inputCommands = cms.untracked.vstring(
'keep *',
'drop l1tTkPrimaryVertexs_L1TkPrimaryVertex__*' # This is to skip this branch causing issue in D88 reco files with older CMSSW <= 12_4_0-pre4
'drop l1tTkPrimaryVertexs_L1TkPrimaryVertex__*' # This is to skip this branch causing issue in D105 reco files with older CMSSW <= 12_4_0-pre4
),
secondaryFileNames = cms.untracked.vstring()
)
33 changes: 21 additions & 12 deletions Validation/HGCalValidation/scripts/testHGCalReco_cfg.py
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
###############################################################################
# Way to use this:
# cmsRun testHGCalDigi_cfg.py geometry=D92 type=DDD data=mu tag=Def
# cmsRun testHGCalDigi_cfg.py geometry=D110 type=DDD data=mu tag=Def
#
# Options for geometry: D88, D92, D93, V17Shift, V18
# Options for geometry: D105, D110, D114, V17Shift, D104
# type: DDD, DD4hep
# data: mu, tt
# tag: Def, Thr, 0Noise
@@ -16,10 +16,10 @@
### SETUP OPTIONS
options = VarParsing.VarParsing('standard')
options.register('geometry',
"D92",
"D110",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"geometry of operations: D88, D92, D93, V17Shift, V18")
"geometry of operations: D105, D110, D114, V17Shift, D104")
options.register('type',
"DDD",
VarParsing.VarParsing.multiplicity.singleton,
@@ -42,25 +42,34 @@

####################################################################
# Use the options
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
if (options.type == "DD4hep"):
from Configuration.ProcessModifiers.dd4hep_cff import dd4hep
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
globalTag = "auto:phase2_realistic_T21"
globalTag = "auto:phase2_realistic_T33"
inFile = "file:step2" + options.type + options.geometry + options.data + ".root"
outFile = "file:step3" + options.type + options.geometry + options.data + ".root"
fileName = "missing" + options.type + options.geometry + options.data + options.tag + ".root"
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
###############################################################################
# Way to use this:
# cmsRun testHGCalSingleMuonPt100_cfg.py geometry=D92 type=DDD
# cmsRun testHGCalSingleMuonPt100_cfg.py geometry=D110 type=DDD
#
# Options for geometry: D88, D92, D93, V17Shift, V18
# Options for geometry: D105, D110, D114, V17Shift, D104
# type: DDD, DD4hep
#
###############################################################################
@@ -14,10 +14,10 @@
### SETUP OPTIONS
options = VarParsing.VarParsing('standard')
options.register('geometry',
"D92",
"D110",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"geometry of operations: D88, D92, D93, V17Shift, V18")
"geometry of operations: D105, D110, D114, V17Shift, D104")
options.register('type',
"DDD",
VarParsing.VarParsing.multiplicity.singleton,
@@ -32,26 +32,35 @@
####################################################################
# Use the options

from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
if (options.type == "DD4hep"):
from Configuration.ProcessModifiers.dd4hep_cff import dd4hep
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"

globalTag = "auto:phase2_realistic_T21"
globalTag = "auto:phase2_realistic_T33"
outFile = "file:step1" + options.type + options.geometry + "mu.root"

print("Geometry file: ", geomFile)
33 changes: 21 additions & 12 deletions Validation/HGCalValidation/scripts/testHGCalSimTTBar_cfg.py
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
###############################################################################
# Way to use this:
# cmsRun testHGCalTTBar_cfg.py geometry=D92 type=DDD
# cmsRun testHGCalTTBar_cfg.py geometry=D110 type=DDD
#
# Options for geometry: D88, D92, D93, V17Shift, V18
# Options for geometry: D105, D110, D114, V17Shift, D104
# type: DDD, DD4hep
#
###############################################################################
@@ -14,10 +14,10 @@
### SETUP OPTIONS
options = VarParsing.VarParsing('standard')
options.register('geometry',
"D92",
"D110",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"geometry of operations: D88, D92, D93, V17Shift, V18")
"geometry of operations: D105, D110, D114, V17Shift, D104")
options.register('type',
"DDD",
VarParsing.VarParsing.multiplicity.singleton,
@@ -32,26 +32,35 @@
####################################################################
# Use the options

from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
if (options.type == "DD4hep"):
from Configuration.ProcessModifiers.dd4hep_cff import dd4hep
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.type + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9,dd4hep)
geomFile = "Configuration.Geometry.Geometry" + options.type +"Extended2026" + options.geometry + "Reco_cff"
else:
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
if (options.geometry == "V17Shift"):
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "V18"):
geomFile = "Geometry.HGCalCommonData.testHGCal" + options.geometry + "Reco_cff"
elif (options.geometry == "D104"):
from Configuration.Eras.Era_Phase2C22I13M9_cff import Phase2C22I13M9
process = cms.Process('SingleMuonSim',Phase2C22I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
else:
from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
process = cms.Process('SingleMuonSim',Phase2C17I13M9)
geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"

globalTag = "auto:phase2_realistic_T21"
globalTag = "auto:phase2_realistic_T33"
outFile = "file:step1" + options.type + options.geometry + "tt.root"

print("Geometry file: ", geomFile)