ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitAna/DataTree/interface/Muon.h
Revision: 1.31
Committed: Fri Jul 17 13:01:13 2009 UTC (15 years, 9 months ago) by loizides
Content type: text/plain
Branch: MAIN
CVS Tags: Mit_013pre1, Mit_012i, Mit_012h, Mit_012g, Mit_012f, Mit_012e, Mit_012d, Mit_012c, Mit_012b, Mit_012a, Mit_012, Mit_011a, Mit_011, Mit_010a, Mit_010
Changes since 1.30: +5 -5 lines
Log Message:
Renamed from TrackerOnly to Tracker.

File Contents

# User Rev Content
1 loizides 1.1 //--------------------------------------------------------------------------------------------------
2 loizides 1.31 // $Id: Muon.h,v 1.30 2009/06/19 07:41:34 loizides Exp $
3 loizides 1.1 //
4     // Muon
5 loizides 1.23 //
6 pharris 1.22 // Classified into 3 Types
7 loizides 1.23 // Standalone : Muon fit from hits in the muon chambers
8     // Tracker : Tracker matched track to muon chambers with segments
9     // Global : Combined track and muon chambers fit
10 loizides 1.1 //
11 loizides 1.23 // In the muon class all three muon classes exist and can be determined by checking
12     // whether a fit track exists for such a Muon (eg HasTrackerTrk() == kTRUE). When
13     // BestTrk function is called the global fit is returned, if that does not
14     // exist the tracker fit is returned upon that failure the standalone muon
15     // is subsequently returned. To consult more see:
16 pharris 1.22 // http://cmsdoc.cern.ch/cms/Physics/muon/MPOG/Notes/MuonReco.pdf
17     //
18     // Quality Varaibles for Selection
19     // Isolation : decomposed to IsoRNNXXXXX
20     // NN = 03,05 to denote R =0.3,0.5 in Isolation Cone
21 loizides 1.23 // XXXXX = SumPt - Sum of Pt of Tracks in Cone (using all Pt Tracks)
22     // XXXXX = EmEt - Sum of Et of Em component of Calo Towers (Veto=0.08)
23     // XXXXX = HadEt - Sum of Et of Had component of Calo Towers (Veto=0.1)
24     // XXXXX = HoEt - Sum of Pt of Ho component of Calo Towers (Veto=0.1)
25     // XXXXX = NTrk - Number of Tracks within Cone (using all Pt Tracks)
26     // XXXXX = NJet - Sum of Jets inside Cone
27     //
28 pharris 1.22 // Muon Quality Variables
29     // Extracted from the Muon Tracking Associator
30 loizides 1.23 // For more se https://twiki.cern.ch/twiki/bin/view/CMS/TrackAssociator
31 pharris 1.22 // Energy Variables : Energy within (5x5 Ecal Crystal/Hcal/Ho) along Track
32 loizides 1.23 // using calo towers
33     // S9 Energy Variables : Energy using reconstructed hits calo towers default
34 pharris 1.22 // Segment Variables :
35     // Segment variables are stored in array of 8:
36 loizides 1.23 // 0-3: Correspond to segments in 4 DT chambers in->outward(3 has no Y res)
37     // 4-7: Correspond to segments in 4 CSC chambers in->out
38 pharris 1.22 // If value undetermined 999999 is returned.
39     // Variables:
40 loizides 1.23 // DX,DY,PullX,PullY: Uncertainties/Pulls of propagated track
41     // with respect to local muon segment
42     // TrackDist,TrackDistErr: Summed Dist/Err (in X,Y) of segment with respect
43     // to propagated track
44     // NSegments : Number of segments in muon using Segment+Track Arbitration
45     // NChambers : Number of muon chambers traversed in propagated track
46     // LastHit : Returns farest (from center) station with a recorded segment
47     // LastStation : Returns farest station using Segment+Track Arbitration
48     //
49 pharris 1.22 // Muon Id Methods: Please see Muon Id Note(as of now unpublished):
50     // https://www.cmsaf.mit.edu/twiki/pub/CmsHep/HepWAnalysis/MuonID-ingo.pdf
51 loizides 1.23 // TMOneStation : Returns bool whether muon passes "Tracker Muon One StationId"
52     // -Matching criteria for one of all segments matched to muon
53     // TMLastStation : Returns bool whether muon passes "LastStation Id" criteria
54     // -Matching criteria for last most segment
55 loizides 1.1 //
56 loizides 1.16 // Authors: J.Bendavid, C.Loizides, C.Paus, P.Harris
57 loizides 1.1 //--------------------------------------------------------------------------------------------------
58 loizides 1.6
59 loizides 1.16 #ifndef MITANA_DATATREE_MUON_H
60     #define MITANA_DATATREE_MUON_H
61 loizides 1.8
62 loizides 1.16 #include "MitAna/DataTree/interface/Track.h"
63 bendavid 1.11 #include "MitAna/DataTree/interface/ChargedParticle.h"
64 pharris 1.15
65     namespace mithep {
66     class Muon : public ChargedParticle {
67 loizides 1.1 public:
68 loizides 1.20 Muon();
69 loizides 1.18
70     enum EClassType {
71     kNone, //no track assigned
72     kGlobal, //"Global"
73 loizides 1.31 kTracker, //"Tracker"
74 loizides 1.18 kSta, //"Standalone"
75     kAny //any "best" of the above
76     };
77 bendavid 1.2
78 loizides 1.23 const Track *BestTrk() const;
79 loizides 1.30 const Track *GlobalTrk() const { return fGlobalTrkRef.Obj(); }
80     const Track *StandaloneTrk() const { return fStaTrkRef.Obj(); }
81     const Track *TrackerTrk() const { return fTrkTrkRef.Obj(); }
82     const Track *Trk() const { return BestTrk(); }
83     Double_t EmEnergy() const { return fEmEnergy; }
84     Double_t EmS9Energy() const { return fEmS9Energy; }
85 loizides 1.20 Double_t GetDX(Int_t iStation) const;
86     Double_t GetDY(Int_t iStation) const;
87     Double_t GetPullX(Int_t iStation) const;
88     Double_t GetPullY(Int_t iStation) const;
89     Double_t GetTrackDist(Int_t iStation) const;
90     Double_t GetTrackDistErr(Int_t iStation) const;
91     Int_t GetNSegments(Int_t iStation) const;
92     Bool_t Has(EClassType type) const;
93 loizides 1.23 Bool_t HasTrk() const;
94 loizides 1.25 Bool_t HasGlobalTrk() const { return fGlobalTrkRef.IsValid(); }
95     Bool_t HasStandaloneTrk() const { return fStaTrkRef.IsValid(); }
96     Bool_t HasTrackerTrk() const { return fTrkTrkRef.IsValid(); }
97 loizides 1.23 Double_t HadEnergy() const { return fHadEnergy; }
98     Double_t HadS9Energy() const { return fHadS9Energy; }
99     Double_t HoEnergy() const { return fHoEnergy; }
100     Double_t HoS9Energy() const { return fHoS9Energy; }
101 loizides 1.20 EClassType Is() const;
102 bendavid 1.27 Bool_t IsGlobalMuon() const { return fIsGlobalMuon; }
103     Bool_t IsTrackerMuon() const { return fIsTrackerMuon; }
104     Bool_t IsStandaloneMuon() const { return fIsStandaloneMuon; }
105     Bool_t IsCaloMuon() const { return fIsCaloMuon; }
106 loizides 1.23 Double_t IsoR03SumPt() const { return fIsoR03SumPt; }
107     Double_t IsoR03EmEt() const { return fIsoR03EmEt; }
108     Double_t IsoR03HadEt() const { return fIsoR03HadEt; }
109     Double_t IsoR03HoEt() const { return fIsoR03HoEt; }
110 loizides 1.28 UShort_t IsoR03NTracks() const { return fIsoR03NTracks; }
111     UShort_t IsoR03NJets() const { return fIsoR03NJets; }
112 loizides 1.23 Double_t IsoR05SumPt() const { return fIsoR05SumPt; }
113     Double_t IsoR05EmEt() const { return fIsoR05EmEt; }
114     Double_t IsoR05HadEt() const { return fIsoR05HadEt; }
115     Double_t IsoR05HoEt() const { return fIsoR05HoEt; }
116 loizides 1.28 UShort_t IsoR05NTracks() const { return fIsoR05NTracks; }
117     UShort_t IsoR05NJets() const { return fIsoR05NJets; }
118     UInt_t NChambers() const { return fNTraversedChambers; }
119     UInt_t NSegments() const { return fStationMask.NBitsSet(); }
120 loizides 1.20 Int_t LastHit() const;
121 loizides 1.18 Int_t LastStation(Double_t iMaxD, Double_t iMaxP) const;
122     Int_t LastStation(Int_t iMax=8) const;
123 loizides 1.30 Int_t ObjId() const { return ObjType()*1000+Is(); }
124 loizides 1.25 EObjType ObjType() const { return kMuon; }
125 loizides 1.29 void Print(Option_t *opt="") const;
126 loizides 1.18 Bool_t PromptTight(EClassType type) const;
127 loizides 1.23 Bool_t StationBit(Int_t bit) const { return fStationMask.TestBit(bit); }
128 loizides 1.28 Bool_t TMLastStation(Double_t iDYMin = 3, Double_t iPYMin = 3,
129     Double_t iDXMin = 3, Double_t iPXMin = 3,UInt_t iN = 2) const;
130     Bool_t TMOneStation(Double_t iDYMin = 3, Double_t iPYMin = 3,
131     Double_t iDXMin = 3, Double_t iPXMin = 3,UInt_t iN = 1) const;
132     void SetGlobalTrk(const Track *t)
133     { fGlobalTrkRef = t; ClearMom(); ClearCharge(); }
134     void SetStandaloneTrk(const Track *t)
135     { fStaTrkRef = t; ClearMom(); ClearCharge(); }
136     void SetTrackerTrk(const Track *t)
137     { fTrkTrkRef = t; ClearMom(); ClearCharge(); }
138 loizides 1.23 void SetDX(Int_t iStation, Double_t iDX);
139     void SetDY(Int_t iStation, Double_t iDY);
140 loizides 1.25 void SetEmEnergy(Double_t EmEnergy) { fEmEnergy = EmEnergy; }
141     void SetEmS9Energy(Double_t EmS9Energy) { fEmS9Energy = EmS9Energy; }
142     void SetHadEnergy(Double_t HadEnergy) { fHadEnergy = HadEnergy; }
143     void SetHadS9Energy(Double_t HadS9Energy) { fHadS9Energy = HadS9Energy; }
144     void SetHoEnergy(Double_t HoEnergy) { fHoEnergy = HoEnergy; }
145     void SetHoS9Energy(Double_t HoS9Energy) { fHoS9Energy = HoS9Energy; }
146 bendavid 1.27 void SetIsGlobalMuon(Bool_t b) { fIsGlobalMuon = b; }
147     void SetIsTrackerMuon(Bool_t b) { fIsTrackerMuon = b; }
148     void SetIsStandaloneMuon(Bool_t b) { fIsStandaloneMuon = b; }
149     void SetIsCaloMuon(Bool_t b) { fIsCaloMuon = b; }
150 loizides 1.25 void SetIsoR03SumPt(Double_t isoR03SumPt) { fIsoR03SumPt = isoR03SumPt; }
151     void SetIsoR03EmEt(Double_t isoR03EmEt) { fIsoR03EmEt = isoR03EmEt; }
152     void SetIsoR03HadEt(Double_t isoR03HadEt) { fIsoR03HadEt = isoR03HadEt; }
153     void SetIsoR03HoEt(Double_t isoR03HoEt) { fIsoR03HoEt = isoR03HoEt; }
154 loizides 1.28 void SetIsoR03NTracks(UShort_t isoR03NTrk) { fIsoR03NTracks = isoR03NTrk; }
155     void SetIsoR03NJets(UShort_t isoR03NJets) { fIsoR03NJets = isoR03NJets; }
156 loizides 1.25 void SetIsoR05SumPt(Double_t isoR05SumPt) { fIsoR05SumPt = isoR05SumPt; }
157     void SetIsoR05EmEt(Double_t isoR05EmEt) { fIsoR05EmEt = isoR05EmEt; }
158     void SetIsoR05HadEt(Double_t isoR05HadEt) { fIsoR05HadEt = isoR05HadEt; }
159     void SetIsoR05HoEt(Double_t isoR05HoEt) { fIsoR05HoEt = isoR05HoEt; }
160 loizides 1.28 void SetIsoR05NTracks(UShort_t isoR05NTrk) { fIsoR05NTracks = isoR05NTrk; }
161     void SetIsoR05NJets(UShort_t isoR05NJets) { fIsoR05NJets = isoR05NJets; }
162     void SetNChambers(UShort_t iNTraCh) { fNTraversedChambers = iNTraCh; }
163 loizides 1.23 void SetNSegments(Int_t iStation, Int_t NSegments);
164 loizides 1.16 void SetPullX(Int_t iStation, Double_t iPullX);
165     void SetPullY(Int_t iStation, Double_t iPullY);
166 loizides 1.25 void SetStationMask(UInt_t iStMask) { fStationMask.SetBits(iStMask); }
167 loizides 1.16 void SetTrackDist(Int_t iStation, Double_t iDist);
168     void SetTrackDistErr(Int_t iStation, Double_t iDistErr);
169 pharris 1.15
170 loizides 1.8 protected:
171 loizides 1.25 Double_t GetMass() const { return 105.658369e-3; }
172    
173     Ref<Track> fGlobalTrkRef; //global combined track reference
174     Ref<Track> fStaTrkRef; //standalone muon track reference
175     Ref<Track> fTrkTrkRef; //tracker track reference
176 loizides 1.28 Double32_t fIsoR03SumPt; //[0,0,14]isolation size R=0.3 sum pt
177     Double32_t fIsoR03EmEt; //[0,0,14]isolation size R=0.3 em trans energy
178     Double32_t fIsoR03HadEt; //[0,0,14]isolation size R=0.3 had trans energy
179     Double32_t fIsoR03HoEt; //[0,0,14]isolation size R=0.3 ho trans energy
180     UShort_t fIsoR03NTracks; //isolation size R=0.3 number of tracks
181     UShort_t fIsoR03NJets; //isolation size R=0.3 number of jets
182     Double32_t fIsoR05SumPt; //[0,0,14]isolation size R=0.5 sum pt
183     Double32_t fIsoR05EmEt; //[0,0,14]isolation size R=0.5 em trans energy
184     Double32_t fIsoR05HadEt; //[0,0,14]isolation size R=0.5 had trans energy
185     Double32_t fIsoR05HoEt; //[0,0,14]isolation size R=0.5 ho trans energy
186     UShort_t fIsoR05NTracks; //isolation size R=0.5 number of tracks
187     UShort_t fIsoR05NJets; //isolation size R=0.5 number of jets
188     Double32_t fEmEnergy; //[0,0,14]energy deposit in ecal
189     Double32_t fHadEnergy; //[0,0,14]energy deposit in hcal
190     Double32_t fHoEnergy; //[0,0,14]energy deposit in outer hcal
191     Double32_t fEmS9Energy; //[0,0,14]energy deposit in 3x3 ecal
192     Double32_t fHadS9Energy; //[0,0,14]energy deposit in 3x3 hcal
193     Double32_t fHoS9Energy; //[0,0,14]energy deposit in 3x3 outer hcal
194     UShort_t fNTraversedChambers; //number of traversed chambers
195 loizides 1.18 BitMask8 fStationMask; //bitmap of station with tracks, 0-3 DT, 4-7 CSCs
196 loizides 1.28 Double32_t fDX[8]; //[0,0,14]uncertainty in x in given muon chamber
197     Double32_t fDY[8]; //[0,0,14]uncertainty in y in given muon chamber
198     Double32_t fPullX[8]; //[0,0,14]pull in x in given muon chamber
199     Double32_t fPullY[8]; //[0,0,14]pull in y in given muon chamber
200     Double32_t fTrackDist[8]; //[0,0,14]dist. to track in trans. plane in muon chamber
201     Double32_t fTrackDistErr[8]; //[0,0,14]error of dist. to track in trans. plane
202 loizides 1.18 Int_t fNSegments[8]; //number of segments in given muon chamber
203 loizides 1.28 Bool_t fIsGlobalMuon; //GlobalMuon algo flag
204     Bool_t fIsTrackerMuon; //TrackerMuon algo flag
205     Bool_t fIsStandaloneMuon; //StandaloneMuon algo flag
206     Bool_t fIsCaloMuon; //CaloMuon algo flag
207 loizides 1.16
208 loizides 1.7 ClassDef(Muon, 1) // Muon class
209 loizides 1.1 };
210 loizides 1.4 }
211 loizides 1.8
212 loizides 1.9 //--------------------------------------------------------------------------------------------------
213 loizides 1.10 inline const mithep::Track *mithep::Muon::BestTrk() const
214     {
215     // Return "best" track.
216    
217 loizides 1.23 if (HasGlobalTrk())
218 loizides 1.10 return GlobalTrk();
219 loizides 1.23 else if (HasTrackerTrk())
220 bendavid 1.17 return TrackerTrk();
221 loizides 1.23 else if (HasStandaloneTrk())
222 loizides 1.10 return StandaloneTrk();
223    
224 loizides 1.20 Error("BestTrk", "No track reference found, returning NULL pointer.");
225 loizides 1.10 return 0;
226     }
227    
228     //--------------------------------------------------------------------------------------------------
229 loizides 1.16 inline Double_t mithep::Muon::GetDX(Int_t iStation) const
230     {
231     // Return uncertainty in x in given chamber.
232    
233     assert(iStation >= 0 && iStation < 8);
234     return fDX[iStation];
235     }
236    
237     //--------------------------------------------------------------------------------------------------
238     inline Double_t mithep::Muon::GetDY(Int_t iStation) const
239     {
240     // Return uncertainty in y in given chamber.
241    
242     assert(iStation >= 0 && iStation < 8);
243     return fDY[iStation];
244     }
245    
246     //--------------------------------------------------------------------------------------------------
247     inline Double_t mithep::Muon::GetPullX(Int_t iStation) const
248     {
249     // Return pull in x in given chamber.
250    
251     assert(iStation >= 0 && iStation < 8);
252     return fPullX[iStation];
253     }
254    
255     //--------------------------------------------------------------------------------------------------
256     inline Double_t mithep::Muon::GetPullY(Int_t iStation) const
257     {
258     // Return pull in y in given chamber.
259    
260     assert(iStation >= 0 && iStation < 8);
261     return fPullY[iStation];
262     }
263    
264     //--------------------------------------------------------------------------------------------------
265     inline Double_t mithep::Muon::GetTrackDist(Int_t iStation) const
266     {
267     // Return track distance in given chamber.
268    
269     assert(iStation >= 0 && iStation < 8);
270     return fTrackDist[iStation];
271     }
272    
273     //--------------------------------------------------------------------------------------------------
274     inline Double_t mithep::Muon::GetTrackDistErr(Int_t iStation) const
275     {
276     // Return error of track distance in given chamber.
277    
278     assert(iStation >= 0 && iStation < 8);
279     return fTrackDistErr[iStation];
280     }
281    
282     //--------------------------------------------------------------------------------------------------
283     inline Int_t mithep::Muon::GetNSegments(Int_t iStation) const
284     {
285     // Return number of segments in chamber.
286    
287     assert(iStation >= 0 && iStation < 8);
288     return fNSegments[iStation];
289     }
290    
291     //--------------------------------------------------------------------------------------------------
292 loizides 1.18 inline Bool_t mithep::Muon::Has(EClassType type) const
293     {
294     // Return true if the muon has a track of given class.
295    
296     switch (type) {
297     case kAny:
298 loizides 1.23 if (HasTrk())
299 loizides 1.18 return kTRUE;
300     break;
301     case kGlobal:
302 loizides 1.23 if (HasGlobalTrk())
303 loizides 1.18 return kTRUE;
304     break;
305 loizides 1.31 case kTracker:
306 loizides 1.23 if (HasTrackerTrk())
307 loizides 1.18 return kTRUE;
308     break;
309     case kSta:
310 loizides 1.23 if (HasStandaloneTrk())
311 loizides 1.18 return kTRUE;
312     break;
313     case kNone:
314 loizides 1.23 if (!HasTrk())
315 loizides 1.18 return kTRUE;
316     break;
317     default:
318     break;
319     }
320    
321     return kFALSE;
322     }
323    
324     //--------------------------------------------------------------------------------------------------
325 loizides 1.23 inline Bool_t mithep::Muon::HasTrk() const
326     {
327     // Return true if the muon has assigned any kind of track.
328    
329     return (HasGlobalTrk() || HasTrackerTrk() || HasStandaloneTrk());
330     }
331    
332     //--------------------------------------------------------------------------------------------------
333 loizides 1.18 inline mithep::Muon::EClassType mithep::Muon::Is() const
334     {
335     // Return the "best" classification of the muon according to the assigned tracks.
336    
337 loizides 1.23 if (HasGlobalTrk())
338 loizides 1.18 return kGlobal;
339 loizides 1.23 else if (HasTrackerTrk())
340 loizides 1.31 return kTracker;
341 loizides 1.23 else if (HasStandaloneTrk())
342 loizides 1.18 return kSta;
343    
344     return kNone;
345     }
346    
347     //--------------------------------------------------------------------------------------------------
348     inline Int_t mithep::Muon::LastHit() const
349     {
350     // Return the last hit, or -1 if no one found.
351    
352     Int_t lId = -1;
353     for (Int_t i0 = 0; i0 < 8; ++i0) {
354     if (GetDX(i0) < 99999 || GetDY(i0) < 99999)
355     lId = i0;
356     }
357     return lId;
358     }
359    
360     //--------------------------------------------------------------------------------------------------
361     inline Int_t mithep::Muon::LastStation(Int_t iMax) const
362     {
363     // Return the last station, or -1 of no one found.
364    
365     Int_t lId = -1;
366     if(TMath::Abs(iMax) > 8)
367     iMax = 8;
368     for (Int_t i0 = 0; i0 < iMax; ++i0) {
369     if (StationBit(i0) && ((lId % 4) < (i0 % 4)))
370     lId = i0;
371     }
372     return lId;
373     }
374    
375     //--------------------------------------------------------------------------------------------------
376     inline Int_t mithep::Muon::LastStation(Double_t iMaxD, Double_t iMaxP) const
377     {
378     // Return the last station for given deviation and relative error, or -1 if no one found.
379    
380     Int_t lId = -1;
381     for (Int_t i0 = 0; i0 < 8; ++i0) {
382     if ((lId % 4) > (i0 % 4))
383     continue;
384     if (GetTrackDist(i0) < iMaxD &&
385     GetTrackDist(i0)/GetTrackDistErr(i0) < iMaxP)
386     lId = i0;
387     }
388     return lId;
389     }
390    
391     //--------------------------------------------------------------------------------------------------
392     inline Bool_t mithep::Muon::PromptTight(EClassType type) const
393     {
394     // Return whether track for given class matches tight quality criteria.
395    
396     const mithep::Track *lTrack = 0;
397     switch (type) {
398     case kAny:
399     lTrack = BestTrk();
400     break;
401     case kGlobal:
402     lTrack = GlobalTrk();
403     break;
404 loizides 1.31 case kTracker:
405 loizides 1.18 lTrack = TrackerTrk();
406     break;
407     case kSta:
408     lTrack = StandaloneTrk();
409     break;
410     default:
411     break;
412     }
413    
414     if (lTrack == 0)
415     return kFALSE;
416     if (lTrack->NHits() < 11)
417     return kFALSE;
418     if (lTrack->Chi2()/lTrack->Ndof() > 10.)
419     return kFALSE;
420     if (lTrack->D0() > 0.2)
421     return kFALSE;
422     if ((LastHit() % 4) == 0)
423     return kFALSE;
424     return kTRUE;
425     }
426    
427     //--------------------------------------------------------------------------------------------------
428     inline Bool_t mithep::Muon::TMOneStation(Double_t iDYMin, Double_t iPYMin,
429 loizides 1.28 Double_t iDXMin, Double_t iPXMin, UInt_t iN) const
430 loizides 1.18 {
431     // Check if the muon is matched to at least one station (chamber).
432    
433     if (NSegments() < iN)
434     return kFALSE; //second last one
435    
436     Bool_t pGoodX = kFALSE;
437 ceballos 1.19 Bool_t pBadY = kFALSE;
438 loizides 1.18 for (Int_t i0 = 0; i0 < 8; ++i0) {
439     if ((TMath::Abs(GetDX(i0)) < iDXMin ||
440     TMath::Abs(GetPullX(i0)) < iPXMin))
441     pGoodX = kTRUE;
442     if (pGoodX &&
443     (TMath::Abs(GetDY(i0)) < iDYMin ||
444     TMath::Abs(GetPullY(i0)) < iPYMin))
445     return kTRUE;
446     if (TMath::Abs(GetDY(i0)) < 999999)
447     pBadY = kTRUE;
448     if (i0 == 3 && pGoodX && !pBadY)
449     return kTRUE;
450     }
451     return kFALSE;
452     }
453    
454     //--------------------------------------------------------------------------------------------------
455     inline Bool_t mithep::Muon::TMLastStation(Double_t iDYMin, Double_t iPYMin,
456 loizides 1.28 Double_t iDXMin, Double_t iPXMin, UInt_t iN) const
457 loizides 1.18 {
458     // Check if the muon is matched to its last station (chamber).
459    
460     if (NSegments() < iN)
461     return kFALSE; //second last one
462    
463     Int_t lLast = LastStation(-3.,-3.); // last required station
464     if (lLast < 0)
465     return kFALSE;
466     if (GetDX(lLast) > 9999.)
467     return kFALSE;
468     lLast = LastStation(); //no requirements
469     if (lLast < 0)
470     return kFALSE;
471     if (!(TMath::Abs(GetDX(lLast)) < iDXMin ||
472     TMath::Abs(GetPullX(lLast)) < iPXMin))
473     return kFALSE;
474     if (lLast == 3)
475     lLast = LastStation(3);
476     if (lLast < 0)
477     return kFALSE;
478     if (!(TMath::Abs(GetDY(lLast)) < iDYMin ||
479     TMath::Abs(GetPullY(lLast)) < iPYMin))
480     return kFALSE;
481     return kTRUE;
482     }
483    
484     //--------------------------------------------------------------------------------------------------
485 loizides 1.16 inline void mithep::Muon::SetDX(Int_t iStation, Double_t iDX)
486     {
487     // Return uncertainty in x in given chamber.
488    
489     assert(iStation >= 0 && iStation < 8);
490     fDX[iStation] = iDX;
491     }
492    
493     //--------------------------------------------------------------------------------------------------
494     inline void mithep::Muon::SetDY(Int_t iStation, Double_t iDY)
495     {
496     // Return uncertainty in y in given chamber.
497    
498     assert(iStation >= 0 && iStation < 8);
499     fDY[iStation] = iDY;
500     }
501    
502     //--------------------------------------------------------------------------------------------------
503     inline void mithep::Muon::SetPullX(Int_t iStation, Double_t iPullX)
504     {
505     // Set pull in x in given chamber.
506    
507     assert(iStation >= 0 && iStation < 8);
508     fPullX[iStation] = iPullX;
509     }
510    
511     //--------------------------------------------------------------------------------------------------
512     inline void mithep::Muon::SetPullY(Int_t iStation, Double_t iPullY)
513     {
514     // Set pull in y in given chamber.
515    
516     assert(iStation >= 0 && iStation < 8);
517     fPullY[iStation] = iPullY;
518     }
519    
520     //--------------------------------------------------------------------------------------------------
521     inline void mithep::Muon::SetTrackDist(Int_t iStation, Double_t iDist)
522 loizides 1.23
523 loizides 1.16 {
524     // Set track distance in given chamber.
525    
526     assert(iStation >= 0 && iStation < 8);
527     fTrackDist[iStation] = iDist;
528     }
529    
530     //--------------------------------------------------------------------------------------------------
531     inline void mithep::Muon::SetTrackDistErr(Int_t iStation, Double_t iDistErr)
532     {
533     // Set error of track distance in given chamber.
534 pharris 1.15
535 loizides 1.16 assert(iStation >= 0 && iStation < 8);
536     fTrackDistErr[iStation] = iDistErr;
537     }
538 pharris 1.15
539 loizides 1.16 //--------------------------------------------------------------------------------------------------
540     inline void mithep::Muon::SetNSegments(Int_t iStation, Int_t NSegments)
541     {
542     // Set number of segments in chamber.
543 pharris 1.15
544 loizides 1.16 assert(iStation >= 0 && iStation < 8);
545     fNSegments[iStation] = NSegments;
546     }
547     #endif
548 loizides 1.18
549    
550    
551    
552