CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Types | Public Member Functions | Static Public Member Functions | Private Types | Private Member Functions | Private Attributes
RectangularEtaPhiTrackingRegion Class Reference

#include <RectangularEtaPhiTrackingRegion.h>

Inheritance diagram for RectangularEtaPhiTrackingRegion:
TrackingRegion

Public Types

typedef
TkTrackingRegionsMargin< float > 		Margin
 
enum  UseMeasurementTracker { UseMeasurementTracker::kNever = -1, UseMeasurementTracker::kForSiStrips = 0, UseMeasurementTracker::kAlways = 1 }
 
- Public Types inherited from TrackingRegion
using ctfHits = ctfseeding::SeedingLayer::Hits
 
typedef
SeedingLayerSetsHits::ConstRecHitPointer 		Hit
 
typedef SeedingLayerSetsHits::Hits Hits
 
typedef PixelRecoRange< float > Range
 

Public Member Functions

virtual HitRZCompatibility * checkRZ (const DetLayer *layer, const Hit &outerHit, const edm::EventSetup &iSetup, const DetLayer *outerlayer=0, float lr=0, float gz=0, float dr=0, float dz=0) const
 
virtual
RectangularEtaPhiTrackingRegion * 		clone () const
 
const Range & etaRange () const
 allowed eta range [eta_min, eta_max] interval More...
 
virtual TrackingRegion::Hits hits (const edm::Event &ev, const edm::EventSetup &es, const SeedingLayerSetsHits::SeedingLayer &layer) const override
 get hits from layer compatible with region constraints More...
 
bool isPrecise () const
 is precise error calculation switched on More...
 
virtual std::string name () const
 
RectangularEtaPhiTrackingRegion & operator= (RectangularEtaPhiTrackingRegion const &)=delete
 
RectangularEtaPhiTrackingRegion & operator= (RectangularEtaPhiTrackingRegion &&)=default
 
const Margin & phiMargin () const
 
virtual std::string print () const
 
 RectangularEtaPhiTrackingRegion (RectangularEtaPhiTrackingRegion const &rh)
 
 RectangularEtaPhiTrackingRegion (RectangularEtaPhiTrackingRegion &&)=default
 
 RectangularEtaPhiTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, UseMeasurementTracker whereToUseMeasurementTracker=UseMeasurementTracker::kNever, bool precise=true, const MeasurementTrackerEvent *measurementTracker=0, bool etaPhiRegion=false)
 
 RectangularEtaPhiTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, Margin etaMargin, Margin phiMargin, UseMeasurementTracker whereToUseMeasurementTracker=UseMeasurementTracker::kNever, bool precise=true, const MeasurementTrackerEvent *measurementTracker=0, bool etaPhiRegion=false)
 
 RectangularEtaPhiTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, Range invPtRange, float rVertex, float zVertex, Margin etaMargin, Margin phiMargin, UseMeasurementTracker whereToUseMeasurementTracker=UseMeasurementTracker::kNever, bool precise=true, const MeasurementTrackerEvent *measurementTracker=0, bool etaPhiRegion=false)
 
- Public Member Functions inherited from TrackingRegion
GlobalVector const & direction () const
 the direction around which region is constructed More...
 
Range invPtRange () const
 inverse pt range More...
 
GlobalPoint const & origin () const
 
float originRBound () const
 bounds the particle vertex in the transverse plane More...
 
float originZBound () const
 bounds the particle vertex in the longitudinal plane More...
 
float phiDirection () const
 
float ptMin () const
 minimal pt of interest More...
 
TrackingRegion * restrictedRegion (const GlobalPoint &originPos, const float &originRBound, const float &originZBound) const
 clone region with new vertex position More...
 
 TrackingRegion (const GlobalVector &direction, const GlobalPoint &originPos, const Range &invPtRange, const float &originRBound, const float &originZBound)
 
GlobalVector const & unitDirection () const
 
virtual ~TrackingRegion ()
 

Static Public Member Functions

static UseMeasurementTracker doubleToUseMeasurementTracker (double value)
 
static UseMeasurementTracker intToUseMeasurementTracker (int value)
 

Private Types

using cacheHitPointer = mayown_ptr< BaseTrackerRecHit >
 
using cacheHits = std::vector< cacheHitPointer >
 

Private Member Functions

HitRZCompatibility * checkRZOld (const DetLayer *layer, const TrackingRecHit *outerHit, const edm::EventSetup &iSetup) const
 
OuterEstimator * estimator (const BarrelDetLayer *layer, const edm::EventSetup &iSetup) const
 
OuterEstimator * estimator (const ForwardDetLayer *layer, const edm::EventSetup &iSetup) const
 
void initEtaRange (const GlobalVector &dir, const Margin &margin)
 
OuterHitPhiPrediction phiWindow (const edm::EventSetup &iSetup) const
 
HitRZConstraint rzConstraint () const
 

Private Attributes

cacheHits cache
 
Range theEtaRange
 
Range theLambdaRange
 
float theMeanLambda
 
const MeasurementTrackerEvent * theMeasurementTracker
 
const UseMeasurementTracker theMeasurementTrackerUsage
 
Margin thePhiMargin
 
bool thePrecise
 
bool theUseEtaPhi
 

Detailed Description

A concrete implementation of TrackingRegion. Apart of vertex constraint from TrackingRegion in this implementation the region of interest is further constrainted in phi and eta around the direction of the region

Definition at line 29 of file RectangularEtaPhiTrackingRegion.h.

Member Typedef Documentation

using RectangularEtaPhiTrackingRegion::cacheHitPointer = mayown_ptr<BaseTrackerRecHit>
private

Definition at line 203 of file RectangularEtaPhiTrackingRegion.h.

using RectangularEtaPhiTrackingRegion::cacheHits = std::vector<cacheHitPointer>
private

Definition at line 204 of file RectangularEtaPhiTrackingRegion.h.

typedef TkTrackingRegionsMargin<float> RectangularEtaPhiTrackingRegion::Margin

Definition at line 64 of file RectangularEtaPhiTrackingRegion.h.

Member Enumeration Documentation

enum RectangularEtaPhiTrackingRegion::UseMeasurementTracker
strong
Enumerator
kNever 
kForSiStrips 
kAlways 

Definition at line 31 of file RectangularEtaPhiTrackingRegion.h.

31  {
32  kNever = -1,
33  kForSiStrips = 0,
34  kAlways = 1
35  };

Constructor & Destructor Documentation

RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion ( RectangularEtaPhiTrackingRegion const &  rh)
inline

