#include <CosmicTrackingRegion.h>

Inheritance diagram for CosmicTrackingRegion:

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

CosmicTrackingRegion *	clone () const

	CosmicTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, float dummy=0.)

	CosmicTrackingRegion (const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, const edm::ParameterSet &extra)

	CosmicTrackingRegion (CosmicTrackingRegion const &rh)

virtual TrackingRegion::ctfHits	hits (const edm::Event &ev, const edm::EventSetup &es, const ctfseeding::SeedingLayer *layer) const

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

std::string	name () const

virtual	~CosmicTrackingRegion ()

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

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

Private Types
using	cacheHitPointer = mayown_ptr< BaseTrackerRecHit >

using	cacheHits = std::vector< cacheHitPointer >

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

Private Attributes
cacheHits	cache

std::string	measurementTrackerName_

Additional Inherited Members
Public Types inherited from TrackingRegion
using	ctfHits = ctfseeding::SeedingLayer::Hits

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

using CosmicTrackingRegion::cacheHitPointer = mayown_ptr<BaseTrackerRecHit>

private

Definition at line 113 of file CosmicTrackingRegion.h.

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

private

Definition at line 114 of file CosmicTrackingRegion.h.

Constructor & Destructor Documentation

virtual CosmicTrackingRegion::~CosmicTrackingRegion ( )

inlinevirtual

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

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

     : TrackingRegionBase( dir, vertexPos, Range( -1/ptMin, 1/ptMin), 
               rVertex, zVertex),
       measurementTrackerName_("")
   { }

CosmicTrackingRegion::CosmicTrackingRegion	(	const GlobalVector &	dir,
		const GlobalPoint &	vertexPos,
		float	ptMin,
		float	rVertex,
		float	zVertex,
		float	deltaEta,
		float	deltaPhi,
		const edm::ParameterSet &	extra
	)

inline

Definition at line 64 of file CosmicTrackingRegion.h.

References edm::ParameterSet::getParameter(), measurementTrackerName_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : TrackingRegionBase( dir, vertexPos, Range( -1/ptMin, 1/ptMin),
               rVertex, zVertex)
   {
     measurementTrackerName_ = extra.getParameter<std::string>("measurementTrackerName");
   }

CosmicTrackingRegion::CosmicTrackingRegion ( CosmicTrackingRegion const & rh )

inline

Definition at line 75 of file CosmicTrackingRegion.h.

                                                         : 
   TrackingRegionBase(rh),
   measurementTrackerName_(rh.measurementTrackerName_){}

Member Function Documentation

virtual 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

inlinevirtual

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

Implements TrackingRegion.

Definition at line 91 of file CosmicTrackingRegion.h.

96 {return 0; }

CosmicTrackingRegion* CosmicTrackingRegion::clone ( void ) const

inlinevirtual

Implements TrackingRegion.

Definition at line 98 of file CosmicTrackingRegion.h.

References CosmicTrackingRegion().

98 { return new CosmicTrackingRegion(*this); }

CosmicTrackingRegion::CosmicTrackingRegion

CosmicTrackingRegion(const GlobalVector &dir, const GlobalPoint &vertexPos, float ptMin, float rVertex, float zVertex, float deltaEta, float deltaPhi, float dummy=0.)

Definition: CosmicTrackingRegion.h:54

TrackingRegion::ctfHits CosmicTrackingRegion::hits	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		const ctfseeding::SeedingLayer *	layer
	)		const

virtual

Definition at line 27 of file CosmicTrackingRegion.cc.

References h, query::result, and tmp.

 {
   TrackingRegion::ctfHits result;
   TrackingRegion::Hits tmp;
   hits_(ev, es, *layer, tmp);
   result.reserve(tmp.size());
   for ( auto h : tmp) result.emplace_back(*h); // not owned
   return result;
 }

TrackingRegion::Hits CosmicTrackingRegion::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 39 of file CosmicTrackingRegion.cc.

References query::result.

 {
   TrackingRegion::Hits result;
   hits_(ev, es, layer, result);
   return result;
 }

template<typename T >

void CosmicTrackingRegion::hits_	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		const T &	layer,
		TrackingRegion::Hits &	result
	)		const

private

