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

typedef ForwardDetRing RBorderFinder::Det

Definition at line 25 of file RBorderFinder.h.

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

Definition at line 26 of file RBorderFinder.h.

Constructor & Destructor Documentation

RBorderFinder::RBorderFinder

(

const std::vector< const Det * > &

utheDets

)

Definition at line 3 of file RBorderFinder.cc.

References beamvalidation::br, alignBH_cfg::disks, Exception, plotBeamSpotDB::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.

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

virtual RBorderFinder::~RBorderFinder ( )

inlinevirtual

Definition at line 30 of file RBorderFinder.h.

{};

Member Function Documentation

int RBorderFinder::binIndex ( int  i ) const

inlineprivate

Definition at line 56 of file RBorderFinder.h.

References SiStripPI::max, and min().

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

bool RBorderFinder::isROverlapping ( ) const

inline

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

Definition at line 36 of file RBorderFinder.h.

References isROverlapping_.

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

{ return isROverlapping_; }

bool RBorderFinder::isRPeriodic ( ) const

inline

Returns true if the Dets are periodic in R.

Definition at line 33 of file RBorderFinder.h.

References isRPeriodic_.

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

{ return isRPeriodic_; }

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

inline

The centers of the Dets.

Definition at line 43 of file RBorderFinder.h.

References theRBins.

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

{ return theRBins; }

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 40 of file RBorderFinder.h.

References theRBorders.

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

{ return theRBorders; }

Member Data Documentation

bool RBorderFinder::isROverlapping_

private

Definition at line 52 of file RBorderFinder.h.

Referenced by isROverlapping(), and RBorderFinder().

bool RBorderFinder::isRPeriodic_

private

Definition at line 51 of file RBorderFinder.h.

Referenced by isRPeriodic(), and RBorderFinder().

int RBorderFinder::theNbins

private

Definition at line 50 of file RBorderFinder.h.

Referenced by RBorderFinder().

std::vector<double> RBorderFinder::theRBins

private

Definition at line 54 of file RBorderFinder.h.

Referenced by RBins(), and RBorderFinder().

std::vector<double> RBorderFinder::theRBorders

private

Definition at line 53 of file RBorderFinder.h.

Referenced by RBorderFinder(), and RBorders().