PhiBorderFinder Class Reference

#include <PhiBorderFinder.h>

Public Types
typedef DetRod	Det
typedef geomsort::ExtractPhi< Det, float >	DetPhi

Public Member Functions
bool	isPhiOverlapping () const
	Returns true if any 2 of the Det overlap in phi. More...
bool	isPhiPeriodic () const
	Returns true if the Dets are periodic in phi. More...
unsigned int	nBins ()
const std::vector< double > &	phiBins () const
	The centers of the Dets. More...
	PhiBorderFinder (const std::vector< const Det * > &utheDets)
const std::vector< double > &	phiBorders () const
virtual	~PhiBorderFinder ()

Private Member Functions
int	binIndex (int i) const
double	positiveRange (double phi) const

Private Attributes
bool	isPhiOverlapping_
bool	isPhiPeriodic_
unsigned int	theNbins
std::vector< double >	thePhiBins
std::vector< double >	thePhiBorders

Detailed Description

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

Author

N. Amapane - INFN Torino

Definition at line 27 of file PhiBorderFinder.h.

Member Typedef Documentation

Det

typedef DetRod PhiBorderFinder::Det

Definition at line 29 of file PhiBorderFinder.h.

DetPhi

typedef geomsort::ExtractPhi<Det, float> PhiBorderFinder::DetPhi

Definition at line 30 of file PhiBorderFinder.h.

Constructor & Destructor Documentation

PhiBorderFinder()

PhiBorderFinder::PhiBorderFinder

(

const std::vector< const Det * > &

utheDets

)

Definition at line 3 of file PhiBorderFinder.cc.

    : theNbins(utheDets.size()), isPhiPeriodic_(false), isPhiOverlapping_(false) {
  std::vector<const Det*> theDets = utheDets;
  precomputed_value_sort(theDets.begin(), theDets.end(), DetPhi());
 
  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|PhiBorderFinder";
 
  double step = 2. * Geom::pi() / theNbins;
 
  LogTrace(metname) << "RecoMuonDetLayers::PhiBorderFinder "
                    << "step w: " << step << " # of bins: " << theNbins;
 
  std::vector<double> spread(theNbins);
  std::vector<std::pair<double, double> > phiEdge;
  phiEdge.reserve(theNbins);
  thePhiBorders.reserve(theNbins);
  thePhiBins.reserve(theNbins);
  for (unsigned int i = 0; i < theNbins; i++) {
    thePhiBins.push_back(theDets[i]->position().phi());
    spread.push_back(theDets[i]->position().phi() - (theDets[0]->position().phi() + i * step));
 
    LogTrace(metname) << "bin: " << i << " phi bin: " << thePhiBins[i] << " spread: " << spread[i];
 
    ConstReferenceCountingPointer<BoundPlane> plane = dynamic_cast<const BoundPlane*>(&theDets[i]->surface());
    if (plane == nullptr) {
      //FIXME
      throw cms::Exception("UnexpectedState") << ("PhiBorderFinder: det surface is not a BoundPlane");
    }
 
    std::vector<GlobalPoint> dc = BoundingBox().corners(*plane);
 
    float phimin(999.), phimax(-999.);
    for (std::vector<GlobalPoint>::const_iterator pt = dc.begin(); pt != dc.end(); pt++) {
      float phi = (*pt).phi();
      //    float z = pt->z();
      if (phi < phimin)
        phimin = phi;
      if (phi > phimax)
        phimax = phi;
    }
    if (phimin * phimax < 0. &&          //Handle pi border:
        phimax - phimin > Geom::pi()) {  //Assume that the Det is on
      //the shortest side
      std::swap(phimin, phimax);
    }
    phiEdge.push_back(std::pair<double, double>(phimin, phimax));
  }
 
  for (unsigned int i = 0; i < theNbins; i++) {
    // Put the two phi values in the [0,2pi] range
    double firstPhi = positiveRange(phiEdge[i].first);
    double secondPhi = positiveRange(phiEdge[binIndex(i - 1)].second);
 
    // Reformat them in the [-pi,pi] range
    Geom::Phi<double> firstEdge(firstPhi);
    Geom::Phi<double> secondEdge(secondPhi);
 
    // Calculate the mean and format the result in the [-pi,pi] range
    // Get their value in order to perform the mean in the correct way
    Geom::Phi<double> mean((firstEdge.value() + secondEdge.value()) / 2.);
 
    // Special case: at +-pi there is a discontinuity
    if (phiEdge[i].first * phiEdge[binIndex(i - 1)].second < 0 && fabs(firstPhi - secondPhi) < Geom::pi())
      mean = Geom::pi() - mean;
 
    thePhiBorders.push_back(mean);
  }
 
  for (unsigned int i = 0; i < theNbins; i++) {
    if (Geom::Phi<double>(phiEdge[i].first) - Geom::Phi<double>(phiEdge[binIndex(i - 1)].second) < 0) {
      isPhiOverlapping_ = true;
      break;
    }
  }
 
  double rms = stat_RMS(spread);
  if (rms < 0.01 * step) {
    isPhiPeriodic_ = true;
  }
 
  //Check that everything is proper
  if (thePhiBorders.size() != theNbins || thePhiBins.size() != theNbins)
    //FIXME
    throw cms::Exception("UnexpectedState") << "PhiBorderFinder: consistency error";
}

