CompositeTECWedge Class Reference

#include <CompositeTECWedge.h>

Inheritance diagram for CompositeTECWedge:

Public Member Functions
const std::vector< const GeomDet * > &	basicComponents () const override
std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &ts, const Propagator &, const MeasurementEstimator &) const override __attribute__((cold))
const std::vector< const GeometricSearchDet * > &	components () const override __attribute__((cold))
	Returns basic components, if any. More...
	CompositeTECWedge (std::vector< const GeomDet * > &innerDets, std::vector< const GeomDet * > &outerDets) __attribute__((cold))
void	groupedCompatibleDetsV (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetGroup > &result) const override __attribute__((hot))
	~CompositeTECWedge () override __attribute__((cold))
Public Member Functions inherited from TECWedge
virtual const BoundDiskSector &	specificSurface () const final
const BoundSurface &	surface () const final
	The surface of the GeometricSearchDet. More...
	TECWedge ()
Public Member Functions inherited from GeometricSearchDet
virtual std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
virtual void	compatibleDetsV (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetWithState > &result) const
	GeometricSearchDet (bool doHaveGroups)
virtual std::vector< DetGroup >	groupedCompatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
bool	hasGroups () const
virtual const Surface::PositionType &	position () const
	Returns position of the surface. More...
virtual	~GeometricSearchDet ()

Private Member Functions
bool	addClosest (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, std::vector< DetGroup > &result) const __attribute__((hot))
SubLayerCrossings	computeCrossings (const TrajectoryStateOnSurface &tsos, PropagationDirection propDir) const __attribute__((hot))
int	findClosestDet (const GlobalPoint &startPos, int sectorId) const
void	searchNeighbors (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, float window, std::vector< DetGroup > &result, bool checkClosest) const __attribute__((hot))
const std::vector< const GeomDet * > &	subWedge (int ind) const

Private Attributes
std::vector< const GeomDet * >	theBackDets
ReferenceCountingPointer< BoundDiskSector >	theBackSector
std::vector< const GeomDet * >	theDets
std::vector< const GeomDet * >	theFrontDets
ReferenceCountingPointer< BoundDiskSector >	theFrontSector

Additional Inherited Members
Public Types inherited from GeometricSearchDet
typedef std::pair< const GeomDet *, TrajectoryStateOnSurface >	DetWithState
typedef BoundSurface::PositionType	PositionType
typedef BoundSurface::RotationType	RotationType
typedef TrajectoryStateOnSurface	TrajectoryState
Protected Attributes inherited from TECWedge
ReferenceCountingPointer< BoundDiskSector >	theDiskSector
Protected Attributes inherited from GeometricSearchDet
bool	haveGroups
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A concrete implementation for TEC layer built out of TECPetals

Definition at line 12 of file CompositeTECWedge.h.

Constructor & Destructor Documentation

CompositeTECWedge::CompositeTECWedge	(	std::vector< const GeomDet * > &	innerDets,
		std::vector< const GeomDet * > &	outerDets
	)

Definition at line 33 of file CompositeTECWedge.cc.

References LogDebug, PV3DBase< T, PVType, FrameType >::perp(), GeometricSearchDet::position(), theBackDets, theBackSector, theDets, TECWedge::theDiskSector, theFrontDets, theFrontSector, and PV3DBase< T, PVType, FrameType >::z().

    : theFrontDets(innerDets.begin(), innerDets.end()), theBackDets(outerDets.begin(), outerDets.end()) {
  theDets.assign(theFrontDets.begin(), theFrontDets.end());
  theDets.insert(theDets.end(), theBackDets.begin(), theBackDets.end());

  //
  std::sort(theFrontDets.begin(), theFrontDets.end(), DetPhiLess());
  std::sort(theBackDets.begin(), theBackDets.end(), DetPhiLess());

  theFrontSector = ForwardDiskSectorBuilderFromDet()(theFrontDets);
  theBackSector = ForwardDiskSectorBuilderFromDet()(theBackDets);
  theDiskSector = ForwardDiskSectorBuilderFromDet()(theDets);

  //--------- DEBUG INFO --------------
  LogDebug("TkDetLayers") << "DEBUG INFO for CompositeTECWedge"
                          << "\n"
                          << "TECWedge z, perp,innerRadius,outerR: " << this->position().z() << " , "
                          << this->position().perp() << " , " << theDiskSector->innerRadius() << " , "
                          << theDiskSector->outerRadius();

  for (vector<const GeomDet*>::const_iterator it = theFrontDets.begin(); it != theFrontDets.end(); it++) {
    LogDebug("TkDetLayers") << "frontDet phi,z,r: " << (*it)->surface().position().phi() << " , "
                            << (*it)->surface().position().z() << " , " << (*it)->surface().position().perp();
  }

  for (vector<const GeomDet*>::const_iterator it = theBackDets.begin(); it != theBackDets.end(); it++) {
    LogDebug("TkDetLayers") << "backDet phi,z,r: " << (*it)->surface().phi() << " , " << (*it)->surface().position().z()
                            << " , " << (*it)->surface().position().perp();
  }
  //-----------------------------------
}

