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

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

private

Definition at line 51 of file MSLayersAtAngle.h.

Constructor & Destructor Documentation

MSLayersAtAngle::MSLayersAtAngle ( )

inline

Definition at line 19 of file MSLayersAtAngle.h.

{ }

MSLayersAtAngle::MSLayersAtAngle

(

const std::vector< MSLayer > &

layers

)

Definition at line 28 of file MSLayersAtAngle.cc.

References init().

  : theLayers(layers)
{ 
  init();
}

Member Function Documentation

const MSLayer * MSLayersAtAngle::findLayer

(

const MSLayer &

layer

)

const

Definition at line 34 of file MSLayersAtAngle.cc.

References spr::find(), and theLayers.

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

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

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

private

Definition at line 181 of file MSLayersAtAngle.cc.

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

{
  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;
}

void MSLayersAtAngle::init ( void  )

private

Definition at line 15 of file MSLayersAtAngle.cc.

References i, prof2calltree::l, and python.multivaluedict::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;
  }
}

void MSLayersAtAngle::print

(

void

)

const

Definition at line 201 of file MSLayersAtAngle.cc.

References gather_cfg::cout, and theLayers.

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

int MSLayersAtAngle::size ( void  ) const

inline

Definition at line 41 of file MSLayersAtAngle.h.

{ return theLayers.size(); }

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

private

Definition at line 160 of file MSLayersAtAngle.cc.

References SimpleLineRZ::cotLine(), cross(), alignCSCRings::e, and f.

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;
}

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

Definition at line 54 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)));
}

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

Definition at line 68 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)));
}

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

Definition at line 88 of file MSLayersAtAngle.cc.

References gather_cfg::cout, doPrint, findLayer(), geometryCSVtoXML::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;
}

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

Definition at line 123 of file MSLayersAtAngle.cc.

References gather_cfg::cout, doPrint, f, indeces, geometryCSVtoXML::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;

}

void MSLayersAtAngle::update

(

const MSLayer &

layer

)

Definition at line 42 of file MSLayersAtAngle.cc.

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

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
  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

std::vector<int> MSLayersAtAngle::indeces

private

Definition at line 46 of file MSLayersAtAngle.h.

Referenced by sumX0D().

std::vector<MSLayer> MSLayersAtAngle::theLayers

private

Definition at line 45 of file MSLayersAtAngle.h.

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