#include <GlobalTrackingRegion.h>

Inheritance diagram for GlobalTrackingRegion:

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 GlobalTrackingRegion *	clone () const

	GlobalTrackingRegion (float ptMin, const GlobalPoint &origin, float originRadius, float originHalfLength, bool precise=false)

	GlobalTrackingRegion (float ptMin=1., float originRadius=0.2, float originHalfLength=22.7, float originZPos=0., bool precise=false)

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

virtual std::string	name () const

virtual std::string	print () const

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

Private Attributes
bool	thePrecise

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

An implementation of the TrackingRegion where the region of interest is global, ie there are no constraints on the allowed direction of particles of interest

Definition at line 13 of file GlobalTrackingRegion.h.

Constructor & Destructor Documentation

GlobalTrackingRegion::GlobalTrackingRegion	(	float	ptMin,
		const GlobalPoint &	origin,
		float	originRadius,
		float	originHalfLength,
		bool	precise = `false`
	)

inline

Construct from minimal track P_t, and origin size and position. The origin region is a cylinder of radius originRadius, half length originHalfLength, positioned at "origin". This class DOES provide the possibility to displace the origin in the transverse plane.

Definition at line 22 of file GlobalTrackingRegion.h.

Referenced by clone().

     :  TrackingRegionBase(GlobalVector( 0, 0, 0), origin,
       Range( -1/ptMin, 1/ptMin), originRadius, originHalfLength),
       thePrecise(precise) { }

GlobalTrackingRegion::GlobalTrackingRegion	(	float	ptMin = `1.`,
		float	originRadius = `0.2`,
		float	originHalfLength = `22.7`,
		float	originZPos = `0.`,
		bool	precise = `false`
	)

inline

Definition at line 29 of file GlobalTrackingRegion.h.

     : TrackingRegionBase(GlobalVector( 0, 0, 0), GlobalPoint( 0, 0, originZPos),
       Range( -1/ptMin, 1/ptMin), originRadius, originHalfLength),
       thePrecise(precise) { }

Member Function Documentation

HitRZCompatibility * GlobalTrackingRegion::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

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

Implements TrackingRegion.

Definition at line 37 of file GlobalTrackingRegion.cc.