Definition at line 49 of file RectangularEtaPhiTrackingRegion.h.

Referenced by clone().

49  :
50  TrackingRegion(rh),
51  theEtaRange(rh.theEtaRange),
52  theLambdaRange(rh.theLambdaRange),
53  thePhiMargin(rh.thePhiMargin),
54  theMeanLambda(rh.theMeanLambda),
55  theMeasurementTrackerUsage(rh.theMeasurementTrackerUsage),
56  thePrecise(rh.thePrecise),
57  theUseEtaPhi(rh.theUseEtaPhi),
58  theMeasurementTracker(rh.theMeasurementTracker) {}
RectangularEtaPhiTrackingRegion::theMeasurementTrackerUsage
const UseMeasurementTracker theMeasurementTrackerUsage
Definition: RectangularEtaPhiTrackingRegion.h:196
TrackingRegion::TrackingRegion
TrackingRegion(const GlobalVector &direction, const GlobalPoint &originPos, const Range &invPtRange, const float &originRBound, const float &originZBound)
Definition: TrackingRegion.h:52
RectangularEtaPhiTrackingRegion::theMeasurementTracker
const MeasurementTrackerEvent * theMeasurementTracker
Definition: RectangularEtaPhiTrackingRegion.h:199
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion ( RectangularEtaPhiTrackingRegion &&  )
default
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion ( const GlobalVector &  dir,
const GlobalPoint &  vertexPos,
float  ptMin,
float  rVertex,
float  zVertex,
float  deltaEta,
float  deltaPhi,
UseMeasurementTracker  whereToUseMeasurementTracker = UseMeasurementTracker::kNever,
bool  precise = true,
const MeasurementTrackerEvent *  measurementTracker = 0,
bool  etaPhiRegion = false 
)
inline

constructor (symmetric eta and phi margins).
dir - the direction around which region is constructed
the initial direction of the momentum of the particle should be in the range
phi of dir +- deltaPhi
eta of dir +- deltaEta

vertexPos - the position of the vertex (origin) of the of the region.
It is a centre of cylinder constraind with rVertex, zVertex. The track of the particle should cross the cylinder
WARNING: in the current implementaion the vertexPos is supposed to be placed on the beam line, i.e. to be of the form (0,0,float)

ptMin - minimal pt of interest
rVertex - radius of the cylinder around beam line where the tracks of interest should point to.
zVertex - half height of the cylinder around the beam line where the tracks of interest should point to.
deltaEta - allowed deviation of the initial direction of particle in eta in respect to direction of the region
deltaPhi - allowed deviation of the initial direction of particle in phi in respect to direction of the region whereToUseMeasurementTracker: 1=everywhere, 0=outside pixles, -1=nowhere

Definition at line 91 of file RectangularEtaPhiTrackingRegion.h.

99  : RectangularEtaPhiTrackingRegion(dir, vertexPos, Range( -1/ptMin, 1/ptMin), rVertex, zVertex,
100  Margin(std::abs(deltaEta), std::abs(deltaEta)),
101  Margin(std::abs(deltaPhi), std::abs(deltaPhi)),
102  whereToUseMeasurementTracker, precise,
103  measurementTracker, etaPhiRegion)
104  {}
TrackingRegion::Range
PixelRecoRange< float > Range
Definition: TrackingRegion.h:45
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
SiPixelRawToDigiRegional_cfi.deltaPhi
tuple deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9
HLTFastRecoForTau_cff.deltaEta
tuple deltaEta
Definition: HLTFastRecoForTau_cff.py:35
TrackingRegion::ptMin
float ptMin() const
minimal pt of interest
Definition: TrackingRegion.h:83
HLTFastRecoForTau_cff.precise
tuple precise
Definition: HLTFastRecoForTau_cff.py:25
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion
RectangularEtaPhiTrackingRegion(RectangularEtaPhiTrackingRegion const &rh)
Definition: RectangularEtaPhiTrackingRegion.h:49
RectangularEtaPhiTrackingRegion::Margin
TkTrackingRegionsMargin< float > Margin
Definition: RectangularEtaPhiTrackingRegion.h:64
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion ( const GlobalVector &  dir,
const GlobalPoint &  vertexPos,
float  ptMin,
float  rVertex,
float  zVertex,
Margin  etaMargin,
Margin  phiMargin,
UseMeasurementTracker  whereToUseMeasurementTracker = UseMeasurementTracker::kNever,
bool  precise = true,
const MeasurementTrackerEvent *  measurementTracker = 0,
bool  etaPhiRegion = false 
)
inline

constructor (asymmetrinc eta and phi margins).
non equal left-right eta and phi bounds around direction are possible. The ranges are defined using Margin. the meaning of other arguments is the same as in the case of symmetring bounds to direction of the region.

Definition at line 112 of file RectangularEtaPhiTrackingRegion.h.

121  : RectangularEtaPhiTrackingRegion(dir, vertexPos, Range( -1/ptMin, 1/ptMin), rVertex, zVertex,
122  etaMargin, phiMargin,
123  whereToUseMeasurementTracker, precise,
124  measurementTracker, etaPhiRegion)
125  {}
TrackingRegion::Range
PixelRecoRange< float > Range
Definition: TrackingRegion.h:45
TrackingRegion::ptMin
float ptMin() const
minimal pt of interest
Definition: TrackingRegion.h:83
HLTFastRecoForTau_cff.precise
tuple precise
Definition: HLTFastRecoForTau_cff.py:25
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion
RectangularEtaPhiTrackingRegion(RectangularEtaPhiTrackingRegion const &rh)
Definition: RectangularEtaPhiTrackingRegion.h:49
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion ( const GlobalVector &  dir,
const GlobalPoint &  vertexPos,
Range  invPtRange,
float  rVertex,
float  zVertex,
Margin  etaMargin,
Margin  phiMargin,
UseMeasurementTracker  whereToUseMeasurementTracker = UseMeasurementTracker::kNever,
bool  precise = true,
const MeasurementTrackerEvent *  measurementTracker = 0,
bool  etaPhiRegion = false 
)
inline

constructor (explicit pt range, asymmetrinc eta and phi margins).
the meaning of other arguments is the same as in the case of symmetring bounds to direction of the region.

Definition at line 131 of file RectangularEtaPhiTrackingRegion.h.

References initEtaRange().

