MSLayersAtAngle Class Reference

#include <MSLayersAtAngle.h>

Public Member Functions
const MSLayer *	findLayer (const MSLayer &layer) const
	MSLayersAtAngle ()
	MSLayersAtAngle (const std::vector< MSLayer > &layers)
void	print () const
int	size () const
float	sumX0D (const PixelRecoPointRZ &pointI, const PixelRecoPointRZ &pointO) const
float	sumX0D (int il, int ol, const PixelRecoPointRZ &pointI, const PixelRecoPointRZ &pointO) const
float	sumX0D (const PixelRecoPointRZ &pointI, const PixelRecoPointRZ &pointM, const PixelRecoPointRZ &pointO) const
float	sumX0D (float zV, int il, int ol, const PixelRecoPointRZ &pointI, const PixelRecoPointRZ &pointO) const
void	update (const MSLayer &layer)

Private Types
typedef std::vector< MSLayer >::const_iterator	LayerItr

Private Member Functions
LayerItr	findLayer (const PixelRecoPointRZ &point, LayerItr i1, LayerItr i2) const
void	init ()
float	sum2RmRn (LayerItr i1, LayerItr i2, float rTarget, const SimpleLineRZ &line) const

Private Attributes
std::vector< int >	indeces
std::vector< MSLayer >	theLayers

Detailed Description

Definition at line 16 of file MSLayersAtAngle.h.

Member Typedef Documentation

LayerItr

typedef std::vector<MSLayer>::const_iterator MSLayersAtAngle::LayerItr

private

Definition at line 43 of file MSLayersAtAngle.h.

Constructor & Destructor Documentation

MSLayersAtAngle() [1/2]

MSLayersAtAngle::MSLayersAtAngle ( )

inline

Definition at line 18 of file MSLayersAtAngle.h.

{}

MSLayersAtAngle() [2/2]

MSLayersAtAngle::MSLayersAtAngle

(

const std::vector< MSLayer > &

layers

)

Definition at line 32 of file MSLayersAtAngle.cc.

References init().

: theLayers(layers) { init(); }

Member Function Documentation

findLayer() [1/2]

const MSLayer * MSLayersAtAngle::findLayer

(

const MSLayer &

layer

)

const

Definition at line 34 of file MSLayersAtAngle.cc.

References spr::find(), pixelTopology::layer, and theLayers.

Referenced by MSLayersKeeperX0AtEta::MSLayersKeeperX0AtEta(), MSLayersKeeperX0Averaged::MSLayersKeeperX0Averaged(), sumX0D(), MSLayer::sumX0D(), and MSLayer::x0().

                                                                    {
  vector<MSLayer>::const_iterator it = find(theLayers.begin(), theLayers.end(), layer);
  return it == theLayers.end() ? nullptr : &(*it);
}

findLayer() [2/2]

MSLayersAtAngle::LayerItr MSLayersAtAngle::findLayer ( const PixelRecoPointRZ &  point, MSLayersAtAngle::LayerItr  ibeg, MSLayersAtAngle::LayerItr  iend  ) const

private

Definition at line 151 of file MSLayersAtAngle.cc.

References ztail::d, kinem::EPSILON, f, and point.

                                                                                         {
  const float BIG = 99999.f;
  const float EPSILON = 1.e-8f;
  LayerItr theIt = ibeg;
  float dist = BIG;
  for (LayerItr it = ibeg; it < iend; it++) {
    float d = it->distance2(point);
    if (d < dist) {
      if (d < EPSILON)
        return it;
      dist = d;
      theIt = it;
    }
  }
  return theIt;
}

init()

void MSLayersAtAngle::init ( void  )

private

Definition at line 16 of file MSLayersAtAngle.cc.

References mps_fire::i, MainPageGenerator::l, and jetUpdater_cfi::sort.

