| Revision: | 1.1 |
| Committed: | Tue Nov 11 23:01:22 2008 UTC (16 years, 5 months ago) by kukartse |
| Content type: | text/plain |
| Branch: | MAIN |
| CVS Tags: | V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15, V00-01-14, V00-01-13, V00-01-12, V00-01-11, V00-01-10, gak031009, gak030509, gak022309, gak021209, gak040209, gak012809, V00-01-09, V00-01-08, V00-01-07, V00-01-06, V00-01-05, V00-01-04, V00-00-07, V00-00-06, V00-00-05, V00-00-04, V00-01-03, V00-00-02, V00-00-01, HEAD |
| Branch point for: | ZMorph-V00-03-01, CMSSW_22X_branch |
| Error occurred while calculating annotation data. | |
| Log Message: | initial creation of a package for LJMET multivariate analysis |
| # | Content |
|---|---|
| 1 | #include "LJMet/MultivariateAnalysis/interface/TMBLorentzVector.h" |
| 2 | #include "TLorentzRotation.h" |
| 3 | |
| 4 | Double_t TMBLorentzVector::operator () (int i) const { |
| 5 | // return i'th component of 4 vector |
| 6 | switch (i) |
| 7 | { |
| 8 | case kX: |
| 9 | case kY: |
| 10 | case kZ: |
| 11 | return Vect ()(i); |
| 12 | case kE:return E (); |
| 13 | default:Error ("operator()()", "bad index (%d) returning 0", i); |
| 14 | } |
| 15 | return 0.; |
| 16 | } |
| 17 | |
| 18 | void TMBLorentzVector::Boost(const TVector3 & b) { |
| 19 | // apply Lorentz boost by vector b |
| 20 | Double_t b2 = b.Mag2 (); |
| 21 | Double_t gamma = 1.0 / TMath::Sqrt (1.0 - b2); |
| 22 | Double_t bp = b.Dot (Vect ()); |
| 23 | Double_t gamma2 = b2 > 0 ? (gamma - 1.0) / b2 : 0.0; |
| 24 | Vect () += (gamma2 * bp + gamma * E ()) * b; |
| 25 | // SetE (gamma * (E () + bp)); |
| 26 | } |
| 27 | |
| 28 | TMBLorentzVector & TMBLorentzVector::operator *= (const TLorentzRotation & m) { |
| 29 | // apply lorentz rotation, see TLorentzRotation::VectorMultiplication() |
| 30 | return *this = m.VectorMultiplication((TLorentzVector)*this); |
| 31 | } |
| 32 | |
| 33 | TMBLorentzVector & TMBLorentzVector::Transform (const TLorentzRotation & m){ |
| 34 | // apply lorentz transformation, see TLorentzRotation::VectorMultiplication() |
| 35 | return *this = m.VectorMultiplication((TLorentzVector)*this); |
| 36 | } |
| 37 | |
| 38 | /// Equivalence operator. True if all components are |
| 39 | /// equivalent within machine precision. |
| 40 | Bool_t TMBLorentzVector::is_equal (const TMBLorentzVector &lv) const |
| 41 | { |
| 42 | return TMBVector3::is_equal (lv.E(), E()) && Vect().is_equal (lv.Vect()); |
| 43 | } |
| 44 | |
| 45 | //ClassImp(TMBLorentzVector) |