Phase2OTBarrelRod Class Reference

#include <Phase2OTBarrelRod.h>

Inheritance diagram for Phase2OTBarrelRod:

Public Types
typedef GenericBinFinderInZ< float, GeomDet >	BinFinderType
Public Types inherited from GeometricSearchDet
typedef std::pair< const GeomDet *, TrajectoryStateOnSurface >	DetWithState
typedef BoundSurface::PositionType	PositionType
typedef BoundSurface::RotationType	RotationType
typedef TrajectoryStateOnSurface	TrajectoryState

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...
void	groupedCompatibleDetsV (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetGroup > &result) const override __attribute__((hot))
	Phase2OTBarrelRod (std::vector< const GeomDet * > &innerDets, std::vector< const GeomDet * > &outerDets, std::vector< const GeomDet * > &innerDetBrothers, std::vector< const GeomDet * > &outerDetBrothers) __attribute__((cold))
	~Phase2OTBarrelRod () override __attribute__((cold))
Public Member Functions inherited from DetRod
virtual const Plane &	specificSurface () const final
	Return the rod surface as a Plane. More...
const BoundSurface &	surface () const final
	The surface of the GeometricSearchDet. More...
	~DetRod () override
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, std::vector< DetGroup > &brotherresult) const __attribute__((hot))
SubLayerCrossings	computeCrossings (const TrajectoryStateOnSurface &tsos, PropagationDirection propDir) const __attribute__((hot))
float	computeWindowSize (const GeomDet *det, const TrajectoryStateOnSurface &tsos, const MeasurementEstimator &est) const __attribute__((hot))
void	searchNeighbors (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, float window, std::vector< DetGroup > &result, std::vector< DetGroup > &brotherresult, bool checkClosest) const __attribute__((hot))
const std::vector< const GeomDet * > &	subRod (int ind) const
const std::vector< const GeomDet * > &	subRodBrothers (int ind) const

Private Attributes
std::vector< const GeomDet * >	theDets
BinFinderType	theInnerBinFinder
std::vector< const GeomDet * >	theInnerDetBrothers
std::vector< const GeomDet * >	theInnerDets
ReferenceCountingPointer< Plane >	theInnerPlane
BinFinderType	theOuterBinFinder
std::vector< const GeomDet * >	theOuterDetBrothers
std::vector< const GeomDet * >	theOuterDets
ReferenceCountingPointer< Plane >	theOuterPlane

Additional Inherited Members
Protected Member Functions inherited from DetRod
void	setPlane (Plane *plane)
	Set the rod's plane. More...
Protected Attributes inherited from GeometricSearchDet
bool	haveGroups
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A concrete implementation for TOB Rod

Definition at line 16 of file Phase2OTBarrelRod.h.

Member Typedef Documentation

typedef GenericBinFinderInZ<float,GeomDet> Phase2OTBarrelRod::BinFinderType

Definition at line 18 of file Phase2OTBarrelRod.h.

Constructor & Destructor Documentation

Phase2OTBarrelRod::Phase2OTBarrelRod	(	std::vector< const GeomDet * > &	innerDets,
		std::vector< const GeomDet * > &	outerDets,
		std::vector< const GeomDet * > &	innerDetBrothers,
		std::vector< const GeomDet * > &	outerDetBrothers
	)

Definition at line 30 of file Phase2OTBarrelRod.cc.

