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.13 by bendavid, Mon Jan 12 10:27:42 2009 UTC vs.
Revision 1.21 by loizides, Wed Mar 18 15:44:32 2009 UTC

# Line 3 | Line 3
3   //
4   // Jet
5   //
6 < // Jet class to hold reconstructed jet information.
6 > // Base jet class to hold reconstructed jet information.
7   //
8 < // Authors: S.Xie, C.Loizides
8 > // Authors: S.Xie, C.Loizides, J.Bendavid
9   //--------------------------------------------------------------------------------------------------
10  
11   #ifndef MITANA_DATATREE_JET_H
12   #define MITANA_DATATREE_JET_H
13  
14 + #include "MitCommon/DataFormats/interface/Vect4M.h"
15   #include "MitAna/DataTree/interface/Particle.h"
15 #include "MitAna/DataTree/interface/Track.h"
16 #include "MitAna/DataTree/interface/CaloTower.h"
16  
17   namespace mithep
18   {
19    class Jet : public Particle
20    {
21      public:
22 <      Jet() : fMaxEInEmTowers(0), fMaxEInHadTowers(0), fEnergyFractionH(0), fEnergyFractionEm(0),
23 <              fHadEnergyInHB(0), fHadEnergyInHO(0), fHadEnergyInHE(0), fHadEnergyInHF(0),
24 <              fEmEnergyInEB(0), fEmEnergyInEE(0), fEmEnergyInHF(0), fTowersArea(0), fN(0),
25 <              fN60(0), fN90(0), fMatchedMCFlavor(0) {}
26 <      Jet(Double_t px, Double_t py, Double_t pz, Double_t e) : fMom(px,py,pz,e),
27 <              fMaxEInEmTowers(0), fMaxEInHadTowers(0), fEnergyFractionH(0), fEnergyFractionEm(0),
28 <              fHadEnergyInHB(0), fHadEnergyInHO(0), fHadEnergyInHE(0), fHadEnergyInHF(0),
29 <              fEmEnergyInEB(0), fEmEnergyInEE(0), fEmEnergyInHF(0), fTowersArea(0), fN(0),
30 <              fN60(0), fN90(0), fMatchedMCFlavor(0) {}
31 <      ~Jet() {}
32 <
33 <      Double_t   Alpha()                       const { return fAlpha;                     }
34 <      Double_t   Beta()                        const { return fBeta;                      }
35 <      Double_t   Charge()                      const { return 0;                          }
36 <      Double_t   CombinedSecondaryVertexBJetTagsDisc()          const
37 <                   { return fCombinedSecondaryVertexBJetTagsDisc;    }
38 <      Double_t   CombinedSecondaryVertexMVABJetTagsDisc()       const
39 <                   { return fCombinedSecondaryVertexMVABJetTagsDisc; }
40 <      Double_t   EmEnergyInEB()                const { return fEmEnergyInEB;              }
41 <      Double_t   EmEnergyInEE()                const { return fEmEnergyInEE;              }
42 <      Double_t   EmEnergyInHF()                const { return fEmEnergyInHF;              }
43 <      Double_t   EnergyFractionH()             const { return fEnergyFractionH;           }
44 <      Double_t   EnergyFractionEm()            const { return fEnergyFractionEm;          }
45 <      Double_t   HadEnergyInHO()               const { return fHadEnergyInHO;             }
46 <      Double_t   HadEnergyInHB()               const { return fHadEnergyInHB;             }
47 <      Double_t   HadEnergyInHF()               const { return fHadEnergyInHF;             }
48 <      Double_t   HadEnergyInHE()               const { return fHadEnergyInHE;             }
49 <      Double_t   ImpactParameterMVABJetTagsDisc()               const
50 <                   { return fImpactParameterMVABJetTagsDisc;         }
51 <      Double_t   JetProbabilityBJetTagsDisc()                   const
52 <                   { return fJetProbabilityBJetTagsDisc;             }
53 <      Double_t   JetBProbabilityBJetTagsDisc()                  const
54 <                   { return fJetBProbabilityBJetTagsDisc;            }
55 <      Int_t      MatchedMCFlavor()             const { return fMatchedMCFlavor;           }
56 <      Double_t   MaxEInEmTowers()              const { return fMaxEInEmTowers;            }
57 <      Double_t   MaxEInHadTowers()             const { return fMaxEInHadTowers;           }
58 <      FourVector Mom()                         const { return fMom;                       }
59 <      UShort_t   N()                           const { return fN;                         }
60 <      UShort_t   N60()                         const { return fN60;                       }
61 <      UShort_t   N90()                         const { return fN90;                       }
62 <      UInt_t     NTowers()                     const { return fTowers.GetEntries();       }
63 <      Double_t   L2RelativeCorrectionScale()   const { return fL2RelativeCorrectionScale; }
64 <      Double_t   L3AbsoluteCorrectionScale()   const { return fL3AbsoluteCorrectionScale; }
65 <      Double_t   L4EMFCorrectionScale()        const { return fL4EMFCorrectionScale;      }
66 <      Double_t   L5FlavorCorrectionScale()     const { return fL5FlavorCorrectionScale;   }
67 <      Double_t   L7PartonCorrectionScale()     const { return fL7PartonCorrectionScale;   }
68 <      EObjType    ObjType()                    const { return kJet;                       }      
69 <      Double_t   SimpleSecondaryVertexBJetTagsDisc()            const
70 <                   { return fSimpleSecondaryVertexBJetTagsDisc;      }
71 <      Double_t   SoftMuonBJetTagsDisc()                         const
72 <                   { return fSoftMuonBJetTagsDisc;                   }
73 <      Double_t   SoftMuonNoIPBJetTagsDisc()                     const
74 <                   { return fSoftMuonNoIPBJetTagsDisc;               }
75 <      Double_t   SoftElectronBJetTagsDisc()                     const
76 <                   { return fSoftElectronBJetTagsDisc;               }
77 <      Double_t   TrackCountingHighEffBJetTagsDisc()             const
78 <                   { return fTrackCountingHighEffBJetTagsDisc;       }
79 <      Double_t   TrackCountingHighPurBJetTagsDisc()             const
80 <                   { return fTrackCountingHighPurBJetTagsDisc;       }
81 <      const CaloTower  
82 <                *Tower(UInt_t i)               const { return fTowers.At(i);           }
83 <      Double_t   TowersArea()                  const { return fTowersArea;             }
84 <      void       AddTower(const CaloTower *tower)  { fTowers.Add(tower);               }
85 <      void       SetAlpha(Double_t val)            { fAlpha            = val;          }
86 <      void       SetBeta(Double_t val)             { fBeta             = val;          }
87 <      void       SetCombinedSecondaryVertexBJetTagsDisc(Double_t d)
88 <                                                   { fCombinedSecondaryVertexBJetTagsDisc = d;     }
89 <      void       SetCombinedSecondaryVertexMVABJetTagsDisc(Double_t d)
90 <                                                   { fCombinedSecondaryVertexMVABJetTagsDisc = d;  }
91 <      void       SetEmEnergyInEB(Double_t val)     { fEmEnergyInEB     = val;          }
92 <      void       SetEmEnergyInEE(Double_t val)     { fEmEnergyInEE     = val;          }
93 <      void       SetEmEnergyInHF(Double_t val)     { fEmEnergyInHF     = val;          }
94 <      void       SetEnergyFractionH(Double_t val)  { fEnergyFractionH  = val;          }
95 <      void       SetEnergyFractionEm(Double_t val) { fEnergyFractionEm = val;          }
96 <      void       SetHadEnergyInHO(Double_t val)    { fHadEnergyInHO    = val;          }
97 <      void       SetHadEnergyInHB(Double_t val)    { fHadEnergyInHB    = val;          }
98 <      void       SetHadEnergyInHF(Double_t val)    { fHadEnergyInHF    = val;          }
99 <      void       SetHadEnergyInHE(Double_t val)    { fHadEnergyInHE    = val;          }
100 <      void       SetImpactParameterMVABJetTagsDisc (Double_t d)
101 <                                                   {   fImpactParameterMVABJetTagsDisc = d;        }
102 <      void       SetJetProbabilityBJetTagsDisc(Double_t d)  { fJetProbabilityBJetTagsDisc = d;     }
103 <      void       SetJetBProbabilityBJetTagsDisc(Double_t d) { fJetBProbabilityBJetTagsDisc = d;    }
104 <      void       SetMom(const FourVector &mom)     { fMom              = mom;          }
105 <      void       SetN(UShort_t n)                  { fN                = n;            }
106 <      void       SetN60(UShort_t n)                { fN60              = n;            }
107 <      void       SetN90(UShort_t n)                { fN90              = n;            }
108 <      void       SetMatchedMCFlavor(Int_t flavor)  { fMatchedMCFlavor = flavor;        }
109 <      void       SetMaxEInEmTowers(Double_t val)   { fMaxEInEmTowers   = val;          }
110 <      void       SetMaxEInHadTowers(Double_t val)  { fMaxEInHadTowers  = val;          }
111 <      void       SetL2RelativeCorrectionScale(Double_t s )   { fL2RelativeCorrectionScale = s;   }
112 <      void       SetL3AbsoluteCorrectionScale(Double_t s )   { fL3AbsoluteCorrectionScale = s;   }
113 <      void       SetL4EMFCorrectionScale(Double_t s )        { fL4EMFCorrectionScale = s;        }
114 <      void       SetL5FlavorCorrectionScale(Double_t s )     { fL5FlavorCorrectionScale = s;     }
115 <      void       SetL7PartonCorrectionScale(Double_t s )     { fL7PartonCorrectionScale = s;     }
116 <      void       SetSimpleSecondaryVertexBJetTagsDisc(Double_t d)
117 <                                                   { fSimpleSecondaryVertexBJetTagsDisc = d;     }
118 <      void       SetSoftMuonBJetTagsDisc(Double_t d )        { fSoftMuonBJetTagsDisc = d;        }
119 <      void       SetSoftMuonNoIPBJetTagsDisc(Double_t d )    { fSoftMuonNoIPBJetTagsDisc = d;    }
120 <      void       SetSoftElectronBJetTagsDisc(Double_t d )    { fSoftElectronBJetTagsDisc = d;    }
121 <      void       SetTowersArea(Double_t val)       { fTowersArea       = val;          }
122 <      void       SetTrackCountingHighEffBJetTagsDisc( Double_t d)
123 <                                                   { fTrackCountingHighEffBJetTagsDisc = d;      }
124 <      void       SetTrackCountingHighPurBJetTagsDisc(Double_t d)
125 <                                                   { fTrackCountingHighPurBJetTagsDisc = d;      }
22 >      enum ECorr {
23 >        L1 = 0,
24 >        L2,
25 >        L3,
26 >        L4,
27 >        L5,
28 >        L6,
29 >        L7,
30 >        Custom
31 >      };
32 >    
33 >      Jet() :
34 >        fAlpha(0), fBeta(0), fMatchedMCFlavor(0),
35 >        fJetProbabilityBJetTagsDisc(0), fJetBProbabilityBJetTagsDisc(0),
36 >        fSimpleSecondaryVertexBJetTagsDisc(0), fCombinedSecondaryVertexBJetTagsDisc(0),
37 >        fCombinedSecondaryVertexMVABJetTagsDisc(0), fImpactParameterMVABJetTagsDisc(0),
38 >        fTrackCountingHighEffBJetTagsDisc(0), fTrackCountingHighPurBJetTagsDisc(0),
39 >        fSoftMuonBJetTagsDisc(0), fSoftMuonNoIPBJetTagsDisc(0),
40 >        fSoftElectronBJetTagsDisc(0), fL2RelativeCorrectionScale(0),
41 >        fL3AbsoluteCorrectionScale(0), fL4EMFCorrectionScale(0),
42 >        fL5FlavorCorrectionScale(0), fL7PartonCorrectionScale(0),
43 >        fCustomCorrectionScale(0) {}
44 >      Jet(Double_t px, Double_t py, Double_t pz, Double_t e) :
45 >        fRawMom(FourVector(px,py,pz,e)),
46 >        fAlpha(0), fBeta(0), fMatchedMCFlavor(0),
47 >        fJetProbabilityBJetTagsDisc(0), fJetBProbabilityBJetTagsDisc(0),
48 >        fSimpleSecondaryVertexBJetTagsDisc(0), fCombinedSecondaryVertexBJetTagsDisc(0),
49 >        fCombinedSecondaryVertexMVABJetTagsDisc(0), fImpactParameterMVABJetTagsDisc(0),
50 >        fTrackCountingHighEffBJetTagsDisc(0), fTrackCountingHighPurBJetTagsDisc(0),
51 >        fSoftMuonBJetTagsDisc(0), fSoftMuonNoIPBJetTagsDisc(0),
52 >        fSoftElectronBJetTagsDisc(0), fL2RelativeCorrectionScale(0),
53 >        fL3AbsoluteCorrectionScale(0), fL4EMFCorrectionScale(0),
54 >        fL5FlavorCorrectionScale(0), fL7PartonCorrectionScale(0),
55 >        fCustomCorrectionScale(0) {}
56 >
57 >      Double_t      Alpha()                       const { return fAlpha;                        }
58 >      Double_t      Beta()                        const { return fBeta;                         }
59 >      Double_t      CombinedSecondaryVertexBJetTagsDisc()          const
60 >                      { return fCombinedSecondaryVertexBJetTagsDisc;    }
61 >      Double_t      CombinedSecondaryVertexMVABJetTagsDisc()       const
62 >                      { return fCombinedSecondaryVertexMVABJetTagsDisc; }
63 >      Double_t      CustomCorrectionScale()       const { return fCustomCorrectionScale;        }
64 >      void          DisableCorrection(ECorr c)          { fCorrections.ClearBit(c); ClearMom(); }
65 >      void          DisableCorrections()                { fCorrections.Clear(); ClearMom();     }
66 >      void          EnableCorrection(ECorr c)           { fCorrections.SetBit(c); ClearMom();   }
67 >      Bool_t        CorrectionActive(ECorr c)     const { return fCorrections.TestBit(c);       }
68 >      Double_t      ImpactParameterMVABJetTagsDisc()               const
69 >                      { return fImpactParameterMVABJetTagsDisc; }
70 >      Double_t      JetProbabilityBJetTagsDisc()                   const
71 >                      { return fJetProbabilityBJetTagsDisc;     }
72 >      Double_t      JetBProbabilityBJetTagsDisc()                  const
73 >                      { return fJetBProbabilityBJetTagsDisc;    }
74 >      Int_t         MatchedMCFlavor()             const { return fMatchedMCFlavor;           }
75 >      virtual
76 >      UInt_t        NConstituents()               const { return 0;                          }
77 >      UInt_t        N()                           const { return NConstituents();            }
78 >      Double_t      L2RelativeCorrectionScale()   const { return fL2RelativeCorrectionScale; }
79 >      Double_t      L3AbsoluteCorrectionScale()   const { return fL3AbsoluteCorrectionScale; }
80 >      Double_t      L4EMFCorrectionScale()        const { return fL4EMFCorrectionScale;      }
81 >      Double_t      L5FlavorCorrectionScale()     const { return fL5FlavorCorrectionScale;   }
82 >      Double_t      L7PartonCorrectionScale()     const { return fL7PartonCorrectionScale;   }
83 >      EObjType      ObjType()                     const { return kJet;                       }      
84 >      Double_t      SimpleSecondaryVertexBJetTagsDisc()            const
85 >                      { return fSimpleSecondaryVertexBJetTagsDisc; }
86 >      Double_t      SoftMuonBJetTagsDisc()                         const
87 >                      { return fSoftMuonBJetTagsDisc;              }
88 >      Double_t      SoftMuonNoIPBJetTagsDisc()                     const
89 >                      { return fSoftMuonNoIPBJetTagsDisc;          }
90 >      Double_t      SoftElectronBJetTagsDisc()                     const
91 >                      { return fSoftElectronBJetTagsDisc;          }
92 >      Double_t      TrackCountingHighEffBJetTagsDisc()             const
93 >                      { return fTrackCountingHighEffBJetTagsDisc;  }
94 >      Double_t      TrackCountingHighPurBJetTagsDisc()             const
95 >                      { return fTrackCountingHighPurBJetTagsDisc;  }
96 >      FourVectorM   RawMom()                      const { return fRawMom.V();                }
97 >      void          SetAlpha(Double_t val)              { fAlpha = val;                      }
98 >      void          SetBeta(Double_t val)               { fBeta  = val;                      }
99 >      void          SetCombinedSecondaryVertexBJetTagsDisc(Double_t d)
100 >                      { fCombinedSecondaryVertexBJetTagsDisc = d;    }
101 >      void          SetCombinedSecondaryVertexMVABJetTagsDisc(Double_t d)
102 >                      { fCombinedSecondaryVertexMVABJetTagsDisc = d; }
103 >      void          SetImpactParameterMVABJetTagsDisc (Double_t d)
104 >                      { fImpactParameterMVABJetTagsDisc = d; }
105 >      void          SetJetProbabilityBJetTagsDisc(Double_t d)  
106 >                      { fJetProbabilityBJetTagsDisc = d;     }
107 >      void          SetJetBProbabilityBJetTagsDisc(Double_t d)
108 >                      { fJetBProbabilityBJetTagsDisc = d;    }
109 >      void          SetRawMom(const FourVectorM &mom) { fRawMom = mom; ClearMom();           }
110 >      void          SetMatchedMCFlavor(Int_t flavor)  { fMatchedMCFlavor = flavor;           }
111 >      void          SetL2RelativeCorrectionScale(Double_t s )  
112 >                      { fL2RelativeCorrectionScale = s; ClearMom(); }
113 >      void          SetL3AbsoluteCorrectionScale(Double_t s )  
114 >                      { fL3AbsoluteCorrectionScale = s; ClearMom(); }
115 >      void          SetL4EMFCorrectionScale(Double_t s )        
116 >                      { fL4EMFCorrectionScale = s; ClearMom();      }
117 >      void          SetL5FlavorCorrectionScale(Double_t s )    
118 >                      { fL5FlavorCorrectionScale = s; ClearMom();   }
119 >      void          SetL7PartonCorrectionScale(Double_t s )    
120 >                      { fL7PartonCorrectionScale = s; ClearMom();   }
121 >      void          SetCustomCorrectionScale(Double_t s)        
122 >                      { fCustomCorrectionScale = s; ClearMom(); }
123 >      void          SetSimpleSecondaryVertexBJetTagsDisc(Double_t d)
124 >                      { fSimpleSecondaryVertexBJetTagsDisc = d;     }
125 >      void          SetSoftMuonBJetTagsDisc(Double_t d)     { fSoftMuonBJetTagsDisc = d;     }
126 >      void          SetSoftMuonNoIPBJetTagsDisc(Double_t d) { fSoftMuonNoIPBJetTagsDisc = d; }
127 >      void          SetSoftElectronBJetTagsDisc(Double_t d) { fSoftElectronBJetTagsDisc = d; }
128 >      void          SetTrackCountingHighEffBJetTagsDisc(Double_t d)
129 >                      { fTrackCountingHighEffBJetTagsDisc = d;      }
130 >      void          SetTrackCountingHighPurBJetTagsDisc(Double_t d)
131 >                      { fTrackCountingHighPurBJetTagsDisc = d;      }
132  
133      protected:
134 <      FourVector fMom;                    //four momentum of jet
135 <      Double_t   fMaxEInEmTowers;         //maximum energy in EM towers
136 <      Double_t   fMaxEInHadTowers;        //maximum energy in HCAL towers
137 <      Double_t   fEnergyFractionH;        //hadronic energy fraction
138 <      Double_t   fEnergyFractionEm;       //electromagnetic energy fraction
139 <      Double_t   fHadEnergyInHB;          //hadronic energy in HB
140 <      Double_t   fHadEnergyInHO;          //hadronic energy in HO
141 <      Double_t   fHadEnergyInHE;          //hadronic energy in HE
142 <      Double_t   fHadEnergyInHF;          //hadronic energy in HF
143 <      Double_t   fEmEnergyInEB;           //electromagnetic energy in EB
144 <      Double_t   fEmEnergyInEE;           //electromagnetic energy in EE
145 <      Double_t   fEmEnergyInHF;           //electromagnetic energy extracted from HF
146 <      Double_t   fTowersArea;             //area of contributing towers
147 <      UShort_t   fN;                      //number of constituents
148 <      UShort_t   fN60;                    //number constituents with 60% of total energy
149 <      UShort_t   fN90;                    //number constituents with 90% of total energy
150 <      Double_t   fAlpha;                  //jet vertex association alpha variable
151 <      Double_t   fBeta;                   //jet vertex association beta variable
152 <      Int_t      fMatchedMCFlavor;        //pdgID of the quark flavor that the jet matched to
153 <      Double_t   fJetProbabilityBJetTagsDisc;              //Discriminants for b-tagging algos
154 <      Double_t   fJetBProbabilityBJetTagsDisc;             //
155 <      Double_t   fSimpleSecondaryVertexBJetTagsDisc;       //
156 <      Double_t   fCombinedSecondaryVertexBJetTagsDisc;     //
157 <      Double_t   fCombinedSecondaryVertexMVABJetTagsDisc;  //
153 <      Double_t   fImpactParameterMVABJetTagsDisc;          //
154 <      Double_t   fTrackCountingHighEffBJetTagsDisc;        //
155 <      Double_t   fTrackCountingHighPurBJetTagsDisc;        //
156 <      Double_t   fSoftMuonBJetTagsDisc;                    //
157 <      Double_t   fSoftMuonNoIPBJetTagsDisc;                //
158 <      Double_t   fSoftElectronBJetTagsDisc;                //
159 <      Double_t   fL2RelativeCorrectionScale;               //
160 <      Double_t   fL3AbsoluteCorrectionScale;               //
161 <      Double_t   fL4EMFCorrectionScale;                    //
162 <      Double_t   fL5FlavorCorrectionScale;                 //
163 <      Double_t   fL7PartonCorrectionScale;                 //
164 <      RefArray<CaloTower,1024> fTowers;                    //calo towers assigned to this jet
134 >      void          GetMom()                   const;
135 >
136 >      Vect4M        fRawMom;                                 //uncorrected four momentum of jet
137 >      Double32_t    fAlpha;                                  //[0,0,14]jet vertex alpha variable
138 >      Double32_t    fBeta;                                   //[0,0,14]jet vertex beta variable
139 >      Int_t         fMatchedMCFlavor;                        //[0,0,14]pdg of matched quark flavor
140 >      Double32_t    fJetProbabilityBJetTagsDisc;             //[0,0,14]discriminants b-tagging algos
141 >      Double32_t    fJetBProbabilityBJetTagsDisc;            //[0,0,14]discriminants b-tagging algos
142 >      Double32_t    fSimpleSecondaryVertexBJetTagsDisc;      //[0,0,14]discriminants b-tagging algos
143 >      Double32_t    fCombinedSecondaryVertexBJetTagsDisc;    //[0,0,14]discriminants b-tagging algos
144 >      Double32_t    fCombinedSecondaryVertexMVABJetTagsDisc; //[0,0,14]discriminants b-tagging algos
145 >      Double32_t    fImpactParameterMVABJetTagsDisc;         //[0,0,14]discriminants b-tagging algos
146 >      Double32_t    fTrackCountingHighEffBJetTagsDisc;       //[0,0,14]discriminants b-tagging algos
147 >      Double32_t    fTrackCountingHighPurBJetTagsDisc;       //[0,0,14]discriminants b-tagging algos
148 >      Double32_t    fSoftMuonBJetTagsDisc;                   //[0,0,14]discriminants b-tagging algos
149 >      Double32_t    fSoftMuonNoIPBJetTagsDisc;               //[0,0,14]discriminants b-tagging algos
150 >      Double32_t    fSoftElectronBJetTagsDisc;               //[0,0,14]discriminants b-tagging algos
151 >      Double32_t    fL2RelativeCorrectionScale;              //[0,0,14]L2 correction scale
152 >      Double32_t    fL3AbsoluteCorrectionScale;              //[0,0,14]L3 correction scale
153 >      Double32_t    fL4EMFCorrectionScale;                   //[0,0,14]L4 correction scale
154 >      Double32_t    fL5FlavorCorrectionScale;                //[0,0,14]L5 correction scale
155 >      Double32_t    fL7PartonCorrectionScale;                //[0,0,14]L7 correction scale
156 >      Double32_t    fCustomCorrectionScale;                  //[0,0,14]custom correction scale
157 >      BitMask8      fCorrections;                            //mask of corrections to be applied
158  
159      ClassDef(Jet, 1) // Jet class
160    };
161   }
162 +
163 + //--------------------------------------------------------------------------------------------------
164 + inline void mithep::Jet::GetMom() const
165 + {
166 +  // Get raw momentum values from stored values and apply all enabled corrections.
167 +
168 +  fCachedMom.SetCoordinates(fRawMom.Pt(),fRawMom.Eta(),fRawMom.Phi(),fRawMom.M());
169 +
170 +  if (CorrectionActive(L2))
171 +    fCachedMom *= fL2RelativeCorrectionScale;
172 +    
173 +  if (CorrectionActive(L3))
174 +    fCachedMom *= fL3AbsoluteCorrectionScale;
175 +  
176 +  if (CorrectionActive(L4))
177 +    fCachedMom *= fL4EMFCorrectionScale;
178 +    
179 +  if (CorrectionActive(L5))
180 +    fCachedMom *= fL5FlavorCorrectionScale;
181 +    
182 +  if (CorrectionActive(L7))
183 +    fCachedMom *= fL7PartonCorrectionScale;
184 +    
185 +  if (CorrectionActive(Custom))
186 +    fCachedMom *= fCustomCorrectionScale;
187 + }
188   #endif

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines