ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/Morgan/interface/TRootMuon.h
(Generate patch)

Comparing UserCode/Morgan/interface/TRootMuon.h (file contents):
Revision 1.2 by lethuill, Wed Nov 19 19:03:16 2008 UTC vs.
Revision 1.10 by lethuill, Fri Apr 17 15:17:38 2009 UTC

# Line 3 | Line 3
3  
4   #include "../interface/TRootParticle.h"
5  
6
6   using namespace std;
7  
8   class TRootMuon : public TRootParticle
9   {
10 <        
10 >
11   public:
12 <        
12 >
13          TRootMuon() :
14 <                        TRootParticle()
15 <                        ,et_em_(0.)
16 <                        ,et_emS9_(0.)
17 <                        ,et_had_(0.)
18 <                        ,et_hadS9_(0.)
19 <                        ,et_ho_(0.)
20 <                        ,et_hoS9_(0.)
21 <                        ,caloCompatibility_(0.)
22 <                        ,isoR03_emEt_(0.)
23 <                        ,isoR03_hadEt_(0.)
24 <                        ,isoR03_hoEt_(0.)
25 <                        ,isoR03_sumPt_(0.)
26 <                        ,isoR03_nTracks_(0)
27 <                        ,isoR03_nJets_(0)
28 <                        ,isoR05_emEt_(0.)
29 <                        ,isoR05_hadEt_(0.)
30 <                        ,isoR05_hoEt_(0.)
31 <                        ,isoR05_sumPt_(0.)
32 <                        ,isoR05_nTracks_(0)
33 <                        ,isoR05_nJets_(0)
34 <                        ,energyValid_(false)
35 <                        ,matchesValid_(false)
36 <                        ,isolationValid_(false)
37 <                        {;}
38 <                        
14 >                TRootParticle()
15 >                ,et_em_(-9999.)
16 >                ,et_emS9_(-9999.)
17 >                ,et_had_(-9999.)
18 >                ,et_hadS9_(-9999.)
19 >                ,et_ho_(-9999.)
20 >                ,et_hoS9_(-9999.)
21 >                ,caloCompatibility_(-9999.)
22 >                ,isoR03_emEt_(-9999.)
23 >                ,isoR03_hadEt_(-9999.)
24 >                ,isoR03_hoEt_(-9999.)
25 >                ,isoR03_sumPt_(-9999.)
26 >                ,isoR03_nTracks_(-9999)
27 >                ,isoR03_nJets_(-9999)
28 >                ,isoR05_emEt_(-9999.)
29 >                ,isoR05_hadEt_(-9999.)
30 >                ,isoR05_hoEt_(-9999.)
31 >                ,isoR05_sumPt_(-9999.)
32 >                ,isoR05_nTracks_(-9999)
33 >                ,isoR05_nJets_(-9999)
34 >                ,energyValid_(false)
35 >                ,matchesValid_(false)
36 >                ,isolationValid_(false)
37 >                ,direction_(-9999)
38 >                ,algo_(-9999)
39 >                ,id_(-9999)
40 >                ,pixelLayersWithMeasurement_(-1)
41 >                ,stripLayersWithMeasurement_(-1)
42 >                ,d0_(-9999.)
43 >                ,d0Error_(-9999.)
44 >                ,dsz_(-9999.)
45 >                ,dszError_(-9999.)
46 >                ,normalizedChi2_(-9999.)
47 >                ,ptError_(-9999.)
48 >                ,etaError_(-9999.)
49 >                ,phiError_(-9999.)
50 >                ,ip3DSignificance_(-9999.)
51 >                {;}
52 >
53          TRootMuon(const TRootMuon& muon) :
54 <                        TRootParticle(muon)
55 <                        ,et_em_(muon.et_em_)
56 <                        ,et_emS9_(muon.et_emS9_)
57 <                        ,et_had_(muon.et_had_)
58 <                        ,et_hadS9_(muon.et_hadS9_)
59 <                        ,et_ho_(muon.et_ho_)
60 <                        ,et_hoS9_(muon.et_hoS9_)
61 <                        ,caloCompatibility_(muon.caloCompatibility_)
62 <                        ,isoR03_emEt_(muon.isoR03_emEt_)
63 <                        ,isoR03_hadEt_(muon.isoR03_hadEt_)
64 <                        ,isoR03_hoEt_(muon.isoR03_hoEt_)
65 <                        ,isoR03_sumPt_(muon.isoR03_sumPt_)
66 <                        ,isoR03_nTracks_(muon.isoR03_nTracks_)
67 <                        ,isoR03_nJets_(muon.isoR03_nJets_)
68 <                        ,isoR05_emEt_(muon.isoR05_emEt_)
69 <                        ,isoR05_hadEt_(muon.isoR05_hadEt_)
70 <                        ,isoR05_hoEt_(muon.isoR05_hoEt_)
71 <                        ,isoR05_sumPt_(muon.isoR05_sumPt_)
72 <                        ,isoR05_nTracks_(muon.isoR05_nTracks_)
73 <                        ,isoR05_nJets_(muon.isoR05_nJets_)
74 <                        ,energyValid_(muon.energyValid_)
75 <                        ,matchesValid_(muon.matchesValid_)
76 <                        ,isolationValid_(muon.isolationValid_)
77 <                        {;}
78 <                                
54 >                TRootParticle(muon)
55 >                ,et_em_(muon.et_em_)
56 >                ,et_emS9_(muon.et_emS9_)
57 >                ,et_had_(muon.et_had_)
58 >                ,et_hadS9_(muon.et_hadS9_)
59 >                ,et_ho_(muon.et_ho_)
60 >                ,et_hoS9_(muon.et_hoS9_)
61 >                ,caloCompatibility_(muon.caloCompatibility_)
62 >                ,isoR03_emEt_(muon.isoR03_emEt_)
63 >                ,isoR03_hadEt_(muon.isoR03_hadEt_)
64 >                ,isoR03_hoEt_(muon.isoR03_hoEt_)
65 >                ,isoR03_sumPt_(muon.isoR03_sumPt_)
66 >                ,isoR03_nTracks_(muon.isoR03_nTracks_)
67 >                ,isoR03_nJets_(muon.isoR03_nJets_)
68 >                ,isoR05_emEt_(muon.isoR05_emEt_)
69 >                ,isoR05_hadEt_(muon.isoR05_hadEt_)
70 >                ,isoR05_hoEt_(muon.isoR05_hoEt_)
71 >                ,isoR05_sumPt_(muon.isoR05_sumPt_)
72 >                ,isoR05_nTracks_(muon.isoR05_nTracks_)
73 >                ,isoR05_nJets_(muon.isoR05_nJets_)
74 >                ,energyValid_(muon.energyValid_)
75 >                ,matchesValid_(muon.matchesValid_)
76 >                ,isolationValid_(muon.isolationValid_)
77 >                ,direction_(muon.direction_)
78 >                ,algo_(muon.algo_)
79 >                ,id_(muon.id_)
80 >                ,pixelLayersWithMeasurement_(muon.pixelLayersWithMeasurement_)
81 >                ,stripLayersWithMeasurement_(muon.stripLayersWithMeasurement_)
82 >                ,d0_(muon.d0_)
83 >                ,d0Error_(muon.d0Error_)
84 >                ,dsz_(muon.dsz_)
85 >                ,dszError_(muon.dszError_)
86 >                ,normalizedChi2_(muon.normalizedChi2_)
87 >                ,ptError_(muon.ptError_)
88 >                ,etaError_(muon.etaError_)
89 >                ,phiError_(muon.phiError_)
90 >                ,ip3DSignificance_(muon.ip3DSignificance_)
91 >                {;}
92 >
93          TRootMuon(Double_t px, Double_t py, Double_t pz, Double_t e) :
94 <                        TRootParticle(px,py,pz,e)
95 <                        ,et_em_(0.)
96 <                        ,et_emS9_(0.)
97 <                        ,et_had_(0.)
98 <                        ,et_hadS9_(0.)
99 <                        ,et_ho_(0.)
100 <                        ,et_hoS9_(0.)
101 <                        ,caloCompatibility_(0.)
102 <                        ,isoR03_emEt_(0.)
103 <                        ,isoR03_hadEt_(0.)
104 <                        ,isoR03_hoEt_(0.)
105 <                        ,isoR03_sumPt_(0.)
106 <                        ,isoR03_nTracks_(0)
107 <                        ,isoR03_nJets_(0)
108 <                        ,isoR05_emEt_(0.)
109 <                        ,isoR05_hadEt_(0.)
110 <                        ,isoR05_hoEt_(0.)
111 <                        ,isoR05_sumPt_(0.)
112 <                        ,isoR05_nTracks_(0)
113 <                        ,isoR05_nJets_(0)
114 <                        ,energyValid_(false)
115 <                        ,matchesValid_(false)
116 <                        ,isolationValid_(false)
117 <                        {;}
118 <        
94 >                TRootParticle(px,py,pz,e)
95 >                ,et_em_(-9999.)
96 >                ,et_emS9_(-9999.)
97 >                ,et_had_(-9999.)
98 >                ,et_hadS9_(-9999.)
99 >                ,et_ho_(-9999.)
100 >                ,et_hoS9_(-9999.)
101 >                ,caloCompatibility_(-9999.)
102 >                ,isoR03_emEt_(-9999.)
103 >                ,isoR03_hadEt_(-9999.)
104 >                ,isoR03_hoEt_(-9999.)
105 >                ,isoR03_sumPt_(-9999.)
106 >                ,isoR03_nTracks_(-9999)
107 >                ,isoR03_nJets_(-9999)
108 >                ,isoR05_emEt_(-9999.)
109 >                ,isoR05_hadEt_(-9999.)
110 >                ,isoR05_hoEt_(-9999.)
111 >                ,isoR05_sumPt_(-9999.)
112 >                ,isoR05_nTracks_(-9999)
113 >                ,isoR05_nJets_(-9999)
114 >                ,energyValid_(false)
115 >                ,matchesValid_(false)
116 >                ,isolationValid_(false)
117 >                ,direction_(-9999)
118 >                ,algo_(-9999)
119 >                ,id_(-9999)
120 >                ,pixelLayersWithMeasurement_(-1)
121 >                ,stripLayersWithMeasurement_(-1)
122 >                ,d0_(-9999.)
123 >                ,d0Error_(-9999.)
124 >                ,dsz_(-9999.)
125 >                ,dszError_(-9999.)
126 >                ,normalizedChi2_(-9999.)
127 >                ,ptError_(-9999.)
128 >                ,etaError_(-9999.)
129 >                ,phiError_(-9999.)
130 >                ,ip3DSignificance_(-9999.)
131 >                {;}
132 >
133          TRootMuon(Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vtx_x, Double_t vtx_y, Double_t vtx_z) :
134 <                        TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z)
135 <                        ,et_em_(0.)
136 <                        ,et_emS9_(0.)
137 <                        ,et_had_(0.)
138 <                        ,et_hadS9_(0.)
139 <                        ,et_ho_(0.)
140 <                        ,et_hoS9_(0.)
141 <                        ,caloCompatibility_(0.)
142 <                        ,isoR03_emEt_(0.)
143 <                        ,isoR03_hadEt_(0.)
144 <                        ,isoR03_hoEt_(0.)
145 <                        ,isoR03_sumPt_(0.)
146 <                        ,isoR03_nTracks_(0)
147 <                        ,isoR03_nJets_(0)
148 <                        ,isoR05_emEt_(0.)
149 <                        ,isoR05_hadEt_(0.)
150 <                        ,isoR05_hoEt_(0.)
151 <                        ,isoR05_sumPt_(0.)
152 <                        ,isoR05_nTracks_(0)
153 <                        ,isoR05_nJets_(0)
154 <                        ,energyValid_(false)
155 <                        ,matchesValid_(false)
156 <                        ,isolationValid_(false)
157 <                        {;}
158 <        
134 >                TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z)
135 >                ,et_em_(-9999.)
136 >                ,et_emS9_(-9999.)
137 >                ,et_had_(-9999.)
138 >                ,et_hadS9_(-9999.)
139 >                ,et_ho_(-9999.)
140 >                ,et_hoS9_(-9999.)
141 >                ,caloCompatibility_(-9999.)
142 >                ,isoR03_emEt_(-9999.)
143 >                ,isoR03_hadEt_(-9999.)
144 >                ,isoR03_hoEt_(-9999.)
145 >                ,isoR03_sumPt_(-9999.)
146 >                ,isoR03_nTracks_(-9999)
147 >                ,isoR03_nJets_(-9999)
148 >                ,isoR05_emEt_(-9999.)
149 >                ,isoR05_hadEt_(-9999.)
150 >                ,isoR05_hoEt_(-9999.)
151 >                ,isoR05_sumPt_(-9999.)
152 >                ,isoR05_nTracks_(-9999)
153 >                ,isoR05_nJets_(-9999)
154 >                ,energyValid_(false)
155 >                ,matchesValid_(false)
156 >                ,isolationValid_(false)
157 >                ,direction_(-9999)
158 >                ,algo_(-9999)
159 >                ,id_(-9999)
160 >                ,pixelLayersWithMeasurement_(-1)
161 >                ,stripLayersWithMeasurement_(-1)
162 >                ,d0_(-9999.)
163 >                ,d0Error_(-9999.)
164 >                ,dsz_(-9999.)
165 >                ,dszError_(-9999.)
166 >                ,normalizedChi2_(-9999.)
167 >                ,ptError_(-9999.)
168 >                ,etaError_(-9999.)
169 >                ,phiError_(-9999.)
170 >                ,ip3DSignificance_(-9999.)
171 >                {;}
172 >
173          TRootMuon(Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vtx_x, Double_t vtx_y, Double_t vtx_z, Int_t type, Float_t charge) :
174 <                        TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z,type,charge)
175 <                        ,et_em_(0.)
176 <                        ,et_emS9_(0.)
177 <                        ,et_had_(0.)
178 <                        ,et_hadS9_(0.)
179 <                        ,et_ho_(0.)
180 <                        ,et_hoS9_(0.)
181 <                        ,caloCompatibility_(0.)
182 <                        ,isoR03_emEt_(0.)
183 <                        ,isoR03_hadEt_(0.)
184 <                        ,isoR03_hoEt_(0.)
185 <                        ,isoR03_sumPt_(0.)
186 <                        ,isoR03_nTracks_(0)
187 <                        ,isoR03_nJets_(0)
188 <                        ,isoR05_emEt_(0.)
189 <                        ,isoR05_hadEt_(0.)
190 <                        ,isoR05_hoEt_(0.)
191 <                        ,isoR05_sumPt_(0.)
192 <                        ,isoR05_nTracks_(0)
193 <                        ,isoR05_nJets_(0)
194 <                        ,energyValid_(false)
195 <                        ,matchesValid_(false)
196 <                        ,isolationValid_(false)
197 <                        {;}
198 <        
174 >                TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z,type,charge)
175 >                ,et_em_(-9999.)
176 >                ,et_emS9_(-9999.)
177 >                ,et_had_(-9999.)
178 >                ,et_hadS9_(-9999.)
179 >                ,et_ho_(-9999.)
180 >                ,et_hoS9_(-9999.)
181 >                ,caloCompatibility_(-9999.)
182 >                ,isoR03_emEt_(-9999.)
183 >                ,isoR03_hadEt_(-9999.)
184 >                ,isoR03_hoEt_(-9999.)
185 >                ,isoR03_sumPt_(-9999.)
186 >                ,isoR03_nTracks_(-9999)
187 >                ,isoR03_nJets_(-9999)
188 >                ,isoR05_emEt_(-9999.)
189 >                ,isoR05_hadEt_(-9999.)
190 >                ,isoR05_hoEt_(-9999.)
191 >                ,isoR05_sumPt_(-9999.)
192 >                ,isoR05_nTracks_(-9999)
193 >                ,isoR05_nJets_(-9999)
194 >                ,energyValid_(false)
195 >                ,matchesValid_(false)
196 >                ,isolationValid_(false)
197 >                ,direction_(-9999)
198 >                ,algo_(-9999)
199 >                ,id_(-9999)
200 >                ,pixelLayersWithMeasurement_(-1)
201 >                ,stripLayersWithMeasurement_(-1)
202 >                ,d0_(-9999.)
203 >                ,d0Error_(-9999.)
204 >                ,dsz_(-9999.)
205 >                ,dszError_(-9999.)
206 >                ,normalizedChi2_(-9999.)
207 >                ,ptError_(-9999.)
208 >                ,etaError_(-9999.)
209 >                ,phiError_(-9999.)
210 >                ,ip3DSignificance_(-9999.)
211 >                {;}
212 >
213          TRootMuon(const TLorentzVector &momentum) :
214 <                        TRootParticle(momentum)
215 <                        ,et_em_(0.)
216 <                        ,et_emS9_(0.)
217 <                        ,et_had_(0.)
218 <                        ,et_hadS9_(0.)
219 <                        ,et_ho_(0.)
220 <                        ,et_hoS9_(0.)
221 <                        ,caloCompatibility_(0.)
222 <                        ,isoR03_emEt_(0.)
223 <                        ,isoR03_hadEt_(0.)
224 <                        ,isoR03_hoEt_(0.)
225 <                        ,isoR03_sumPt_(0.)
226 <                        ,isoR03_nTracks_(0)
227 <                        ,isoR03_nJets_(0)
228 <                        ,isoR05_emEt_(0.)
229 <                        ,isoR05_hadEt_(0.)
230 <                        ,isoR05_hoEt_(0.)
231 <                        ,isoR05_sumPt_(0.)
232 <                        ,isoR05_nTracks_(0)
233 <                        ,isoR05_nJets_(0)
234 <                        ,energyValid_(false)
235 <                        ,matchesValid_(false)
236 <                        ,isolationValid_(false)
237 <                        {;}
238 <        
214 >                TRootParticle(momentum)
215 >                ,et_em_(-9999.)
216 >                ,et_emS9_(-9999.)
217 >                ,et_had_(-9999.)
218 >                ,et_hadS9_(-9999.)
219 >                ,et_ho_(-9999.)
220 >                ,et_hoS9_(-9999.)
221 >                ,caloCompatibility_(-9999.)
222 >                ,isoR03_emEt_(-9999.)
223 >                ,isoR03_hadEt_(-9999.)
224 >                ,isoR03_hoEt_(-9999.)
225 >                ,isoR03_sumPt_(-9999.)
226 >                ,isoR03_nTracks_(-9999)
227 >                ,isoR03_nJets_(-9999)
228 >                ,isoR05_emEt_(-9999.)
229 >                ,isoR05_hadEt_(-9999.)
230 >                ,isoR05_hoEt_(-9999.)
231 >                ,isoR05_sumPt_(-9999.)
232 >                ,isoR05_nTracks_(-9999)
233 >                ,isoR05_nJets_(-9999)
234 >                ,energyValid_(false)
235 >                ,matchesValid_(false)
236 >                ,isolationValid_(false)
237 >                ,direction_(-9999)
238 >                ,algo_(-9999)
239 >                ,id_(-9999)
240 >                ,pixelLayersWithMeasurement_(-1)
241 >                ,stripLayersWithMeasurement_(-1)
242 >                ,d0_(-9999.)
243 >                ,d0Error_(-9999.)
244 >                ,dsz_(-9999.)
245 >                ,dszError_(-9999.)
246 >                ,normalizedChi2_(-9999.)
247 >                ,ptError_(-9999.)
248 >                ,etaError_(-9999.)
249 >                ,phiError_(-9999.)
250 >                ,ip3DSignificance_(-9999.)
251 >                {;}
252 >
253          TRootMuon(const TLorentzVector &momentum, const TVector3 &vertex, Int_t type, Float_t charge) :
254 <                        TRootParticle(momentum, vertex, type, charge)
255 <                        ,et_em_(0.)
256 <                        ,et_emS9_(0.)
257 <                        ,et_had_(0.)
258 <                        ,et_hadS9_(0.)
259 <                        ,et_ho_(0.)
260 <                        ,et_hoS9_(0.)
261 <                        ,caloCompatibility_(0.)
262 <                        ,isoR03_emEt_(0.)
263 <                        ,isoR03_hadEt_(0.)
264 <                        ,isoR03_hoEt_(0.)
265 <                        ,isoR03_sumPt_(0.)
266 <                        ,isoR03_nTracks_(0)
267 <                        ,isoR03_nJets_(0)
268 <                        ,isoR05_emEt_(0.)
269 <                        ,isoR05_hadEt_(0.)
270 <                        ,isoR05_hoEt_(0.)
271 <                        ,isoR05_sumPt_(0.)
272 <                        ,isoR05_nTracks_(0)
273 <                        ,isoR05_nJets_(0)
274 <                        ,energyValid_(false)
275 <                        ,matchesValid_(false)
276 <                        ,isolationValid_(false)
277 <                        {;}
278 <        
254 >                TRootParticle(momentum, vertex, type, charge)
255 >                ,et_em_(-9999.)
256 >                ,et_emS9_(-9999.)
257 >                ,et_had_(-9999.)
258 >                ,et_hadS9_(-9999.)
259 >                ,et_ho_(-9999.)
260 >                ,et_hoS9_(-9999.)
261 >                ,caloCompatibility_(-9999.)
262 >                ,isoR03_emEt_(-9999.)
263 >                ,isoR03_hadEt_(-9999.)
264 >                ,isoR03_hoEt_(-9999.)
265 >                ,isoR03_sumPt_(-9999.)
266 >                ,isoR03_nTracks_(-9999)
267 >                ,isoR03_nJets_(-9999)
268 >                ,isoR05_emEt_(-9999.)
269 >                ,isoR05_hadEt_(-9999.)
270 >                ,isoR05_hoEt_(-9999.)
271 >                ,isoR05_sumPt_(-9999.)
272 >                ,isoR05_nTracks_(-9999)
273 >                ,isoR05_nJets_(-9999)
274 >                ,energyValid_(false)
275 >                ,matchesValid_(false)
276 >                ,isolationValid_(false)
277 >                ,direction_(-9999)
278 >                ,algo_(-9999)
279 >                ,id_(-9999)
280 >                ,pixelLayersWithMeasurement_(-1)
281 >                ,stripLayersWithMeasurement_(-1)
282 >                ,d0_(-9999.)
283 >                ,d0Error_(-9999.)
284 >                ,dsz_(-9999.)
285 >                ,dszError_(-9999.)
286 >                ,normalizedChi2_(-9999.)
287 >                ,ptError_(-9999.)
288 >                ,etaError_(-9999.)
289 >                ,phiError_(-9999.)
290 >                ,ip3DSignificance_(-9999.)
291 >                {;}
292 >
293          ~TRootMuon() {;}
294 <        
295 <        
294 >
295 >
296          Float_t et_em() const { return et_em_;}
297          Float_t et_emS9() const { return et_emS9_;}
298          Float_t et_had() const { return et_had_;}
# Line 222 | Line 319 | public:
319          Bool_t matchesValid() const { return matchesValid_;}
320          Bool_t isolationValid() const { return isolationValid_;}
321          
322 <        
322 >        Int_t direction() const { return direction_;}
323 >        Int_t algo() const { return algo_;}
324 >
325 >        Bool_t isGlobalMuon() const { return algo_ & 2; }
326 >        Bool_t isTrackerMuon() const { return algo_ & 4; }
327 >        Bool_t isStandAloneMuon() const { return algo_ & 8; }
328 >        Bool_t isCaloMuon() const { return algo_ & 16; }
329 >
330 >        Int_t id() const { return id_;}
331 >        Bool_t idTrackerMuonArbitrated() const { return id_ & 1; }
332 >        Bool_t idAllArbitrated() const { return id_ & 2; }
333 >        Bool_t idGlobalMuonPromptTight() const { return id_ & 4; }
334 >        Bool_t idTMLastStationLoose() const { return id_ & 8; }
335 >        Bool_t idTMLastStationTight() const { return id_ & 16; }
336 >        Bool_t idTM2DCompatibilityLoose() const { return id_ & 32; }
337 >        Bool_t idTM2DCompatibilityTight() const { return id_ & 64; }
338 >
339 >        Int_t pixelLayersWithMeasurement() const { return pixelLayersWithMeasurement_; }
340 >        Int_t stripLayersWithMeasurement() const { return stripLayersWithMeasurement_; }
341 >        Float_t d0() const { return d0_ ;}
342 >        Float_t d0Error() const { return d0Error_ ;}
343 >        Float_t dsz() const { return dsz_ ;}
344 >        Float_t dszError() const { return dszError_ ;}
345 >        Float_t normalizedChi2() const { return normalizedChi2_ ;}
346 >        Float_t ptError() const { return ptError_ ;}
347 >        Float_t etaError() const { return etaError_ ;}
348 >        Float_t phiError() const { return phiError_ ;}
349 >        Float_t ip3DSignificance() const { return ip3DSignificance_ ;}
350 >
351 >
352 >        //TObject* genMuon() const { return genMuon_.GetObject() ;}
353 >        virtual TString typeName() const { return "TRootMuon"; }
354 >
355 >
356          void setCaloEnergy(Float_t et_em, Float_t et_emS9, Float_t et_had, Float_t et_hadS9, Float_t et_ho, Float_t et_hoS9, Float_t caloCompatibility)
357          {
358                  et_em_ = et_em;
# Line 254 | Line 384 | public:
384                  isoR05_nJets_ = isoR05_nJets;
385          }
386          
387 +        void setValidity(Bool_t energyValid, Bool_t matchesValid, Bool_t isolationValid)
388 +        {
389 +                energyValid_ = energyValid;
390 +                matchesValid_ = matchesValid;
391 +                isolationValid_ = isolationValid;
392 +        }
393 +
394 +        void setDirection(Int_t direction) { direction_ = direction; }
395 +        void setAlgo(Int_t algo) { algo_ = algo; }
396 +        void setID(Int_t id) { id_ = id; }
397 +        void setID(
398 +                Int_t trackerMuonArbitrated
399 +                ,Int_t allArbitrated
400 +                ,Int_t globalMuonPromptTight
401 +                ,Int_t tmLastStationLoose
402 +                ,Int_t tmLastStationTight
403 +                ,Int_t tm2DCompatibilityLoose
404 +                ,Int_t tm2DCompatibilityTight
405 +        )
406 +        { id_ = trackerMuonArbitrated*1 + allArbitrated*2 + globalMuonPromptTight*4 + tmLastStationLoose*8 + tmLastStationTight*16 + tm2DCompatibilityLoose*32 + tm2DCompatibilityTight*64; }
407 +
408 +        void setPixelLayersWithMeasurement(Int_t pixelLayersWithMeasurement) { pixelLayersWithMeasurement_ = pixelLayersWithMeasurement; }
409 +        void setStripLayersWithMeasurement(Int_t stripLayersWithMeasurement) { stripLayersWithMeasurement_ = stripLayersWithMeasurement; }
410 +        void setD0(Float_t d0) { d0_ = d0; }
411 +        void setD0Error(Float_t d0Error) { d0Error_ = d0Error; }
412 +        void setDsz(Float_t dsz) { dsz_ = dsz; }
413 +        void setDszError(Float_t dszError) { dszError_ = dszError; }
414 +        void setNormalizedChi2(Float_t normalizedChi2) { normalizedChi2_ = normalizedChi2; }
415 +        void setPtError(Float_t ptError) { ptError_ = ptError; }
416 +        void setEtaError(Float_t etaError) { etaError_ = etaError; }
417 +        void setPhiError(Float_t phiError) { phiError_ = phiError; }
418 +        void setIP3DSignificance(Float_t ip3DSignificance) { ip3DSignificance_ = ip3DSignificance; }
419 +
420 +
421          friend std::ostream& operator<< (std::ostream& stream, const TRootMuon& muon) {
422 <                stream << "TRootMuon - Charge=" << muon.charge() << " (Et,eta,phi)=("<< muon.Et() <<","<< muon.Eta() <<","<< muon.Phi() << ")  vertex(x,y,z)=("<< muon.vx() <<","<< muon.vy() <<","<< muon.vz() << ")"
423 <                                << "  caloCompatibility="<<  muon.caloCompatibility_ << "  validity(energy, matches, isolation)=(" << muon.energyValid_ <<","<< muon.matchesValid_ <<","<< muon.isolationValid_ << ")" << endl
422 >                stream << "TRootMuon - Charge=" << muon.charge() << " (Et,eta,phi)=("<< muon.Et() <<","<< muon.Eta() <<","<< muon.Phi() << ")  vertex(x,y,z)=("<< muon.vx() <<","<< muon.vy() <<","<< muon.vz() << ")" << endl
423 >                                << "            Type(G,T,S,C)=(" << muon.isGlobalMuon() << ","  << muon.isTrackerMuon() << ","  << muon.isStandAloneMuon() << ","  << muon.isCaloMuon() << ") "
424 >                                << "  ID=(" << muon.idTrackerMuonArbitrated() << ","  << muon.idAllArbitrated() << ","  << muon.idGlobalMuonPromptTight() << ","  << muon.idTMLastStationLoose()
425 >                                << ","  << muon.idTMLastStationTight() << ","  << muon.idTM2DCompatibilityLoose() << ","  << muon.idTM2DCompatibilityTight() << ")" << "  Direction=" << muon.direction_
426 >                                << "  caloCompatibility="<<  muon.caloCompatibility_ << "  validity(energy,matches,isolation)=(" << muon.energyValid_ <<","<< muon.matchesValid_ <<","<< muon.isolationValid_ << ")" << endl
427                                  << "            Et_em=" << muon.et_em_ << " Et_emS9=" << muon.et_emS9_ << "  Et_had=" << muon.et_had_ << " Et_hadS9=" << muon.et_hadS9_ << "  Et_ho=" << muon.et_ho_ << "  Et_hoS9=" << muon.et_hoS9_ << endl
428                                  << "            isolation cone 0.3: Et_em=" << muon.isoR03_emEt_ << "  Et_had=" << muon.isoR03_hadEt_ <<"  Et_ho=" << muon.isoR03_hoEt_  << "  nTracks=" << muon.isoR03_nTracks_ <<"  nJets=" << muon.isoR03_nJets_ << endl
429 <                                << "            isolation cone 0.5: Et_em=" << muon.isoR05_emEt_ << "  Et_had=" << muon.isoR05_hadEt_ <<"  Et_ho=" << muon.isoR05_hoEt_  << "  nTracks=" << muon.isoR05_nTracks_ <<"  nJets=" << muon.isoR05_nJets_;
429 >                                << "            isolation cone 0.5: Et_em=" << muon.isoR05_emEt_ << "  Et_had=" << muon.isoR05_hadEt_ <<"  Et_ho=" << muon.isoR05_hoEt_  << "  nTracks=" << muon.isoR05_nTracks_ <<"  nJets=" << muon.isoR05_nJets_ << endl
430 >                                << "            IP 3D Significance=" << muon.ip3DSignificance_;
431                  return stream;
432          };
433  
434          
435          private:
436          
437 +        // Variables from reco::Muon
438          Float_t et_em_;
439          Float_t et_emS9_;
440          Float_t et_had_;
# Line 291 | Line 460 | public:
460          Bool_t energyValid_;
461          Bool_t matchesValid_;
462          Bool_t isolationValid_;
463 +
464 +        Int_t direction_;       //  OutsideIn = -1, Undefined = 0, InsideOut = 1
465 +        Int_t algo_; // binary => GlobalMuon=00010 , TrackerMuon=00100 , StandAloneMuon=01000 , CaloMuon=10000
466 +        // MuonId coded in binary word id_ ==> TrackerMuonArbitrated=0000001 , AllArbitrated=0000010 , GlobalMuonPromptTight=0000100 ,
467 +        // TMLastStationLoose=0001000 , TMLastStationTight=0010000 , TM2DCompatibilityLoose=0100000 , TM2DCompatibilityTight=1000000
468 +        Int_t id_;
469 +
470          
471 <        ClassDef (TRootMuon,1);
471 >        // Variables from reco::GsfTrack
472 >
473 >        /* cf: http://cmslxr.fnal.gov/lxr/source/DataFormats/TrackReco/interface/TrackBase.h
474 >        For tracks reconstructed in the CMS Tracker, the reference position is the point of closest approach to the centre
475 >        of CMS. For muons, this is not necessarily true. Parameters associated to the 5D curvilinear covariance matrix:
476 >                qoverp = q / abs(p) = signed inverse of momentum [1/GeV]
477 >                lambda = pi/2 - polar angle at the given point
478 >                phi = azimuth angle at the given point
479 >                dxy = -vx*sin(phi) + vy*cos(phi) [cm]
480 >                dsz = vz*cos(lambda) - (vx*cos(phi)+vy*sin(phi))*sin(lambda) [cm]
481 >
482 >        Geometrically, dxy is the signed distance in the XY plane between the straight line passing through (vx,vy) with
483 >        azimuthal angle phi and the point (0,0). The dsz parameter is the signed distance in the SZ plane between the straight
484 >        line passing through (vx,vy,vz) with angles (phi, lambda) and the point (s=0,z=0). The S axis is defined by the projection
485 >        of the straight line onto the XY plane. The convention is to assign the S coordinate for (vx,vy) as the value
486 >        vx*cos(phi)+vy*sin(phi). This value is zero when (vx,vy) is the point of minimum transverse distance to (0,0).
487 >
488 >        Note that dxy and dsz provide sensible estimates of the distance from the true particle trajectory to (0,0,0) ONLY
489 >        in two cases:
490 >                - When (vx,vy,vz) already correspond to the point of minimum transverse distance to (0,0,0) or is close to it
491 >                (so that the differences between considering the exact trajectory or a straight line in this range  are negligible).
492 >                This is usually true for Tracker tracks.
493 >                - When the track has infinite or extremely high momentum */
494 >
495 >        Int_t pixelLayersWithMeasurement_; // Number of pixel layers with at least one valid hit
496 >        Int_t stripLayersWithMeasurement_; // Number of strip layers with at least one valid hit
497 >
498 >        Float_t d0_;             // d0=-dxy
499 >     Float_t d0Error_;        // error on d0_
500 >     Float_t dsz_;            // dsz parameter
501 >     Float_t dszError_;       // error on dsz_
502 >     Float_t normalizedChi2_; // chi-squared divided by n.d.o.f. of track fit
503 >
504 >     Float_t ptError_; // needed ?  ptError()
505 >     Float_t etaError_; // needed ?  etaError()
506 >     Float_t phiError_; // needed ?  phiError()
507 >
508 >     Float_t ip3DSignificance_;
509 >
510 >        ClassDef (TRootMuon,3);
511   };
512  
513   #endif

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines