CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
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 18 of file OuterDetCompatibility.h.

Member Typedef Documentation

typedef PixelRecoRange<float> OuterDetCompatibility::Range

Definition at line 20 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 22 of file OuterDetCompatibility.h.

26  : theLayer(layer), barrel(true),
27  hitDetPhiRange(phiRange), hitDetRRange(rRange), hitDetZRange(zRange) { }
OuterDetCompatibility::zRange
const Range & zRange() const
Definition: OuterDetCompatibility.h:47
OuterDetCompatibility::hitDetPhiRange
OuterHitPhiPrediction::Range hitDetPhiRange
Definition: OuterDetCompatibility.h:60
OuterDetCompatibility::rRange
const Range & rRange() const
Definition: OuterDetCompatibility.h:46
OuterDetCompatibility::OuterDetCompatibility ( const ForwardDetLayer *  layer,
const OuterHitPhiPrediction::Range &  phiRange,
const Range &  rRange,
const Range &  zRange 
)
inline

Definition at line 29 of file OuterDetCompatibility.h.

33  : theLayer(layer), barrel(false),
34  hitDetPhiRange(phiRange), hitDetRRange(rRange), hitDetZRange(zRange) { }
OuterDetCompatibility::zRange
const Range & zRange() const
Definition: OuterDetCompatibility.h:47
OuterDetCompatibility::hitDetPhiRange
OuterHitPhiPrediction::Range hitDetPhiRange
Definition: OuterDetCompatibility.h:60
OuterDetCompatibility::rRange
const Range & rRange() const
Definition: OuterDetCompatibility.h:46

Member Function Documentation

GlobalPoint OuterDetCompatibility::center ( ) const

Definition at line 34 of file OuterDetCompatibility.cc.

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

35 {
36  float phi = hitDetPhiRange.mean();
37  float r = hitDetRRange.mean();
38  return GlobalPoint( r*cos(phi), r*sin(phi), hitDetZRange.mean() );
39 }
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
OuterDetCompatibility::hitDetPhiRange
OuterHitPhiPrediction::Range hitDetPhiRange
Definition: OuterDetCompatibility.h:60
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
PixelRecoRange::mean
T mean() const
Definition: PixelRecoRange.h:32
bool OuterDetCompatibility::checkPhi ( const OuterHitPhiPrediction::Range &  detPhiRange) const
private

Definition at line 20 of file OuterDetCompatibility.cc.

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

22 { return rangesIntersect(detPhiRange, hitDetPhiRange,
23  [](auto x, auto y){ return Geom::phiLess(x,y); }); }
OuterDetCompatibility::hitDetPhiRange
OuterHitPhiPrediction::Range hitDetPhiRange
Definition: OuterDetCompatibility.h:60
rangesIntersect
bool rangesIntersect(const Range &a, const Range &b)
Definition: rangesIntersect.h:14
Geom::phiLess
bool phiLess(float phi1, float phi2)
Definition: VectorUtil.h:23
bool OuterDetCompatibility::checkR ( const Range &  detRRange) const
private

Definition at line 25 of file OuterDetCompatibility.cc.

References rangesIntersect().

27 { return rangesIntersect(detRRange, hitDetRRange); }
rangesIntersect
bool rangesIntersect(const Range &a, const Range &b)
Definition: rangesIntersect.h:14
bool OuterDetCompatibility::checkZ ( const Range &  detZRange) const
private

Definition at line 29 of file OuterDetCompatibility.cc.

References rangesIntersect().

31 { return rangesIntersect(detZRange, hitDetZRange); }
rangesIntersect
bool rangesIntersect(const Range &a, const Range &b)
Definition: rangesIntersect.h:14
double OuterDetCompatibility::loc_dist ( double  radius,
double  ts_phi,
double  range_phi,
double  cosGamma 
) const
private

Definition at line 114 of file OuterDetCompatibility.cc.

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