References binIndex(), BoundingBox::corners(), Exception, first, mps_fire::i, isPhiOverlapping_, isPhiPeriodic_, LogTrace, SiStripPI::mean, metname, phi, phimax, phimin, Geom::pi(), position, positiveRange(), precomputed_value_sort(), DiDispStaMuonMonitor_cfi::pt, SiStripPI::rms, edm::second(), stat_RMS(), AlCaHLTBitMon_QueryRunRegistry::string, std::swap(), theNbins, thePhiBins, thePhiBorders, and Geom::Phi< T1, Range >::value().

~PhiBorderFinder()

virtual PhiBorderFinder::~PhiBorderFinder ( )

inlinevirtual

Definition at line 34 of file PhiBorderFinder.h.

{};

Member Function Documentation

binIndex()

int PhiBorderFinder::binIndex ( int  i ) const

inlineprivate

Definition at line 63 of file PhiBorderFinder.h.

                            {
    int ind = i % (int)theNbins;
    return (ind < 0) ? ind + theNbins : ind;
  }

References mps_fire::i, createfilelist::int, and theNbins.

Referenced by PhiBorderFinder().

isPhiOverlapping()

bool PhiBorderFinder::isPhiOverlapping ( ) const

inline

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

Definition at line 42 of file PhiBorderFinder.h.

{ return isPhiOverlapping_; }

References isPhiOverlapping_.

Referenced by MTDTrayBarrelLayer::MTDTrayBarrelLayer(), and MuRodBarrelLayer::MuRodBarrelLayer().

isPhiPeriodic()

bool PhiBorderFinder::isPhiPeriodic ( ) const

inline

Returns true if the Dets are periodic in phi.

Definition at line 39 of file PhiBorderFinder.h.

{ return isPhiPeriodic_; }

References isPhiPeriodic_.

Referenced by MTDTrayBarrelLayer::MTDTrayBarrelLayer(), and MuRodBarrelLayer::MuRodBarrelLayer().

nBins()

unsigned int PhiBorderFinder::nBins ( )

inline

Definition at line 36 of file PhiBorderFinder.h.

{ return theNbins; }

References theNbins.

phiBins()

const std::vector<double>& PhiBorderFinder::phiBins ( ) const

inline

The centers of the Dets.

Definition at line 49 of file PhiBorderFinder.h.

{ return thePhiBins; }

References thePhiBins.

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

phiBorders()

const std::vector<double>& PhiBorderFinder::phiBorders ( ) 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 46 of file PhiBorderFinder.h.

{ return thePhiBorders; }

References thePhiBorders.

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

positiveRange()

double PhiBorderFinder::positiveRange ( double  phi ) const

inlineprivate

Definition at line 61 of file PhiBorderFinder.h.

{ return (phi > 0) ? phi : phi + Geom::twoPi(); }

References phi, and Geom::twoPi().

Referenced by PhiBorderFinder().

Member Data Documentation

isPhiOverlapping_

bool PhiBorderFinder::isPhiOverlapping_

private

Definition at line 57 of file PhiBorderFinder.h.

Referenced by isPhiOverlapping(), and PhiBorderFinder().

isPhiPeriodic_

bool PhiBorderFinder::isPhiPeriodic_

private

Definition at line 56 of file PhiBorderFinder.h.

Referenced by isPhiPeriodic(), and PhiBorderFinder().

theNbins

unsigned int PhiBorderFinder::theNbins

private

Definition at line 55 of file PhiBorderFinder.h.

Referenced by binIndex(), nBins(), and PhiBorderFinder().

thePhiBins

std::vector<double> PhiBorderFinder::thePhiBins

private

Definition at line 59 of file PhiBorderFinder.h.

Referenced by phiBins(), and PhiBorderFinder().

thePhiBorders

std::vector<double> PhiBorderFinder::thePhiBorders

private

Definition at line 58 of file PhiBorderFinder.h.

Referenced by PhiBorderFinder(), and phiBorders().