141  : TrackingRegionBase( dir, vertexPos, invPtRange, rVertex, zVertex),
142  thePhiMargin( phiMargin), theMeasurementTrackerUsage(whereToUseMeasurementTracker), thePrecise(precise),theUseEtaPhi(etaPhiRegion),
143  theMeasurementTracker(measurementTracker)
144  { initEtaRange(dir, etaMargin); }
RectangularEtaPhiTrackingRegion::theMeasurementTrackerUsage
const UseMeasurementTracker theMeasurementTrackerUsage
Definition: RectangularEtaPhiTrackingRegion.h:196
RectangularEtaPhiTrackingRegion::theMeasurementTracker
const MeasurementTrackerEvent * theMeasurementTracker
Definition: RectangularEtaPhiTrackingRegion.h:199
RectangularEtaPhiTrackingRegion::initEtaRange
void initEtaRange(const GlobalVector &dir, const Margin &margin)
Definition: RectangularEtaPhiTrackingRegion.cc:46
TrackingRegion::invPtRange
Range invPtRange() const
inverse pt range
Definition: TrackingRegion.h:86
HLTFastRecoForTau_cff.precise
tuple precise
Definition: HLTFastRecoForTau_cff.py:25
TrackingRegionBase
TrackingRegion TrackingRegionBase
Definition: TrackingRegion.h:139

Member Function Documentation

virtual HitRZCompatibility* RectangularEtaPhiTrackingRegion::checkRZ ( const DetLayer *  layer,
const Hit &  outerHit,
const edm::EventSetup &  iSetup,
const DetLayer *  outerlayer = 0,
float  lr = 0,
float  gz = 0,
float  dr = 0,
float  dz = 0 
) const
inlinevirtual

utility to check eta/theta hit compatibility with region constraints and outer hit constraint

Implements TrackingRegion.

Definition at line 162 of file RectangularEtaPhiTrackingRegion.h.

References checkRZOld().

167  { return checkRZOld(layer,outerHit->hit(),iSetup); }
RectangularEtaPhiTrackingRegion::checkRZOld
HitRZCompatibility * checkRZOld(const DetLayer *layer, const TrackingRecHit *outerHit, const edm::EventSetup &iSetup) const
Definition: RectangularEtaPhiTrackingRegion.cc:54
HitRZCompatibility * RectangularEtaPhiTrackingRegion::checkRZOld ( const DetLayer *  layer,
const TrackingRecHit *  outerHit,
const edm::EventSetup &  iSetup 
) const
private

Definition at line 54 of file RectangularEtaPhiTrackingRegion.cc.