116 {
117  double sinDphi = sin(ts_phi - range_phi);
118  double cosDphi = sqrt(1-sinDphi*sinDphi);
119  double sinGamma = sqrt(1-cosGamma*cosGamma);
120  double sinBeta = fabs(cosDphi*cosGamma - sinDphi* sinGamma);
121  return radius * fabs(sinDphi) / sinBeta;
122 }
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
TCMET_cfi.radius
radius
Definition: TCMET_cfi.py:32
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 42 of file OuterDetCompatibility.cc.

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

44 {
45  float x_loc = 0.;
46  float y_loc = 0.;
47  if(barrel) {
48  double radius = ts.perp();
49  GlobalVector planeNorm = plane.normalVector();
50  GlobalVector tsDir = GlobalVector( ts.x(), ts.y(),0. ).unit();
51  double ts_phi = tsDir.phi();
52  if (! hitDetPhiRange.inside(ts_phi) ) {
53  while (ts_phi >= hitDetPhiRange.max() ) ts_phi -= 2*M_PI;
54  while (ts_phi < hitDetPhiRange.min() ) ts_phi += 2*M_PI;
55  if (!hitDetPhiRange.inside(ts_phi)) return MeasurementEstimator::Local2DVector(0.,0.);
56  }
57  double cosGamma = tsDir.dot(planeNorm);
58 
59  double dx1 = loc_dist( radius, ts_phi, hitDetPhiRange.min(), cosGamma);
60  double dx2 = loc_dist( radius, ts_phi, hitDetPhiRange.max(), cosGamma);
61 
62  double ts_z = ts.z();
63  double dy1 = ts_z - hitDetZRange.min();
64  double dy2 = hitDetZRange.max() - ts_z;
65 
66  x_loc = max(dx1,dx2);
67  y_loc = max(dy1,dy2);
68 
69  // debug only
70 /*
71  double r1 = dx1 * fabs(cosGamma) / sin(ts_phi-hitDetPhiRange.min());
72  double r2 = dx2 * fabs(cosGamma) / sin(hitDetPhiRange.max()-ts_phi);
73  GlobalPoint p1( r1* cos(hitDetPhiRange.min()), r1 * sin(hitDetPhiRange.min()), hitDetZRange.min());
74  GlobalPoint p2( r2* cos(hitDetPhiRange.max()), r2 * sin(hitDetPhiRange.max()), hitDetZRange.min());
75  GlobalPoint p3( r1* cos(hitDetPhiRange.min()), r1 * sin(hitDetPhiRange.min()), hitDetZRange.max());
76  GlobalPoint p4( r2* cos(hitDetPhiRange.max()), r2 * sin(hitDetPhiRange.max()), hitDetZRange.max());
77  cout << " Local1: " << plane.toLocal(ts-p1) << endl;
78  cout << " Local2: " << plane.toLocal(ts-p2) << endl;
79  cout << " Local3: " << plane.toLocal(ts-p3) << endl;
80  cout << " Local4: " << plane.toLocal(ts-p4) << endl;
81 */
82  }
83  else {
84  LocalPoint ts_loc = plane.toLocal(ts);
85  GlobalVector planeNorm = plane.normalVector();
86 
87  double x_glob[4], y_glob[4], z_glob[4];
88  x_glob[0] = hitDetRRange.min()*cos(hitDetPhiRange.min());
89  y_glob[0] = hitDetRRange.min()*sin(hitDetPhiRange.min());
90  x_glob[1] = hitDetRRange.max()*cos(hitDetPhiRange.min());
91  y_glob[1] = hitDetRRange.max()*sin(hitDetPhiRange.min());
92  x_glob[2] = hitDetRRange.min()*cos(hitDetPhiRange.max());
93  y_glob[2] = hitDetRRange.min()*sin(hitDetPhiRange.max());
94  x_glob[3] = hitDetRRange.max()*cos(hitDetPhiRange.max());
95  y_glob[3] = hitDetRRange.max()*sin(hitDetPhiRange.max());
96 
97  for (int idx = 0; idx < 4; idx++) {
98  double dx_glob = x_glob[idx] - ts.x();
99  double dy_glob = y_glob[idx] - ts.y();
100  double dz_glob = -(dx_glob * planeNorm.x() + dy_glob*planeNorm.y()) / planeNorm.z();
101  z_glob[idx] = dz_glob + ts.z();
102  }
103 
104  for (int idx=0; idx <4; idx++) {
105  LocalPoint lp = plane.toLocal( GlobalPoint( x_glob[idx], y_glob[idx], z_glob[idx]));
106  x_loc = max(x_loc, fabs(lp.x()-ts_loc.x()));
107  y_loc = max(y_loc, fabs(lp.y()-ts_loc.y()));
108  }
109  }
110  MeasurementEstimator::Local2DVector distance(x_loc,y_loc);
111  return distance;
112 }
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
PixelRecoRange::max
T max() const
Definition: PixelRecoRange.h:31
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
OuterDetCompatibility::loc_dist
double loc_dist(double radius, double ts_phi, double range_phi, double cosGamma) const
Definition: OuterDetCompatibility.cc:114
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
Vector3DBase::dot
PreciseFloatType< T, U >::Type dot(const Vector3DBase< U, FrameTag > &v) const
Definition: Vector3DBase.h:107
PixelRecoRange::min
T min() const
Definition: PixelRecoRange.h:30
OuterDetCompatibility::hitDetPhiRange
OuterHitPhiPrediction::Range hitDetPhiRange
Definition: OuterDetCompatibility.h:60
SoftLeptonByDistance_cfi.distance
distance
Definition: SoftLeptonByDistance_cfi.py:6
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
TCMET_cfi.radius
radius
Definition: TCMET_cfi.py:32
training_settings.idx
idx
Definition: training_settings.py:16
PixelRecoRange::inside
bool inside(const T &value) const
Definition: PixelRecoRange.h:39
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
bool OuterDetCompatibility::operator() ( const BoundPlane &  plane) const