CompositeTECWedge::~CompositeTECWedge ( )

override

Definition at line 65 of file CompositeTECWedge.cc.

{}

Member Function Documentation

bool CompositeTECWedge::addClosest ( const TrajectoryStateOnSurface &  tsos, const Propagator &  prop, const MeasurementEstimator &  est, const SubLayerCrossing &  crossing, std::vector< DetGroup > &  result  ) const

private

Definition at line 162 of file CompositeTECWedge.cc.

References CompatibleDetToGroupAdder::add(), SubLayerCrossing::closestDetIndex(), LogDebug, mps_fire::result, SubLayerCrossing::subLayerIndex(), and subWedge().

Referenced by basicComponents(), and groupedCompatibleDetsV().

                                                                   {
  const vector<const GeomDet*>& sWedge(subWedge(crossing.subLayerIndex()));

  LogDebug("TkDetLayers") << "in CompositeTECWedge,adding GeomDet at r,z,phi: ("
                          << sWedge[crossing.closestDetIndex()]->position().perp() << ","
                          << sWedge[crossing.closestDetIndex()]->position().z() << ","
                          << sWedge[crossing.closestDetIndex()]->position().phi() << ")";

  return CompatibleDetToGroupAdder().add(*sWedge[crossing.closestDetIndex()], tsos, prop, est, result);
}

const std::vector<const GeomDet*>& CompositeTECWedge::basicComponents ( ) const

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 20 of file CompositeTECWedge.h.

References addClosest(), compatible(), components(), computeCrossings(), findClosestDet(), groupedCompatibleDetsV(), mps_fire::result, searchNeighbors(), theDets, and svgfig::window().

{ return theDets; }

pair< bool, TrajectoryStateOnSurface > CompositeTECWedge::compatible ( const TrajectoryStateOnSurface &  ts, const Propagator &  , const MeasurementEstimator &    ) const

overridevirtual

tests the geometrical compatibility of the Det with the predicted state. The FreeTrajectoryState argument is propagated to the Det surface using the Propagator argument. The resulting TrajectoryStateOnSurface is tested for compatibility with the surface bounds. If compatible, a std::pair< true, propagatedState> is returned. If the propagation fails, or if the state is not compatible, a std::pair< false, propagatedState> is returned.

Implements GeometricSearchDet.

Definition at line 71 of file CompositeTECWedge.cc.

Referenced by basicComponents().

                                                                                                      {
  edm::LogError("TkDetLayers") << "temporary dummy implementation of CompositeTECWedge::compatible()!!";
  return pair<bool, TrajectoryStateOnSurface>();
}

const vector< const GeometricSearchDet * > & CompositeTECWedge::components ( ) const

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 67 of file CompositeTECWedge.cc.

Referenced by basicComponents().

                                                                             {
  throw DetLayerException("CompositeTECWedge doesn't have GeometricSearchDet components");
}

SubLayerCrossings CompositeTECWedge::computeCrossings ( const TrajectoryStateOnSurface &  tsos, PropagationDirection  propDir  ) const

private

