16 #pragma implementation
32 }
else if (
ee < w.
ee ) {
59 double limit = std::fabs(
pp.
dot(w.
pp));
64 return (delta <= limit );
70 if ( (wdw > 0) && (delta < wdw) ) {
71 return std::sqrt (delta/wdw);
72 }
else if ( (wdw == 0) && (delta == 0) ) {
87 double tTotal = (ee + w.
ee);
89 double vTotal2 = vTotal.
mag2();
91 if ( vTotal2 >= tTotal*tTotal ) {
100 return (isNear(w, epsilon));
105 double tRecip = 1./tTotal;
116 double b2 = vTotal2*tRecip*tRecip;
118 double ggamma = std::sqrt(1./(1.-b2));
119 double boostDotV1 = bboost.
dot(
pp);
120 double gm1_b2 = (ggamma-1)/b2;
123 ggamma * (ee + boostDotV1) );
125 double boostDotV2 = bboost.dot(w.
pp);
127 ggamma * (w.
ee + boostDotV2) );
129 return (w1.
isNear(w2, epsilon));
135 double tTotal = (
ee + w.
ee);
137 double vTotal2 = vTotal.
mag2();
139 if ( vTotal2 >= tTotal*tTotal ) {
151 if ( vTotal2 == 0 ) {
157 double tRecip = 1./tTotal;
168 double b2 = vTotal2*tRecip*tRecip;
173 double ggamma = std::sqrt(1./(1.-b2));
174 double boostDotV1 = bboost.
dot(
pp);
175 double gm1_b2 = (ggamma-1)/b2;
178 ggamma * (ee + boostDotV1) );
180 double boostDotV2 = bboost.dot(w.
pp);
182 ggamma * (w.
ee + boostDotV2) );
200 return std::sqrt ( a*a + b*b );
245 return (x1 < 1) ? x1 : 1;
251 double twoT2 = 2*
ee*
ee;
std::vector< ExP01TrackerHit * > a
double deltaR(const HepLorentzVector &v) const
int compare(const HepLorentzVector &w) const
bool operator<(const HepLorentzVector &w) const
Float_t x1[n_points_granero]
int compare(const Hep3Vector &v) const
double deltaPhi(const Hep3Vector &v2) const
double dot(const Hep3Vector &) const
bool operator<=(const HepLorentzVector &w) const
bool isParallel(const HepLorentzVector &w, double epsilon=tolerance) const
bool operator>(const HepLorentzVector &w) const
double euclideanNorm2() const
bool isNearCM(const HepLorentzVector &w, double epsilon=tolerance) const
double howLightlike() const
double epsilon(double density, double temperature)
bool operator>=(const HepLorentzVector &w) const
double howParallel(const HepLorentzVector &w) const
bool isNear(const HepLorentzVector &w, double epsilon=tolerance) const
double howNearCM(const HepLorentzVector &w) const
double howNear(const HepLorentzVector &w) const
double euclideanNorm() const