References mps_fire::i, LogDebug, DetRod::setPlane(), theDets, theInnerBinFinder, theInnerDetBrothers, theInnerDets, theInnerPlane, theOuterBinFinder, theOuterDetBrothers, theOuterDets, and theOuterPlane.

                                                                  :
  DetRod(true),
  theInnerDets(innerDets),theOuterDets(outerDets),theInnerDetBrothers(innerDetBrothers),theOuterDetBrothers(outerDetBrothers)
{
  theDets.assign(theInnerDets.begin(),theInnerDets.end());
  theDets.insert(theDets.end(),theOuterDets.begin(),theOuterDets.end());
  theDets.insert(theDets.end(),theInnerDetBrothers.begin(),theInnerDetBrothers.end());
  theDets.insert(theDets.end(),theOuterDetBrothers.begin(),theOuterDetBrothers.end());


  RodPlaneBuilderFromDet planeBuilder;
  setPlane( planeBuilder( theDets));
  theInnerPlane = planeBuilder( theInnerDets);
  theOuterPlane = planeBuilder( theOuterDets);

  // modules be already sorted in z

  theInnerBinFinder = BinFinderType(theInnerDets.begin(), theInnerDets.end());
  theOuterBinFinder = BinFinderType(theOuterDets.begin(), theOuterDets.end());


 
#ifdef EDM_ML_DEBUG
  LogDebug("TkDetLayers") << "==== DEBUG Phase2OTBarrelRod =====" ; 
  for (vector<const GeomDet*>::const_iterator i=theInnerDets.begin();
       i != theInnerDets.end(); i++){
    LogDebug("TkDetLayers") << "inner Phase2OTBarrelRod's Det pos z,perp,eta,phi: " 
                << (**i).position().z() << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta() << " , " 
                << (**i).position().phi() ;
  }
  
  for (vector<const GeomDet*>::const_iterator i=theInnerDetBrothers.begin();
       i != theInnerDetBrothers.end(); i++){
    LogDebug("TkDetLayers") << "inner Phase2OTBarrelRod's Det Brother pos z,perp,eta,phi: " 
                << (**i).position().z() << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta() << " , " 
                << (**i).position().phi() ;
  }

  for (vector<const GeomDet*>::const_iterator i=theOuterDets.begin();
       i != theOuterDets.end(); i++){
    LogDebug("TkDetLayers") << "outer Phase2OTBarrelRod's Det pos z,perp,eta,phi: " 
                << (**i).position().z() << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta() << " , " 
                << (**i).position().phi() ;
  }
  
  for (vector<const GeomDet*>::const_iterator i=theOuterDetBrothers.begin();
       i != theOuterDetBrothers.end(); i++){
    LogDebug("TkDetLayers") << "outer Phase2OTBarrelRod's Det Brother pos z,perp,eta,phi: " 
                << (**i).position().z() << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta() << " , " 
                << (**i).position().phi() ;
  }
  LogDebug("TkDetLayers") << "==== end DEBUG Phase2OTBarrelRod =====" ; 
#endif  


}

Phase2OTBarrelRod::~Phase2OTBarrelRod ( )

override

Definition at line 98 of file Phase2OTBarrelRod.cc.

                                     {
  
} 

Member Function Documentation

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

private

Definition at line 229 of file Phase2OTBarrelRod.cc.

References CompatibleDetToGroupAdder::add(), SubLayerCrossing::closestDetIndex(), mps_fire::result, SubLayerCrossing::subLayerIndex(), subRod(), and subRodBrothers().

Referenced by basicComponents(), and groupedCompatibleDetsV().

{

  const vector<const GeomDet*>& sRod( subRod( crossing.subLayerIndex()));
  bool firstgroup = CompatibleDetToGroupAdder::add( *sRod[crossing.closestDetIndex()], 
                            tsos, prop, est, result);
  // it assumes that the closestDetIndex is ok also for the brother detectors: the crossing is NOT recomputed
  const vector<const GeomDet*>& sRodBrothers( subRodBrothers( crossing.subLayerIndex()));
  bool brothergroup = CompatibleDetToGroupAdder::add( *sRodBrothers[crossing.closestDetIndex()], 
                              tsos, prop, est, brotherresult);

  return firstgroup || brothergroup;
}

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

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 28 of file Phase2OTBarrelRod.h.

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

{return theDets;}

pair< bool, TrajectoryStateOnSurface > Phase2OTBarrelRod::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 109 of file Phase2OTBarrelRod.cc.

Referenced by basicComponents().

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

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

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 104 of file Phase2OTBarrelRod.cc.

Referenced by basicComponents().

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

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

private

Definition at line 190 of file Phase2OTBarrelRod.cc.

References funct::abs(), GenericBinFinderInZ< T, G >::binIndex(), GenericBinFinderInZ< T, G >::binPosition(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), rho, theInnerBinFinder, theInnerPlane, theOuterBinFinder, theOuterPlane, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by basicComponents(), and groupedCompatibleDetsV().