Referenced by MSLayersAtAngle(), and update().

                           {
  std::sort(theLayers.begin(), theLayers.end());
  int i = -1;
  indeces.clear();
  for (auto const& l : theLayers) {
    ++i;
    int sq = l.seqNum();
    if (sq < 0)
      continue;
    if (sq >= int(indeces.size()))
      indeces.resize(sq + 1, -1);
    indeces[sq] = i;
  }
}

print()

void MSLayersAtAngle::print

(

void

)

const

Definition at line 171 of file MSLayersAtAngle.cc.

References gather_cfg::cout, and theLayers.

                                  {
  for (LayerItr it = theLayers.begin(); it != theLayers.end(); it++)
    cout << *it << endl;
}

size()

int MSLayersAtAngle::size ( void  ) const

inline

Definition at line 33 of file MSLayersAtAngle.h.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return theLayers.size(); }

sum2RmRn()

float MSLayersAtAngle::sum2RmRn ( MSLayersAtAngle::LayerItr  i1, MSLayersAtAngle::LayerItr  i2, float  rTarget, const SimpleLineRZ &  line  ) const

private

Definition at line 131 of file MSLayersAtAngle.cc.

References cross(), l1ctLayer1_cff::dr, MillePedeFileConverter_cfg::e, f, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, and mps_splice::line.

Referenced by sumX0D().

                                                                {
  float sum2 = 0.f;
  float cotTh = line.cotLine();
  for (LayerItr it = i1; it < i2; it++) {
    std::pair<PixelRecoPointRZ, bool> cross = it->crossing(line);
    if (cross.second) {
      float x0 = it->x0(cotTh);
      float dr = rTarget - cross.first.r();
      if (x0 > 1.e-5f)
        dr *= 1.f + 0.038f * unsafe_logf<2>(x0);
      sum2 += x0 * dr * dr;
    }
    //  cout << *it << " crossing: "<<cross.second<<endl;
  }
  return sum2;
}

sumX0D() [1/4]

float MSLayersAtAngle::sumX0D	(	const PixelRecoPointRZ &	pointI,
		const PixelRecoPointRZ &	pointO
	)		const

Definition at line 51 of file MSLayersAtAngle.cc.

References findLayer(), PixelRecoPointRZ::r(), mathSSE::sqrt(), sum2RmRn(), and theLayers.

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

                                                                                                  {
  LayerItr iO = findLayer(pointO, theLayers.begin(), theLayers.end());
  //  cout << "outer Layer: "<<*iO<<endl;
  LayerItr iI = findLayer(pointI, theLayers.begin(), iO);
  //  cout << "inner Layer: "<<*iI<<endl;

  return sqrt(sum2RmRn(iI, iO, pointO.r(), SimpleLineRZ(pointI, pointO)));
}

sumX0D() [2/4]

float MSLayersAtAngle::sumX0D	(	int	il,
		int	ol,
		const PixelRecoPointRZ &	pointI,
		const PixelRecoPointRZ &	pointO
	)		const

Definition at line 60 of file MSLayersAtAngle.cc.

References indeces, PixelRecoPointRZ::r(), mathSSE::sqrt(), sum2RmRn(), sumX0D(), and theLayers.

                                                                                                                  {
  // if layer not at this angle (WHY???) revert to slow comp
  if (il >= int(indeces.size()) || ol >= int(indeces.size()) || indeces[il] < 0 || indeces[ol] < 0)
    return sumX0D(pointI, pointO);

  LayerItr iI = theLayers.begin() + indeces[il];
  LayerItr iO = theLayers.begin() + indeces[ol];

  return sqrt(sum2RmRn(iI, iO, pointO.r(), SimpleLineRZ(pointI, pointO)));
}

sumX0D() [3/4]

float MSLayersAtAngle::sumX0D	(	const PixelRecoPointRZ &	pointI,
		const PixelRecoPointRZ &	pointM,
		const PixelRecoPointRZ &	pointO
	)		const

Definition at line 75 of file MSLayersAtAngle.cc.