References funct::abs(), GeomDetEnumerators::barrel, SimpleLineRZ::cotLine(), f, TrackingRecHit::geographicalId(), edm::EventSetup::get(), GeomDetEnumerators::isBarrel(), TrackingRecHit::localPosition(), DetLayer::location(), bookConverter::max, min(), SurfaceOrientation::outer, PtMinSelector_cfg::ptMin, funct::sqr(), mathSSE::sqrt(), patCandidatesForDimuonsSequences_cff::tracker, vtxMean(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PixelRecoPointRZ::z(), PV3DBase< T, PVType, FrameType >::z(), and PixelRecoLineRZ::zAtR().

Referenced by checkRZ().

55 {
56  edm::ESHandle<TrackerGeometry> tracker;
57  iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
58 
59  bool isBarrel = (layer->location() == GeomDetEnumerators::barrel);
60  GlobalPoint ohit = tracker->idToDet(outerHit->geographicalId())->surface().toGlobal(outerHit->localPosition());
61  float outerred_r = sqrt( sqr(ohit.x()-origin().x())+sqr(ohit.y()-origin().y()) );
62  //PixelRecoPointRZ outer(ohit.perp(), ohit.z());
63  PixelRecoPointRZ outer(outerred_r, ohit.z());
64 
65  float zMinOrigin = origin().z() - originZBound();
66  float zMaxOrigin = origin().z() + originZBound();
67 
68  if (!thePrecise) {
69  float vcotMin = (outer.z() > zMaxOrigin) ?
70  (outer.z()-zMaxOrigin)/(outer.r()+originRBound())
71  : (outer.z()-zMaxOrigin)/(outer.r()-originRBound());
72  float vcotMax = (outer.z() > zMinOrigin) ?
73  (outer.z()-zMinOrigin)/(outer.r()-originRBound())
74  : (outer.z()-zMinOrigin)/(outer.r()+originRBound());
75  float cotRight = std::max(vcotMin,theLambdaRange.min());
76  float cotLeft = std::min(vcotMax, theLambdaRange.max());
77  return new HitEtaCheck( isBarrel, outer, cotLeft, cotRight);
78  }
79  float hitZErr = 0.;
80  float hitRErr = 0.;
81 
82  PixelRecoPointRZ outerL, outerR;
83  if (layer->location() == GeomDetEnumerators::barrel) {
84  outerL = PixelRecoPointRZ(outer.r(), outer.z()-hitZErr);
85  outerR = PixelRecoPointRZ(outer.r(), outer.z()+hitZErr);
86  } else if (outer.z() > 0) {
87  outerL = PixelRecoPointRZ(outer.r()+hitRErr, outer.z());
88  outerR = PixelRecoPointRZ(outer.r()-hitRErr, outer.z());
89  } else {
90  outerL = PixelRecoPointRZ(outer.r()-hitRErr, outer.z());
91  outerR = PixelRecoPointRZ(outer.r()+hitRErr, outer.z());
92  }
93  //CHECK
94  MultipleScatteringParametrisation oSigma(layer,iSetup);
95  float cotThetaOuter = theMeanLambda;
96  float sinThetaOuter = 1/std::sqrt(1+sqr(cotThetaOuter));
97  float outerZscatt = 3.f*oSigma(ptMin(),cotThetaOuter) / sinThetaOuter;
98 
99  PixelRecoLineRZ boundL(outerL, theLambdaRange.max());
100  PixelRecoLineRZ boundR(outerR, theLambdaRange.min());
101  float zMinLine = boundL.zAtR(0.)-outerZscatt;
102  float zMaxLine = boundR.zAtR(0.)+outerZscatt;
103  PixelRecoPointRZ vtxL(0.,max(zMinLine, zMinOrigin));
104  PixelRecoPointRZ vtxR(0.,min(zMaxLine, zMaxOrigin));
105  PixelRecoPointRZ vtxMean(0.,(vtxL.z()+vtxR.z())*0.5f);
106  //CHECK
107  MultipleScatteringParametrisation iSigma(layer,iSetup);
108  float innerScatt = 3.f * iSigma(ptMin(),vtxMean, outer);
109 
110  SimpleLineRZ leftLine( vtxL, outerL);
111  SimpleLineRZ rightLine( vtxR, outerR);
112 
113  HitRZConstraint rzConstraint(leftLine, rightLine);
114  float cotTheta = std::abs(leftLine.cotLine()+rightLine.cotLine())*0.5f;
115 
116 // float bendR = longitudinalBendingCorrection(outer.r(),ptMin());
117 
118  // std::cout << "RectangularEtaPhiTrackingRegion " << outer.r()<<','<< outer.z() << " " << innerScatt << " " << cotTheta << " " << hitZErr << std::endl;
119 
120  if (isBarrel) {
121  float sinTheta = 1/std::sqrt(1+sqr(cotTheta));
122  float corrZ = innerScatt/sinTheta + hitZErr;
123  return new HitZCheck(rzConstraint, HitZCheck::Margin(corrZ,corrZ));
124  } else {
125  float cosTheta = 1/std::sqrt(1+sqr(1/cotTheta));
126  float corrR = innerScatt/cosTheta + hitRErr;
127  return new HitRCheck( rzConstraint, HitRCheck::Margin(corrR,corrR));
128  }
129 }
TrackingRegion::originRBound
float originRBound() const
bounds the particle vertex in the transverse plane
Definition: TrackingRegion.h:77
TrackingRegion::origin
GlobalPoint const & origin() const
Definition: TrackingRegion.h:74
PixelRecoRange::max
T max() const
Definition: PixelRecoRange.h:31
DetLayer::location
virtual Location location() const =0
Which part of the detector (barrel, endcap)
GeomDetEnumerators::isBarrel
bool isBarrel(GeomDetEnumerators::SubDetector m)
Definition: GeomDetEnumerators.cc:32
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
PixelRecoRange::min
T min() const
Definition: PixelRecoRange.h:30
PixelRecoLineRZ::zAtR
float zAtR(float r) const
Definition: PixelRecoLineRZ.h:41
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
vtxMean
static GlobalPoint vtxMean(const GlobalPoint &p1, const GlobalError &e1, const GlobalPoint &p2, const GlobalError &e2)
Definition: GhostTrackVertexFinder.cc:119
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
min
T min(T a, T b)
Definition: MathUtil.h:58
TrackingRegion::originZBound
float originZBound() const
bounds the particle vertex in the longitudinal plane
Definition: TrackingRegion.h:80
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
TrackingRegion::ptMin
float ptMin() const
minimal pt of interest
Definition: TrackingRegion.h:83
funct::sqr
Square< F >::type sqr(const F &f)
Definition: Square.h:13
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:116
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
TrackingRecHit::localPosition
virtual LocalPoint localPosition() const =0
virtual RectangularEtaPhiTrackingRegion* RectangularEtaPhiTrackingRegion::clone ( void  ) const
inlinevirtual

Implements TrackingRegion.

Definition at line 169 of file RectangularEtaPhiTrackingRegion.h.

References RectangularEtaPhiTrackingRegion().

169  {
170  return new RectangularEtaPhiTrackingRegion(*this);
171  }
RectangularEtaPhiTrackingRegion::RectangularEtaPhiTrackingRegion
RectangularEtaPhiTrackingRegion(RectangularEtaPhiTrackingRegion const &rh)
Definition: RectangularEtaPhiTrackingRegion.h:49
static UseMeasurementTracker RectangularEtaPhiTrackingRegion::doubleToUseMeasurementTracker ( double  value)
inlinestatic

Definition at line 42 of file RectangularEtaPhiTrackingRegion.h.

References kAlways, kForSiStrips, and kNever.

Referenced by MuonTrackingRegionBuilder::build(), CandidateSeededTrackingRegionsProducer::CandidateSeededTrackingRegionsProducer(), L3MumuTrackingRegion::L3MumuTrackingRegion(), PointSeededTrackingRegionsProducer::PointSeededTrackingRegionsProducer(), TauRegionalPixelSeedGenerator::TauRegionalPixelSeedGenerator(), and TrackingRegionsFromBeamSpotAndL2Tau::TrackingRegionsFromBeamSpotAndL2Tau().

42  {
43  // mimic the old behaviour
44  if(value > 0.5) return UseMeasurementTracker::kAlways;
45  if(value > -0.5) return UseMeasurementTracker::kForSiStrips;
46  return UseMeasurementTracker::kNever;
47  }
OuterEstimator * RectangularEtaPhiTrackingRegion::estimator ( const BarrelDetLayer *  layer,
const edm::EventSetup &  iSetup 
) const
private

Definition at line 132 of file RectangularEtaPhiTrackingRegion.cc.

References funct::abs(), Surface::bounds(), PixelRecoRange< T >::empty(), reco::helper::VirtualJetProducerHelper::intersection(), PixelRecoRange< T >::intersection(), Bounds::length(), PixelRecoUtilities::longitudinalBendingCorrection(), PixelRecoRange< T >::max(), bookConverter::max, PixelRecoRange< T >::mean(), PixelRecoRange< T >::min(), min(), GeometricSearchDet::position(), PtMinSelector_cfg::ptMin, CosmicsPD_Skims::radius, HitZCheck::range(), OuterHitPhiPrediction::setTolerance(), HitZCheck::setTolerance(), BarrelDetLayer::specificSurface(), funct::sqr(), mathSSE::sqrt(), BarrelDetLayer::surface(), Bounds::thickness(), w2, and PV3DBase< T, PVType, FrameType >::z().

133 {
134 
135  // det dimensions
136  float halfLength = 0.5f*layer->surface().bounds().length();
137  float halfThickness = 0.5f*layer->surface().bounds().thickness();
138  float z0 = layer->position().z();
139  float radius = layer->specificSurface().radius();
140 
141  // det ranges
142  Range detRWindow (radius-halfThickness, radius+halfThickness);
143  Range detZWindow(z0-halfLength,z0+halfLength);
144 
145  // z prediction, skip if not intersection
146  HitZCheck zPrediction(rzConstraint());
147  Range hitZWindow = zPrediction.range(detRWindow.min()).
148  intersection(detZWindow);
149  if (hitZWindow.empty()) return 0;
150 
151  // phi prediction
152  OuterHitPhiPrediction phiPrediction = phiWindow(iSetup);
153 
154  //
155  // optional corrections for tolerance (mult.scatt, error, bending)
156  //
157  OuterHitPhiPrediction::Range phiRange;
158  if (thePrecise) {
159  float cotTheta = (hitZWindow.mean()-origin().z()) / radius;
160  float sinTheta = 1/std::sqrt(1+sqr(cotTheta));
161  MultipleScatteringParametrisation msSigma(layer,iSetup);
162  float scatt = 3.f * msSigma(ptMin(), cotTheta);
163  float bendR = longitudinalBendingCorrection(radius,ptMin(),iSetup);
164 
165  float hitErrRPhi = 0.;
166  float hitErrZ = 0.;
167  float corrPhi = (scatt+ hitErrRPhi)/radius;
168  float corrZ = scatt/sinTheta + bendR*std::abs(cotTheta) + hitErrZ;
169 
170  phiPrediction.setTolerance(OuterHitPhiPrediction::Margin(corrPhi,corrPhi));
171  zPrediction.setTolerance(HitZCheck::Margin(corrZ,corrZ));
172 
173  //
174  // hit ranges in det
175  //
176  OuterHitPhiPrediction::Range phi1 = phiPrediction(detRWindow.min());
177  OuterHitPhiPrediction::Range phi2 = phiPrediction(detRWindow.max());
178  phiRange = Range( std::min(phi1.min(),phi2.min()), std::max(phi1.max(),phi2.max()));
179  Range w1 = zPrediction.range(detRWindow.min());
180  Range w2 = zPrediction.range(detRWindow.max());
181  hitZWindow = Range(std::min(w1.min(),w2.min()), std::max(w1.max(),w2.max())).intersection(detZWindow);
182  }
183  else {
184  phiRange = phiPrediction(detRWindow.mean());
185  }
186 
187  return new OuterEstimator(
188  OuterDetCompatibility( layer, phiRange, detRWindow, hitZWindow),
189  OuterHitCompatibility( phiPrediction, zPrediction ),
190  iSetup);
191 }
Bounds::length
virtual float length() const =0
w2
common ppss p3p6s2 common epss epspn46 common const1 w2
Definition: inclppp.h:1
TrackingRegion::origin
GlobalPoint const & origin() const
Definition: TrackingRegion.h:74
PixelRecoRange::max
T max() const
Definition: PixelRecoRange.h:31
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:128
PixelRecoRange::min
T min() const
Definition: PixelRecoRange.h:30
Bounds::thickness
virtual float thickness() const =0
TrackingRegion::Range
PixelRecoRange< float > Range
Definition: TrackingRegion.h:45
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
min
T min(T a, T b)
Definition: MathUtil.h:58
RectangularEtaPhiTrackingRegion::phiWindow
OuterHitPhiPrediction phiWindow(const edm::EventSetup &iSetup) const
Definition: RectangularEtaPhiTrackingRegion.cc:252
BarrelDetLayer::surface
virtual const BoundSurface & surface() const
GeometricSearchDet interface.
Definition: BarrelDetLayer.h:32
PixelRecoUtilities::longitudinalBendingCorrection
double longitudinalBendingCorrection(double radius, double pt, const edm::EventSetup &iSetup)
Definition: PixelRecoUtilities.h:51
GeometricSearchDet::position
virtual const Surface::PositionType & position() const
Returns position of the surface.
Definition: GeometricSearchDet.h:33
BarrelDetLayer::specificSurface
virtual const BoundCylinder & specificSurface() const
Extension of the interface.
Definition: BarrelDetLayer.h:43
TrackingRegion::ptMin
float ptMin() const
minimal pt of interest
Definition: TrackingRegion.h:83
PixelRecoRange::intersection
PixelRecoRange< T > intersection(const PixelRecoRange< T > &r) const
Definition: PixelRecoRange.h:44
funct::sqr
Square< F >::type sqr(const F &f)
Definition: Square.h:13
OuterHitPhiPrediction::setTolerance
void setTolerance(const Margin &tolerance)
Definition: OuterHitPhiPrediction.h:25
CosmicsPD_Skims.radius
tuple radius
Definition: CosmicsPD_Skims.py:135
OuterEstimator * RectangularEtaPhiTrackingRegion::estimator ( const ForwardDetLayer *  layer,
const edm::EventSetup &  iSetup 
) const
private

Definition at line 194 of file RectangularEtaPhiTrackingRegion.cc.

References funct::abs(), Surface::bounds(), PixelRecoRange< T >::empty(), reco::helper::VirtualJetProducerHelper::intersection(), PixelRecoRange< T >::intersection(), PixelRecoUtilities::longitudinalBendingCorrection(), PixelRecoRange< T >::max(), PixelRecoRange< T >::mean(), PixelRecoRange< T >::min(), GeometricSearchDet::position(), PtMinSelector_cfg::ptMin, OuterHitPhiPrediction::setTolerance(), ForwardDetLayer::specificSurface(), funct::sqr(), mathSSE::sqrt(), ForwardDetLayer::surface(), Bounds::thickness(), w2, and PV3DBase< T, PVType, FrameType >::z().

195 {
196 
197  // det dimensions, ranges
198  float halfThickness = 0.5f*layer->surface().bounds().thickness();
199  float zLayer = layer->position().z() ;
200  Range detZWindow( zLayer-halfThickness, zLayer+halfThickness);
201  Range detRWindow( layer->specificSurface().innerRadius(),
202  layer->specificSurface().outerRadius());
203 
204  // r prediction, skip if not intersection
205  HitRCheck rPrediction(rzConstraint());
206  Range hitRWindow = rPrediction.range(zLayer).intersection(detRWindow);
207  if (hitRWindow.empty()) return 0;
208 
209  // phi prediction
210  OuterHitPhiPrediction phiPrediction = phiWindow(iSetup);
211  OuterHitPhiPrediction::Range phiRange = phiPrediction(detRWindow.max());
212 
213  //
214  // optional corrections for tolerance (mult.scatt, error, bending)
215  //
216  if (thePrecise) {
217  float cotTheta = (detZWindow.mean()-origin().z())/hitRWindow.mean();
218  float cosTheta = cotTheta/std::sqrt(1+sqr(cotTheta));
219  MultipleScatteringParametrisation msSigma(layer,iSetup);
220  float scatt = 3.f * msSigma(ptMin(),cotTheta);
221  float bendR = longitudinalBendingCorrection(hitRWindow.max(),ptMin(),iSetup);
222  float hitErrRPhi = 0.;
223  float hitErrR = 0.;
224  float corrPhi = (scatt+hitErrRPhi)/detRWindow.min();
225  float corrR = scatt/std::abs(cosTheta) + bendR + hitErrR;
226 
227  phiPrediction.setTolerance(OuterHitPhiPrediction::Margin(corrPhi,corrPhi));
228  rPrediction.setTolerance(HitRCheck::Margin(corrR,corrR));
229 
230  //
231  // hit ranges in det
232  //
233  Range w1,w2;
234  if (zLayer > 0) {
235  w1 = rPrediction.range(detZWindow.min());
236  w2 = rPrediction.range(detZWindow.max());
237  } else {
238  w1 = rPrediction.range(detZWindow.max());
239  w2 = rPrediction.range(detZWindow.min());
240  }
241  hitRWindow = Range(w1.min(),w2.max()).intersection(detRWindow);
242  }
243 
244  return new OuterEstimator(
245  OuterDetCompatibility( layer, phiRange, hitRWindow, detZWindow),
246  OuterHitCompatibility( phiPrediction, rPrediction),iSetup );
247 }
w2
common ppss p3p6s2 common epss epspn46 common const1 w2
Definition: inclppp.h:1
TrackingRegion::origin
GlobalPoint const & origin() const
Definition: TrackingRegion.h:74
ForwardDetLayer::surface
virtual const BoundSurface & surface() const
The surface of the GeometricSearchDet.
Definition: ForwardDetLayer.h:31
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:128
Bounds::thickness
virtual float thickness() const =0
TrackingRegion::Range
PixelRecoRange< float > Range
Definition: TrackingRegion.h:45
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
ForwardDetLayer::specificSurface
virtual const BoundDisk & specificSurface() const
Definition: ForwardDetLayer.h:41
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
RectangularEtaPhiTrackingRegion::phiWindow
OuterHitPhiPrediction phiWindow(const edm::EventSetup &iSetup) const
Definition: RectangularEtaPhiTrackingRegion.cc:252
PixelRecoUtilities::longitudinalBendingCorrection
double longitudinalBendingCorrection(double radius, double pt, const edm::EventSetup &iSetup)
Definition: PixelRecoUtilities.h:51
GeometricSearchDet::position
virtual const Surface::PositionType & position() const
Returns position of the surface.
Definition: GeometricSearchDet.h:33
TrackingRegion::ptMin
float ptMin() const
minimal pt of interest
Definition: TrackingRegion.h:83
PixelRecoRange::intersection
PixelRecoRange< T > intersection(const PixelRecoRange< T > &r) const
Definition: PixelRecoRange.h:44
funct::sqr
Square< F >::type sqr(const F &f)
Definition: Square.h:13
OuterHitPhiPrediction::setTolerance
void setTolerance(const Margin &tolerance)
Definition: OuterHitPhiPrediction.h:25
const Range& RectangularEtaPhiTrackingRegion::etaRange ( ) const
inline

allowed eta range [eta_min, eta_max] interval

Definition at line 148 of file RectangularEtaPhiTrackingRegion.h.

References theEtaRange.

Referenced by FastTSGFromL2Muon::clean(), FastTSGFromIOHit::clean(), TrackerSeedCleaner::clean(), GlobalTrajectoryBuilderBase::defineRegionOfInterest(), and SeedFromConsecutiveHitsCreator::makeSeed().

148 { return theEtaRange; }
TrackingRegion::Hits RectangularEtaPhiTrackingRegion::hits ( const edm::Event &  ev,
const edm::EventSetup &  es,
const SeedingLayerSetsHits::SeedingLayer &  layer 
) const
overridevirtual

get hits from layer compatible with region constraints

Implements TrackingRegion.

Definition at line 286 of file RectangularEtaPhiTrackingRegion.cc.

References alongMomentum, GeomDetEnumerators::barrel, newFWLiteAna::build, funct::cos(), SeedingLayerSetsHits::SeedingLayer::detLayer(), dir, GeomDetEnumerators::endcap, HLT_25ns14e33_v1_cff::estimator, f, edm::EventSetup::get(), h, OuterEstimator::hitCompatibility(), SeedingLayerSetsHits::SeedingLayer::hits(), DetLayer::location(), LogDebug, eostools::move(), phi, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, edm::ESHandle< class >::product(), LayerMeasurements::recHits(), query::result, makeMuonMisalignmentScenario::rot, funct::sin(), and DetLayer::subDetector().

289  {
290  TrackingRegion::Hits result;
291 
292  //ESTIMATOR
293 
294  const DetLayer * detLayer = layer.detLayer();
295  OuterEstimator * est = 0;
296 
297  bool measurementMethod = false;
298  if(theMeasurementTrackerUsage == UseMeasurementTracker::kAlways) measurementMethod = true;
299  else if(theMeasurementTrackerUsage == UseMeasurementTracker::kForSiStrips &&
300  !(detLayer->subDetector() == GeomDetEnumerators::PixelBarrel ||
301  detLayer->subDetector() == GeomDetEnumerators::PixelEndcap) ) measurementMethod = true;
302 
303  if(measurementMethod) {
304  edm::ESHandle<MagneticField> field;
305  es.get<IdealMagneticFieldRecord>().get(field);
306  const MagneticField * magField = field.product();
307 
308  const GlobalPoint vtx = origin();
309  GlobalVector dir = direction();
310 
311  if (detLayer->subDetector() == GeomDetEnumerators::PixelBarrel || (!theUseEtaPhi && detLayer->location() == GeomDetEnumerators::barrel)){
312  const BarrelDetLayer& bl = dynamic_cast<const BarrelDetLayer&>(*detLayer);
313  est = estimator(&bl,es);
314  } else if (detLayer->subDetector() == GeomDetEnumerators::PixelEndcap || (!theUseEtaPhi && detLayer->location() == GeomDetEnumerators::endcap)) {
315  const ForwardDetLayer& fl = dynamic_cast<const ForwardDetLayer&>(*detLayer);
316  est = estimator(&fl,es);
317  }
318 
319  EtaPhiMeasurementEstimator etaPhiEstimator ((theEtaRange.second-theEtaRange.first)*0.5f,
320  (thePhiMargin.left()+thePhiMargin.right())*0.5f);
321  MeasurementEstimator * findDetAndHits = &etaPhiEstimator;
322  if (est){
323  LogDebug("RectangularEtaPhiTrackingRegion")<<"use pixel specific estimator.";
324  findDetAndHits = est;
325  }
326  else{
327  LogDebug("RectangularEtaPhiTrackingRegion")<<"use generic etat phi estimator.";
328  }
329 
330  // TSOS
331  float phi = phiDirection();
332  // std::cout << "dir " << direction().x()/direction().perp() <<','<< direction().y()/direction().perp() << " " << sin(phi) <<','<<cos(phi)<< std::endl;
333  Surface::RotationType rot( sin(phi), -cos(phi), 0,
334  0, 0, -1,
335  cos(phi), sin(phi), 0);
336 
337  Plane::PlanePointer surface = Plane::build(GlobalPoint(0.,0.,0.), rot);
338  //TrajectoryStateOnSurface tsos(lpar, *surface, magField);
339 
340  FreeTrajectoryState fts( GlobalTrajectoryParameters(vtx, dir, 1, magField) );
341  TrajectoryStateOnSurface tsos(fts, *surface);
342 
343  // propagator
344  StraightLinePropagator prop( magField, alongMomentum);
345 
346  LayerMeasurements lm(theMeasurementTracker->measurementTracker(), *theMeasurementTracker);
347 
348  LayerMeasurements:: SimpleHitContainer hits;
349  lm.recHits(hits,*detLayer, tsos, prop, *findDetAndHits);
350  /*
351  { // old code
352  vector<TrajectoryMeasurement> meas = lm.measurements(*detLayer, tsos, prop, *findDetAndHits);
353  auto n=0UL;
354  for (auto const & im : meas)
355  if(im.recHit()->isValid()) ++n;
356  assert(n==hits.size());
357  // std::cout << "old/new " << n <<'/'<<hits.size() << std::endl;
358  }
359  */
360 
361  result.reserve(hits.size());
362  for (auto h : hits) {
363  cache.emplace_back(h);
364  result.emplace_back(h);
365  }
366 
367  LogDebug("RectangularEtaPhiTrackingRegion")<<" found "<< hits.size()<<" minus one measurements on layer: "<<detLayer->subDetector();
368  // std::cout << "RectangularEtaPhiTrackingRegion" <<" found "<< meas.size()<<" minus one measurements on layer: "<<detLayer->subDetector() << std::endl;
369 
370  } else {
371  //
372  // temporary solution
373  //
374  if (detLayer->location() == GeomDetEnumerators::barrel) {
375  const BarrelDetLayer& bl = dynamic_cast<const BarrelDetLayer&>(*detLayer);
376  est = estimator(&bl,es);
377  } else {
378  const ForwardDetLayer& fl = dynamic_cast<const ForwardDetLayer&>(*detLayer);
379  est = estimator(&fl,es);
380  }
381  if (!est) return result;
382 
383  auto layerHits = layer.hits();
384  result.reserve(layerHits.size());
385  for (auto && ih : layerHits) {
386  if ( est->hitCompatibility()(*ih) ) {
387  result.emplace_back( std::move(ih) );
388  }
389  }
390  }
391 
392  // std::cout << "RectangularEtaPhiTrackingRegion hits " << result.size() << std::endl;
393  delete est;
394 
395  return result;
396 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
RectangularEtaPhiTrackingRegion::estimator
OuterEstimator * estimator(const BarrelDetLayer *layer, const edm::EventSetup &iSetup) const
Definition: RectangularEtaPhiTrackingRegion.cc:132
TrackingRegion::origin
GlobalPoint const & origin() const
Definition: TrackingRegion.h:74
DetLayer::location
virtual Location location() const =0
Which part of the detector (barrel, endcap)
OuterEstimator::hitCompatibility
const OuterHitCompatibility & hitCompatibility() const
Definition: OuterEstimator.h:62
LayerMeasurements::SimpleHitContainer
std::vector< BaseTrackerRecHit * > SimpleHitContainer
Definition: LayerMeasurements.h:24
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
DetLayer::subDetector
virtual SubDetector subDetector() const =0
The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap)
RectangularEtaPhiTrackingRegion::theMeasurementTrackerUsage
const UseMeasurementTracker theMeasurementTrackerUsage
Definition: RectangularEtaPhiTrackingRegion.h:196
TrackingRegion::direction
GlobalVector const & direction() const
the direction around which region is constructed
Definition: TrackingRegion.h:65
RectangularEtaPhiTrackingRegion::theMeasurementTracker
const MeasurementTrackerEvent * theMeasurementTracker
Definition: RectangularEtaPhiTrackingRegion.h:199
RectangularEtaPhiTrackingRegion::hits
virtual TrackingRegion::Hits hits(const edm::Event &ev, const edm::EventSetup &es, const SeedingLayerSetsHits::SeedingLayer &layer) const override
get hits from layer compatible with region constraints
Definition: RectangularEtaPhiTrackingRegion.cc:286
query.result
tuple result
Definition: query.py:137
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
eostools.move
def move
Definition: eostools.py:508
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
TrackingRegion::phiDirection
float phiDirection() const
Definition: TrackingRegion.h:68
TrackingRegion::Hits
SeedingLayerSetsHits::Hits Hits
Definition: TrackingRegion.h:47
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
SeedingLayerSetsHits::SeedingLayer::detLayer
const DetLayer * detLayer() const
Definition: SeedingLayerSetsHits.h:60
dir
dbl *** dir
Definition: mlp_gen.cc:35
phi
Definition: DDAxes.h:10
void RectangularEtaPhiTrackingRegion::initEtaRange ( const GlobalVector &  dir,
const Margin &  margin 
)
private

Definition at line 46 of file RectangularEtaPhiTrackingRegion.cc.

References PV3DBase< T, PVType, FrameType >::eta(), eta(), TkTrackingRegionsMargin< T >::left(), and TkTrackingRegionsMargin< T >::right().

Referenced by RectangularEtaPhiTrackingRegion().

46  {
47  float eta = dir.eta();
48  theEtaRange = Range(eta-margin.left(), eta+margin.right());
49  theLambdaRange=Range(std::sinh(theEtaRange.min()),std::sinh(theEtaRange.max()));
50  theMeanLambda = std::sinh(theEtaRange.mean());
51 }
PixelRecoRange::max
T max() const
Definition: PixelRecoRange.h:31
PixelRecoRange::min
T min() const
Definition: PixelRecoRange.h:30
eta
T eta() const
Definition: Basic3DVectorLD.h:177
TrackingRegion::Range
PixelRecoRange< float > Range
Definition: TrackingRegion.h:45
margin
Double_t margin
Definition: trackSplitPlot.h:45
PixelRecoRange::mean
T mean() const
Definition: PixelRecoRange.h:32
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
static UseMeasurementTracker RectangularEtaPhiTrackingRegion::intToUseMeasurementTracker ( int  value)
inlinestatic

Definition at line 37 of file RectangularEtaPhiTrackingRegion.h.

References assert(), and relativeConstraints::value.

Referenced by SeedFilter::SeedFilter().

37  {
38  assert(value >= -1 && value <= 1);
39  return static_cast<UseMeasurementTracker>(value);
40  }
assert
assert(m_qm.get())
bool RectangularEtaPhiTrackingRegion::isPrecise ( ) const
inline

is precise error calculation switched on

Definition at line 155 of file RectangularEtaPhiTrackingRegion.h.

References thePrecise.

155 { return thePrecise; }
virtual std::string RectangularEtaPhiTrackingRegion::name ( void  ) const
inlinevirtual

Reimplemented from TrackingRegion.

Definition at line 173 of file RectangularEtaPhiTrackingRegion.h.

Referenced by cuy.divideElement::__init__(), cuy.plotElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), and config.CFG::__str__().