Definition at line 49 of file CosmicTrackingRegion.cc.

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

 {
 
   //get and name collections
   //++++++++++++++++++++++++
 
  
   //detector layer
   const DetLayer * detLayer = layer.detLayer();
   LogDebug("CosmicTrackingRegion") << "Looking at hits on subdet/layer " << layer.name();
   EtaPhiMeasurementEstimator est(0.3,0.3);
 
   //magnetic field
   edm::ESHandle<MagneticField> field;
   es.get<IdealMagneticFieldRecord>().get(field);
   const MagneticField * magField = field.product();
 
   //region
   const GlobalPoint vtx = origin();
   GlobalVector dir = direction();
   LogDebug("CosmicTrackingRegion") <<"The initial region characteristics are:" << "\n"
                    <<" Origin    = " << origin() << "\n"
                    <<" Direction = " << direction() << "\n" 
                    <<" Eta = " << origin().eta()  << "\n" 
                    <<" Phi = " << origin().phi();
      
   //trajectory state on surface
   float phi = dir.phi();
   Surface::RotationType rot( sin(phi), -cos(phi),           0,
                              0,                0,          -1,
                              cos(phi),  sin(phi),           0);
 
   Plane::PlanePointer surface = Plane::build(vtx, rot);
   FreeTrajectoryState fts( GlobalTrajectoryParameters(vtx, dir, 1, magField) );
   TrajectoryStateOnSurface tsos(fts, *surface);
   LogDebug("CosmicTrackingRegion") 
     << "The state used to find measurement with the measurement tracker is:\n" << tsos;
 
   //propagator
   AnalyticalPropagator prop( magField, alongMomentum);
 
   //propagation verification (debug)
   //++++++++++++++++++++++++++++++++
 
   //creation of the state
   TrajectoryStateOnSurface stateOnLayer = prop.propagate( *tsos.freeState(),
                               detLayer->surface());
   
   //verification of the state
   if (stateOnLayer.isValid()){
     LogDebug("CosmicTrackingRegion") << "The initial state propagates to the layer surface: \n" << stateOnLayer
                      << "R   = " << stateOnLayer.globalPosition().perp() << "\n"
                      << "Eta = " << stateOnLayer.globalPosition().eta() << "\n"
                      << "Phi = " << stateOnLayer.globalPosition().phi();
 
   }
   else{
     LogDebug("CosmicTrackingRegion") << "The initial state does not propagate to the layer surface.";
   }
 
   //number of compatible dets
   typedef DetLayer::DetWithState DetWithState;
   vector<DetWithState> compatDets = detLayer->compatibleDets(tsos, prop, est);
   LogDebug("CosmicTrackingRegion") << "Compatible dets = " << compatDets.size();
   
 
   //get hits
   //++++++++
 
   //measurement tracker (find hits)
   edm::ESHandle<MeasurementTracker> measurementTrackerESH;
   es.get<CkfComponentsRecord>().get(measurementTrackerName_,measurementTrackerESH);
   const MeasurementTracker * measurementTracker = measurementTrackerESH.product(); 
   edm::Handle<MeasurementTrackerEvent> mte;
   ev.getByLabel(edm::InputTag("MeasurementTrackerEvent"), mte);
   LayerMeasurements lm(*measurementTracker, *mte);
   vector<TrajectoryMeasurement> meas = lm.measurements(*detLayer, tsos, prop, est);
   LogDebug("CosmicTrackingRegion") << "Number of Trajectory measurements = " << meas.size()
                    <<" but the last one is always an invalid hit, by construction.";
 
   //trajectory measurement
 
   // std::cout <<"CRegion b " << cache.size() << std::endl;
 
   // waiting for a migration at LayerMeasurements level and at seed builder level
   for (auto const & im : meas) {
     if(!im.recHit()->isValid()) continue;
     assert(!trackerHitRTTI::isUndef(*im.recHit()->hit()));
     auto ptrHit = (BaseTrackerRecHit *)(im.recHit()->hit()->clone());
     cache.emplace_back(ptrHit);
     result.emplace_back(ptrHit);
   }
 
   // std::cout <<"CRegion a " << cache.size() << std::endl;
 
 }

std::string CosmicTrackingRegion::name ( void ) const

inlinevirtual

Reimplemented from TrackingRegion.

Definition at line 100 of file CosmicTrackingRegion.h.

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

100 { return "CosmicTrackingRegion"; }

Member Data Documentation

cacheHits CosmicTrackingRegion::cache

mutableprivate

Definition at line 119 of file CosmicTrackingRegion.h.

std::string CosmicTrackingRegion::measurementTrackerName_

private

Definition at line 111 of file CosmicTrackingRegion.h.

Referenced by CosmicTrackingRegion().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation