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

typedef DetRod PhiBorderFinder::Det

Definition at line 30 of file PhiBorderFinder.h.

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

Definition at line 31 of file PhiBorderFinder.h.

Constructor & Destructor Documentation

PhiBorderFinder::PhiBorderFinder ( const std::vector< const Det * > &  utheDets )

inline

Definition at line 34 of file PhiBorderFinder.h.

References binIndex(), BoundingBox::corners(), edm::hlt::Exception, first, i, isPhiOverlapping_, isPhiPeriodic_, LogTrace, timingPdfMaker::mean, metname, phi, phimax, phimin, Geom::pi(), position, positiveRange(), precomputed_value_sort(), RecoTauCleanerPlugins::pt, plotscripts::rms(), edm::second(), stat_RMS(), relval_parameters_module::step, 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==0) {
    //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 128 of file PhiBorderFinder.h.

{};

Member Function Documentation

int PhiBorderFinder::binIndex ( int  i ) const

inlineprivate

Definition at line 160 of file PhiBorderFinder.h.

References theNbins.

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

References isPhiOverlapping_.

Referenced by MuRodBarrelLayer::MuRodBarrelLayer().

{ return isPhiOverlapping_; }

bool PhiBorderFinder::isPhiPeriodic ( ) const

inline

Returns true if the Dets are periodic in phi.

Definition at line 133 of file PhiBorderFinder.h.

References isPhiPeriodic_.

Referenced by MuRodBarrelLayer::MuRodBarrelLayer().

{ return isPhiPeriodic_; }

unsigned int PhiBorderFinder::nBins ( )

inline

Definition at line 130 of file PhiBorderFinder.h.

References theNbins.

{return theNbins;}

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

inline

The centers of the Dets.

Definition at line 143 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 140 of file PhiBorderFinder.h.

References thePhiBorders.

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

{ return thePhiBorders; }

double PhiBorderFinder::positiveRange ( double  phi ) const

inlineprivate

Definition at line 155 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 151 of file PhiBorderFinder.h.

Referenced by isPhiOverlapping(), and PhiBorderFinder().

bool PhiBorderFinder::isPhiPeriodic_

private

Definition at line 150 of file PhiBorderFinder.h.

Referenced by isPhiPeriodic(), and PhiBorderFinder().

unsigned int PhiBorderFinder::theNbins

private

Definition at line 149 of file PhiBorderFinder.h.

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

std::vector<double> PhiBorderFinder::thePhiBins

private

Definition at line 153 of file PhiBorderFinder.h.

Referenced by phiBins(), and PhiBorderFinder().

std::vector<double> PhiBorderFinder::thePhiBorders

private

Definition at line 152 of file PhiBorderFinder.h.

Referenced by PhiBorderFinder(), and phiBorders().