173 { return "RectangularEtaPhiTrackingRegion"; }
RectangularEtaPhiTrackingRegion& RectangularEtaPhiTrackingRegion::operator= ( RectangularEtaPhiTrackingRegion const &  )
delete
RectangularEtaPhiTrackingRegion& RectangularEtaPhiTrackingRegion::operator= ( RectangularEtaPhiTrackingRegion &&  )
default
const Margin& RectangularEtaPhiTrackingRegion::phiMargin ( ) const
inline

defined phi range around phi0, margin is [phi_left,phi_right]. region is defined in a range: [phi0-phi_left, phi0+phi_right]

Definition at line 152 of file RectangularEtaPhiTrackingRegion.h.

References thePhiMargin.

Referenced by FastTSGFromL2Muon::clean(), FastTSGFromIOHit::clean(), TrackerSeedCleaner::clean(), GlobalTrajectoryBuilderBase::defineRegionOfInterest(), and SeedFromConsecutiveHitsCreator::makeSeed().

152 { return thePhiMargin; }
OuterHitPhiPrediction RectangularEtaPhiTrackingRegion::phiWindow ( const edm::EventSetup &  iSetup) const
private

Definition at line 252 of file RectangularEtaPhiTrackingRegion.cc.

References PixelRecoUtilities::curvature(), bookConverter::max, and min().

