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

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

Definition at line 25 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, first, mps_fire::i, isROverlapping_, isRPeriodic_, 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";
 }