Definition at line 7 of file OuterDetCompatibility.cc.

References Reference_intrackfit_cff::barrel.

8 {
9  if (barrel) {
10  if (!checkPhi(plane.phiSpan())) return false;
11  if (!checkZ(plane.zSpan())) return false;
12  } else {
13  if (!checkPhi(plane.phiSpan())) return false;
14  if (!checkR(plane.rSpan())) return false;
15  }
16  return true;
17 }
OuterDetCompatibility::checkZ
bool checkZ(const Range &detZRange) const
Definition: OuterDetCompatibility.cc:29
OuterDetCompatibility::checkPhi
bool checkPhi(const OuterHitPhiPrediction::Range &detPhiRange) const
Definition: OuterDetCompatibility.cc:20
OuterDetCompatibility::checkR
bool checkR(const Range &detRRange) const
Definition: OuterDetCompatibility.cc:25
const OuterHitPhiPrediction::Range& OuterDetCompatibility::phiRange ( ) const
inline

Definition at line 45 of file OuterDetCompatibility.h.

45 {return hitDetPhiRange;}
OuterDetCompatibility::hitDetPhiRange
OuterHitPhiPrediction::Range hitDetPhiRange
Definition: OuterDetCompatibility.h:60
const Range& OuterDetCompatibility::rRange ( ) const
inline

Definition at line 46 of file OuterDetCompatibility.h.

46 { return hitDetRRange; }
const Range& OuterDetCompatibility::zRange ( ) const
inline

Definition at line 47 of file OuterDetCompatibility.h.

References TCMET_cfi::radius.

47 { return hitDetZRange; }

Member Data Documentation

bool OuterDetCompatibility::barrel
private

Definition at line 59 of file OuterDetCompatibility.h.

OuterHitPhiPrediction::Range OuterDetCompatibility::hitDetPhiRange
private

Definition at line 60 of file OuterDetCompatibility.h.

Range OuterDetCompatibility::hitDetRRange
private

Definition at line 61 of file OuterDetCompatibility.h.

Range OuterDetCompatibility::hitDetZRange
private

Definition at line 61 of file OuterDetCompatibility.h.

const DetLayer* OuterDetCompatibility::theLayer
private

Definition at line 58 of file OuterDetCompatibility.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11