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.8 by loizides, Wed Sep 10 03:33:26 2008 UTC vs.
Revision 1.14 by bendavid, Tue Feb 17 14:23:41 2009 UTC

#	Line 3 | Line 3
3		//
4		// Jet
5		//
6	<	// Details to be worked out...
6	>	// Jet class to hold reconstructed jet information.
7		//
8	<	// Authors: C.Loizides
8	>	// Authors: S.Xie, C.Loizides
9		//--------------------------------------------------------------------------------------------------
10
11		#ifndef MITANA_DATATREE_JET_H
#	Line 31 | Line 31 | namespace mithep
31		fN60(0), fN90(0), fMatchedMCFlavor(0) {}
32		~Jet() {}
33
34	<	FourVector            Mom()                             const { return fMom; }
35	<	Double_t              Charge()                          const { return 0; }
36	<	Double_t              MaxEInEmTowers()                  const { return fMaxEInEmTowers; }
37	<	Double_t              MaxEInHadTowers()                 const { return fMaxEInHadTowers; }
38	<	Double_t              HadEnergyInHO()                   const { return fHadEnergyInHO; }
39	<	Double_t              HadEnergyInHB()                   const { return fHadEnergyInHB; }
40	<	Double_t              HadEnergyInHF()                   const { return fHadEnergyInHF; }
41	<	Double_t              HadEnergyInHE()                   const { return fHadEnergyInHE; }
42	<	Double_t              EmEnergyInEB()                    const { return fEmEnergyInEB; }
43	<	Double_t              EmEnergyInEE()                    const { return fEmEnergyInEE; }
44	<	Double_t              EmEnergyInHF()                    const { return fEmEnergyInHF; }
45	<	Double_t              EnergyFractionH()                 const { return fEnergyFractionH; }
46	<	Double_t              EnergyFractionEm()                const { return fEnergyFractionEm; }
47	<	const CaloTower      *Tower(UInt_t i)                   const { return fTowers.At(i); }
48	<	Double_t              TowersArea()                      const { return fTowersArea; }
49	<	UShort_t              N()                               const { return fN; }
50	<	UShort_t              N60()                             const { return fN60; }
51	<	UShort_t              N90()                             const { return fN90; }
52	<	UInt_t                NTowers()                         const { return fTowers.GetEntries(); }
53	<	Int_t                 MatchedMCFlavor()                 const { return fMatchedMCFlavor; }
54	<	Double_t JetProbabilityBJetTagsDisc()                   const
55	<	{ return fJetProbabilityBJetTagsDisc;             }
56	<	Double_t JetBProbabilityBJetTagsDisc()                  const
57	<	{ return fJetBProbabilityBJetTagsDisc;            }
58	<	Double_t SimpleSecondaryVertexBJetTagsDisc()            const
59	<	{ return fSimpleSecondaryVertexBJetTagsDisc;      }
60	<	Double_t CombinedSecondaryVertexBJetTagsDisc()          const
61	<	{ return fCombinedSecondaryVertexBJetTagsDisc;    }
62	<	Double_t CombinedSecondaryVertexMVABJetTagsDisc()       const
63	<	{ return fCombinedSecondaryVertexMVABJetTagsDisc; }
64	<	Double_t ImpactParameterMVABJetTagsDisc()               const
65	<	{ return fImpactParameterMVABJetTagsDisc;         }
66	<	Double_t TrackCountingHighEffBJetTagsDisc()             const
67	<	{ return fTrackCountingHighEffBJetTagsDisc;       }
68	<	Double_t TrackCountingHighPurBJetTagsDisc()             const
69	<	{ return fTrackCountingHighPurBJetTagsDisc;       }
70	<	Double_t SoftMuonBJetTagsDisc()                         const
71	<	{ return fSoftMuonBJetTagsDisc;                   }
72	<	Double_t SoftMuonNoIPBJetTagsDisc()                     const
73	<	{ return fSoftMuonNoIPBJetTagsDisc;               }
74	<	Double_t SoftElectronBJetTagsDisc()                     const
75	<	{ return fSoftElectronBJetTagsDisc;               }
76	<	void                  AddTower(CaloTower *tower)         { fTowers.Add(tower);      }
77	<	void                  SetMaxEInEmTowers(Double_t val)    { fMaxEInEmTowers   = val;   }
78	<	void                  SetMaxEInHadTowers(Double_t val)   { fMaxEInHadTowers  = val;   }
79	<	void                  SetHadEnergyInHO(Double_t val)     { fHadEnergyInHO    = val;   }
80	<	void                  SetHadEnergyInHB(Double_t val)     { fHadEnergyInHB    = val;   }
81	<	void                  SetHadEnergyInHF(Double_t val)     { fHadEnergyInHF    = val;   }
82	<	void                  SetHadEnergyInHE(Double_t val)     { fHadEnergyInHE    = val;   }
83	<	void                  SetEmEnergyInEB(Double_t val)      { fEmEnergyInEB     = val;   }
84	<	void                  SetEmEnergyInEE(Double_t val)      { fEmEnergyInEE     = val;   }
85	<	void                  SetEmEnergyInHF(Double_t val)      { fEmEnergyInHF     = val;   }
86	<	void                  SetEnergyFractionH(Double_t val)   { fEnergyFractionH  = val;   }
87	<	void                  SetEnergyFractionEm(Double_t val)  { fEnergyFractionEm = val;   }
88	<	void                  SetTowersArea(Double_t val)        { fTowersArea       = val;   }
89	<	void                  SetN(UShort_t n)                   { fN = n;                    }
90	<	void                  SetN60(UShort_t n)                 { fN60 = n;                  }
91	<	void                  SetN90(UShort_t n)                 { fN90 = n;                  }
92	<	void                  SetMatchedMCFlavor(Int_t flavor)   { fMatchedMCFlavor = flavor; }
93	<	void       SetJetProbabilityBJetTagsDisc ( Double_t d )  { fJetProbabilityBJetTagsDisc = d;  }
94	<	void       SetJetBProbabilityBJetTagsDisc ( Double_t d ) { fJetBProbabilityBJetTagsDisc = d; }
95	<	void       SetSimpleSecondaryVertexBJetTagsDisc ( Double_t d )
96	<	{ fSimpleSecondaryVertexBJetTagsDisc = d;       }
97	<	void       SetCombinedSecondaryVertexBJetTagsDisc ( Double_t d )
34	>	Double_t   Alpha()                       const { return fAlpha;                     }
35	>	Double_t   Beta()                        const { return fBeta;                      }
36	>	Double_t   Charge()                      const { return 0;                          }
37	>	Double_t   CombinedSecondaryVertexBJetTagsDisc()          const
38	>	{ return fCombinedSecondaryVertexBJetTagsDisc;    }
39	>	Double_t   CombinedSecondaryVertexMVABJetTagsDisc()       const
40	>	{ return fCombinedSecondaryVertexMVABJetTagsDisc; }
41	>	Double_t   EmEnergyInEB()                const { return fEmEnergyInEB;              }
42	>	Double_t   EmEnergyInEE()                const { return fEmEnergyInEE;              }
43	>	Double_t   EmEnergyInHF()                const { return fEmEnergyInHF;              }
44	>	Double_t   EnergyFractionH()             const { return fEnergyFractionH;           }
45	>	Double_t   EnergyFractionEm()            const { return fEnergyFractionEm;          }
46	>	Double_t   HadEnergyInHO()               const { return fHadEnergyInHO;             }
47	>	Double_t   HadEnergyInHB()               const { return fHadEnergyInHB;             }
48	>	Double_t   HadEnergyInHF()               const { return fHadEnergyInHF;             }
49	>	Double_t   HadEnergyInHE()               const { return fHadEnergyInHE;             }
50	>	Bool_t     HasTower(const CaloTower *t)  const { return fTowers.HasObject(t);       }
51	>	Double_t   ImpactParameterMVABJetTagsDisc()               const
52	>	{ return fImpactParameterMVABJetTagsDisc;         }
53	>	Double_t   JetProbabilityBJetTagsDisc()                   const
54	>	{ return fJetProbabilityBJetTagsDisc;             }
55	>	Double_t   JetBProbabilityBJetTagsDisc()                  const
56	>	{ return fJetBProbabilityBJetTagsDisc;            }
57	>	Int_t      MatchedMCFlavor()             const { return fMatchedMCFlavor;           }
58	>	Double_t   MaxEInEmTowers()              const { return fMaxEInEmTowers;            }
59	>	Double_t   MaxEInHadTowers()             const { return fMaxEInHadTowers;           }
60	>	FourVector Mom()                         const { return fMom;                       }
61	>	UShort_t   N()                           const { return fN;                         }
62	>	UShort_t   N60()                         const { return fN60;                       }
63	>	UShort_t   N90()                         const { return fN90;                       }
64	>	UInt_t     NTowers()                     const { return fTowers.GetEntries();       }
65	>	Double_t   L2RelativeCorrectionScale()   const { return fL2RelativeCorrectionScale; }
66	>	Double_t   L3AbsoluteCorrectionScale()   const { return fL3AbsoluteCorrectionScale; }
67	>	Double_t   L4EMFCorrectionScale()        const { return fL4EMFCorrectionScale;      }
68	>	Double_t   L5FlavorCorrectionScale()     const { return fL5FlavorCorrectionScale;   }
69	>	Double_t   L7PartonCorrectionScale()     const { return fL7PartonCorrectionScale;   }
70	>	EObjType    ObjType()                    const { return kJet;                       }
71	>	Double_t   SimpleSecondaryVertexBJetTagsDisc()            const
72	>	{ return fSimpleSecondaryVertexBJetTagsDisc;      }
73	>	Double_t   SoftMuonBJetTagsDisc()                         const
74	>	{ return fSoftMuonBJetTagsDisc;                   }
75	>	Double_t   SoftMuonNoIPBJetTagsDisc()                     const
76	>	{ return fSoftMuonNoIPBJetTagsDisc;               }
77	>	Double_t   SoftElectronBJetTagsDisc()                     const
78	>	{ return fSoftElectronBJetTagsDisc;               }
79	>	Double_t   TrackCountingHighEffBJetTagsDisc()             const
80	>	{ return fTrackCountingHighEffBJetTagsDisc;       }
81	>	Double_t   TrackCountingHighPurBJetTagsDisc()             const
82	>	{ return fTrackCountingHighPurBJetTagsDisc;       }
83	>	const CaloTower
84	>	*Tower(UInt_t i)               const { return fTowers.At(i);           }
85	>	Double_t   TowersArea()                  const { return fTowersArea;             }
86	>	void       AddTower(const CaloTower *tower)  { fTowers.Add(tower);               }
87	>	void       SetAlpha(Double_t val)            { fAlpha            = val;          }
88	>	void       SetBeta(Double_t val)             { fBeta             = val;          }
89	>	void       SetCombinedSecondaryVertexBJetTagsDisc(Double_t d)
90		{ fCombinedSecondaryVertexBJetTagsDisc = d;     }
91	<	void       SetCombinedSecondaryVertexMVABJetTagsDisc ( Double_t d )
91	>	void       SetCombinedSecondaryVertexMVABJetTagsDisc(Double_t d)
92		{ fCombinedSecondaryVertexMVABJetTagsDisc = d;  }
93	<	void       SetImpactParameterMVABJetTagsDisc ( Double_t d )
93	>	void       SetEmEnergyInEB(Double_t val)     { fEmEnergyInEB     = val;          }
94	>	void       SetEmEnergyInEE(Double_t val)     { fEmEnergyInEE     = val;          }
95	>	void       SetEmEnergyInHF(Double_t val)     { fEmEnergyInHF     = val;          }
96	>	void       SetEnergyFractionH(Double_t val)  { fEnergyFractionH  = val;          }
97	>	void       SetEnergyFractionEm(Double_t val) { fEnergyFractionEm = val;          }
98	>	void       SetHadEnergyInHO(Double_t val)    { fHadEnergyInHO    = val;          }
99	>	void       SetHadEnergyInHB(Double_t val)    { fHadEnergyInHB    = val;          }
100	>	void       SetHadEnergyInHF(Double_t val)    { fHadEnergyInHF    = val;          }
101	>	void       SetHadEnergyInHE(Double_t val)    { fHadEnergyInHE    = val;          }
102	>	void       SetImpactParameterMVABJetTagsDisc (Double_t d)
103		{   fImpactParameterMVABJetTagsDisc = d;        }
104	<	void       SetTrackCountingHighEffBJetTagsDisc ( Double_t d )
105	<	{ fTrackCountingHighEffBJetTagsDisc = d;        }
106	<	void       SetTrackCountingHighPurBJetTagsDisc ( Double_t d )
107	<	{ fTrackCountingHighPurBJetTagsDisc = d;        }
108	<	void       SetSoftMuonBJetTagsDisc ( Double_t d )        { fSoftMuonBJetTagsDisc = d;        }
109	<	void       SetSoftMuonNoIPBJetTagsDisc ( Double_t d )    { fSoftMuonNoIPBJetTagsDisc = d;    }
110	<	void       SetSoftElectronBJetTagsDisc ( Double_t d )    { fSoftElectronBJetTagsDisc = d;    }
111	<
104	>	void       SetJetProbabilityBJetTagsDisc(Double_t d)  { fJetProbabilityBJetTagsDisc = d;     }
105	>	void       SetJetBProbabilityBJetTagsDisc(Double_t d) { fJetBProbabilityBJetTagsDisc = d;    }
106	>	void       SetMom(const FourVector &mom)     { fMom              = mom;          }
107	>	void       SetN(UShort_t n)                  { fN                = n;            }
108	>	void       SetN60(UShort_t n)                { fN60              = n;            }
109	>	void       SetN90(UShort_t n)                { fN90              = n;            }
110	>	void       SetMatchedMCFlavor(Int_t flavor)  { fMatchedMCFlavor = flavor;        }
111	>	void       SetMaxEInEmTowers(Double_t val)   { fMaxEInEmTowers   = val;          }
112	>	void       SetMaxEInHadTowers(Double_t val)  { fMaxEInHadTowers  = val;          }
113	>	void       SetL2RelativeCorrectionScale(Double_t s )   { fL2RelativeCorrectionScale = s;   }
114	>	void       SetL3AbsoluteCorrectionScale(Double_t s )   { fL3AbsoluteCorrectionScale = s;   }
115	>	void       SetL4EMFCorrectionScale(Double_t s )        { fL4EMFCorrectionScale = s;        }
116	>	void       SetL5FlavorCorrectionScale(Double_t s )     { fL5FlavorCorrectionScale = s;     }
117	>	void       SetL7PartonCorrectionScale(Double_t s )     { fL7PartonCorrectionScale = s;     }
118	>	void       SetSimpleSecondaryVertexBJetTagsDisc(Double_t d)
119	>	{ fSimpleSecondaryVertexBJetTagsDisc = d;     }
120	>	void       SetSoftMuonBJetTagsDisc(Double_t d )        { fSoftMuonBJetTagsDisc = d;        }
121	>	void       SetSoftMuonNoIPBJetTagsDisc(Double_t d )    { fSoftMuonNoIPBJetTagsDisc = d;    }
122	>	void       SetSoftElectronBJetTagsDisc(Double_t d )    { fSoftElectronBJetTagsDisc = d;    }
123	>	void       SetTowersArea(Double_t val)       { fTowersArea       = val;          }
124	>	void       SetTrackCountingHighEffBJetTagsDisc( Double_t d)
125	>	{ fTrackCountingHighEffBJetTagsDisc = d;      }
126	>	void       SetTrackCountingHighPurBJetTagsDisc(Double_t d)
127	>	{ fTrackCountingHighPurBJetTagsDisc = d;      }
128	>
129		protected:
130		FourVector fMom;                    //four momentum of jet
131		Double_t   fMaxEInEmTowers;         //maximum energy in EM towers
#	Line 125 | Line 143 | namespace mithep
143		UShort_t   fN;                      //number of constituents
144		UShort_t   fN60;                    //number constituents with 60% of total energy
145		UShort_t   fN90;                    //number constituents with 90% of total energy
146	+	Double_t   fAlpha;                  //jet vertex association alpha variable
147	+	Double_t   fBeta;                   //jet vertex association beta variable
148		Int_t      fMatchedMCFlavor;        //pdgID of the quark flavor that the jet matched to
149	<	Double_t   fJetProbabilityBJetTagsDisc;              //Discriminants for various b-tagging
150	<	Double_t   fJetBProbabilityBJetTagsDisc;             //algorithms
151	<	Double_t   fSimpleSecondaryVertexBJetTagsDisc;
152	<	Double_t   fCombinedSecondaryVertexBJetTagsDisc;
153	<	Double_t   fCombinedSecondaryVertexMVABJetTagsDisc;
154	<	Double_t   fImpactParameterMVABJetTagsDisc;
155	<	Double_t   fTrackCountingHighEffBJetTagsDisc;
156	<	Double_t   fTrackCountingHighPurBJetTagsDisc;
157	<	Double_t   fSoftMuonBJetTagsDisc;
158	<	Double_t   fSoftMuonNoIPBJetTagsDisc;
159	<	Double_t   fSoftElectronBJetTagsDisc;
160	<
161	<	RefArray<CaloTower> fTowers;        //||CaloTowers in this jet
149	>	Double_t   fJetProbabilityBJetTagsDisc;              //Discriminants for b-tagging algos
150	>	Double_t   fJetBProbabilityBJetTagsDisc;             //
151	>	Double_t   fSimpleSecondaryVertexBJetTagsDisc;       //
152	>	Double_t   fCombinedSecondaryVertexBJetTagsDisc;     //
153	>	Double_t   fCombinedSecondaryVertexMVABJetTagsDisc;  //
154	>	Double_t   fImpactParameterMVABJetTagsDisc;          //
155	>	Double_t   fTrackCountingHighEffBJetTagsDisc;        //
156	>	Double_t   fTrackCountingHighPurBJetTagsDisc;        //
157	>	Double_t   fSoftMuonBJetTagsDisc;                    //
158	>	Double_t   fSoftMuonNoIPBJetTagsDisc;                //
159	>	Double_t   fSoftElectronBJetTagsDisc;                //
160	>	Double_t   fL2RelativeCorrectionScale;               //
161	>	Double_t   fL3AbsoluteCorrectionScale;               //
162	>	Double_t   fL4EMFCorrectionScale;                    //
163	>	Double_t   fL5FlavorCorrectionScale;                 //
164	>	Double_t   fL7PartonCorrectionScale;                 //
165	>	RefArray<CaloTower,1024> fTowers;                    //calo towers assigned to this jet
166
167		ClassDef(Jet, 1) // Jet class
168		};

Diff Legend

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