{
  GlobalPoint startPos( startingState.globalPosition());
  GlobalVector startDir( startingState.globalMomentum());
  double rho( startingState.transverseCurvature());

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


  std::pair<bool,double> outerPath = crossing.pathLength( *theOuterPlane);
  if (!outerPath.first) return SubLayerCrossings();
  GlobalPoint gOuterPoint( crossing.position(outerPath.second));

  std::pair<bool,double> innerPath = crossing.pathLength( *theInnerPlane);
  if (!innerPath.first) return SubLayerCrossings();
  GlobalPoint gInnerPoint( crossing.position(innerPath.second));


  int innerIndex = theInnerBinFinder.binIndex(gInnerPoint.z());
  float innerDist = std::abs( theInnerBinFinder.binPosition(innerIndex) - gInnerPoint.z());
  SubLayerCrossing innerSLC( 0, innerIndex, gInnerPoint);

  int outerIndex = theOuterBinFinder.binIndex(gOuterPoint.z());
  float outerDist = std::abs( theOuterBinFinder.binPosition(outerIndex) - gOuterPoint.z());
  SubLayerCrossing outerSLC( 1, outerIndex, gOuterPoint);

  if (innerDist < outerDist) {
    return SubLayerCrossings( innerSLC, outerSLC, 0);
  }
  else {
    return SubLayerCrossings( outerSLC, innerSLC, 1);
  } 
}

float Phase2OTBarrelRod::computeWindowSize ( const GeomDet *  det, const TrajectoryStateOnSurface &  tsos, const MeasurementEstimator &  est  ) const

private

Definition at line 249 of file Phase2OTBarrelRod.cc.

References funct::abs(), Surface::bounds(), constexpr, Bounds::length(), MeasurementEstimator::maximalLocalDisplacement(), muon::overlap(), GeomDet::surface(), GloballyPositioned< T >::toLocal(), svgfig::window(), y, and PV3DBase< T, PVType, FrameType >::y().

Referenced by basicComponents(), and groupedCompatibleDetsV().

{
  return
    est.maximalLocalDisplacement(tsos, det->surface()).y();
}

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

overridevirtual

Reimplemented from GeometricSearchDet.

Definition at line 120 of file Phase2OTBarrelRod.cc.

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

Referenced by basicComponents().

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

  std::vector<DetGroup> closestResult;
  std::vector<DetGroup> closestBrotherResult;
  addClosest( tsos, prop, est, crossings.closest(), closestResult, closestBrotherResult);
  if (closestResult.empty()){
    std::vector<DetGroup> nextResult;
    std::vector<DetGroup> nextBrotherResult;
    addClosest( tsos, prop, est, crossings.other(), nextResult, nextBrotherResult);
    if(nextResult.empty())    return;

    DetGroupElement nextGel( nextResult.front().front());  
    int crossingSide = LayerCrossingSide().barrelSide( nextGel.trajectoryState(), prop);
    std::vector<DetGroup> closestCompleteResult;
    DetGroupMerger::orderAndMergeTwoLevels(std::move(closestResult),std::move(closestBrotherResult),closestCompleteResult,
                       0, crossingSide);
    std::vector<DetGroup> nextCompleteResult;
    DetGroupMerger::orderAndMergeTwoLevels(std::move(nextResult),std::move(nextBrotherResult),nextCompleteResult,
                       0, crossingSide);

    DetGroupMerger::orderAndMergeTwoLevels( std::move(closestCompleteResult), std::move(nextCompleteResult), result, 
                       crossings.closestIndex(), crossingSide);   
  } else {
  
    DetGroupElement closestGel( closestResult.front().front());
    int crossingSide = LayerCrossingSide().barrelSide( closestGel.trajectoryState(), prop);
    float window = computeWindowSize( closestGel.det(), closestGel.trajectoryState(), est);

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

    std::vector<DetGroup> closestCompleteResult;
    DetGroupMerger::orderAndMergeTwoLevels(std::move(closestResult),std::move(closestBrotherResult),closestCompleteResult,
                       0, crossingSide);

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

    std::vector<DetGroup> nextCompleteResult;
    DetGroupMerger::orderAndMergeTwoLevels(std::move(nextResult),std::move(nextBrotherResult),nextCompleteResult,
                       0, crossingSide);

    DetGroupMerger::orderAndMergeTwoLevels( std::move(closestCompleteResult), std::move(nextCompleteResult), result, 
                        crossings.closestIndex(), crossingSide);
  }

  //due to propagator problems, when we add single pt sub modules, we should order them in r (barrel)
  sort(result.begin(),result.end(),DetGroupElementPerpLess());
  for (auto&  grp : result) {
    if ( grp.empty() )  continue;
    LogTrace("TkDetLayers") <<"New group in Phase2OTBarrelRod made by : ";
    for (auto const & det : grp) {
      LogTrace("TkDetLayers") <<" geom det at r: " << det.det()->position().perp() <<" id:" << det.det()->geographicalId().rawId()
                              <<" tsos at:" << det.trajectoryState().globalPosition();
    }
  }

}

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

