CMS 3D CMS Logo

zSelection_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 zSelectionLoose = cms.PSet(
4  cut = cms.string("charge = 0 & daughter(0).pt > 15 & daughter(1).pt > 15 & abs(daughter(0).eta)<2.4 & abs(daughter(1).eta)<2.4 & mass > 0"),
5  isoCut = cms.double(1000.),
6  ptThreshold = cms.untracked.double(1.5),
7  etEcalThreshold = cms.untracked.double(0.2),
8  etHcalThreshold = cms.untracked.double(0.5),
9  deltaRVetoTrk = cms.untracked.double(0.015),
10  deltaRTrk = cms.untracked.double(0.3),
11  deltaREcal = cms.untracked.double(0.25),
12  deltaRHcal = cms.untracked.double(0.25),
13  alpha = cms.untracked.double(0.),
14  beta = cms.untracked.double(-0.75),
15  relativeIsolation = cms.bool(False)
16 
17 # For standard isolation (I_Tkr<3GeV) choose this configuration:
18 # isoCut = cms.double(3.),
19 # ptThreshold = cms.untracked.double(1.5),
20 # etEcalThreshold = cms.untracked.double(0.2),
21 # etHcalThreshold = cms.untracked.double(0.5),
22 # deltaRVetoTrk = cms.untracked.double(0.015),
23 # deltaRTrk = cms.untracked.double(0.3),
24 # deltaREcal = cms.untracked.double(0.25),
25 # deltaRHcal = cms.untracked.double(0.25),
26 # alpha = cms.untracked.double(0.),
27 # beta = cms.untracked.double(-0.75),
28 # relativeIsolation = cms.bool(False)
29  )
30 
31 
32 ##### I = alpha /2 (( 1 + beta) HCal + (1 - beta) Ecal ) + (1 - alpha)Trk
33 
34 ####### combined isolation
35 #zSelection = cms.PSet(
36 # cut = cms.string("charge = 0 & daughter(0).pt > 20. & daughter(1).pt > 20. & abs(daughter(0).eta)<2.1 & abs(daughter(1).eta)<2.1 & mass > 0"),
37 # isoCut = cms.double(.45), ### with alpha = 2/3 and beta =0, so 0.45 is equivalent to 0.15......
38 # ptThreshold = cms.untracked.double(0.),
39 # etEcalThreshold = cms.untracked.double(0.),
40 # etHcalThreshold = cms.untracked.double(0.),
41 # deltaRVetoTrk = cms.untracked.double(0.01),
42 # deltaRTrk = cms.untracked.double(0.3),
43 # deltaREcal = cms.untracked.double(0.3),
44 # deltaRHcal = cms.untracked.double(0.3),
45 # alpha = cms.untracked.double(0.666667),
46 # beta = cms.untracked.double(0.0),
47 # relativeIsolation = cms.bool(True)
48 # )
49 
50 
51 #### tracker isolation
52 zSelection = cms.PSet(
53  cut = cms.string("charge = 0 & daughter(0).pt > 20. & daughter(1).pt > 20. & abs(daughter(0).eta)<2.1 & abs(daughter(1).eta)<2.1 & mass > 0"),
54  isoCut = cms.double(3.00),
55  ptThreshold = cms.untracked.double(0.),
56  etEcalThreshold = cms.untracked.double(0.),
57  etHcalThreshold = cms.untracked.double(0.),
58  deltaRVetoTrk = cms.untracked.double(0.01),
59  deltaRTrk = cms.untracked.double(0.3),
60  deltaREcal = cms.untracked.double(0.3),
61  deltaRHcal = cms.untracked.double(0.3),
62  alpha = cms.untracked.double(0.),
63  beta = cms.untracked.double(0.0),
64  relativeIsolation = cms.bool(False)
65  )
66 
67 
68 
69 
70 ### region A: |eta|<2.1, region B: 2.1< |eta| <2.4
71 
72 zSelectionABLoose = cms.PSet(
73  cut = cms.string("charge = 0 & daughter(0).pt > 15 & daughter(1).pt > 15 & ( (abs(daughter(0).eta)<2.1 & 2.1< abs(daughter(1).eta)<2.4 ) || (abs(daughter(1).eta)<2.1 & 2.1< abs(daughter(0).eta)<2.4 ) ) & mass > 0"),
74  isoCut = cms.double(1000.),
75  ptThreshold = cms.untracked.double(1.5),
76  etEcalThreshold = cms.untracked.double(0.2),
77  etHcalThreshold = cms.untracked.double(0.5),
78  deltaRVetoTrk = cms.untracked.double(0.015),
79  deltaRTrk = cms.untracked.double(0.3),
80  deltaREcal = cms.untracked.double(0.25),
81  deltaRHcal = cms.untracked.double(0.25),
82  alpha = cms.untracked.double(0.),
83  beta = cms.untracked.double(-0.75),
84  relativeIsolation = cms.bool(False)
85  )
86 
87 
88 zSelectionAB = cms.PSet(
89  cut = cms.string("charge = 0 & daughter(0).pt > 20. & daughter(1).pt > 20. & ( (abs(daughter(0).eta)<2.1 & 2.1< abs(daughter(1).eta)<2.4 ) || (abs(daughter(1).eta)<2.1 & 2.1< abs(daughter(0).eta)<2.4 ) ) & mass > 0"),
90  isoCut = cms.double(1000.),
91  ptThreshold = cms.untracked.double(1.5),
92  etEcalThreshold = cms.untracked.double(0.2),
93  etHcalThreshold = cms.untracked.double(0.5),
94  deltaRVetoTrk = cms.untracked.double(0.015),
95  deltaRTrk = cms.untracked.double(0.3),
96  deltaREcal = cms.untracked.double(0.25),
97  deltaRHcal = cms.untracked.double(0.25),
98  alpha = cms.untracked.double(0.),
99  beta = cms.untracked.double(-0.75),
100  relativeIsolation = cms.bool(False)
101  )
102 
103 
104 zSelectionBBLoose = cms.PSet(
105  cut = cms.string("charge = 0 & daughter(0).pt > 15 & daughter(1).pt > 15 & ( 2.1< abs(daughter(0).eta)<2.4 & 2.1< abs(daughter(1).eta)<2.4 ) & mass > 0"),
106  isoCut = cms.double(1000.),
107  ptThreshold = cms.untracked.double(1.5),
108  etEcalThreshold = cms.untracked.double(0.2),
109  etHcalThreshold = cms.untracked.double(0.5),
110  deltaRVetoTrk = cms.untracked.double(0.015),
111  deltaRTrk = cms.untracked.double(0.3),
112  deltaREcal = cms.untracked.double(0.25),
113  deltaRHcal = cms.untracked.double(0.25),
114  alpha = cms.untracked.double(0.),
115  beta = cms.untracked.double(-0.75),
116  relativeIsolation = cms.bool(False)
117  )
118 
119 
120 zSelectionBB = cms.PSet(
121  cut = cms.string("charge = 0 & daughter(0).pt > 20 & daughter(1).pt > 20 & ( 2.1< abs(daughter(0).eta)<2.4 & 2.1< abs(daughter(1).eta)<2.4 ) & mass > 0"),
122  isoCut = cms.double(1000.),
123  ptThreshold = cms.untracked.double(1.5),
124  etEcalThreshold = cms.untracked.double(0.2),
125  etHcalThreshold = cms.untracked.double(0.5),
126  deltaRVetoTrk = cms.untracked.double(0.015),
127  deltaRTrk = cms.untracked.double(0.3),
128  deltaREcal = cms.untracked.double(0.25),
129  deltaRHcal = cms.untracked.double(0.25),
130  alpha = cms.untracked.double(0.),
131  beta = cms.untracked.double(-0.75),
132  relativeIsolation = cms.bool(False)
133  )
134 
135 
136 
137 
138 goodZTight = cms.EDFilter(
139  "ZToMuMuIsolatedIDSelector",
140  zSelection,
141  src = cms.InputTag("goodZ"),
142  filter = cms.bool(True)
143 )
144