ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/Morgan/interface/TRootMuon.h
Revision: 1.3
Committed: Tue Dec 2 17:55:01 2008 UTC (16 years, 5 months ago) by lethuill
Content type: text/plain
Branch: MAIN
Changes since 1.2: +72 -4 lines
Log Message:
more MuonID variables

File Contents

# User Rev Content
1 mlethuil 1.1 #ifndef TRootMuon_h
2     #define TRootMuon_h
3    
4     #include "../interface/TRootParticle.h"
5    
6    
7     using namespace std;
8    
9     class TRootMuon : public TRootParticle
10     {
11    
12     public:
13    
14 lethuill 1.2 TRootMuon() :
15     TRootParticle()
16     ,et_em_(0.)
17     ,et_emS9_(0.)
18     ,et_had_(0.)
19     ,et_hadS9_(0.)
20     ,et_ho_(0.)
21     ,et_hoS9_(0.)
22     ,caloCompatibility_(0.)
23     ,isoR03_emEt_(0.)
24     ,isoR03_hadEt_(0.)
25     ,isoR03_hoEt_(0.)
26     ,isoR03_sumPt_(0.)
27     ,isoR03_nTracks_(0)
28     ,isoR03_nJets_(0)
29     ,isoR05_emEt_(0.)
30     ,isoR05_hadEt_(0.)
31     ,isoR05_hoEt_(0.)
32     ,isoR05_sumPt_(0.)
33     ,isoR05_nTracks_(0)
34     ,isoR05_nJets_(0)
35     ,energyValid_(false)
36     ,matchesValid_(false)
37     ,isolationValid_(false)
38 lethuill 1.3 ,direction_(0)
39     ,algo_(0)
40     ,id_(0)
41 lethuill 1.2 {;}
42 mlethuil 1.1
43 lethuill 1.2 TRootMuon(const TRootMuon& muon) :
44     TRootParticle(muon)
45     ,et_em_(muon.et_em_)
46     ,et_emS9_(muon.et_emS9_)
47     ,et_had_(muon.et_had_)
48     ,et_hadS9_(muon.et_hadS9_)
49     ,et_ho_(muon.et_ho_)
50     ,et_hoS9_(muon.et_hoS9_)
51     ,caloCompatibility_(muon.caloCompatibility_)
52     ,isoR03_emEt_(muon.isoR03_emEt_)
53     ,isoR03_hadEt_(muon.isoR03_hadEt_)
54     ,isoR03_hoEt_(muon.isoR03_hoEt_)
55     ,isoR03_sumPt_(muon.isoR03_sumPt_)
56     ,isoR03_nTracks_(muon.isoR03_nTracks_)
57     ,isoR03_nJets_(muon.isoR03_nJets_)
58     ,isoR05_emEt_(muon.isoR05_emEt_)
59     ,isoR05_hadEt_(muon.isoR05_hadEt_)
60     ,isoR05_hoEt_(muon.isoR05_hoEt_)
61     ,isoR05_sumPt_(muon.isoR05_sumPt_)
62     ,isoR05_nTracks_(muon.isoR05_nTracks_)
63     ,isoR05_nJets_(muon.isoR05_nJets_)
64     ,energyValid_(muon.energyValid_)
65     ,matchesValid_(muon.matchesValid_)
66     ,isolationValid_(muon.isolationValid_)
67 lethuill 1.3 ,direction_(muon.direction_)
68     ,algo_(muon.algo_)
69     ,id_(muon.id_)
70 lethuill 1.2 {;}
71 mlethuil 1.1
72 lethuill 1.2 TRootMuon(Double_t px, Double_t py, Double_t pz, Double_t e) :
73     TRootParticle(px,py,pz,e)
74     ,et_em_(0.)
75     ,et_emS9_(0.)
76     ,et_had_(0.)
77     ,et_hadS9_(0.)
78     ,et_ho_(0.)
79     ,et_hoS9_(0.)
80     ,caloCompatibility_(0.)
81     ,isoR03_emEt_(0.)
82     ,isoR03_hadEt_(0.)
83     ,isoR03_hoEt_(0.)
84     ,isoR03_sumPt_(0.)
85     ,isoR03_nTracks_(0)
86     ,isoR03_nJets_(0)
87     ,isoR05_emEt_(0.)
88     ,isoR05_hadEt_(0.)
89     ,isoR05_hoEt_(0.)
90     ,isoR05_sumPt_(0.)
91     ,isoR05_nTracks_(0)
92     ,isoR05_nJets_(0)
93     ,energyValid_(false)
94     ,matchesValid_(false)
95     ,isolationValid_(false)
96 lethuill 1.3 ,direction_(0)
97     ,algo_(0)
98     ,id_(0)
99 lethuill 1.2 {;}
100    
101     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) :
102     TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z)
103     ,et_em_(0.)
104     ,et_emS9_(0.)
105     ,et_had_(0.)
106     ,et_hadS9_(0.)
107     ,et_ho_(0.)
108     ,et_hoS9_(0.)
109     ,caloCompatibility_(0.)
110     ,isoR03_emEt_(0.)
111     ,isoR03_hadEt_(0.)
112     ,isoR03_hoEt_(0.)
113     ,isoR03_sumPt_(0.)
114     ,isoR03_nTracks_(0)
115     ,isoR03_nJets_(0)
116     ,isoR05_emEt_(0.)
117     ,isoR05_hadEt_(0.)
118     ,isoR05_hoEt_(0.)
119     ,isoR05_sumPt_(0.)
120     ,isoR05_nTracks_(0)
121     ,isoR05_nJets_(0)
122     ,energyValid_(false)
123     ,matchesValid_(false)
124     ,isolationValid_(false)
125 lethuill 1.3 ,direction_(0)
126     ,algo_(0)
127     ,id_(0)
128 lethuill 1.2 {;}
129    
130     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) :
131     TRootParticle(px,py,pz,e,vtx_x,vtx_y,vtx_z,type,charge)
132     ,et_em_(0.)
133     ,et_emS9_(0.)
134     ,et_had_(0.)
135     ,et_hadS9_(0.)
136     ,et_ho_(0.)
137     ,et_hoS9_(0.)
138     ,caloCompatibility_(0.)
139     ,isoR03_emEt_(0.)
140     ,isoR03_hadEt_(0.)
141     ,isoR03_hoEt_(0.)
142     ,isoR03_sumPt_(0.)
143     ,isoR03_nTracks_(0)
144     ,isoR03_nJets_(0)
145     ,isoR05_emEt_(0.)
146     ,isoR05_hadEt_(0.)
147     ,isoR05_hoEt_(0.)
148     ,isoR05_sumPt_(0.)
149     ,isoR05_nTracks_(0)
150     ,isoR05_nJets_(0)
151     ,energyValid_(false)
152     ,matchesValid_(false)
153     ,isolationValid_(false)
154 lethuill 1.3 ,direction_(0)
155     ,algo_(0)
156     ,id_(0)
157 lethuill 1.2 {;}
158    
159     TRootMuon(const TLorentzVector &momentum) :
160     TRootParticle(momentum)
161     ,et_em_(0.)
162     ,et_emS9_(0.)
163     ,et_had_(0.)
164     ,et_hadS9_(0.)
165     ,et_ho_(0.)
166     ,et_hoS9_(0.)
167     ,caloCompatibility_(0.)
168     ,isoR03_emEt_(0.)
169     ,isoR03_hadEt_(0.)
170     ,isoR03_hoEt_(0.)
171     ,isoR03_sumPt_(0.)
172     ,isoR03_nTracks_(0)
173     ,isoR03_nJets_(0)
174     ,isoR05_emEt_(0.)
175     ,isoR05_hadEt_(0.)
176     ,isoR05_hoEt_(0.)
177     ,isoR05_sumPt_(0.)
178     ,isoR05_nTracks_(0)
179     ,isoR05_nJets_(0)
180     ,energyValid_(false)
181     ,matchesValid_(false)
182     ,isolationValid_(false)
183 lethuill 1.3 ,direction_(0)
184     ,algo_(0)
185     ,id_(0)
186 lethuill 1.2 {;}
187 mlethuil 1.1
188 lethuill 1.2 TRootMuon(const TLorentzVector &momentum, const TVector3 &vertex, Int_t type, Float_t charge) :
189     TRootParticle(momentum, vertex, type, charge)
190     ,et_em_(0.)
191     ,et_emS9_(0.)
192     ,et_had_(0.)
193     ,et_hadS9_(0.)
194     ,et_ho_(0.)
195     ,et_hoS9_(0.)
196     ,caloCompatibility_(0.)
197     ,isoR03_emEt_(0.)
198     ,isoR03_hadEt_(0.)
199     ,isoR03_hoEt_(0.)
200     ,isoR03_sumPt_(0.)
201     ,isoR03_nTracks_(0)
202     ,isoR03_nJets_(0)
203     ,isoR05_emEt_(0.)
204     ,isoR05_hadEt_(0.)
205     ,isoR05_hoEt_(0.)
206     ,isoR05_sumPt_(0.)
207     ,isoR05_nTracks_(0)
208     ,isoR05_nJets_(0)
209     ,energyValid_(false)
210     ,matchesValid_(false)
211     ,isolationValid_(false)
212 lethuill 1.3 ,direction_(0)
213     ,algo_(0)
214     ,id_(0)
215 lethuill 1.2 {;}
216 mlethuil 1.1
217     ~TRootMuon() {;}
218    
219    
220 lethuill 1.2 Float_t et_em() const { return et_em_;}
221     Float_t et_emS9() const { return et_emS9_;}
222     Float_t et_had() const { return et_had_;}
223     Float_t et_hadS9() const { return et_hadS9_;}
224     Float_t et_ho() const { return et_ho_;}
225     Float_t et_hoS9() const { return et_hoS9_;}
226     Float_t caloCompatibility() const { return caloCompatibility_;}
227    
228     Float_t isoR03_emEt() const { return isoR03_emEt_;}
229     Float_t isoR03_hadEt() const { return isoR03_hadEt_;}
230     Float_t isoR03_hoEt() const { return isoR03_hoEt_;}
231     Float_t isoR03_sumPt() const { return isoR03_sumPt_;}
232     Int_t isoR03_nTracks() const { return isoR03_nTracks_;}
233     Int_t isoR03_nJets() const { return isoR03_nJets_;}
234    
235     Float_t isoR05_emEt() const { return isoR05_emEt_;}
236     Float_t isoR05_hadEt() const { return isoR05_hadEt_;}
237     Float_t isoR05_hoEt() const { return isoR05_hoEt_;}
238     Float_t isoR05_sumPt() const { return isoR05_sumPt_;}
239     Int_t isoR05_nTracks() const { return isoR05_nTracks_;}
240     Int_t isoR05_nJets() const { return isoR05_nJets_;}
241    
242     Bool_t energyValid() const { return energyValid_;}
243     Bool_t matchesValid() const { return matchesValid_;}
244     Bool_t isolationValid() const { return isolationValid_;}
245 mlethuil 1.1
246 lethuill 1.3 Int_t direction() const { return direction_;}
247     Int_t algo() const { return algo_;}
248    
249     Bool_t isGlobalMuon() const { return algo_ & 2; }
250     Bool_t isTrackerMuon() const { return algo_ & 4; }
251     Bool_t isStandAloneMuon() const { return algo_ & 8; }
252     Bool_t isCaloMuon() const { return algo_ & 16; }
253    
254     Int_t id() const { return id_;}
255     Bool_t idTrackerMuonArbitrated() const { return id_ & 1; }
256     Bool_t idAllArbitrated() const { return id_ & 2; }
257     Bool_t idGlobalMuonPromptTight() const { return id_ & 4; }
258     Bool_t idTMLastStationLoose() const { return id_ & 8; }
259     Bool_t idTMLastStationTight() const { return id_ & 16; }
260     Bool_t idTM2DCompatibilityLoose() const { return id_ & 32; }
261     Bool_t idTM2DCompatibilityTight() const { return id_ & 64; }
262    
263 mlethuil 1.1 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)
264     {
265     et_em_ = et_em;
266     et_emS9_ = et_emS9;
267     et_had_ = et_had;
268     et_hadS9_ = et_hadS9;
269     et_ho_ = et_ho;
270     et_hoS9_ = et_hoS9;
271     caloCompatibility_ = caloCompatibility;
272     }
273    
274     void setIsoR03(Float_t isoR03_emEt, Float_t isoR03_hadEt, Float_t isoR03_hoEt, Float_t isoR03_sumPt, Int_t isoR03_nTracks, Int_t isoR03_nJets)
275     {
276     isoR03_emEt_ = isoR03_emEt;
277     isoR03_hadEt_ = isoR03_hadEt;
278     isoR03_hoEt_ = isoR03_hoEt;
279     isoR03_sumPt_ = isoR03_sumPt;
280     isoR03_nTracks_ = isoR03_nTracks;
281     isoR03_nJets_ = isoR03_nJets;
282     }
283    
284     void setIsoR05(Float_t isoR05_emEt, Float_t isoR05_hadEt, Float_t isoR05_hoEt, Float_t isoR05_sumPt, Int_t isoR05_nTracks, Int_t isoR05_nJets)
285     {
286     isoR05_emEt_ = isoR05_emEt;
287     isoR05_hadEt_ = isoR05_hadEt;
288     isoR05_hoEt_ = isoR05_hoEt;
289     isoR05_sumPt_ = isoR05_sumPt;
290     isoR05_nTracks_ = isoR05_nTracks;
291     isoR05_nJets_ = isoR05_nJets;
292     }
293    
294 lethuill 1.3 void setValidity(Bool_t energyValid, Bool_t matchesValid, Bool_t isolationValid)
295     {
296     energyValid_ = energyValid;
297     matchesValid_ = matchesValid;
298     isolationValid_ = isolationValid;
299     }
300    
301     void setDirection(Int_t direction) { direction_ = direction; }
302     void setAlgo(Int_t algo) { algo_ = algo; }
303     void setID(Int_t id) { id_ = id; }
304     void setID(
305     Int_t trackerMuonArbitrated
306     ,Int_t allArbitrated
307     ,Int_t globalMuonPromptTight
308     ,Int_t tmLastStationLoose
309     ,Int_t tmLastStationTight
310     ,Int_t tm2DCompatibilityLoose
311     ,Int_t tm2DCompatibilityTight
312     )
313     { id_ = trackerMuonArbitrated*1 + allArbitrated*2 + globalMuonPromptTight*4 + tmLastStationLoose*8 + tmLastStationTight*16 + tm2DCompatibilityLoose*32 + tm2DCompatibilityTight*64; }
314    
315    
316 mlethuil 1.1 friend std::ostream& operator<< (std::ostream& stream, const TRootMuon& muon) {
317 lethuill 1.3 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
318     << " Type(G,T,S,C)=(" << muon.isGlobalMuon() << "," << muon.isTrackerMuon() << "," << muon.isStandAloneMuon() << "," << muon.isCaloMuon() << ") "
319     << " ID=(" << muon.idTrackerMuonArbitrated() << "," << muon.idAllArbitrated() << "," << muon.idGlobalMuonPromptTight() << "," << muon.idTMLastStationLoose()
320     << "," << muon.idTMLastStationTight() << "," << muon.idTM2DCompatibilityLoose() << "," << muon.idTM2DCompatibilityTight() << ")" << " Direction=" << muon.direction_
321     << " caloCompatibility="<< muon.caloCompatibility_ << " validity(energy,matches,isolation)=(" << muon.energyValid_ <<","<< muon.matchesValid_ <<","<< muon.isolationValid_ << ")" << endl
322 mlethuil 1.1 << " 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
323     << " 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
324     << " 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_;
325     return stream;
326     };
327    
328    
329     private:
330    
331     Float_t et_em_;
332     Float_t et_emS9_;
333     Float_t et_had_;
334     Float_t et_hadS9_;
335     Float_t et_ho_;
336     Float_t et_hoS9_;
337     Float_t caloCompatibility_;
338    
339     Float_t isoR03_emEt_;
340     Float_t isoR03_hadEt_;
341     Float_t isoR03_hoEt_;
342     Float_t isoR03_sumPt_;
343     Int_t isoR03_nTracks_;
344     Int_t isoR03_nJets_;
345    
346     Float_t isoR05_emEt_;
347     Float_t isoR05_hadEt_;
348     Float_t isoR05_hoEt_;
349     Float_t isoR05_sumPt_;
350     Int_t isoR05_nTracks_;
351     Int_t isoR05_nJets_;
352    
353     Bool_t energyValid_;
354     Bool_t matchesValid_;
355     Bool_t isolationValid_;
356 lethuill 1.3
357     Int_t direction_; // OutsideIn = -1, Undefined = 0, InsideOut = 1
358     Int_t algo_; // binary => GlobalMuon=00010 , TrackerMuon=00100 , StandAloneMuon=01000 , CaloMuon=10000
359     // MuonId coded in binary word id_ ==> TrackerMuonArbitrated=0000001 , AllArbitrated=0000010 , GlobalMuonPromptTight=0000100 ,
360     // TMLastStationLoose=0001000 , TMLastStationTight=0010000 , TM2DCompatibilityLoose=0100000 , TM2DCompatibilityTight=1000000
361     Int_t id_;
362    
363 mlethuil 1.1 ClassDef (TRootMuon,1);
364     };
365    
366     #endif
367    
368