Definition at line 110 of file CompositeTECWedge.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, findClosestDet(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), LogDebug, HelixForwardPlaneCrossing::pathLength(), HelixForwardPlaneCrossing::position(), rho, theBackDets, theBackSector, theFrontDets, theFrontSector, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by basicComponents(), and groupedCompatibleDetsV().

                                                                                          {
  HelixPlaneCrossing::PositionType startPos(startingState.globalPosition());
  HelixPlaneCrossing::DirectionType startDir(startingState.globalMomentum());

  auto rho = startingState.transverseCurvature();

  HelixForwardPlaneCrossing crossing(startPos, startDir, rho, propDir);

  pair<bool, double> frontPath = crossing.pathLength(*theFrontSector);
  if (!frontPath.first)
    return SubLayerCrossings();

  pair<bool, double> backPath = crossing.pathLength(*theBackSector);
  if (!backPath.first)
    return SubLayerCrossings();

  GlobalPoint gFrontPoint(crossing.position(frontPath.second));
  GlobalPoint gBackPoint(crossing.position(backPath.second));

  LogDebug("TkDetLayers") << "in TECWedge,front crossing r,z,phi: (" << gFrontPoint.perp() << "," << gFrontPoint.z()
                          << "," << gFrontPoint.phi() << ")";

  LogDebug("TkDetLayers") << "in TECWedge,back crossing r,z,phi: (" << gBackPoint.perp() << "," << gBackPoint.z() << ","
                          << gBackPoint.phi() << ")" << endl;

  int frontIndex = findClosestDet(gFrontPoint, 0);
  SubLayerCrossing frontSLC(0, frontIndex, gFrontPoint);

  int backIndex = findClosestDet(gBackPoint, 1);
  SubLayerCrossing backSLC(1, backIndex, gBackPoint);

  auto frontDist = std::abs(Geom::deltaPhi(gFrontPoint.barePhi(), theFrontDets[frontIndex]->surface().phi()));
  /*
    float frontDist = theFrontDets[frontIndex]->surface().phi()  - gFrontPoint.phi(); 
    frontDist *= Geom::phiLess( theFrontDets[frontIndex]->surface().phi(),gFrontPoint.barePhi()) ? -1. : 1.; 
    if (frontDist < 0.) { frontDist += 2.*Geom::pi();}
  */
  auto backDist = std::abs(Geom::deltaPhi(gBackPoint.barePhi(), theBackDets[backIndex]->surface().phi()));
  /*
    float backDist = theBackDets[backIndex]->surface().phi()  - gBackPoint.phi(); 
    backDist  *= Geom::phiLess( theBackDets[backIndex]->surface().phi(),gBackPoint.barePhi()) ? -1. : 1.;
    if ( backDist < 0.) { backDist  += 2.*Geom::pi();}
  */

  if (frontDist < backDist) {
    return SubLayerCrossings(frontSLC, backSLC, 0);
  } else {
    return SubLayerCrossings(backSLC, frontSLC, 1);
  }
}

int CompositeTECWedge::findClosestDet ( const GlobalPoint &  startPos, int  sectorId  ) const

private

Definition at line 217 of file CompositeTECWedge.cc.

References funct::abs(), mps_fire::i, theBackDets, theFrontDets, toLocal(), and x.

Referenced by basicComponents(), and computeCrossings().

                                                                                     {
  vector<const GeomDet*> const& myDets = sectorId == 0 ? theFrontDets : theBackDets;

  int close = 0;
  auto closeDist = std::abs((myDets.front()->toLocal(startPos)).x());
  for (unsigned int i = 1; i < myDets.size(); i++) {
    auto dist = std::abs((myDets[i]->toLocal(startPos)).x());
    if (dist < closeDist) {
      close = i;
      closeDist = dist;
    }
  }
  return close;
}

void CompositeTECWedge::groupedCompatibleDetsV ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop, const MeasurementEstimator &  est, std::vector< DetGroup > &  result  ) const

overridevirtual

Reimplemented from GeometricSearchDet.

Definition at line 78 of file CompositeTECWedge.cc.

References addClosest(), SubLayerCrossings::closest(), SubLayerCrossings::closestIndex(), computeCrossings(), tkDetUtil::computeWindowSize(), LayerCrossingSide::endcapSide(), SubLayerCrossings::isValid(), LogDebug, eostools::move(), DetGroupMerger::orderAndMergeTwoLevels(), SubLayerCrossings::other(), Propagator::propagationDirection(), searchNeighbors(), and svgfig::window().

