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.

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().

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

~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.

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

Referenced by PhiBorderFinder().

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

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.

References isPhiOverlapping_.

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

{ return isPhiOverlapping_; }

isPhiPeriodic()

bool PhiBorderFinder::isPhiPeriodic ( ) const

inline

Returns true if the Dets are periodic in phi.

Definition at line 39 of file PhiBorderFinder.h.

References isPhiPeriodic_.

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

{ return isPhiPeriodic_; }

nBins()

unsigned int PhiBorderFinder::nBins ( )

inline

Definition at line 36 of file PhiBorderFinder.h.

References theNbins.

{ return theNbins; }

phiBins()

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

inline

The centers of the Dets.

Definition at line 49 of file PhiBorderFinder.h.

References thePhiBins.

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

{ return thePhiBins; }

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.

References thePhiBorders.

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

{ return thePhiBorders; }

positiveRange()

double PhiBorderFinder::positiveRange ( double  phi ) const

inlineprivate

Definition at line 61 of file PhiBorderFinder.h.

References phi, and Geom::twoPi().

Referenced by PhiBorderFinder().

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

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().