CMS 3D CMS Logo

List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
CosmicTrackingRegion Class Reference

#include <CosmicTrackingRegion.h>

Inheritance diagram for CosmicTrackingRegion:
TrackingRegion

Public Member Functions

HitRZCompatibilitycheckRZ (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 override
 
CosmicTrackingRegionclone () const override
 
 CosmicTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, float dummy=0., const MeasurementTrackerEvent *measurementTracker=0)
 
 CosmicTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, const edm::ParameterSet &extra, const MeasurementTrackerEvent *measurementTracker=0)
 
 CosmicTrackingRegion (CosmicTrackingRegion const &rh)
 
TrackingRegion::Hits hits (const edm::EventSetup &es, const SeedingLayerSetsHits::SeedingLayer &layer) const override
 get hits from layer compatible with region constraints More...
 
std::string name () const override
 
 ~CosmicTrackingRegion () override
 
- 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
 
virtual std::string print () const
 
float ptMin () const
 minimal pt of interest More...
 
TrackingRegionrestrictedRegion (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 ()
 

Private Types

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

Private Member Functions

template<typename T >
void hits_ (const edm::EventSetup &es, const T &layer, TrackingRegion::Hits &result) const
 

Private Attributes

cacheHits cache
 
std::string measurementTrackerName_
 
const MeasurementTrackerEventtheMeasurementTracker_
 

Additional Inherited Members

- Public Types inherited from TrackingRegion
typedef SeedingLayerSetsHits::ConstRecHitPointer Hit
 
typedef SeedingLayerSetsHits::Hits Hits
 
typedef PixelRecoRange< float > Range
 

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 24 of file CosmicTrackingRegion.h.

Member Typedef Documentation

Definition at line 111 of file CosmicTrackingRegion.h.

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

Definition at line 112 of file CosmicTrackingRegion.h.

Constructor & Destructor Documentation

CosmicTrackingRegion::~CosmicTrackingRegion ( )
inlineoverride

Definition at line 28 of file CosmicTrackingRegion.h.

28 { }
CosmicTrackingRegion::CosmicTrackingRegion ( const GlobalVector dir,
const GlobalPoint vertexPos,
float  ptMin,
float  rVertex,
float  zVertex,
float  deltaEta,
float  deltaPhi,
float  dummy = 0.,
const MeasurementTrackerEvent measurementTracker = 0 
)
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

Definition at line 54 of file CosmicTrackingRegion.h.

Referenced by clone().

CosmicTrackingRegion::CosmicTrackingRegion ( const GlobalVector dir,
const GlobalPoint vertexPos,
float  ptMin,
float  rVertex,
float  zVertex,
float  deltaEta,
float  deltaPhi,
const edm::ParameterSet extra,
const MeasurementTrackerEvent measurementTracker = 0 
)
inline
CosmicTrackingRegion::CosmicTrackingRegion ( CosmicTrackingRegion const &  rh)
inline

Definition at line 79 of file CosmicTrackingRegion.h.

References hits().

79  :
81  theMeasurementTracker_(rh.theMeasurementTracker_),
82  measurementTrackerName_(rh.measurementTrackerName_){}
const MeasurementTrackerEvent * theMeasurementTracker_
TrackingRegion TrackingRegionBase

Member Function Documentation

HitRZCompatibility* CosmicTrackingRegion::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
inlineoverridevirtual

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

Implements TrackingRegion.

Definition at line 89 of file CosmicTrackingRegion.h.

94  {return nullptr; }
CosmicTrackingRegion* CosmicTrackingRegion::clone ( void  ) const
inlineoverridevirtual

Implements TrackingRegion.

Definition at line 96 of file CosmicTrackingRegion.h.

References CosmicTrackingRegion().

96 { return new CosmicTrackingRegion(*this); }
CosmicTrackingRegion(const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, float dummy=0., const MeasurementTrackerEvent *measurementTracker=0)
TrackingRegion::Hits CosmicTrackingRegion::hits ( const edm::EventSetup es,
const SeedingLayerSetsHits::SeedingLayer layer 
) const
overridevirtual

get hits from layer compatible with region constraints

Implements TrackingRegion.

Definition at line 25 of file CosmicTrackingRegion.cc.

References mps_fire::result.

Referenced by CosmicTrackingRegion().

27 {
29  hits_(es, layer, result);
30  return result;
31 }
void hits_(const edm::EventSetup &es, const T &layer, TrackingRegion::Hits &result) const
SeedingLayerSetsHits::Hits Hits
template<typename T >
void CosmicTrackingRegion::hits_ ( const edm::EventSetup es,
const T layer,
TrackingRegion::Hits result 
) const
private

Definition at line 34 of file CosmicTrackingRegion.cc.

References alongMomentum, AnalyticalPropagator_cfi::AnalyticalPropagator, Plane::build(), utilities::cache(), GeometricSearchDet::compatibleDets(), funct::cos(), dir, PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::freeState(), edm::EventSetup::get(), TrajectoryStateOnSurface::globalPosition(), trackerHitRTTI::isUndef(), TrajectoryStateOnSurface::isValid(), LogDebug, LayerMeasurements::measurements(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), edm::ESHandle< T >::product(), makeMuonMisalignmentScenario::rot, funct::sin(), GeometricSearchDet::surface(), and extraflags_cff::vtx.

Referenced by name().

36 {
37 
38  //get and name collections
39  //++++++++++++++++++++++++
40 
41 
42  //detector layer
43  const DetLayer * detLayer = layer.detLayer();
44  LogDebug("CosmicTrackingRegion") << "Looking at hits on subdet/layer " << layer.name();
45  EtaPhiMeasurementEstimator est(0.3,0.3);
46 
47  //magnetic field
49  es.get<IdealMagneticFieldRecord>().get(field);
50  const MagneticField * magField = field.product();
51 
52  //region
53  const GlobalPoint vtx = origin();
55  LogDebug("CosmicTrackingRegion") <<"The initial region characteristics are:" << "\n"
56  <<" Origin = " << origin() << "\n"
57  <<" Direction = " << direction() << "\n"
58  <<" Eta = " << origin().eta() << "\n"
59  <<" Phi = " << origin().phi();
60 
61  //trajectory state on surface
62  float phi = dir.phi();
63  Surface::RotationType rot( sin(phi), -cos(phi), 0,
64  0, 0, -1,
65  cos(phi), sin(phi), 0);
66 
67  Plane::PlanePointer surface = Plane::build(vtx, rot);
68  FreeTrajectoryState fts( GlobalTrajectoryParameters(vtx, dir, 1, magField) );
69  TrajectoryStateOnSurface tsos(fts, *surface);
70  LogDebug("CosmicTrackingRegion")
71  << "The state used to find measurement with the measurement tracker is:\n" << tsos;
72 
73  //propagator
74  AnalyticalPropagator prop( magField, alongMomentum);
75 
76  //propagation verification (debug)
77  //++++++++++++++++++++++++++++++++
78 
79  //creation of the state
80  TrajectoryStateOnSurface stateOnLayer = prop.propagate( *tsos.freeState(),
81  detLayer->surface());
82 
83  //verification of the state
84  if (stateOnLayer.isValid()){
85  LogDebug("CosmicTrackingRegion") << "The initial state propagates to the layer surface: \n" << stateOnLayer
86  << "R = " << stateOnLayer.globalPosition().perp() << "\n"
87  << "Eta = " << stateOnLayer.globalPosition().eta() << "\n"
88  << "Phi = " << stateOnLayer.globalPosition().phi();
89 
90  }
91  else{
92  LogDebug("CosmicTrackingRegion") << "The initial state does not propagate to the layer surface.";
93  }
94 
95  //number of compatible dets
97  vector<DetWithState> compatDets = detLayer->compatibleDets(tsos, prop, est);
98  LogDebug("CosmicTrackingRegion") << "Compatible dets = " << compatDets.size();
99 
100 
101  //get hits
102  //++++++++
103 
104  //measurement tracker (find hits)
106  vector<TrajectoryMeasurement> meas = lm.measurements(*detLayer, tsos, prop, est);
107  LogDebug("CosmicTrackingRegion") << "Number of Trajectory measurements = " << meas.size()
108  <<" but the last one is always an invalid hit, by construction.";
109 
110  //trajectory measurement
111 
112  // std::cout <<"CRegion b " << cache.size() << std::endl;
113 
114  // waiting for a migration at LayerMeasurements level and at seed builder level
115  for (auto const & im : meas) {
116  if(!im.recHit()->isValid()) continue;
117  assert(!trackerHitRTTI::isUndef(*im.recHit()->hit()));
118  auto ptrHit = (BaseTrackerRecHit *)(im.recHit()->hit()->clone());
119  cache.emplace_back(ptrHit);
120  result.emplace_back(ptrHit);
121  }
122 
123  // std::cout <<"CRegion a " << cache.size() << std::endl;
124 
125 }
#define LogDebug(id)
std::vector< TrajectoryMeasurement > measurements(const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
T perp() const
Definition: PV3DBase.h:72
GlobalPoint const & origin() const
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
GlobalPoint globalPosition() const
virtual std::vector< DetWithState > compatibleDets(const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
GlobalVector const & direction() const
the direction around which region is constructed
static PlanePointer build(Args &&...args)
Definition: Plane.h:33
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
virtual const BoundSurface & surface() const =0
The surface of the GeometricSearchDet.
std::pair< const GeomDet *, TrajectoryStateOnSurface > DetWithState
T eta() const
Definition: PV3DBase.h:76
T get() const
Definition: EventSetup.h:71
std::pair< const GeomDet *, TrajectoryStateOnSurface > DetWithState
bool isUndef(TrackingRecHit const &hit)
const MeasurementTrackerEvent * theMeasurementTracker_
dbl *** dir
Definition: mlp_gen.cc:35
T const * product() const
Definition: ESHandle.h:86
std::string CosmicTrackingRegion::name ( void  ) const
inlineoverridevirtual

Reimplemented from TrackingRegion.

Definition at line 98 of file CosmicTrackingRegion.h.

References hits_(), and mps_fire::result.

Referenced by config.CFG::__str__(), and validation.Sample::digest().

98 { return "CosmicTrackingRegion"; }

Member Data Documentation

cacheHits CosmicTrackingRegion::cache
mutableprivate

Definition at line 117 of file CosmicTrackingRegion.h.

std::string CosmicTrackingRegion::measurementTrackerName_
private

Definition at line 109 of file CosmicTrackingRegion.h.

Referenced by CosmicTrackingRegion().

const MeasurementTrackerEvent* CosmicTrackingRegion::theMeasurementTracker_
private

Definition at line 108 of file CosmicTrackingRegion.h.