253 {
254  float phi0 = phiDirection();
255  return OuterHitPhiPrediction(
256  OuterHitPhiPrediction::Range( phi0-thePhiMargin.left(),
257  phi0+thePhiMargin.right()),
258  OuterHitPhiPrediction::Range( curvature(invPtRange().min(),iSetup),
259  curvature(invPtRange().max(),iSetup)),
260  originRBound());
261 }
TrackingRegion::originRBound
float originRBound() const
bounds the particle vertex in the transverse plane
Definition: TrackingRegion.h:77
PixelRecoUtilities::curvature
T curvature(T InversePt, const edm::EventSetup &iSetup)
Definition: PixelRecoUtilities.h:41
min
T min(T a, T b)
Definition: MathUtil.h:58
TrackingRegion::phiDirection
float phiDirection() const
Definition: TrackingRegion.h:68
TrackingRegion::invPtRange
Range invPtRange() const
inverse pt range
Definition: TrackingRegion.h:86
std::string RectangularEtaPhiTrackingRegion::print ( void  ) const
virtual

Reimplemented from TrackingRegion.

Definition at line 398 of file RectangularEtaPhiTrackingRegion.cc.

References TrackingRegion::print().

398  {
399  std::ostringstream str;
400  str << TrackingRegionBase::print()
401  <<" eta: "<<theEtaRange<<" phi:"<<thePhiMargin
402  << "precise: "<<thePrecise;
403  return str.str();
404 }
TrackingRegion::print
virtual std::string print() const
Definition: TrackingRegion.h:117
HitRZConstraint RectangularEtaPhiTrackingRegion::rzConstraint ( ) const
private

