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

Comparing UserCode/MitAna/DataTree/interface/Jet.h (file contents):
Revision 1.22 by bendavid, Tue Jul 7 08:30:28 2009 UTC vs.
Revision 1.25 by sixie, Fri Mar 26 14:31:31 2010 UTC

# Line 60 | Line 60 | namespace mithep
60                        { return fCombinedSecondaryVertexBJetTagsDisc;    }
61        Double_t      CombinedSecondaryVertexMVABJetTagsDisc()       const
62                        { return fCombinedSecondaryVertexMVABJetTagsDisc; }
63 +      Double_t      CombinedCorrectionScale()     const;
64        Double_t      CustomCorrectionScale()       const { return fCustomCorrectionScale;        }
65        void          DisableCorrection(ECorr c)          { fCorrections.ClearBit(c); ClearMom(); }
66        void          DisableCorrections()                { fCorrections.Clear(); ClearMom();     }
# Line 70 | Line 71 | namespace mithep
71        Double_t      JetBProbabilityBJetTagsDisc()                  const
72                        { return fJetBProbabilityBJetTagsDisc;    }
73        Int_t         MatchedMCFlavor()             const { return fMatchedMCFlavor;           }
74 +      virtual Jet  *MakeCopy()                    const { return new Jet(*this);             }
75        virtual
76        UInt_t        NConstituents()               const { return 0;                          }
77        UInt_t        N()                           const { return NConstituents();            }
# Line 96 | Line 98 | namespace mithep
98        Double_t      TrackCountingHighPurBJetTagsDisc()             const
99                        { return fTrackCountingHighPurBJetTagsDisc;  }
100        FourVectorM   RawMom()                      const { return fRawMom.V();                }
101 +      Double_t      SigmaEta()                    const { return fSigmaEta;                  }
102 +      Double_t      SigmaPhi()                    const { return fSigmaPhi;                  }
103        void          SetAlpha(Double_t val)              { fAlpha = val;                      }
104        void          SetBeta(Double_t val)               { fBeta  = val;                      }
105 +      void          SetSigmaEta(Double_t val)           { fSigmaEta = val;                   }
106 +      void          SetSigmaPhi(Double_t val)           { fSigmaPhi = val;                   }
107        void          SetCombinedSecondaryVertexBJetTagsDisc(Double_t d)
108                        { fCombinedSecondaryVertexBJetTagsDisc = d;    }
109        void          SetCombinedSecondaryVertexMVABJetTagsDisc(Double_t d)
# Line 139 | Line 145 | namespace mithep
145        Double32_t    fAlpha;                                  //[0,0,14]jet vertex alpha variable
146        Double32_t    fBeta;                                   //[0,0,14]jet vertex beta variable
147        Int_t         fMatchedMCFlavor;                        //[0,0,14]pdg of matched quark flavor
148 +      Double32_t    fSigmaEta;                               //[0,0,14]sqrt(etaetaMoment)
149 +      Double32_t    fSigmaPhi;                               //[0,0,14]sqrt(phiphiMoment)
150        Double32_t    fJetProbabilityBJetTagsDisc;             //[0,0,14]discriminants b-tagging algos
151        Double32_t    fJetBProbabilityBJetTagsDisc;            //[0,0,14]discriminants b-tagging algos
152        Double32_t    fSimpleSecondaryVertexBJetTagsDisc;      //[0,0,14]discriminants b-tagging algos
# Line 159 | Line 167 | namespace mithep
167        Double32_t    fCustomCorrectionScale;                  //[0,0,14]custom correction scale
168        BitMask8      fCorrections;                            //mask of corrections to be applied
169  
170 <    ClassDef(Jet, 2) // Jet class
170 >    ClassDef(Jet, 3) // Jet class
171    };
172   }
173  
174   //--------------------------------------------------------------------------------------------------
175 < inline void mithep::Jet::GetMom() const
175 > inline Double_t mithep::Jet::CombinedCorrectionScale() const
176   {
177 <  // Get raw momentum values from stored values and apply all enabled corrections.
178 <
171 <  fCachedMom.SetCoordinates(fRawMom.Pt(),fRawMom.Eta(),fRawMom.Phi(),fRawMom.M());
177 >  // compute combined correction scale from all enabled corrections
178 >  Double_t scale = 1.0;
179  
180    if (CorrectionActive(L2))
181 <    fCachedMom *= fL2RelativeCorrectionScale;
181 >    scale *= fL2RelativeCorrectionScale;
182      
183    if (CorrectionActive(L3))
184 <    fCachedMom *= fL3AbsoluteCorrectionScale;
184 >    scale *= fL3AbsoluteCorrectionScale;
185    
186    if (CorrectionActive(L4))
187 <    fCachedMom *= fL4EMFCorrectionScale;
187 >    scale *= fL4EMFCorrectionScale;
188      
189    if (CorrectionActive(L5))
190 <    fCachedMom *= fL5FlavorCorrectionScale;
190 >    scale *= fL5FlavorCorrectionScale;
191      
192    if (CorrectionActive(L7))
193 <    fCachedMom *= fL7PartonCorrectionScale;
193 >    scale *= fL7PartonCorrectionScale;
194      
195    if (CorrectionActive(Custom))
196 <    fCachedMom *= fCustomCorrectionScale;
196 >    scale *= fCustomCorrectionScale;
197 >
198 >  return scale;
199 > }
200 >
201 > //--------------------------------------------------------------------------------------------------
202 > inline void mithep::Jet::GetMom() const
203 > {
204 >  // Get raw momentum values from stored values and apply all enabled corrections.
205 >
206 >  fCachedMom.SetCoordinates(fRawMom.Pt(),fRawMom.Eta(),fRawMom.Phi(),fRawMom.M());
207 >
208 >  fCachedMom *= CombinedCorrectionScale();
209 >
210   }
211   #endif

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines