-
Notifications
You must be signed in to change notification settings - Fork 3
/
muDpgNtuples_cfg.py
executable file
·157 lines (125 loc) · 6.3 KB
/
muDpgNtuples_cfg.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
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
import FWCore.ParameterSet.Config as cms
import FWCore.ParameterSet.VarParsing as VarParsing
from Configuration.StandardSequences.Eras import eras
from Configuration.AlCa.GlobalTag import GlobalTag
import os
import subprocess
import sys
options = VarParsing.VarParsing()
options.register('globalTag',
'124X_dataRun3_Prompt_v4',
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"Global Tag")
options.register('nEvents',
-1, #to run on a sub-sample
#-1, #default value
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.int,
"Maximum number of processed events")
options.register('isMC',
False, #default value
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.bool,
"Dataset is MC")
options.register('SaveVFATStatus',
True, #default value
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.bool,
"Save VFAT status info from unpacker")
options.register('GE21',
True, #default value
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.bool,
"Save GE21 data in the ntuples")
options.register('inputFolder',
#/eos/cms/store/
#"/eos/cms/store/data/Run2022D/Muon/RAW-RECO/ZMu-PromptReco-v2/000/357/734/00000/",
"/eos/cms/store/data/Run2022D/Muon/RAW/v1/000/357/542/00000/",
#"/eos/cms/store/group/dpg_gem/comm_gem/reRECO/SingleMuon/GEM-reRECO-GEM-only__Run2022B-ZMu-PromptReco-v1__RAW-RECO/220721_151149/0000/",
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"EOS folder with input files")
options.register('secondaryInputFolder',
'', #default value
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"EOS folder with input files for secondary files")
# options.register('fileNumber',
# "1", #default value
# VarParsing.VarParsing.multiplicity.singleton,
# VarParsing.VarParsing.varType.string,
# "FileNumber to be processed")
options.register('ntupleName',
'MuDPGNtuple', #default value
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"Name for output ntuple")
options.parseArguments()
process = cms.Process("MUNTUPLES",eras.Run3)#Run2_2018)
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
#SkipEvent = cms.untracked.vstring('ProductNotFound')
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True))
process.MessageLogger.cerr.FwkReport.reportEvery = 100
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(options.nEvents))
#process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
process.GlobalTag.globaltag = cms.string(options.globalTag)
process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(),
secondaryFileNames = cms.untracked.vstring()
)
if "eos/cms" in options.inputFolder:
#files = subprocess.check_output(['xrdfs', 'root://eoscms.cern.ch/', 'ls', options.inputFolder]) ## Did work with CMSSW 11XX, not anymore w CMSSW 12
files = os.listdir(options.inputFolder)
#process.source.fileNames = ["file:"+options.inputFolder + f for f in files if "07a64f0e-25eb-40b6-b2a6-e8971a4e0ce8.root" in f]
process.source.fileNames = ["root://xrootd-cms.infn.it//store/data/Run2022E/Muon/RAW-RECO/ZMu-PromptReco-v1/000/360/019/00001/4933b7ab-622b-4b80-87ee-8165f976a332.root"]
elif "/xrd/" in options.inputFolder:
files = subprocess.check_output(['xrdfs', 'root://cms-xrdr.sdfarm.kr/', 'ls', options.inputFolder])
process.source.fileNames = ["root://cms-xrdr.sdfarm.kr//" +f for f in files.split()]
else:
files = subprocess.check_output(['ls', options.inputFolder])
process.source.fileNames = ["file://" + options.inputFolder + "/" + f for f in files.split()]
if options.secondaryInputFolder != "" :
files = subprocess.check_output(["ls", options.secondaryInputFolder])
process.source.secondaryFileNames = ["file://" + options.secondaryInputFolder + "/" + f for f in files.split()]
process.TFileService = cms.Service('TFileService',
fileName = cms.string(options.ntupleName+".root")
)
process.load('Configuration/StandardSequences/GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load('RecoLocalMuon.GEMRecHit.gemRecHits_cfi')
process.load("TrackingTools/TransientTrack/TransientTrackBuilder_cfi")
process.load('TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorAny_cfi')
process.load('TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorAlong_cfi')
process.load('TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorOpposite_cfi')
process.load('Configuration.StandardSequences.RawToDigi_Data_cff')
process.load('MuDPGAnalysis.MuonDPGNtuples.muNtupleProducer_cfi')
process.muNtupleProducer.isMC = cms.bool(options.isMC)
process.muNtupleProducer.SaveVFATStatus = cms.bool(options.SaveVFATStatus)
if options.SaveVFATStatus:
processDigi = True
if options.GE21:
processRechit = True
else:
processRechit = False
else:
processDigi = False
processRechit = False
if options.GE21:
processDigi = True
processRechit = True
if processRechit:
process.gemRecHits.ge21Off = cms.bool(False)
process.p = cms.Path(
process.muonGEMDigis *
process.gemRecHits *
process.muNtupleProducer)
elif processDigi and not processRechit:
process.p = cms.Path(
process.muonGEMDigis *
process.muNtupleProducer)
else:
process.p = cms.Path(
process.muNtupleProducer)