Definition at line 265 of file RectangularEtaPhiTrackingRegion.cc.

References PixelRecoPointRZ::z().

265  {
266  HitRZConstraint::Point pLeft,pRight;
267  float zMin = origin().z() - originZBound();
268  float zMax = origin().z() + originZBound();
269  float rMin = -originRBound();
270  float rMax = originRBound();
271  if(theEtaRange.max() > 0) {
272  pRight = HitRZConstraint::Point(rMin,zMax);
273  } else {
274  pRight = HitRZConstraint::Point(rMax,zMax);
275  }
276  if (theEtaRange.min() > 0.) {
277  pLeft = HitRZConstraint::Point(rMax, zMin);
278  } else {
279  pLeft = HitRZConstraint::Point(rMin, zMin);
280  }
281  return HitRZConstraint(pLeft, theLambdaRange.min(),
282  pRight,theLambdaRange.max()
283  );
284 }
TrackingRegion::originRBound
float originRBound() const
bounds the particle vertex in the transverse plane
Definition: TrackingRegion.h:77
TrackingRegion::origin
GlobalPoint const & origin() const
Definition: TrackingRegion.h:74
PixelRecoRange::max
T max() const
Definition: PixelRecoRange.h:31
PixelRecoRange::min
T min() const
Definition: PixelRecoRange.h:30
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
TrackingRegion::originZBound
float originZBound() const
bounds the particle vertex in the longitudinal plane
Definition: TrackingRegion.h:80
HitRZConstraint::Point
SimpleLineRZ::Point Point
Definition: HitRZConstraint.h:11

