ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitAna/DataTree/interface/Jet.h

Comparing UserCode/MitAna/DataTree/interface/Jet.h (file contents):
Revision 1.17 by bendavid, Thu Feb 26 20:28:59 2009 UTC vs.
Revision 1.25 by sixie, Fri Mar 26 14:31:31 2010 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		{
#	Line 31 | Line 30 | namespace mithep
30		Custom
31		};
32
33	<	Jet() : fMaxEInEmTowers(0), fMaxEInHadTowers(0), fEnergyFractionH(0), fEnergyFractionEm(0),
34	<	fHadEnergyInHB(0), fHadEnergyInHO(0), fHadEnergyInHE(0), fHadEnergyInHF(0),
35	<	fEmEnergyInEB(0), fEmEnergyInEE(0), fEmEnergyInHF(0), fTowersArea(0), fN(0),
36	<	fN60(0), fN90(0), fAlpha(0), fBeta(0), fMatchedMCFlavor(0),
37	<	fJetProbabilityBJetTagsDisc(0), fJetBProbabilityBJetTagsDisc(0),
38	<	fSimpleSecondaryVertexBJetTagsDisc(0), fCombinedSecondaryVertexBJetTagsDisc(0),
39	<	fCombinedSecondaryVertexMVABJetTagsDisc(0), fImpactParameterMVABJetTagsDisc(0),
40	<	fTrackCountingHighEffBJetTagsDisc(0), fTrackCountingHighPurBJetTagsDisc(0),
41	<	fSoftMuonBJetTagsDisc(0), fSoftMuonNoIPBJetTagsDisc(0),
42	<	fSoftElectronBJetTagsDisc(0), fL2RelativeCorrectionScale(0),
43	<	fL3AbsoluteCorrectionScale(0), fL4EMFCorrectionScale(0),
44	<	fL5FlavorCorrectionScale(0), fL7PartonCorrectionScale(0)
45	<	{}
46	<	Jet(Double_t px, Double_t py, Double_t pz, Double_t e) : fRawMom(FourVector(px,py,pz,e)),
47	<	fMaxEInEmTowers(0), fMaxEInHadTowers(0), fEnergyFractionH(0), fEnergyFractionEm(0),
48	<	fHadEnergyInHB(0), fHadEnergyInHO(0), fHadEnergyInHE(0), fHadEnergyInHF(0),
49	<	fEmEnergyInEB(0), fEmEnergyInEE(0), fEmEnergyInHF(0), fTowersArea(0), fN(0),
50	<	fN60(0), fN90(0), fAlpha(0), fBeta(0), fMatchedMCFlavor(0),
51	<	fJetProbabilityBJetTagsDisc(0), fJetBProbabilityBJetTagsDisc(0),
52	<	fSimpleSecondaryVertexBJetTagsDisc(0), fCombinedSecondaryVertexBJetTagsDisc(0),
53	<	fCombinedSecondaryVertexMVABJetTagsDisc(0), fImpactParameterMVABJetTagsDisc(0),
54	<	fTrackCountingHighEffBJetTagsDisc(0), fTrackCountingHighPurBJetTagsDisc(0),
55	<	fSoftMuonBJetTagsDisc(0), fSoftMuonNoIPBJetTagsDisc(0),
56	<	fSoftElectronBJetTagsDisc(0), fL2RelativeCorrectionScale(0),
57	<	fL3AbsoluteCorrectionScale(0), fL4EMFCorrectionScale(0),
58	<	fL5FlavorCorrectionScale(0), fL7PartonCorrectionScale(0)
59	<	{}
60	<
61	<	Double_t   Alpha()                       const { return fAlpha;                     }
62	<	Double_t   Beta()                        const { return fBeta;                      }
63	<	Double_t   CombinedSecondaryVertexBJetTagsDisc()          const
64	<	{ return fCombinedSecondaryVertexBJetTagsDisc;    }
65	<	Double_t   CombinedSecondaryVertexMVABJetTagsDisc()       const
66	<	{ return fCombinedSecondaryVertexMVABJetTagsDisc; }
67	<	Double_t   CustomCorrectionScale()       const { return fCustomCorrectionScale;     }
68	<	void       DisableCorrection(ECorr c)          { fCorrections.ClearBit(c); ClearMom(); }
69	<	void       DisableCorrections()                { fCorrections.Clear(); ClearMom();  }
70	<	Double_t   EmEnergyInEB()                const { return fEmEnergyInEB;              }
71	<	Double_t   EmEnergyInEE()                const { return fEmEnergyInEE;              }
72	<	Double_t   EmEnergyInHF()                const { return fEmEnergyInHF;              }
73	<	void       EnableCorrection(ECorr c)           { fCorrections.SetBit(c); ClearMom();}
74	<	Bool_t     CorrectionActive(ECorr c)     const { return fCorrections.TestBit(c);  }
75	<	Double_t   EnergyFractionH()             const { return fEnergyFractionH;           }
76	<	Double_t   EnergyFractionEm()            const { return fEnergyFractionEm;          }
77	<	Double_t   HadEnergyInHO()               const { return fHadEnergyInHO;             }
78	<	Double_t   HadEnergyInHB()               const { return fHadEnergyInHB;             }
79	<	Double_t   HadEnergyInHF()               const { return fHadEnergyInHF;             }
80	<	Double_t   HadEnergyInHE()               const { return fHadEnergyInHE;             }
81	<	Bool_t     HasTower(const CaloTower *t)  const { return fTowers.HasObject(t);       }
82	<	Double_t   ImpactParameterMVABJetTagsDisc()               const
83	<	{ return fImpactParameterMVABJetTagsDisc;         }
84	<	Double_t   JetProbabilityBJetTagsDisc()                   const
85	<	{ return fJetProbabilityBJetTagsDisc;             }
86	<	Double_t   JetBProbabilityBJetTagsDisc()                  const
87	<	{ return fJetBProbabilityBJetTagsDisc;            }
88	<	Int_t      MatchedMCFlavor()             const { return fMatchedMCFlavor;           }
89	<	Double_t   MaxEInEmTowers()              const { return fMaxEInEmTowers;            }
90	<	Double_t   MaxEInHadTowers()             const { return fMaxEInHadTowers;           }
91	<	UShort_t   N()                           const { return fN;                         }
92	<	UShort_t   N60()                         const { return fN60;                       }
93	<	UShort_t   N90()                         const { return fN90;                       }
94	<	UInt_t     NTowers()                     const { return fTowers.GetEntries();       }
95	<	Double_t   L2RelativeCorrectionScale()   const { return fL2RelativeCorrectionScale; }
96	<	Double_t   L3AbsoluteCorrectionScale()   const { return fL3AbsoluteCorrectionScale; }
97	<	Double_t   L4EMFCorrectionScale()        const { return fL4EMFCorrectionScale;      }
98	<	Double_t   L5FlavorCorrectionScale()     const { return fL5FlavorCorrectionScale;   }
99	<	Double_t   L7PartonCorrectionScale()     const { return fL7PartonCorrectionScale;   }
100	<	EObjType   ObjType()                     const { return kJet;                       }
101	<	Double_t   SimpleSecondaryVertexBJetTagsDisc()            const
102	<	{ return fSimpleSecondaryVertexBJetTagsDisc; }
103	<	Double_t   SoftMuonBJetTagsDisc()                         const
104	<	{ return fSoftMuonBJetTagsDisc;              }
105	<	Double_t   SoftMuonNoIPBJetTagsDisc()                     const
106	<	{ return fSoftMuonNoIPBJetTagsDisc;          }
107	<	Double_t   SoftElectronBJetTagsDisc()                     const
108	<	{ return fSoftElectronBJetTagsDisc;          }
109	<	Double_t   TrackCountingHighEffBJetTagsDisc()             const
110	<	{ return fTrackCountingHighEffBJetTagsDisc;  }
111	<	Double_t   TrackCountingHighPurBJetTagsDisc()             const
112	<	{ return fTrackCountingHighPurBJetTagsDisc;  }
113	<	const CaloTower
114	<	*Tower(UInt_t i)               const { return fTowers.At(i);           }
115	<	Double_t   TowersArea()                  const { return fTowersArea;             }
116	<	void       AddTower(const CaloTower *tower)  { fTowers.Add(tower);               }
117	<	const FourVectorM &RawMom()              const { return fRawMom;                 }
118	<	void       SetAlpha(Double_t val)            { fAlpha            = val;          }
119	<	void       SetBeta(Double_t val)             { fBeta             = val;          }
120	<	void       SetCombinedSecondaryVertexBJetTagsDisc(Double_t d)
121	<	{ fCombinedSecondaryVertexBJetTagsDisc = d;    }
122	<	void       SetCombinedSecondaryVertexMVABJetTagsDisc(Double_t d)
123	<	{ fCombinedSecondaryVertexMVABJetTagsDisc = d; }
124	<	void       SetEmEnergyInEB(Double_t val)     { fEmEnergyInEB     = val;          }
125	<	void       SetEmEnergyInEE(Double_t val)     { fEmEnergyInEE     = val;          }
126	<	void       SetEmEnergyInHF(Double_t val)     { fEmEnergyInHF     = val;          }
127	<	void       SetEnergyFractionH(Double_t val)  { fEnergyFractionH  = val;          }
128	<	void       SetEnergyFractionEm(Double_t val) { fEnergyFractionEm = val;          }
129	<	void       SetHadEnergyInHO(Double_t val)    { fHadEnergyInHO    = val;          }
130	<	void       SetHadEnergyInHB(Double_t val)    { fHadEnergyInHB    = val;          }
131	<	void       SetHadEnergyInHF(Double_t val)    { fHadEnergyInHF    = val;          }
132	<	void       SetHadEnergyInHE(Double_t val)    { fHadEnergyInHE    = val;          }
133	<	void       SetImpactParameterMVABJetTagsDisc (Double_t d)
134	<	{ fImpactParameterMVABJetTagsDisc = d;        }
135	<	void       SetJetProbabilityBJetTagsDisc(Double_t d)  { fJetProbabilityBJetTagsDisc = d;   }
136	<	void       SetJetBProbabilityBJetTagsDisc(Double_t d) { fJetBProbabilityBJetTagsDisc = d;  }
137	<	void       SetRawMom(const FourVectorM &mom) { fRawMom = mom; ClearMom(); }
138	<	void       SetN(UShort_t n)                  { fN                = n;            }
139	<	void       SetN60(UShort_t n)                { fN60              = n;            }
141	<	void       SetN90(UShort_t n)                { fN90              = n;            }
142	<	void       SetMatchedMCFlavor(Int_t flavor)  { fMatchedMCFlavor = flavor;        }
143	<	void       SetMaxEInEmTowers(Double_t val)   { fMaxEInEmTowers   = val;          }
144	<	void       SetMaxEInHadTowers(Double_t val)  { fMaxEInHadTowers  = val;          }
145	<	void       SetL2RelativeCorrectionScale(Double_t s )   { fL2RelativeCorrectionScale = s; ClearMom(); }
146	<	void       SetL3AbsoluteCorrectionScale(Double_t s )   { fL3AbsoluteCorrectionScale = s; ClearMom(); }
147	<	void       SetL4EMFCorrectionScale(Double_t s )        { fL4EMFCorrectionScale = s; ClearMom();    }
148	<	void       SetL5FlavorCorrectionScale(Double_t s )     { fL5FlavorCorrectionScale = s; ClearMom(); }
149	<	void       SetL7PartonCorrectionScale(Double_t s )     { fL7PartonCorrectionScale = s; ClearMom(); }
150	<	void       SetCustomCorrectionScale(Double_t s)        { fCustomCorrectionScale = s; ClearMom(); }
151	<	void       SetSimpleSecondaryVertexBJetTagsDisc(Double_t d)
152	<	{ fSimpleSecondaryVertexBJetTagsDisc = d;     }
153	<	void       SetSoftMuonBJetTagsDisc(Double_t d)         { fSoftMuonBJetTagsDisc = d;        }
154	<	void       SetSoftMuonNoIPBJetTagsDisc(Double_t d)     { fSoftMuonNoIPBJetTagsDisc = d;    }
155	<	void       SetSoftElectronBJetTagsDisc(Double_t d)     { fSoftElectronBJetTagsDisc = d;    }
156	<	void       SetTowersArea(Double_t val)       { fTowersArea       = val;          }
157	<	void       SetTrackCountingHighEffBJetTagsDisc(Double_t d)
158	<	{ fTrackCountingHighEffBJetTagsDisc = d;      }
159	<	void       SetTrackCountingHighPurBJetTagsDisc(Double_t d)
160	<	{ fTrackCountingHighPurBJetTagsDisc = d;      }
33	>	Jet() :
34	>	fAlpha(0), fBeta(0), fMatchedMCFlavor(0),
35	>	fJetProbabilityBJetTagsDisc(0), fJetBProbabilityBJetTagsDisc(0),
36	>	fSimpleSecondaryVertexBJetTagsDisc(0), fCombinedSecondaryVertexBJetTagsDisc(0),
37	>	fCombinedSecondaryVertexMVABJetTagsDisc(0), fTrackCountingHighEffBJetTagsDisc(0),
38	>	fTrackCountingHighPurBJetTagsDisc(0), fSoftMuonBJetTagsDisc(0),
39	>	fSoftMuonByIP3dBJetTagsDisc(0), fSoftMuonByPtBJetTagsDisc(0),
40	>	fSoftElectronByIP3dBJetTagsDisc(0), fSoftElectronByPtBJetTagsDisc(0),
41	>	fL2RelativeCorrectionScale(0), fL3AbsoluteCorrectionScale(0),
42	>	fL4EMFCorrectionScale(0), 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), fTrackCountingHighEffBJetTagsDisc(0),
50	>	fTrackCountingHighPurBJetTagsDisc(0), fSoftMuonBJetTagsDisc(0),
51	>	fSoftMuonByIP3dBJetTagsDisc(0), fSoftMuonByPtBJetTagsDisc(0),
52	>	fSoftElectronByIP3dBJetTagsDisc(0), fSoftElectronByPtBJetTagsDisc(0),
53	>	fL2RelativeCorrectionScale(0), fL3AbsoluteCorrectionScale(0),
54	>	fL4EMFCorrectionScale(0), 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      CombinedCorrectionScale()     const;
64	>	Double_t      CustomCorrectionScale()       const { return fCustomCorrectionScale;        }
65	>	void          DisableCorrection(ECorr c)          { fCorrections.ClearBit(c); ClearMom(); }
66	>	void          DisableCorrections()                { fCorrections.Clear(); ClearMom();     }
67	>	void          EnableCorrection(ECorr c)           { fCorrections.SetBit(c); ClearMom();   }
68	>	Bool_t        CorrectionActive(ECorr c)     const { return fCorrections.TestBit(c);       }
69	>	Double_t      JetProbabilityBJetTagsDisc()                   const
70	>	{ return fJetProbabilityBJetTagsDisc;     }
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();            }
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      SoftMuonByIP3dBJetTagsDisc()                     const
89	>	{ return fSoftMuonByIP3dBJetTagsDisc;          }
90	>	Double_t      SoftMuonByPtBJetTagsDisc()                     const
91	>	{ return fSoftMuonByPtBJetTagsDisc;          }
92	>	Double_t      SoftElectronByIP3dBJetTagsDisc()                     const
93	>	{ return fSoftElectronByIP3dBJetTagsDisc;          }
94	>	Double_t      SoftElectronByPtBJetTagsDisc()                     const
95	>	{ return fSoftElectronByPtBJetTagsDisc;          }
96	>	Double_t      TrackCountingHighEffBJetTagsDisc()             const
97	>	{ return fTrackCountingHighEffBJetTagsDisc;  }
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)
110	>	{ fCombinedSecondaryVertexMVABJetTagsDisc = d; }
111	>	void          SetJetProbabilityBJetTagsDisc(Double_t d)
112	>	{ fJetProbabilityBJetTagsDisc = d;     }
113	>	void          SetJetBProbabilityBJetTagsDisc(Double_t d)
114	>	{ fJetBProbabilityBJetTagsDisc = d;    }
115	>	void          SetRawMom(const FourVectorM &mom) { fRawMom = mom; ClearMom();           }
116	>	void          SetMatchedMCFlavor(Int_t flavor)  { fMatchedMCFlavor = flavor;           }
117	>	void          SetL2RelativeCorrectionScale(Double_t s )
118	>	{ fL2RelativeCorrectionScale = s; ClearMom(); }
119	>	void          SetL3AbsoluteCorrectionScale(Double_t s )
120	>	{ fL3AbsoluteCorrectionScale = s; ClearMom(); }
121	>	void          SetL4EMFCorrectionScale(Double_t s )
122	>	{ fL4EMFCorrectionScale = s; ClearMom();      }
123	>	void          SetL5FlavorCorrectionScale(Double_t s )
124	>	{ fL5FlavorCorrectionScale = s; ClearMom();   }
125	>	void          SetL7PartonCorrectionScale(Double_t s )
126	>	{ fL7PartonCorrectionScale = s; ClearMom();   }
127	>	void          SetCustomCorrectionScale(Double_t s)
128	>	{ fCustomCorrectionScale = s; ClearMom(); }
129	>	void          SetSimpleSecondaryVertexBJetTagsDisc(Double_t d)
130	>	{ fSimpleSecondaryVertexBJetTagsDisc = d;     }
131	>	void          SetSoftMuonBJetTagsDisc(Double_t d)     { fSoftMuonBJetTagsDisc = d;     }
132	>	void          SetSoftMuonByIP3dBJetTagsDisc(Double_t d) { fSoftMuonByIP3dBJetTagsDisc = d; }
133	>	void          SetSoftMuonByPtBJetTagsDisc(Double_t d) { fSoftMuonByPtBJetTagsDisc = d; }
134	>	void          SetSoftElectronByIP3dBJetTagsDisc(Double_t d) { fSoftElectronByIP3dBJetTagsDisc = d; }
135	>	void          SetSoftElectronByPtBJetTagsDisc(Double_t d) { fSoftElectronByPtBJetTagsDisc = d; }
136	>	void          SetTrackCountingHighEffBJetTagsDisc(Double_t d)
137	>	{ fTrackCountingHighEffBJetTagsDisc = d;      }
138	>	void          SetTrackCountingHighPurBJetTagsDisc(Double_t d)
139	>	{ fTrackCountingHighPurBJetTagsDisc = d;      }
140
141		protected:
142	<	void       GetMom()                      const;
142	>	void          GetMom()                   const;
143
144	<	FourVectorM32 fRawMom;              //uncorrected four momentum of jet
145	<	Double_t   fMaxEInEmTowers;         //maximum energy in EM towers
146	<	Double_t   fMaxEInHadTowers;        //maximum energy in HCAL towers
147	<	Double_t   fEnergyFractionH;        //hadronic energy fraction
148	<	Double_t   fEnergyFractionEm;       //electromagnetic energy fraction
149	<	Double_t   fHadEnergyInHB;          //hadronic energy in HB
150	<	Double_t   fHadEnergyInHO;          //hadronic energy in HO
151	<	Double_t   fHadEnergyInHE;          //hadronic energy in HE
152	<	Double_t   fHadEnergyInHF;          //hadronic energy in HF
153	<	Double_t   fEmEnergyInEB;           //electromagnetic energy in EB
154	<	Double_t   fEmEnergyInEE;           //electromagnetic energy in EE
155	<	Double_t   fEmEnergyInHF;           //electromagnetic energy extracted from HF
156	<	Double_t   fTowersArea;             //area of contributing towers
157	<	UShort_t   fN;                      //number of constituents
158	<	UShort_t   fN60;                    //number constituents with 60% of total energy
159	<	UShort_t   fN90;                    //number constituents with 90% of total energy
160	<	Double_t   fAlpha;                  //jet vertex association alpha variable
161	<	Double_t   fBeta;                   //jet vertex association beta variable
162	<	Int_t      fMatchedMCFlavor;        //pdgID of the quark flavor that the jet matched to
163	<	Double_t   fJetProbabilityBJetTagsDisc;              //Discriminants for b-tagging algos
164	<	Double_t   fJetBProbabilityBJetTagsDisc;             //
165	<	Double_t   fSimpleSecondaryVertexBJetTagsDisc;       //
166	<	Double_t   fCombinedSecondaryVertexBJetTagsDisc;     //
167	<	Double_t   fCombinedSecondaryVertexMVABJetTagsDisc;  //
168	<	Double_t   fImpactParameterMVABJetTagsDisc;          //
190	<	Double_t   fTrackCountingHighEffBJetTagsDisc;        //
191	<	Double_t   fTrackCountingHighPurBJetTagsDisc;        //
192	<	Double_t   fSoftMuonBJetTagsDisc;                    //
193	<	Double_t   fSoftMuonNoIPBJetTagsDisc;                //
194	<	Double_t   fSoftElectronBJetTagsDisc;                //
195	<	Double_t   fL2RelativeCorrectionScale;           //
196	<	Double_t   fL3AbsoluteCorrectionScale;           //
197	<	Double_t   fL4EMFCorrectionScale;                //
198	<	Double_t   fL5FlavorCorrectionScale;             //
199	<	Double_t   fL7PartonCorrectionScale;             //
200	<	Double_t   fCustomCorrectionScale;               //
201	<	BitMask8   fCorrections;                         //mask of which corrections should be applied
202	<	RefArray<CaloTower>      fTowers;                //calo towers assigned to this jet
144	>	Vect4M        fRawMom;                                 //uncorrected four momentum of jet
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
153	>	Double32_t    fCombinedSecondaryVertexBJetTagsDisc;    //[0,0,14]discriminants b-tagging algos
154	>	Double32_t    fCombinedSecondaryVertexMVABJetTagsDisc; //[0,0,14]discriminants b-tagging algos
155	>	Double32_t    fTrackCountingHighEffBJetTagsDisc;       //[0,0,14]discriminants b-tagging algos
156	>	Double32_t    fTrackCountingHighPurBJetTagsDisc;       //[0,0,14]discriminants b-tagging algos
157	>	Double32_t    fSoftMuonBJetTagsDisc;                   //[0,0,14]discriminants b-tagging algos
158	>	Double32_t    fSoftMuonByIP3dBJetTagsDisc;             //[0,0,14]discriminants b-tagging algos
159	>	Double32_t    fSoftMuonByPtBJetTagsDisc;               //[0,0,14]discriminants b-tagging algos
160	>	Double32_t    fSoftElectronByIP3dBJetTagsDisc;         //[0,0,14]discriminants b-tagging algos
161	>	Double32_t    fSoftElectronByPtBJetTagsDisc;           //[0,0,14]discriminants b-tagging algos
162	>	Double32_t    fL2RelativeCorrectionScale;              //[0,0,14]L2 correction scale
163	>	Double32_t    fL3AbsoluteCorrectionScale;              //[0,0,14]L3 correction scale
164	>	Double32_t    fL4EMFCorrectionScale;                   //[0,0,14]L4 correction scale
165	>	Double32_t    fL5FlavorCorrectionScale;                //[0,0,14]L5 correction scale
166	>	Double32_t    fL7PartonCorrectionScale;                //[0,0,14]L7 correction scale
167	>	Double32_t    fCustomCorrectionScale;                  //[0,0,14]custom correction scale
168	>	BitMask8      fCorrections;                            //mask of corrections to be applied
169
170	<	ClassDef(Jet, 1) // 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 momentum values from stored values.
177	>	// compute combined correction scale from all enabled corrections
178	>	Double_t scale = 1.0;
179
213	–	fCachedMom = fRawMom;
214	–
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