RBorderFinder Class Reference

#include <RBorderFinder.h>

Public Types
typedef ForwardDetRing	Det
typedef geomsort::ExtractR< Det, float >	DetR

Public Member Functions
bool	isROverlapping () const
	Returns true if any 2 of the Det overlap in R. More...
bool	isRPeriodic () const
	Returns true if the Dets are periodic in R. More...
std::vector< double >	RBins () const
	The centers of the Dets. More...
	RBorderFinder (const std::vector< const Det * > &utheDets)
std::vector< double >	RBorders () const
virtual	~RBorderFinder ()

Private Member Functions
int	binIndex (int i) const

Private Attributes
bool	isROverlapping_
bool	isRPeriodic_
int	theNbins
std::vector< double >	theRBins
std::vector< double >	theRBorders

Detailed Description

Find the R binning of a list of detector according to several definitions.

Author

N. Amapane - INFN Torino

Definition at line 22 of file RBorderFinder.h.

Member Typedef Documentation

Det

typedef ForwardDetRing RBorderFinder::Det

Definition at line 24 of file RBorderFinder.h.

DetR

typedef geomsort::ExtractR<Det, float> RBorderFinder::DetR

Definition at line 25 of file RBorderFinder.h.

Constructor & Destructor Documentation

RBorderFinder()

RBorderFinder::RBorderFinder

(

const std::vector< const Det * > &

utheDets

)

Definition at line 3 of file RBorderFinder.cc.

    : theNbins(utheDets.size()), isRPeriodic_(false), isROverlapping_(false) {
  std::vector<const Det*> theDets = utheDets;
  precomputed_value_sort(theDets.begin(), theDets.end(), DetR());
 
  std::vector<ConstReferenceCountingPointer<BoundDisk> > disks(theNbins);
  for (int i = 0; i < theNbins; i++) {
    disks[i] = dynamic_cast<const BoundDisk*>(&(theDets[i]->surface()));
    if (disks[i] == nullptr) {
      throw cms::Exception("UnexpectedState") << "RBorderFinder: implemented for BoundDisks only";
    }
  }
 
  if (theNbins == 1) {    // Trivial case
    isRPeriodic_ = true;  // meaningless in this case
    theRBorders.push_back(disks.front()->innerRadius());
    theRBins.push_back((disks.front()->outerRadius() + disks.front()->innerRadius()));
    //       std::cout << "RBorderFinder:  theNbins " << theNbins << std::endl
    //      << " C: " << theRBins[0]
    //      << " Border: " << theRBorders[0] << std::endl;
  } else {  // More than 1 bin
    double step = (disks.back()->innerRadius() - disks.front()->innerRadius()) / (theNbins - 1);
    std::vector<double> spread;
    std::vector<std::pair<double, double> > REdge;
    REdge.reserve(theNbins);
    theRBorders.reserve(theNbins);
    theRBins.reserve(theNbins);
    spread.reserve(theNbins);
 
    for (int i = 0; i < theNbins; i++) {
      theRBins.push_back((disks[i]->outerRadius() + disks[i]->innerRadius()) / 2.);
      spread.push_back(theRBins.back() - (theRBins[0] + i * step));
      REdge.push_back(std::pair<double, double>(disks[i]->innerRadius(), disks[i]->outerRadius()));
    }
 
    theRBorders.push_back(REdge[0].first);
    for (int i = 1; i < theNbins; i++) {
      // Average borders of previous and next bins
      double br = (REdge[(i - 1)].second + REdge[i].first) / 2.;
      theRBorders.push_back(br);
    }
 
    for (int i = 1; i < theNbins; i++) {
      if (REdge[i].first - REdge[i - 1].second < 0) {
        isROverlapping_ = true;
        break;
      }
    }
 
    double rms = stat_RMS(spread);
    if (rms < 0.01 * step) {
      isRPeriodic_ = true;
    }
  }
 
  //Check that everything is proper
  if ((int)theRBorders.size() != theNbins || (int)theRBins.size() != theNbins)
    throw cms::Exception("UnexpectedState") << "RBorderFinder consistency error";
}

References beamvalidation::br, alignBH_cfg::disks, Exception, first, mps_fire::i, interestingEgammaIsoDetIdsSequence_cff::innerRadius, createfilelist::int, isROverlapping_, isRPeriodic_, interestingEgammaIsoDetIdsSequence_cff::outerRadius, precomputed_value_sort(), SiStripPI::rms, edm::second(), stat_RMS(), theNbins, theRBins, and theRBorders.

~RBorderFinder()

virtual RBorderFinder::~RBorderFinder ( )

inlinevirtual

Definition at line 29 of file RBorderFinder.h.

{};

Member Function Documentation

binIndex()

int RBorderFinder::binIndex ( int  i ) const

inlineprivate

Definition at line 54 of file RBorderFinder.h.

{ return std::min(std::max(i, 0), theNbins - 1); }

References mps_fire::i, SiStripPI::max, min(), and theNbins.

isROverlapping()

bool RBorderFinder::isROverlapping ( ) const

inline

Returns true if any 2 of the Det overlap in R.

Definition at line 35 of file RBorderFinder.h.

{ return isROverlapping_; }

References isROverlapping_.

Referenced by MTDRingForwardLayer::MTDRingForwardLayer(), and MuRingForwardLayer::MuRingForwardLayer().

isRPeriodic()

bool RBorderFinder::isRPeriodic ( ) const

inline

Returns true if the Dets are periodic in R.

Definition at line 32 of file RBorderFinder.h.

{ return isRPeriodic_; }

References isRPeriodic_.

Referenced by MTDRingForwardLayer::MTDRingForwardLayer(), and MuRingForwardLayer::MuRingForwardLayer().

RBins()

std::vector<double> RBorderFinder::RBins ( ) const

inline

The centers of the Dets.

Definition at line 42 of file RBorderFinder.h.

{ return theRBins; }

References theRBins.

Referenced by GeneralBinFinderInR< T >::GeneralBinFinderInR().

RBorders()

std::vector<double> RBorderFinder::RBorders ( ) const

inline

The borders, defined for each det as the middle between its lower edge and the previous Det's upper edge.

Definition at line 39 of file RBorderFinder.h.

{ return theRBorders; }

References theRBorders.

Referenced by GeneralBinFinderInR< T >::GeneralBinFinderInR().

Member Data Documentation

isROverlapping_

bool RBorderFinder::isROverlapping_

private

Definition at line 50 of file RBorderFinder.h.

Referenced by isROverlapping(), and RBorderFinder().

isRPeriodic_

bool RBorderFinder::isRPeriodic_

private

Definition at line 49 of file RBorderFinder.h.

Referenced by isRPeriodic(), and RBorderFinder().

theNbins

int RBorderFinder::theNbins

private

Definition at line 48 of file RBorderFinder.h.

Referenced by binIndex(), and RBorderFinder().

theRBins

std::vector<double> RBorderFinder::theRBins

private

Definition at line 52 of file RBorderFinder.h.

Referenced by RBins(), and RBorderFinder().

theRBorders

std::vector<double> RBorderFinder::theRBorders

private

Definition at line 51 of file RBorderFinder.h.

Referenced by RBorderFinder(), and RBorders().