private

Definition at line 289 of file Phase2OTBarrelRod.cc.

References PVValHelper::add(), SubLayerCrossing::closestDetIndex(), muon::overlap(), SubLayerCrossing::position(), SubLayerCrossing::subLayerIndex(), subRod(), subRodBrothers(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by basicComponents(), and groupedCompatibleDetsV().

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

  const vector<const GeomDet*>& sRod( subRod( crossing.subLayerIndex()));
  const vector<const GeomDet*>& sBrotherRod( subRodBrothers( 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 (gCrossingPos.z() < sRod[closestIndex]->surface().position().z()) {
      posStartIndex = closestIndex;
    }
    else {
      negStartIndex = closestIndex;
    }
  }

  typedef CompatibleDetToGroupAdder Adder;
  for (int idet=negStartIndex; idet >= 0; idet--) {
    if (!overlap( gCrossingPos, *sRod[idet], window)) break;
    if (!Adder::add( *sRod[idet], tsos, prop, est, result)) break;
    // If the two above checks are passed also the brother module will be added with no further checks
    Adder::add( *sBrotherRod[idet], tsos, prop, est, brotherresult);
  }
  for (int idet=posStartIndex; idet < static_cast<int>(sRod.size()); idet++) {
    if (!overlap( gCrossingPos, *sRod[idet], window)) break;
    if (!Adder::add( *sRod[idet], tsos, prop, est, result)) break;
    // If the two above checks are passed also the brother module will be added with no further checks
    Adder::add( *sBrotherRod[idet], tsos, prop, est, brotherresult);
  }
}

const std::vector<const GeomDet*>& Phase2OTBarrelRod::subRod ( int  ind ) const

inlineprivate

Definition at line 70 of file Phase2OTBarrelRod.h.

References theInnerDets, and theOuterDets.

Referenced by addClosest(), and searchNeighbors().

                                                          {
    return (ind==0 ? theInnerDets : theOuterDets);
  }

const std::vector<const GeomDet*>& Phase2OTBarrelRod::subRodBrothers ( int  ind ) const

inlineprivate

Definition at line 74 of file Phase2OTBarrelRod.h.

References theInnerDetBrothers, and theOuterDetBrothers.

Referenced by addClosest(), and searchNeighbors().

                                                                  {
    return (ind==0 ? theInnerDetBrothers : theOuterDetBrothers);
  }

Member Data Documentation

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

private

Definition at line 80 of file Phase2OTBarrelRod.h.

Referenced by basicComponents(), and Phase2OTBarrelRod().

BinFinderType Phase2OTBarrelRod::theInnerBinFinder

private

Definition at line 89 of file Phase2OTBarrelRod.h.

Referenced by computeCrossings(), and Phase2OTBarrelRod().

std::vector<const GeomDet*> Phase2OTBarrelRod::theInnerDetBrothers

private

Definition at line 83 of file Phase2OTBarrelRod.h.

Referenced by Phase2OTBarrelRod(), and subRodBrothers().

std::vector<const GeomDet*> Phase2OTBarrelRod::theInnerDets

private

Definition at line 81 of file Phase2OTBarrelRod.h.

Referenced by Phase2OTBarrelRod(), and subRod().

ReferenceCountingPointer<Plane> Phase2OTBarrelRod::theInnerPlane

private

Definition at line 86 of file Phase2OTBarrelRod.h.

Referenced by computeCrossings(), and Phase2OTBarrelRod().

BinFinderType Phase2OTBarrelRod::theOuterBinFinder

private

Definition at line 90 of file Phase2OTBarrelRod.h.

Referenced by computeCrossings(), and Phase2OTBarrelRod().

std::vector<const GeomDet*> Phase2OTBarrelRod::theOuterDetBrothers

private

Definition at line 84 of file Phase2OTBarrelRod.h.

Referenced by Phase2OTBarrelRod(), and subRodBrothers().

std::vector<const GeomDet*> Phase2OTBarrelRod::theOuterDets

private

Definition at line 82 of file Phase2OTBarrelRod.h.

Referenced by Phase2OTBarrelRod(), and subRod().

ReferenceCountingPointer<Plane> Phase2OTBarrelRod::theOuterPlane

private

Definition at line 87 of file Phase2OTBarrelRod.h.

Referenced by computeCrossings(), and Phase2OTBarrelRod().