Member Data Documentation

cacheHits RectangularEtaPhiTrackingRegion::cache
mutableprivate

Definition at line 212 of file RectangularEtaPhiTrackingRegion.h.

Range RectangularEtaPhiTrackingRegion::theEtaRange
private

Definition at line 192 of file RectangularEtaPhiTrackingRegion.h.

Referenced by etaRange().

Range RectangularEtaPhiTrackingRegion::theLambdaRange
private

Definition at line 193 of file RectangularEtaPhiTrackingRegion.h.

float RectangularEtaPhiTrackingRegion::theMeanLambda
private

Definition at line 195 of file RectangularEtaPhiTrackingRegion.h.

const MeasurementTrackerEvent* RectangularEtaPhiTrackingRegion::theMeasurementTracker
private

Definition at line 199 of file RectangularEtaPhiTrackingRegion.h.

const UseMeasurementTracker RectangularEtaPhiTrackingRegion::theMeasurementTrackerUsage
private

Definition at line 196 of file RectangularEtaPhiTrackingRegion.h.

Margin RectangularEtaPhiTrackingRegion::thePhiMargin
private

Definition at line 194 of file RectangularEtaPhiTrackingRegion.h.

Referenced by phiMargin().

bool RectangularEtaPhiTrackingRegion::thePrecise
private

Definition at line 197 of file RectangularEtaPhiTrackingRegion.h.

Referenced by isPrecise().

bool RectangularEtaPhiTrackingRegion::theUseEtaPhi
private

Definition at line 198 of file RectangularEtaPhiTrackingRegion.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5