MSLayer Class Reference

#include <MSLayer.h>

Classes
struct	DataX0

Public Types
typedef PixelRecoRange< float >	Range

Public Member Functions
std::pair< PixelRecoPointRZ, bool >	crossing (const PixelRecoLineRZ &line) const
std::pair< PixelRecoPointRZ, bool >	crossing (const SimpleLineRZ &line) const
float	distance2 (const PixelRecoPointRZ &point) const
const GeomDetEnumerators::Location &	face () const
float	halfThickness () const
	MSLayer (const DetLayer *layer, const DataX0 &dataX0=DataX0(nullptr))
	MSLayer ()
	MSLayer (GeomDetEnumerators::Location part, float position, Range range, float halfThickness=0., const DataX0 &dataX0=DataX0(nullptr))
bool	operator< (const MSLayer &o) const
bool	operator== (const MSLayer &o) const
float	position () const
const Range &	range () const
int	seqNum () const
float	sumX0D (float cotTheta) const
float	x0 (float cotTheta) const

Private Attributes
GeomDetEnumerators::Location	theFace
float	theHalfThickness
float	thePosition
Range	theRange
int	theSeqNum
DataX0	theX0Data

Friends
class	MSLayersKeeper
std::ostream &	operator<< (std::ostream &s, const MSLayer &l)

Detailed Description

Definition at line 13 of file MSLayer.h.

Member Typedef Documentation

Range

typedef PixelRecoRange<float> MSLayer::Range

Definition at line 15 of file MSLayer.h.

Constructor & Destructor Documentation

MSLayer() [1/3]

MSLayer::MSLayer	(	const DetLayer *	layer,
		const DataX0 &	dataX0 = DataX0(nullptr)
	)

Definition at line 35 of file MSLayer.cc.

References Reference_intrackfit_cff::barrel, Surface::bounds(), gather_cfg::cout, makeMuonMisalignmentScenario::endcap, pixelTopology::layer, Bounds::length(), GeometricSearchDet::position(), ForwardDetLayer::specificSurface(), BarrelDetLayer::specificSurface(), BarrelDetLayer::surface(), theFace, theHalfThickness, thePosition, theRange, and PV3DBase< T, PVType, FrameType >::z().

    : theFace(layer->location()), theSeqNum(layer->seqNum()), theX0Data(dataX0) {
  const BarrelDetLayer* bl;
  const ForwardDetLayer* fl;
  theHalfThickness = layer->surface().bounds().thickness() / 2;

  switch (theFace) {
    case barrel:
      bl = static_cast<const BarrelDetLayer*>(layer);
      thePosition = bl->specificSurface().radius();
      theRange = Range(-bl->surface().bounds().length() / 2, bl->surface().bounds().length() / 2);
      break;
    case endcap:
      fl = static_cast<const ForwardDetLayer*>(layer);
      thePosition = fl->position().z();
      theRange = Range(fl->specificSurface().innerRadius(), fl->specificSurface().outerRadius());
      break;
    default:
      // should throw or simimal
      cout << " ** MSLayer ** unknown part - will not work!" << endl;
      break;
  }
}

MSLayer() [2/3]

MSLayer::MSLayer ( )

inline

Definition at line 38 of file MSLayer.h.

{}

MSLayer() [3/3]

MSLayer::MSLayer	(	GeomDetEnumerators::Location	part,
		float	position,
		Range	range,
		float	halfThickness = 0.,
		const DataX0 &	dataX0 = DataX0(nullptr)
	)

Definition at line 59 of file MSLayer.cc.

    : theFace(part),
      thePosition(position),
      theRange(range),
      theHalfThickness(halfThickness),
      theSeqNum(-1),
      theX0Data(dataX0) {}

Member Function Documentation

crossing() [1/2]

pair< PixelRecoPointRZ, bool > MSLayer::crossing

(

const PixelRecoLineRZ &

line

)

const

Definition at line 84 of file MSLayer.cc.

References Reference_intrackfit_cff::barrel, HLT_2022v15_cff::eps, mps_splice::line, PixelRecoRange< T >::max(), PixelRecoRange< T >::min(), std::swap(), theFace, thePosition, theRange, and z.

