StartingLayerFinder Class Reference

#include <StartingLayerFinder.h>

Public Member Functions
const std::vector< const ForwardDetLayer * >	firstNegPixelFwdLayer () const
const BarrelDetLayer *	firstPixelBarrelLayer () const
const std::vector< const ForwardDetLayer * >	firstPosPixelFwdLayer () const
const Propagator *	propagator () const
	StartingLayerFinder (const Propagator *aPropagator, const MeasurementTracker *tracker)
std::vector< const DetLayer * >	startingLayers (const FTS &aFts, float dr, float dz) const
std::vector< const DetLayer * >	startingLayers (const TrajectorySeed &aSeed) const
	~StartingLayerFinder ()

Private Types
typedef FreeTrajectoryState	FTS
typedef std::pair< float, float >	Range
typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
void	checkPixelLayers () const
bool	rangesIntersect (const Range &a, const Range &b) const

Private Attributes
std::vector< const ForwardDetLayer * >	theFirstNegPixelFwdLayer
const BarrelDetLayer *	theFirstPixelBarrelLayer
std::vector< const ForwardDetLayer * >	theFirstPosPixelFwdLayer
const MeasurementTracker *	theMeasurementTracker
bool	thePixelLayersValid
const Propagator *	thePropagator

Detailed Description

Definition at line 27 of file StartingLayerFinder.h.

Member Typedef Documentation

FTS

typedef FreeTrajectoryState StartingLayerFinder::FTS

private

Definition at line 29 of file StartingLayerFinder.h.

Range

typedef std::pair<float, float> StartingLayerFinder::Range

private

Definition at line 31 of file StartingLayerFinder.h.

TSOS

typedef TrajectoryStateOnSurface StartingLayerFinder::TSOS

private

Definition at line 30 of file StartingLayerFinder.h.

Constructor & Destructor Documentation

StartingLayerFinder()

StartingLayerFinder::StartingLayerFinder ( const Propagator *  aPropagator, const MeasurementTracker *  tracker  )

inline

Definition at line 34 of file StartingLayerFinder.h.

      :
 
        thePropagator(aPropagator),
        theMeasurementTracker(tracker),
        thePixelLayersValid(false),
        theFirstPixelBarrelLayer(nullptr),
        theFirstNegPixelFwdLayer(0),
        theFirstPosPixelFwdLayer(0) {}

~StartingLayerFinder()

StartingLayerFinder::~StartingLayerFinder ( )

inline

Definition at line 44 of file StartingLayerFinder.h.

{}

Member Function Documentation

checkPixelLayers()

void StartingLayerFinder::checkPixelLayers ( ) const

private

Definition at line 122 of file StartingLayerFinder.cc.

                                                 {
  if (!thePixelLayersValid) {
    const GeometricSearchTracker* theGeometricSearchTracker = theMeasurementTracker->geometricSearchTracker();
 
    theFirstPixelBarrelLayer = theGeometricSearchTracker->pixelBarrelLayers().front();
    theFirstNegPixelFwdLayer = theGeometricSearchTracker->negPixelForwardLayers();
    theFirstPosPixelFwdLayer = theGeometricSearchTracker->posPixelForwardLayers();
    thePixelLayersValid = true;
  }
}

References GeometricSearchTracker::negPixelForwardLayers(), GeometricSearchTracker::pixelBarrelLayers(), and GeometricSearchTracker::posPixelForwardLayers().

firstNegPixelFwdLayer()

const vector< const ForwardDetLayer * > StartingLayerFinder::firstNegPixelFwdLayer

(

)

const

Definition at line 112 of file StartingLayerFinder.cc.

                                                                                      {
  checkPixelLayers();
  return theFirstNegPixelFwdLayer;
}

firstPixelBarrelLayer()

const BarrelDetLayer * StartingLayerFinder::firstPixelBarrelLayer

(

)

const

Definition at line 107 of file StartingLayerFinder.cc.

                                                                       {
  checkPixelLayers();
  return theFirstPixelBarrelLayer;
}

firstPosPixelFwdLayer()

const vector< const ForwardDetLayer * > StartingLayerFinder::firstPosPixelFwdLayer

(

)

const

Definition at line 117 of file StartingLayerFinder.cc.

                                                                                      {
  checkPixelLayers();
  return theFirstPosPixelFwdLayer;
}

propagator()

const Propagator* StartingLayerFinder::propagator ( ) const

inline

Definition at line 54 of file StartingLayerFinder.h.

{ return thePropagator; }

References thePropagator.

rangesIntersect()

bool StartingLayerFinder::rangesIntersect ( const Range &  a, const Range &  b  ) const

inlineprivate

Definition at line 66 of file StartingLayerFinder.h.

                                                                    {
    if (a.first > b.second || b.first > a.second)
      return false;
    else
      return true;
  }

References a, and b.

startingLayers() [1/2]

vector< const DetLayer * > StartingLayerFinder::startingLayers	(	const FTS &	aFts,
		float	dr,
		float	dz
	)		const

