OuterDetCompatibility Class Reference

#include <OuterDetCompatibility.h>

Public Types
typedef PixelRecoRange< float >	Range

Public Member Functions
GlobalPoint	center () const
MeasurementEstimator::Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const
MeasurementEstimator::Local2DVector	maximalLocalDisplacement (const GlobalPoint &ts, const BoundPlane &plane) const
bool	operator() (const BoundPlane &plane) const
	OuterDetCompatibility (const BarrelDetLayer *layer, const OuterHitPhiPrediction::Range &phiRange, const Range &rRange, const Range &zRange)
	OuterDetCompatibility (const ForwardDetLayer *layer, const OuterHitPhiPrediction::Range &phiRange, const Range &rRange, const Range &zRange)
const OuterHitPhiPrediction::Range &	phiRange () const
const Range &	rRange () const
const Range &	zRange () const

Private Member Functions
bool	checkPhi (const OuterHitPhiPrediction::Range &detPhiRange) const
bool	checkR (const Range &detRRange) const
bool	checkZ (const Range &detZRange) const
double	loc_dist (double radius, double ts_phi, double range_phi, double cosGamma) const

Private Attributes
bool	barrel
OuterHitPhiPrediction::Range	hitDetPhiRange
Range	hitDetRRange
Range	hitDetZRange
const DetLayer *	theLayer

Detailed Description

check det compatibility by comparistion of det BoundPlane ranges with phi,r,z ranges (given at construction).

Definition at line 16 of file OuterDetCompatibility.h.

Member Typedef Documentation

typedef PixelRecoRange<float> OuterDetCompatibility::Range

Definition at line 18 of file OuterDetCompatibility.h.

Constructor & Destructor Documentation

OuterDetCompatibility::OuterDetCompatibility ( const BarrelDetLayer *  layer, const OuterHitPhiPrediction::Range &  phiRange, const Range &  rRange, const Range &  zRange  )

inline

Definition at line 20 of file OuterDetCompatibility.h.

      : theLayer(layer), barrel(true), hitDetPhiRange(phiRange), hitDetRRange(rRange), hitDetZRange(zRange) {}

OuterDetCompatibility::OuterDetCompatibility ( const ForwardDetLayer *  layer, const OuterHitPhiPrediction::Range &  phiRange, const Range &  rRange, const Range &  zRange  )

inline

Definition at line 26 of file OuterDetCompatibility.h.

      : theLayer(layer), barrel(false), hitDetPhiRange(phiRange), hitDetRRange(rRange), hitDetZRange(zRange) {}

Member Function Documentation

GlobalPoint OuterDetCompatibility::center

(

)

const

Definition at line 30 of file OuterDetCompatibility.cc.

References funct::cos(), alignCSCRings::r, and funct::sin().

                                                {
  float phi = hitDetPhiRange.mean();
  float r = hitDetRRange.mean();
  return GlobalPoint(r * cos(phi), r * sin(phi), hitDetZRange.mean());
}

bool OuterDetCompatibility::checkPhi ( const OuterHitPhiPrediction::Range &  detPhiRange ) const

private

Definition at line 22 of file OuterDetCompatibility.cc.

References Geom::phiLess(), and rangesIntersect().

                                                                                        {
  return rangesIntersect(detPhiRange, hitDetPhiRange, [](auto x, auto y) { return Geom::phiLess(x, y); });
}

bool OuterDetCompatibility::checkR ( const Range &  detRRange ) const

private

Definition at line 26 of file OuterDetCompatibility.cc.

References rangesIntersect().

{ return rangesIntersect(detRRange, hitDetRRange); }

bool OuterDetCompatibility::checkZ ( const Range &  detZRange ) const

private

Definition at line 28 of file OuterDetCompatibility.cc.

References rangesIntersect().

{ return rangesIntersect(detZRange, hitDetZRange); }

double OuterDetCompatibility::loc_dist ( double  radius, double  ts_phi, double  range_phi, double  cosGamma  ) const

private

Definition at line 109 of file OuterDetCompatibility.cc.

References funct::sin(), and mathSSE::sqrt().

                                                                                                            {
  double sinDphi = sin(ts_phi - range_phi);
  double cosDphi = sqrt(1 - sinDphi * sinDphi);
  double sinGamma = sqrt(1 - cosGamma * cosGamma);
  double sinBeta = fabs(cosDphi * cosGamma - sinDphi * sinGamma);
  return radius * fabs(sinDphi) / sinBeta;
}

MeasurementEstimator::Local2DVector OuterDetCompatibility::maximalLocalDisplacement	(	const TrajectoryStateOnSurface &	ts,
		const BoundPlane &	plane
	)		const

MeasurementEstimator::Local2DVector OuterDetCompatibility::maximalLocalDisplacement	(	const GlobalPoint &	ts,
		const BoundPlane &	plane
	)		const

Definition at line 36 of file OuterDetCompatibility.cc.

