CMS 3D CMS Logo

CSCCrossGap.cc
Go to the documentation of this file.
3 #include <cmath>
4 
5 #include <iostream>
7  : theBeta2(0.), theGamma(1.), loggam(0.), theGap(gap), clusters(), electronsInClusters(), steps(), elosses() {
8  logGamma(mass, mom);
9  LogTrace("CSCCrossGap") << "CSCCrossGap: simhit \n"
10  << "mass = " << mass << "GeV/c2, momentum = " << mom << " GeV/c, gap length = " << length()
11  << " cm \n";
12 }
13 
14 double CSCCrossGap::logGamma(double mass, float mom) {
15  theGamma = sqrt((mom / mass) * (mom / mass) + 1.);
16  theBeta2 = 1. - 1. / (theGamma * theGamma);
17  LogTrace("CSCCrossGap") << "gamma = " << theGamma << ", beta2 = " << theBeta2
18  << ", beta*gamma = " << sqrt(theGamma * theGamma - 1.);
19  ;
20 
21  // The lowest value in table (=theGammaBins[0]) is ln(1.1)=0.0953102
22  // (Compensate later if lower)
23  loggam = log(std::max(1.1, theGamma)); // F-P literal IS double by default!
24  LogTrace("CSCCrossGap") << "logGamma = " << loggam;
25 
26  return loggam;
27 }
#define LogTrace(id)
T sqrt(T t)
Definition: SSEVec.h:19
double loggam
Definition: CSCCrossGap.h:59
double theBeta2
Definition: CSCCrossGap.h:57
double theGamma
Definition: CSCCrossGap.h:58
double logGamma()
Definition: CSCCrossGap.h:49
CSCCrossGap(double mass, float mom, LocalVector gap)
Definition: CSCCrossGap.cc:6
float length() const
Definition: CSCCrossGap.h:54