Referenced by MultipleScatteringParametrisation::operator()().

                                                                                {
  const float eps = 1.e-5;
  bool inLayer = true;
  float value = (theFace == barrel) ? line.zAtR(thePosition) : line.rAtZ(thePosition);
  if (value > theRange.max()) {
    value = theRange.max() - eps;
    inLayer = false;
  } else if (value < theRange.min()) {
    value = theRange.min() + eps;
    inLayer = false;
  }
  float z = thePosition;
  if (theFace == barrel)
    std::swap(z, value);  // if barrel value is z
  return make_pair(PixelRecoPointRZ(value, z), inLayer);
}

crossing() [2/2]

pair< PixelRecoPointRZ, bool > MSLayer::crossing

(

const SimpleLineRZ &

line

)

const

Definition at line 100 of file MSLayer.cc.

References Reference_intrackfit_cff::barrel, HLT_2022v15_cff::eps, mps_splice::line, PixelRecoRange< T >::max(), PixelRecoRange< T >::min(), std::swap(), theFace, thePosition, theRange, and z.

                                                                             {
  const float eps = 1.e-5;
  bool inLayer = true;
  float value = (theFace == barrel) ? line.zAtR(thePosition) : line.rAtZ(thePosition);
  if (value > theRange.max()) {
    value = theRange.max() - eps;
    inLayer = false;
  } else if (value < theRange.min()) {
    value = theRange.min() + eps;
    inLayer = false;
  }
  float z = thePosition;
  if (theFace == barrel)
    std::swap(z, value);  // if barrel value is z
  return make_pair(PixelRecoPointRZ(value, z), inLayer);
}

distance2()

float MSLayer::distance2

(

const PixelRecoPointRZ &

point

)

const

Definition at line 118 of file MSLayer.cc.

References funct::abs(), Reference_intrackfit_cff::barrel, PixelRecoRange< T >::inside(), PixelRecoRange< T >::max(), PixelRecoRange< T >::min(), point, sqr(), theFace, theHalfThickness, thePosition, theRange, and findQualityFiles::v.

                                                            {
  float u = (theFace == barrel) ? point.r() : point.z();
  float v = (theFace == barrel) ? point.z() : point.r();

  float du = std::abs(u - thePosition);
  if (theRange.inside(v))
    return (du < theHalfThickness) ? 0.f : du * du;

  float dv = (v > theRange.max()) ? v - theRange.max() : theRange.min() - v;
  return sqr(du) + sqr(dv);
}

face()

const GeomDetEnumerators::Location& MSLayer::face ( ) const

inline

Definition at line 52 of file MSLayer.h.

References theFace.

{ return theFace; }

halfThickness()

float MSLayer::halfThickness ( ) const

inline

Definition at line 54 of file MSLayer.h.

References theHalfThickness.

{ return theHalfThickness; }

operator<()

bool MSLayer::operator<

(

const MSLayer &

o

)

const

Definition at line 72 of file MSLayer.cc.

References funct::abs(), Reference_intrackfit_cff::barrel, makeMuonMisalignmentScenario::endcap, PixelRecoRange< T >::max(), EcalTangentSkim_cfg::o, range(), theFace, and thePosition.

                                              {
  if (theFace == barrel && o.theFace == barrel)
    return thePosition < o.thePosition;
  else if (theFace == barrel && o.theFace == endcap)
    return thePosition < o.range().max();
  else if (theFace == endcap && o.theFace == endcap)
    return std::abs(thePosition) < std::abs(o.thePosition);
  else
    return range().max() < o.thePosition;
}

operator==()

bool MSLayer::operator==

(

const MSLayer &

o

)

const

Definition at line 68 of file MSLayer.cc.

References funct::abs(), f, EcalTangentSkim_cfg::o, theFace, and thePosition.

                                               {
  return theFace == o.theFace && std::abs(thePosition - o.thePosition) < 1.e-3f;
}

position()

float MSLayer::position ( ) const

inline

Definition at line 53 of file MSLayer.h.

References thePosition.

{ return thePosition; }

range()

const Range& MSLayer::range ( ) const

inline

Definition at line 50 of file MSLayer.h.

References theRange.

Referenced by operator<().

{ return theRange; }

seqNum()

int MSLayer::seqNum ( ) const

inline

Definition at line 47 of file MSLayer.h.

References theSeqNum.