References Reference_intrackfit_cff::barrel, funct::cos(), HLT_2018_cff::distance, Vector3DBase< T, FrameTag >::dot(), training_settings::idx, M_PI, SiStripPI::max, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), CosmicsPD_Skims::radius, funct::sin(), Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                   {
  float x_loc = 0.;
  float y_loc = 0.;
  if (barrel) {
    double radius = ts.perp();
    GlobalVector planeNorm = plane.normalVector();
    GlobalVector tsDir = GlobalVector(ts.x(), ts.y(), 0.).unit();
    double ts_phi = tsDir.phi();
    if (!hitDetPhiRange.inside(ts_phi)) {
      while (ts_phi >= hitDetPhiRange.max())
        ts_phi -= 2 * M_PI;
      while (ts_phi < hitDetPhiRange.min())
        ts_phi += 2 * M_PI;
      if (!hitDetPhiRange.inside(ts_phi))
        return MeasurementEstimator::Local2DVector(0., 0.);
    }
    double cosGamma = tsDir.dot(planeNorm);

    double dx1 = loc_dist(radius, ts_phi, hitDetPhiRange.min(), cosGamma);
    double dx2 = loc_dist(radius, ts_phi, hitDetPhiRange.max(), cosGamma);

    double ts_z = ts.z();
    double dy1 = ts_z - hitDetZRange.min();
    double dy2 = hitDetZRange.max() - ts_z;

    x_loc = max(dx1, dx2);
    y_loc = max(dy1, dy2);

    // debug only
    /*
    double r1 = dx1 * fabs(cosGamma) / sin(ts_phi-hitDetPhiRange.min());
    double r2 = dx2 * fabs(cosGamma) / sin(hitDetPhiRange.max()-ts_phi);
    GlobalPoint p1( r1* cos(hitDetPhiRange.min()), r1 * sin(hitDetPhiRange.min()), hitDetZRange.min());
    GlobalPoint p2( r2* cos(hitDetPhiRange.max()), r2 * sin(hitDetPhiRange.max()), hitDetZRange.min());
    GlobalPoint p3( r1* cos(hitDetPhiRange.min()), r1 * sin(hitDetPhiRange.min()), hitDetZRange.max());
    GlobalPoint p4( r2* cos(hitDetPhiRange.max()), r2 * sin(hitDetPhiRange.max()), hitDetZRange.max());
    cout << " Local1: " << plane.toLocal(ts-p1) << endl;
    cout << " Local2: " << plane.toLocal(ts-p2) << endl;
    cout << " Local3: " << plane.toLocal(ts-p3) << endl;
    cout << " Local4: " << plane.toLocal(ts-p4) << endl;
*/
  } else {
    LocalPoint ts_loc = plane.toLocal(ts);
    GlobalVector planeNorm = plane.normalVector();

    double x_glob[4], y_glob[4], z_glob[4];
    x_glob[0] = hitDetRRange.min() * cos(hitDetPhiRange.min());
    y_glob[0] = hitDetRRange.min() * sin(hitDetPhiRange.min());
    x_glob[1] = hitDetRRange.max() * cos(hitDetPhiRange.min());
    y_glob[1] = hitDetRRange.max() * sin(hitDetPhiRange.min());
    x_glob[2] = hitDetRRange.min() * cos(hitDetPhiRange.max());
    y_glob[2] = hitDetRRange.min() * sin(hitDetPhiRange.max());
    x_glob[3] = hitDetRRange.max() * cos(hitDetPhiRange.max());
    y_glob[3] = hitDetRRange.max() * sin(hitDetPhiRange.max());

    for (int idx = 0; idx < 4; idx++) {
      double dx_glob = x_glob[idx] - ts.x();
      double dy_glob = y_glob[idx] - ts.y();
      double dz_glob = -(dx_glob * planeNorm.x() + dy_glob * planeNorm.y()) / planeNorm.z();
      z_glob[idx] = dz_glob + ts.z();
    }

    for (int idx = 0; idx < 4; idx++) {
      LocalPoint lp = plane.toLocal(GlobalPoint(x_glob[idx], y_glob[idx], z_glob[idx]));
      x_loc = max(x_loc, fabs(lp.x() - ts_loc.x()));
      y_loc = max(y_loc, fabs(lp.y() - ts_loc.y()));
    }
  }
  MeasurementEstimator::Local2DVector distance(x_loc, y_loc);
  return distance;
}

bool OuterDetCompatibility::operator()

(

const BoundPlane &

plane

)

const

Definition at line 7 of file OuterDetCompatibility.cc.

References Reference_intrackfit_cff::barrel.

                                                                    {
  if (barrel) {
    if (!checkPhi(plane.phiSpan()))
      return false;
    if (!checkZ(plane.zSpan()))
      return false;
  } else {
    if (!checkPhi(plane.phiSpan()))
      return false;
    if (!checkR(plane.rSpan()))
      return false;
  }
  return true;
}

const OuterHitPhiPrediction::Range& OuterDetCompatibility::phiRange ( ) const

inline

Definition at line 40 of file OuterDetCompatibility.h.

{ return hitDetPhiRange; }

const Range& OuterDetCompatibility::rRange ( ) const

inline

Definition at line 41 of file OuterDetCompatibility.h.

{ return hitDetRRange; }

const Range& OuterDetCompatibility::zRange ( ) const

inline

Definition at line 42 of file OuterDetCompatibility.h.

References CosmicsPD_Skims::radius.

{ return hitDetZRange; }

Member Data Documentation

bool OuterDetCompatibility::barrel

private

Definition at line 53 of file OuterDetCompatibility.h.

OuterHitPhiPrediction::Range OuterDetCompatibility::hitDetPhiRange

private

Definition at line 54 of file OuterDetCompatibility.h.

Range OuterDetCompatibility::hitDetRRange

private

Definition at line 55 of file OuterDetCompatibility.h.

Range OuterDetCompatibility::hitDetZRange

private

Definition at line 55 of file OuterDetCompatibility.h.

const DetLayer* OuterDetCompatibility::theLayer

private

Definition at line 52 of file OuterDetCompatibility.h.