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 28 of file PhiBorderFinder.h.

Member Typedef Documentation

typedef DetRod PhiBorderFinder::Det

Definition at line 31 of file PhiBorderFinder.h.

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

Definition at line 32 of file PhiBorderFinder.h.

Constructor & Destructor Documentation

PhiBorderFinder::PhiBorderFinder

(

const std::vector< const Det * > &

utheDets

)

Definition at line 3 of file PhiBorderFinder.cc.

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

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

virtual PhiBorderFinder::~PhiBorderFinder ( )

inlinevirtual

Definition at line 37 of file PhiBorderFinder.h.

{};

Member Function Documentation

int PhiBorderFinder::binIndex ( int  i ) const

inlineprivate

Definition at line 69 of file PhiBorderFinder.h.

References createfilelist::int.

Referenced by PhiBorderFinder().

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

bool PhiBorderFinder::isPhiOverlapping ( ) const

inline

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

Definition at line 45 of file PhiBorderFinder.h.

References isPhiOverlapping_.

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

{ return isPhiOverlapping_; }

bool PhiBorderFinder::isPhiPeriodic ( ) const

inline

Returns true if the Dets are periodic in phi.

Definition at line 42 of file PhiBorderFinder.h.

References isPhiPeriodic_.

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

{ return isPhiPeriodic_; }

unsigned int PhiBorderFinder::nBins ( )

inline

Definition at line 39 of file PhiBorderFinder.h.

References theNbins.

{return theNbins;}

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

inline

The centers of the Dets.

Definition at line 52 of file PhiBorderFinder.h.

References thePhiBins.

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

{ return thePhiBins; }

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 49 of file PhiBorderFinder.h.

References thePhiBorders.

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

{ return thePhiBorders; }

double PhiBorderFinder::positiveRange ( double  phi ) const

inlineprivate

Definition at line 64 of file PhiBorderFinder.h.

References Geom::twoPi().

Referenced by PhiBorderFinder().

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

Member Data Documentation

bool PhiBorderFinder::isPhiOverlapping_

private

Definition at line 60 of file PhiBorderFinder.h.

Referenced by isPhiOverlapping(), and PhiBorderFinder().

bool PhiBorderFinder::isPhiPeriodic_

private

Definition at line 59 of file PhiBorderFinder.h.

Referenced by isPhiPeriodic(), and PhiBorderFinder().

unsigned int PhiBorderFinder::theNbins

private

Definition at line 58 of file PhiBorderFinder.h.

Referenced by nBins(), and PhiBorderFinder().

std::vector<double> PhiBorderFinder::thePhiBins

private

Definition at line 62 of file PhiBorderFinder.h.

Referenced by phiBins(), and PhiBorderFinder().

std::vector<double> PhiBorderFinder::thePhiBorders

private

Definition at line 61 of file PhiBorderFinder.h.

Referenced by PhiBorderFinder(), and phiBorders().