Referenced by MultipleScatteringParametrisation::operator()().

{ return theSeqNum; }

sumX0D()

float MSLayer::sumX0D

(

float

cotTheta

)

const

Definition at line 144 of file MSLayer.cc.

References MSLayer::DataX0::allLayers, Reference_intrackfit_cff::barrel, MSLayer::DataX0::cotTheta, makeMuonMisalignmentScenario::endcap, f, MSLayersAtAngle::findLayer(), MSLayer::DataX0::hasFSlope, MSLayer::DataX0::hasX0, GeomDetEnumerators::invalidLoc, MSLayersKeeper::layers(), LIKELY, MSLayer::DataX0::slopeSumX0D, mathSSE::sqrt(), MSLayer::DataX0::sumX0D, sumX0D(), theFace, and theX0Data.

Referenced by MSLayersKeeperX0Averaged::MSLayersKeeperX0Averaged(), MultipleScatteringParametrisation::operator()(), and sumX0D().

                                          {
  if LIKELY (theX0Data.hasX0) {
    switch (theFace) {
      case barrel:
        return theX0Data.sumX0D *
               std::sqrt(std::sqrt((1.f + cotTheta * cotTheta) / (1.f + theX0Data.cotTheta * theX0Data.cotTheta)));
      case endcap:
        return (theX0Data.hasFSlope)
                   ? theX0Data.sumX0D + theX0Data.slopeSumX0D * (1.f / cotTheta - 1.f / theX0Data.cotTheta)
                   : theX0Data.sumX0D;
      case invalidLoc:
        break;  // make gcc happy
    }
  } else if (theX0Data.allLayers) {
    const MSLayer* dataLayer = theX0Data.allLayers->layers(cotTheta).findLayer(*this);
    if (dataLayer)
      return dataLayer->sumX0D(cotTheta);
  }
  return 0.;
}

x0()

float MSLayer::x0

(

float

cotTheta

)

const

Definition at line 131 of file MSLayer.cc.

References funct::abs(), MSLayer::DataX0::allLayers, Reference_intrackfit_cff::barrel, f, MSLayersAtAngle::findLayer(), MSLayer::DataX0::hasX0, MSLayersKeeper::layers(), LIKELY, mathSSE::sqrt(), theFace, theX0Data, MSLayer::DataX0::x0, and x0().

Referenced by MultipleScatteringParametrisation::operator()(), and x0().

                                      {
  if LIKELY (theX0Data.hasX0) {
    float OverSinTheta = std::sqrt(1.f + cotTheta * cotTheta);
    return (theFace == barrel) ? theX0Data.x0 * OverSinTheta : theX0Data.x0 * OverSinTheta / std::abs(cotTheta);
  } else if (theX0Data.allLayers) {
    const MSLayer* dataLayer = theX0Data.allLayers->layers(cotTheta).findLayer(*this);
    if (dataLayer)
      return dataLayer->x0(cotTheta);
  }
  return 0.;
}

Friends And Related Function Documentation

MSLayersKeeper

friend class MSLayersKeeper

friend

Definition at line 76 of file MSLayer.h.

operator<<

std::ostream& operator<< ( std::ostream &  s, const MSLayer &  l  )

friend

Member Data Documentation

theFace

GeomDetEnumerators::Location MSLayer::theFace

private

Definition at line 68 of file MSLayer.h.

Referenced by crossing(), distance2(), face(), MSLayer(), operator<(), operator==(), sumX0D(), and x0().

theHalfThickness

float MSLayer::theHalfThickness

private

Definition at line 71 of file MSLayer.h.

Referenced by distance2(), halfThickness(), and MSLayer().

thePosition

float MSLayer::thePosition

private

Definition at line 69 of file MSLayer.h.

Referenced by crossing(), distance2(), MSLayer(), operator<(), operator==(), and position().

theRange

Range MSLayer::theRange

private

Definition at line 70 of file MSLayer.h.

Referenced by crossing(), distance2(), MSLayer(), and range().

theSeqNum

int MSLayer::theSeqNum

private

Definition at line 72 of file MSLayer.h.

Referenced by seqNum().

theX0Data

DataX0 MSLayer::theX0Data

private

Definition at line 74 of file MSLayer.h.

Referenced by sumX0D(), and x0().