-
Notifications
You must be signed in to change notification settings - Fork 0
/
ntProd.py
87 lines (73 loc) · 3.41 KB
/
ntProd.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
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
process.load("Configuration.StandardSequences.MagneticField_cff")
#process.load("TrackPropagation.SteppingHelixPropagator.SteppingHelixPropagatorAny_cfi")
process.load("Configuration.StandardSequences.Geometry_cff")
process.load("Geometry.CommonDetUnit.globalTrackingGeometry_cfi")
process.load("Geometry.CaloEventSetup.CaloGeometry_cff")
process.source = cms.Source("PoolSource",
#debugFlag = cms.untracked.bool(False),
#debugVebosity = cms.untracked.uint32(10),
#skipEvents = cms.untracked.uint32(1683),
fileNames = cms.untracked.vstring(
#'dcache:/pnfs/cms/WAX/resilient/sckao/PAT361/tt_pat356_1_1.root'
'file:/uscms_data/d2/sckao/tt_pat387D6.root'
), duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
)
# replace the source files from a file list
#import TopPhysics.TtAnalysis.ttPATfile_list as fileList
#process.source.fileNames = fileList.fileNames
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1)
)
process.MessageLogger = cms.Service("MessageLogger")
process.ttNtp = cms.EDAnalyzer("TtNtupleProd",
# General Setup
debug = cms.untracked.bool(False),
isData = cms.untracked.bool(False),
eventId = cms.untracked.int32(0),
#rootFileName = cms.untracked.string('/uscms_data/d2/sckao/Ntp361/tt_361_0jpfA.root'),
rootFileName = cms.untracked.string('tt_387_0j.root'),
numberOfJets = cms.untracked.int32(0),
# primary vertex selection
pvSource = cms.InputTag('offlinePrimaryVertices'),
pvNDOF = cms.untracked.double(4),
pvMaxZ = cms.untracked.double(15.),
pvMaxRho = cms.untracked.double(2.),
beamSpotSource = cms.InputTag('offlineBeamSpot'),
# trigger Source
trigSource = cms.InputTag('patTriggerEvent'),
trigTag = cms.untracked.string('HLT_Mu9'),
# Jet/MET Setup
#btag = cms.untracked.bool(False),
bTagCut = cms.untracked.double(5),
bTagAlgo = cms.untracked.string('trackCountingHighEffBJetTags'),
## ( pT, eta, JES emF )
jetSetup = cms.vdouble( 20, 2.4, 1.00, 0.01 ),
jetSource = cms.InputTag("selectedPatJetsPF"),
metSource = cms.InputTag("patMETsPF"),
#jetSource = cms.InputTag("selectedPatJets"),
#metSource = cms.InputTag("patMETs"),
#caloSource = cms.InputTag("towerMaker"),
recoMetSource = cms.InputTag("tcMet"),
#genJetSource = cms.InputTag("ak5GenJets"),
# Muon Setup
## ( pT, eta, Iso, nHits, chi2, d0(Bsp) )
muSetup = cms.vdouble( 15, 2.1, 0.15, 11, 10., 0.02 ),
muonSource = cms.InputTag("selectedPatMuonsPF"),
recoMuons = cms.untracked.string('muons'),
# e/gamma veto Setup
## ( pT, eta, Iso, H/E, E/P )
eleSetup = cms.vdouble(15, 2.5, 0.2, 0.02, 0.8 ),
#eleSetup = cms.vdouble(15, 2.5, 0.2, 0.02, 0.8 ),
electronSource = cms.InputTag("selectedPatElectronsPF"),
photonSource = cms.InputTag("selectedPatPhotons"),
# other setup
genParticles = cms.InputTag("genParticles"),
)
process.p = cms.Path( process.ttNtp )
#process.ttAna.TrackAssociatorParameters.useEcal = False
#process.ttAna.TrackAssociatorParameters.useHcal = False
#process.ttAna.TrackAssociatorParameters.useCalo = True
#process.ttAna.TrackAssociatorParameters.useHO = False
#process.ttAna.TrackAssociatorParameters.useMuon = False