Definition at line 18 of file StartingLayerFinder.cc.

                                                                                                     {
  vector<const DetLayer*> mylayers;
  mylayers.reserve(3);
 
  FTS fastFts(aFts.parameters());
 
  //barrel pixel
  TSOS pTsos = propagator()->propagate(fastFts, firstPixelBarrelLayer()->surface());
 
  if (pTsos.isValid()) {
    Range barrZRange(
        firstPixelBarrelLayer()->position().z() - 0.5 * (firstPixelBarrelLayer()->surface().bounds().length()),
        firstPixelBarrelLayer()->position().z() + 0.5 * (firstPixelBarrelLayer()->surface().bounds().length()));
    Range trajZRange(pTsos.globalPosition().z() - dz, pTsos.globalPosition().z() + dz);
 
    if (rangesIntersect(trajZRange, barrZRange)) {
      mylayers.push_back(firstPixelBarrelLayer());
    }
  }
 
  //negative fwd pixel
 
  for (auto infwd : firstPosPixelFwdLayer()) {
    pTsos = propagator()->propagate(fastFts, infwd->surface());
    if (pTsos.isValid()) {
      Range nfwdRRange(infwd->specificSurface().innerRadius(), infwd->specificSurface().outerRadius());
      Range trajRRange(pTsos.globalPosition().perp() - dr, pTsos.globalPosition().perp() + dr);
      if (rangesIntersect(trajRRange, nfwdRRange)) {
        mylayers.push_back(infwd);
      }
    }
  }
 
  //positive fwd pixel
  for (auto ipfwd : firstPosPixelFwdLayer()) {
    pTsos = propagator()->propagate(fastFts, ipfwd->surface());
    if (pTsos.isValid()) {
      Range pfwdRRange(ipfwd->specificSurface().innerRadius(), ipfwd->specificSurface().outerRadius());
      Range trajRRange(pTsos.globalPosition().perp() - dr, pTsos.globalPosition().perp() + dr);
      if (rangesIntersect(trajRRange, pfwdRRange)) {
        mylayers.push_back(ipfwd);
      }
    }
  }
 
  return mylayers;
}

References flavorHistoryFilter_cfi::dr, PVValHelper::dz, FreeTrajectoryState::parameters(), position, TrackCandidateProducer_cfi::propagator, and rangesIntersect().

startingLayers() [2/2]

vector< const DetLayer * > StartingLayerFinder::startingLayers

(

const TrajectorySeed &

aSeed

)

const

Definition at line 66 of file StartingLayerFinder.cc.

                                                                                             {
  float dr = 0., dz = 0.;
 
  if (propagator()->propagationDirection() != aSeed.direction())
    return vector<const DetLayer*>();
 
  if (aSeed.nHits() != 2)
    return vector<const DetLayer*>();
 
  TrackingRecHitCollection::const_iterator firstHit = aSeed.recHits().first;
  const TrackingRecHit* recHit1 = &(*firstHit);
  const DetLayer* hit1Layer = theMeasurementTracker->geometricSearchTracker()->detLayer(recHit1->geographicalId());
 
  TrackingRecHitCollection::const_iterator secondHit = aSeed.recHits().second;
  const TrackingRecHit* recHit2 = &(*secondHit);
  const DetLayer* hit2Layer = theMeasurementTracker->geometricSearchTracker()->detLayer(recHit2->geographicalId());
 
  GeomDetEnumerators::Location p1 = hit1Layer->location();
  GeomDetEnumerators::Location p2 = hit2Layer->location();
 
  if (p1 == GeomDetEnumerators::barrel && p2 == GeomDetEnumerators::barrel) {
    dr = 0.1;
    dz = 5.;
  } else if (p1 == GeomDetEnumerators::endcap && p2 == GeomDetEnumerators::endcap) {
    dr = 5.;
    dz = 0.1;
  } else {
    dr = 0.1;
    dz = 0.1;
  }
 
  const GeomDet* gdet = theMeasurementTracker->geomTracker()->idToDet(DetId(aSeed.startingState().detId()));
 
  TrajectoryStateOnSurface tsos = trajectoryStateTransform::transientState(
      aSeed.startingState(), &(gdet->surface()), thePropagator->magneticField());
 
  const FreeTrajectoryState* fts = tsos.freeTrajectoryState();
 
  return startingLayers(*fts, dr, dz);
}

References GeomDetEnumerators::barrel, PTrajectoryStateOnDet::detId(), TrajectorySeed::direction(), flavorHistoryFilter_cfi::dr, PVValHelper::dz, GeomDetEnumerators::endcap, TrajectoryStateOnSurface::freeTrajectoryState(), TrackingRecHit::geographicalId(), DetLayer::location(), TrajectorySeed::nHits(), p1, p2, TrackCandidateProducer_cfi::propagator, TrajectorySeed::recHits(), TrajectorySeed::startingState(), GeomDet::surface(), and trajectoryStateTransform::transientState().

Member Data Documentation

theFirstNegPixelFwdLayer

std::vector<const ForwardDetLayer*> StartingLayerFinder::theFirstNegPixelFwdLayer

mutableprivate

Definition at line 61 of file StartingLayerFinder.h.

theFirstPixelBarrelLayer

const BarrelDetLayer* StartingLayerFinder::theFirstPixelBarrelLayer

mutableprivate

Definition at line 60 of file StartingLayerFinder.h.

theFirstPosPixelFwdLayer

std::vector<const ForwardDetLayer*> StartingLayerFinder::theFirstPosPixelFwdLayer

mutableprivate

Definition at line 62 of file StartingLayerFinder.h.

theMeasurementTracker

const MeasurementTracker* StartingLayerFinder::theMeasurementTracker

private

Definition at line 58 of file StartingLayerFinder.h.

thePixelLayersValid

bool StartingLayerFinder::thePixelLayersValid

mutableprivate

Definition at line 59 of file StartingLayerFinder.h.

thePropagator

const Propagator* StartingLayerFinder::thePropagator

private

Definition at line 57 of file StartingLayerFinder.h.

Referenced by propagator().