Referenced by basicComponents().

                                                                                  {
  SubLayerCrossings crossings = computeCrossings(tsos, prop.propagationDirection());
  if (!crossings.isValid())
    return;

  std::vector<DetGroup> closestResult;
  addClosest(tsos, prop, est, crossings.closest(), closestResult);

  LogDebug("TkDetLayers") << "in CompositeTECWedge::groupedCompatibleDets,closestResult.size(): "
                          << closestResult.size();

  if (closestResult.empty())
    return;

  DetGroupElement closestGel(closestResult.front().front());
  auto window = tkDetUtil::computeWindowSize(closestGel.det(), closestGel.trajectoryState(), est);

  searchNeighbors(tsos, prop, est, crossings.closest(), window, closestResult, false);

  std::vector<DetGroup> nextResult;
  searchNeighbors(tsos, prop, est, crossings.other(), window, nextResult, true);

  int crossingSide = LayerCrossingSide().endcapSide(closestGel.trajectoryState(), prop);
  DetGroupMerger::orderAndMergeTwoLevels(
      std::move(closestResult), std::move(nextResult), result, crossings.closestIndex(), crossingSide);
}

void CompositeTECWedge::searchNeighbors ( const TrajectoryStateOnSurface &  tsos, const Propagator &  prop, const MeasurementEstimator &  est, const SubLayerCrossing &  crossing, float  window, std::vector< DetGroup > &  result, bool  checkClosest  ) const

private

Definition at line 177 of file CompositeTECWedge.cc.

References PVValHelper::add(), PV3DBase< T, PVType, FrameType >::barePhi(), SubLayerCrossing::closestDetIndex(), tkDetUtil::overlapInPhi(), Geom::phiLess(), SubLayerCrossing::position(), SubLayerCrossing::subLayerIndex(), and subWedge().

Referenced by basicComponents(), and groupedCompatibleDetsV().

                                                                 {
  const GlobalPoint& gCrossingPos = crossing.position();

  const vector<const GeomDet*>& sWedge(subWedge(crossing.subLayerIndex()));

  int closestIndex = crossing.closestDetIndex();
  int negStartIndex = closestIndex - 1;
  int posStartIndex = closestIndex + 1;

  if (checkClosest) {  // must decide if the closest is on the neg or pos side
    if (Geom::phiLess(gCrossingPos.barePhi(), sWedge[closestIndex]->surface().phi())) {
      posStartIndex = closestIndex;
    } else {
      negStartIndex = closestIndex;
    }
  }

  typedef CompatibleDetToGroupAdder Adder;
  for (int idet = negStartIndex; idet >= 0; idet--) {
    //if(idet <0 || idet>=sWedge.size()) {edm::LogInfo(TkDetLayers) << "==== warning! gone out vector bounds.idet: " << idet ;break;}
    if (!tkDetUtil::overlapInPhi(gCrossingPos, *sWedge[idet], window))
      break;
    if (!Adder::add(*sWedge[idet], tsos, prop, est, result))
      break;
  }
  for (int idet = posStartIndex; idet < static_cast<int>(sWedge.size()); idet++) {
    //if(idet <0 || idet>=sWedge.size()) {edm::LogInfo(TkDetLayers) << "==== warning! gone out vector bounds.idet: " << idet ;break;}
    if (!tkDetUtil::overlapInPhi(gCrossingPos, *sWedge[idet], window))
      break;
    if (!Adder::add(*sWedge[idet], tsos, prop, est, result))
      break;
  }
}

const std::vector<const GeomDet*>& CompositeTECWedge::subWedge ( int  ind ) const

inlineprivate

Definition at line 55 of file CompositeTECWedge.h.

References theBackDets, and theFrontDets.

Referenced by addClosest(), and searchNeighbors().

{ return (ind == 0 ? theFrontDets : theBackDets); }

Member Data Documentation

std::vector<const GeomDet*> CompositeTECWedge::theBackDets

private

Definition at line 59 of file CompositeTECWedge.h.

Referenced by CompositeTECWedge(), computeCrossings(), findClosestDet(), and subWedge().

ReferenceCountingPointer<BoundDiskSector> CompositeTECWedge::theBackSector

private

Definition at line 63 of file CompositeTECWedge.h.

Referenced by CompositeTECWedge(), and computeCrossings().

std::vector<const GeomDet*> CompositeTECWedge::theDets

private

Definition at line 60 of file CompositeTECWedge.h.

Referenced by basicComponents(), and CompositeTECWedge().

std::vector<const GeomDet*> CompositeTECWedge::theFrontDets

private

Definition at line 58 of file CompositeTECWedge.h.

Referenced by CompositeTECWedge(), computeCrossings(), findClosestDet(), and subWedge().

ReferenceCountingPointer<BoundDiskSector> CompositeTECWedge::theFrontSector

private

Definition at line 62 of file CompositeTECWedge.h.

Referenced by CompositeTECWedge(), and computeCrossings().