References gather_cfg::cout, doPrint, findLayer(), mps_splice::line, PixelRecoPointRZ::r(), funct::sqr(), mathSSE::sqrt(), sum2RmRn(), theLayers, and PixelRecoPointRZ::z().

                                                                    {
  LayerItr iO = findLayer(pointO, theLayers.begin(), theLayers.end());
  LayerItr iI = findLayer(pointI, theLayers.begin(), iO);
  LayerItr iM = findLayer(pointM, iI, iO);

  float drOI = pointO.r() - pointI.r();
  float drMO = pointO.r() - pointM.r();
  float drMI = pointM.r() - pointI.r();

  SimpleLineRZ line(pointI, pointO);
  float sum2I = sum2RmRn(iI + 1, iM, pointI.r(), line);
  float sum2O = sum2RmRn(iM, iO, pointO.r(), line);

  float sum = std::sqrt(sum2I * sqr(drMO) + sum2O * sqr(drMI)) / drOI;

  // if (iI!=theLayers.begin() )
  // doPrint = ((*iM).seqNum()<0 || (*iO).seqNum()<0);
  if (doPrint)
    std::cout << "old " << (*iI).seqNum() << " " << iI - theLayers.begin() << ", " << (*iM).seqNum() << " "
              << iM - theLayers.begin() << ", " << (*iO).seqNum() << " " << iO - theLayers.begin() << "  " << pointI.r()
              << " " << pointI.z() << " " << sum << std::endl;

  return sum;
}

sumX0D() [4/4]

float MSLayersAtAngle::sumX0D	(	float	zV,
		int	il,
		int	ol,
		const PixelRecoPointRZ &	pointI,
		const PixelRecoPointRZ &	pointO
	)		const

Definition at line 102 of file MSLayersAtAngle.cc.

References gather_cfg::cout, doPrint, f, indeces, mps_splice::line, PixelRecoPointRZ::r(), funct::sqr(), mathSSE::sqrt(), sum2RmRn(), sumX0D(), and theLayers.

                                                                                                    {
  PixelRecoPointRZ pointV(0.f, zV);

  // if layer not at this angle (WHY???) revert to slow comp
  if (il >= int(indeces.size()) || ol >= int(indeces.size()) || indeces[il] < 0 || indeces[ol] < 0)
    return sumX0D(pointV, pointI, pointO);

  LayerItr iI = theLayers.begin() + indeces[il];
  LayerItr iO = theLayers.begin() + indeces[ol];

  float drOI = pointO.r();
  float drMO = pointO.r() - pointI.r();
  float drMI = pointI.r();

  SimpleLineRZ line(pointV, pointO);
  float sum2I = sum2RmRn(theLayers.begin() + 1, iI, pointV.r(), line);
  float sum2O = sum2RmRn(iI, iO, pointO.r(), line);

  float sum = std::sqrt(sum2I * sqr(drMO) + sum2O * sqr(drMI)) / drOI;

  if (doPrint)
    std::cout << "new " << il << " " << (*iI).seqNum() << " " << iI - theLayers.begin() << ", " << ol << " "
              << (*iO).seqNum() << " " << iO - theLayers.begin() << "  " << zV << " " << sum << std::endl;

  return sum;
}

update()

void MSLayersAtAngle::update

(

const MSLayer &

layer

)

Definition at line 40 of file MSLayersAtAngle.cc.

References spr::find(), init(), pixelTopology::layer, and theLayers.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                 {
  vector<MSLayer>::iterator it = find(theLayers.begin(), theLayers.end(), layer);
  if (it == theLayers.end()) {
    theLayers.push_back(layer);
    init();
  } else {
    *it = layer;
  }
}

Member Data Documentation

indeces

std::vector<int> MSLayersAtAngle::indeces

private

Definition at line 38 of file MSLayersAtAngle.h.

Referenced by sumX0D().

theLayers

std::vector<MSLayer> MSLayersAtAngle::theLayers

private

Definition at line 37 of file MSLayersAtAngle.h.

Referenced by findLayer(), print(), sumX0D(), and update().