References constexpr, PixelRecoLineRZ::cotLine(), SimpleLineRZ::cotLine(), GeomDetEnumerators::isBarrel(), DetLayer::isBarrel(), nSigmaPhi, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, PtMinSelector_cfg::ptMin, DetLayer::seqNum(), funct::sqr(), mathSSE::sqrt(), DetLayer::subDetector(), unlikely, vtxMean(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

 {
 
   bool isBarrel = layer->isBarrel();
   bool isPixel = (layer->subDetector() == PixelBarrel || layer->subDetector() == PixelEndcap);
   
   if unlikely(!outerlayer) {
       GlobalPoint ohit =  outerHit->globalPosition();
     lr = sqrt( sqr(ohit.x()-origin().x())+sqr(ohit.y()-origin().y()) );
     gz = ohit.z();
     dr = outerHit->errorGlobalR();
     dz = outerHit->errorGlobalZ();
     }
   
 
   PixelRecoPointRZ outerred(lr, gz);
 
 
   PixelRecoPointRZ vtxR = (gz > origin().z()+originZBound()) ?
       PixelRecoPointRZ(-originRBound(), origin().z()+originZBound())
     : PixelRecoPointRZ( originRBound(), origin().z()+originZBound());
   PixelRecoPointRZ vtxL = (gz< origin().z()-originZBound()) ?
       PixelRecoPointRZ(-originRBound(), origin().z()-originZBound())
     : PixelRecoPointRZ( originRBound(), origin().z()-originZBound()); 
 
   if unlikely((!thePrecise) &&(isPixel )) {
     double VcotMin = PixelRecoLineRZ( vtxR, outerred).cotLine();
     double VcotMax = PixelRecoLineRZ( vtxL, outerred).cotLine();
     return new HitEtaCheck(isBarrel, outerred, VcotMax, VcotMin);
   }
   
   constexpr float nSigmaPhi = 3.;
  
   dr *= nSigmaPhi;
   dz *= nSigmaPhi;
   PixelRecoPointRZ  outerL, outerR;
 
   if (layer->isBarrel()) {
     outerL = PixelRecoPointRZ(lr, gz-dz);
     outerR = PixelRecoPointRZ(lr, gz+dz);
   } 
   else if (gz > 0) {
     outerL = PixelRecoPointRZ(lr+dr, gz);
     outerR = PixelRecoPointRZ(lr-dr, gz);
   } 
   else {
     outerL = PixelRecoPointRZ(lr-dr, gz);
     outerR = PixelRecoPointRZ(lr+dr, gz);
   }
   
   MultipleScatteringParametrisation iSigma(layer,iSetup);
   PixelRecoPointRZ vtxMean(0.,origin().z());
 
   /*
   float innerScatt=0;
   if (outerlayer) {
     innerScatt = 3.f * iSigma( ptMin(), vtxMean, outerred);
     float anew = 3.f * iSigma( ptMin(), vtxMean, outerred, outerlayer->seqNum());
     if (std::abs( (innerScatt-anew)/innerScatt) > .05)  
     std::cout << "MS old/new in " << outerlayer->seqNum() << " " << layer->seqNum()
               << ": " << innerScatt <<  " / " <<   anew
              << std::endl;
   } else
   innerScatt = 3.f * iSigma( ptMin(), vtxMean, outerred);
   */
 
 
   float innerScatt = 3.f * ( outerlayer ?
      iSigma( ptMin(), vtxMean, outerred, outerlayer->seqNum())
     :  iSigma( ptMin(), vtxMean, outerred) ) ;
 
 
   //
   //
   //
   SimpleLineRZ leftLine( vtxL, outerL);
   SimpleLineRZ rightLine( vtxR, outerR);
   HitRZConstraint rzConstraint(leftLine, rightLine);
   float cotTheta = SimpleLineRZ(vtxMean,outerred).cotLine();
 
   if (isBarrel) {
     float sinTheta = 1/std::sqrt(1+sqr(cotTheta));
     float corrZ = innerScatt/sinTheta + dz;
     return new HitZCheck(rzConstraint, HitZCheck::Margin(corrZ,corrZ));
   } else {
     float cosTheta = 1/std::sqrt(1+sqr(1/cotTheta));
     float corrR = innerScatt/cosTheta + dr;
     return new HitRCheck( rzConstraint, HitRCheck::Margin(corrR,corrR));
   }
 }

virtual GlobalTrackingRegion* GlobalTrackingRegion::clone ( void ) const

inlinevirtual

Implements TrackingRegion.

Definition at line 48 of file GlobalTrackingRegion.h.

References GlobalTrackingRegion().

                                               { 
     return new GlobalTrackingRegion(*this);
   }

TrackingRegion::Hits GlobalTrackingRegion::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 27 of file GlobalTrackingRegion.cc.

References SeedingLayerSetsHits::SeedingLayer::hits().

Referenced by SimpleCosmicBONSeeder::triplets().

                                                            {
   return layer.hits();
 }

virtual std::string GlobalTrackingRegion::name ( void ) const

inlinevirtual

Reimplemented from TrackingRegion.

Definition at line 52 of file GlobalTrackingRegion.h.

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

52 { return "GlobalTrackingRegion"; }

std::string GlobalTrackingRegion::print ( void ) const

virtual

Reimplemented from TrackingRegion.

Definition at line 21 of file GlobalTrackingRegion.cc.

References TrackingRegion::print().

                                             {
   std::ostringstream str;
   str << TrackingRegionBase::print() << "precise: "<<thePrecise;
   return str.str(); 
 }

Member Data Documentation

bool GlobalTrackingRegion::thePrecise

private

Definition at line 56 of file GlobalTrackingRegion.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation