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

Comparing UserCode/MitPhysics/Utils/interface/MetTools.h (file contents):
Revision 1.8 by sixie, Sun Jan 22 17:17:29 2012 UTC vs.
Revision 1.9 by ceballos, Sat Apr 28 19:10:01 2012 UTC

# Line 64 | Line 64 | namespace mithep {
64      Met  GetMinimumTrackMet                       (const Met *UncorrectedMet);
65      Met  GetCorrectedMet()         { return fCorrectedMet; }
66      Met  GetCorrectedTrackMet()    { return fCorrectedTrackMet; }
67 +    Met  GetCHSMet()               { return fCHSMet; }
68 +    Met  GetNHSMet()               { return fNHSMet; }
69      FourVectorM  Recoil()          { return fRecoil; }
70      FourVectorM  ChargedRecoil()   { return fChargedRecoil; }
71  
# Line 77 | Line 79 | namespace mithep {
79    private:
80      Met fCorrectedMet;
81      Met fCorrectedTrackMet;
82 +    Met fCHSMet;
83 +    Met fNHSMet;
84      FourVectorM fRecoil;
85      FourVectorM fChargedRecoil;
86      
# Line 87 | Line 91 | namespace mithep {
91    double MetTools::GetProjectedMet(const V *fV, const Met *UncorrectedMet) {
92      double projectedMet = UncorrectedMet->Pt();
93      double minDPhi = 999;
90    int index = -1;
94      for (UInt_t m = 0; m < fV->GetEntries(); ++m) {
95        if (fabs(MathUtils::DeltaPhi(UncorrectedMet->Phi(), fV->At(m)->Phi())) < minDPhi) {
96          minDPhi = fabs(MathUtils::DeltaPhi(UncorrectedMet->Phi(), fV->At(m)->Phi()));
94        index = m;
97        }
98      }
99      if (minDPhi < TMath::Pi()/2.) return projectedMet = projectedMet * sin(minDPhi);
# Line 102 | Line 104 | namespace mithep {
104    double MetTools::GetProjectedMet(const V *fV) {
105      double projectedMet = fCorrectedMet.Pt();
106      double minDPhi = 999;
105    int index = -1;
107      for (UInt_t m = 0; m < fV->GetEntries(); ++m) {
108        if (fabs(MathUtils::DeltaPhi(fCorrectedMet.Phi(), fV->At(m)->Phi())) < minDPhi) {
109          minDPhi = fabs(MathUtils::DeltaPhi(fCorrectedMet.Phi(), fV->At(m)->Phi()));
109        index = m;
110        }
111      }
112      if (minDPhi < TMath::Pi()/2.) return projectedMet = projectedMet * sin(minDPhi);
# Line 117 | Line 117 | namespace mithep {
117    double MetTools::GetProjectedTrackMet(const V *fV) {
118      double projectedMet = fCorrectedTrackMet.Pt();
119      double minDPhi = 999;
120    int index = -1;
120      for (UInt_t m = 0; m < fV->GetEntries(); ++m) {
121        if (fabs(MathUtils::DeltaPhi(fCorrectedTrackMet.Phi(), fV->At(m)->Phi())) < minDPhi) {
122          minDPhi = fabs(MathUtils::DeltaPhi(fCorrectedTrackMet.Phi(), fV->At(m)->Phi()));
124        index = m;
123        }
124      }
125      if (minDPhi < TMath::Pi()/2.) return projectedMet = projectedMet * sin(minDPhi);

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines