10 using namespace GeomDetEnumerators;
12 template <
class T>
T sqr(
T t) {
return t*
t;}
17 s <<
" face: "<<l.
face()
19 <<
" range:"<<l.
range()<<
", "
27 else if (d.
allLayers) s <<
"x0 by MSLayersKeeper";
28 else s <<
"empty DataX0";
34 : theFace(layer->location()), theSeqNum(layer->seqNum()), theX0Data(dataX0)
54 cout <<
" ** MSLayer ** unknown part - will not work!" <<endl;
62 thePosition(position),
64 theHalfThickness(halfThickness),
91 const float eps = 1.e-5;
106 const float eps = 1.e-5;
140 float OverSinTheta =
std::sqrt(1.
f+cotTheta*cotTheta);
146 if (dataLayer)
return dataLayer->
x0(cotTheta);
172 if (dataLayer)
return dataLayer->
sumX0D(cotTheta);
virtual const BoundSurface & surface() const =0
The surface of the GeometricSearchDet.
const MSLayersKeeper * allLayers
virtual float length() const =0
float sumX0D(float cotTheta) const
const Bounds & bounds() const
float x0(float cotTheta) const
std::ostream & operator<<(std::ostream &out, const ALILine &li)
virtual float thickness() const =0
float zAtR(float r) const
const Range & range() const
const GeomDetEnumerators::Location & face() const
void swap(edm::DataFrameContainer &lhs, edm::DataFrameContainer &rhs)
virtual const BoundDisk & specificSurface() const
bool operator<(const MSLayer &o) const
Abs< T >::type abs(const T &t)
bool operator==(const MSLayer &o) const
virtual const MSLayersAtAngle & layers(float cotTheta) const =0
std::pair< PixelRecoPointRZ, bool > crossing(const PixelRecoLineRZ &line) const
virtual const BoundSurface & surface() const
GeometricSearchDet interface.
virtual const Surface::PositionType & position() const
Returns position of the surface.
bool inside(const T &value) const
virtual const BoundCylinder & specificSurface() const
Extension of the interface.
float rAtZ(float z) const
Square< F >::type sqr(const F &f)
static int position[264][3]
float rAtZ(float z) const
float distance2(const PixelRecoPointRZ &point) const
float zAtR(float r) const
GeomDetEnumerators::Location theFace
PixelRecoRange< float > Range
const MSLayer * findLayer(const MSLayer &layer) const
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point