|
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
|
import FWCore.ParameterSet.Config as cms
puppiCentral = cms.VPSet(
cms.PSet(
algoId = cms.int32(5), #0 is default Puppi
useCharged = cms.bool(True),
applyLowPUCorr = cms.bool(True),
combOpt = cms.int32(0),
cone = cms.double(0.4),
rmsPtMin = cms.double(0.1),
rmsScaleFactor = cms.double(1.0)
)
)
puppiForward = cms.VPSet(
cms.PSet(
algoId = cms.int32(5), #0 is default Puppi
useCharged = cms.bool(False),
applyLowPUCorr = cms.bool(True),
combOpt = cms.int32(0),
cone = cms.double(0.4),
rmsPtMin = cms.double(0.5),
rmsScaleFactor = cms.double(1.0)
)
)
# from cfipython/, generated by PuppiProducer::fillDescriptions
from CommonTools.ParticleFlow.pfNoPileUpJME_cff import primaryVertexAssociationJME
import CommonTools.PileupAlgos.PuppiProducer_cfi as _mod
puppi = _mod.PuppiProducer.clone(
UseDeltaZCutForPileup = False,
DeltaZCut = primaryVertexAssociationJME.assignment.maxDzForPrimaryAssignment,
EtaMinUseDeltaZ = primaryVertexAssociationJME.assignment.EtaMinUseDz,
PtMaxCharged = primaryVertexAssociationJME.assignment.PtMaxCharged,
NumOfPUVtxsForCharged = primaryVertexAssociationJME.assignment.NumOfPUVtxsForCharged,
DeltaZCutForChargedFromPUVtxs = primaryVertexAssociationJME.assignment.DzCutForChargedFromPUVtxs,
PtMaxNeutralsStartSlope = 20.,
PtMaxPhotons = 20.,
UseFromPV2Recovery = True,
PtMinForFromPV2Recovery = 4.,
clonePackedCands = False, # should only be set to True for MiniAOD
algos = {
0: dict(
puppiAlgos = puppiCentral
)
} #only one PSet is defined in reference
)
puppi.algos.append(cms.PSet(
etaMin = cms.vdouble( 2.5, 3.0),
etaMax = cms.vdouble( 3.0, 10.0),
ptMin = cms.vdouble( 0.0, 0.0),
MinNeutralPt = cms.vdouble( 1.7, 2.0),
MinNeutralPtSlope = cms.vdouble(0.08, 0.08),
RMSEtaSF = cms.vdouble(1.20, 0.95),
MedEtaSF = cms.vdouble(0.90, 0.75),
EtaMaxExtrap = cms.double( 2.0),
puppiAlgos = puppiForward
)
)
# cms.PSet(
# etaMin = cms.double(3.0),
# etaMax = cms.double(10.0),
# ptMin = cms.double(0.0),
# MinNeutralPt = cms.double(2.0),
# MinNeutralPtSlope = cms.double(0.07),
# # RMSEtaSF = cms.double(1.18),
# # MedEtaSF = cms.double(0.4397),
# RMSEtaSF = cms.double(1.10),
# MedEtaSF = cms.double(0.90),
# EtaMaxExtrap = cms.double(2.0),
# puppiAlgos = puppiForward
# )
from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
phase2_common.toModify(
puppi,
EtaMinUseDeltaZ = 4.0,
DeltaZCut = 0.1,
algos = cms.VPSet(
cms.PSet(
etaMin = cms.vdouble(0., 2.5),
etaMax = cms.vdouble(2.5, 3.5),
ptMin = cms.vdouble(0., 0.), #Normally 0
MinNeutralPt = cms.vdouble(0.2, 0.2),
MinNeutralPtSlope = cms.vdouble(0.015, 0.030),
RMSEtaSF = cms.vdouble(1.0, 1.0),
MedEtaSF = cms.vdouble(1.0, 1.0),
EtaMaxExtrap = cms.double(2.0),
puppiAlgos = puppiCentral
), cms.PSet(
etaMin = cms.vdouble( 3.5),
etaMax = cms.vdouble(10.0),
ptMin = cms.vdouble( 0.), #Normally 0
MinNeutralPt = cms.vdouble( 2.0),
MinNeutralPtSlope = cms.vdouble(0.08),
RMSEtaSF = cms.vdouble(1.0 ),
MedEtaSF = cms.vdouble(0.75),
EtaMaxExtrap = cms.double( 2.0),
puppiAlgos = puppiForward
)
)
)
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
pp_on_AA.toModify(puppi, algos = [])
puppiNoLep = puppi.clone(
puppiNoLep = True
)
|