MTDRingForwardDoubleLayer Class Reference

#include <MTDRingForwardDoubleLayer.h>

Inheritance diagram for MTDRingForwardDoubleLayer:

Public Member Functions
const MTDRingForwardLayer *	backLayer () const
const std::vector< const GeomDet * > &	basicComponents () const override
std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &, const Propagator &, const MeasurementEstimator &) const override
std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const override
const std::vector< const GeometricSearchDet * > &	components () const override
const MTDRingForwardLayer *	frontLayer () const
std::vector< DetGroup >	groupedCompatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const override
bool	isCrack (const GlobalPoint &gp) const
bool	isInsideOut (const TrajectoryStateOnSurface &tsos) const
	MTDRingForwardDoubleLayer (const std::vector< const ForwardDetRing * > &frontRings, const std::vector< const ForwardDetRing * > &backRings)
	Constructor, takes ownership of pointers. More...
virtual const std::vector < const ForwardDetRing * > &	rings () const
	Return the vector of rings. More...
void	selfTest () const
SubDetector	subDetector () const override
	~MTDRingForwardDoubleLayer () override

Protected Member Functions
BoundDisk *	computeSurface () override

Private Attributes
MTDRingForwardLayer	theBackLayer
std::vector< const GeomDet * >	theBasicComponents
std::vector< const GeometricSearchDet * >	theComponents
MTDRingForwardLayer	theFrontLayer
std::vector< const ForwardDetRing * >	theRings

Detailed Description

A plane composed two layers of disks. The Endcap Timing Layer.

Author

L. Gray

Definition at line 19 of file MTDRingForwardDoubleLayer.h.

Constructor & Destructor Documentation

MTDRingForwardDoubleLayer::MTDRingForwardDoubleLayer	(	const std::vector< const ForwardDetRing * > &	frontRings,
		const std::vector< const ForwardDetRing * > &	backRings
	)

Constructor, takes ownership of pointers.

Definition at line 24 of file MTDRingForwardDoubleLayer.cc.

References basicComponents(), computeSurface(), LogTrace, selfTest(), theBasicComponents, theComponents, and theRings.

    : RingedForwardLayer(true),
      theFrontLayer(frontRings),
      theBackLayer(backRings),
      theRings(frontRings),  // add back later
      theComponents(),
      theBasicComponents() {
  theRings.insert(theRings.end(), backRings.begin(), backRings.end());
  theComponents = std::vector<const GeometricSearchDet*>(theRings.begin(), theRings.end());

  // Cache chamber pointers (the basic components_)
  // and find extension in R and Z
  for (vector<const ForwardDetRing*>::const_iterator it = theRings.begin(); it != theRings.end(); it++) {
    vector<const GeomDet*> tmp2 = (*it)->basicComponents();
    theBasicComponents.insert(theBasicComponents.end(), tmp2.begin(), tmp2.end());
  }

  setSurface(computeSurface());

  LogTrace("MTDDetLayers") << "Constructing MTDRingForwardDoubleLayer: " << basicComponents().size() << " Dets "
                           << theRings.size() << " Rings "
                           << " Z: " << specificSurface().position().z() << " R1: " << specificSurface().innerRadius()
                           << " R2: " << specificSurface().outerRadius();

  selfTest();
}

MTDRingForwardDoubleLayer::~MTDRingForwardDoubleLayer ( )

inlineoverride

Definition at line 25 of file MTDRingForwardDoubleLayer.h.

{}

Member Function Documentation

const MTDRingForwardLayer* MTDRingForwardDoubleLayer::backLayer ( ) const

inline

Definition at line 59 of file MTDRingForwardDoubleLayer.h.

References theBackLayer.

{ return &theBackLayer; }

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

inlineoverride

Definition at line 29 of file MTDRingForwardDoubleLayer.h.

References theBasicComponents.

Referenced by MTDRingForwardDoubleLayer().

{ return theBasicComponents; }

std::pair< bool, TrajectoryStateOnSurface > MTDRingForwardDoubleLayer::compatible ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop, const MeasurementEstimator &  est  ) const

override

Definition at line 75 of file MTDRingForwardDoubleLayer.cc.

References HLT_FULL_cff::deltaR, submitPVValidationJobs::err, TrajectoryStateOnSurface::hasError(), SimpleDiskBounds::inside(), insideOut, isInsideOut(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), LogTrace, HLT_FULL_cff::nSigma, LocalTrajectoryError::positionError(), Propagator::propagate(), mathSSE::sqrt(), funct::tan(), theBackLayer, theFrontLayer, PV3DBase< T, PVType, FrameType >::theta(), createJobs::tmp, LocalError::xx(), and LocalError::yy().

Referenced by compatibleDets().

                                                                                                                  {
  // mostly copied from ForwardDetLayer, except propagates to closest surface,
  // not to center

  bool insideOut = isInsideOut(startingState);
  const MTDRingForwardLayer& closerLayer = (insideOut) ? theFrontLayer : theBackLayer;
  LogTrace("MTDDetLayers") << "MTDRingForwardDoubleLayer::compatible is assuming inside-out direction: " << insideOut;

  TrajectoryStateOnSurface myState = prop.propagate(startingState, closerLayer.specificSurface());
  if (!myState.isValid())
    return make_pair(false, myState);

  // take into account the thickness of the layer
  float deltaR = surface().bounds().thickness() / 2. * fabs(tan(myState.localDirection().theta()));

  // take into account the error on the predicted state
  const float nSigma = 3.;
  if (myState.hasError()) {
    LocalError err = myState.localError().positionError();
    // ignore correlation for the moment...
    deltaR += nSigma * sqrt(err.xx() + err.yy());
  }

  float zPos = (zmax() + zmin()) / 2.;
  SimpleDiskBounds tmp(rmin() - deltaR, rmax() + deltaR, zmin() - zPos, zmax() - zPos);

  return make_pair(tmp.inside(myState.localPosition()), myState);
}

vector< GeometricSearchDet::DetWithState > MTDRingForwardDoubleLayer::compatibleDets ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop, const MeasurementEstimator &  est  ) const

override

Definition at line 105 of file MTDRingForwardDoubleLayer.cc.

References compatible(), groupedCompatibleDets(), LogTrace, and mps_fire::result.

                                                                                                                  {
  vector<DetWithState> result;
  pair<bool, TrajectoryStateOnSurface> compat = compatible(startingState, prop, est);

  if (!compat.first) {
    LogTrace("MTDDetLayers") << "     MTDRingForwardDoubleLayer::compatibleDets: not compatible"
                             << " (should not have been selected!)";
    return result;
  }

  TrajectoryStateOnSurface& tsos = compat.second;

  // standard implementation of compatibleDets() for class which have
  // groupedCompatibleDets implemented.
  // This code should be moved in a common place intead of being
  // copied many times.
  vector<DetGroup> vectorGroups = groupedCompatibleDets(tsos, prop, est);
  for (const auto& thisDG : vectorGroups) {
    for (const auto& thisDGE : thisDG) {
      result.emplace_back(DetWithState(thisDGE.det(), thisDGE.trajectoryState()));
    }
  }
  return result;
}

const std::vector<const GeometricSearchDet*>& MTDRingForwardDoubleLayer::components ( ) const

inlineoverride

Definition at line 31 of file MTDRingForwardDoubleLayer.h.

References theComponents.

{ return theComponents; }

BoundDisk * MTDRingForwardDoubleLayer::computeSurface ( )

overrideprotected

Definition at line 52 of file MTDRingForwardDoubleLayer.cc.

References SiStripPI::max, min(), makeMuonMisalignmentScenario::rot, theBackLayer, and theFrontLayer.

Referenced by MTDRingForwardDoubleLayer().

                                                     {
  const BoundDisk& frontDisk = theFrontLayer.specificSurface();
  const BoundDisk& backDisk = theBackLayer.specificSurface();

  float rmin = min(frontDisk.innerRadius(), backDisk.innerRadius());
  float rmax = max(frontDisk.outerRadius(), backDisk.outerRadius());
  float zmin = frontDisk.position().z();
  float halfThickness = frontDisk.bounds().thickness() / 2.;
  zmin = (zmin > 0) ? zmin - halfThickness : zmin + halfThickness;
  float zmax = backDisk.position().z();
  halfThickness = backDisk.bounds().thickness() / 2.;
  zmax = (zmax > 0) ? zmax + halfThickness : zmax - halfThickness;
  float zPos = (zmax + zmin) / 2.;
  PositionType pos(0., 0., zPos);
  RotationType rot;

  return new BoundDisk(pos, rot, new SimpleDiskBounds(rmin, rmax, zmin - zPos, zmax - zPos));
}

const MTDRingForwardLayer* MTDRingForwardDoubleLayer::frontLayer ( ) const

inline

Definition at line 58 of file MTDRingForwardDoubleLayer.h.

References theFrontLayer.

{ return &theFrontLayer; }

vector< DetGroup > MTDRingForwardDoubleLayer::groupedCompatibleDets ( const TrajectoryStateOnSurface &  startingState, const Propagator &  prop, const MeasurementEstimator &  est  ) const

override

Definition at line 131 of file MTDRingForwardDoubleLayer.cc.

References MTDRingForwardLayer::compatibleDets(), LogTrace, mps_fire::result, theBackLayer, and theFrontLayer.

Referenced by compatibleDets().

                                                                                                         {
  vector<GeometricSearchDet::DetWithState> detWithStates1, detWithStates2;

  LogTrace("MTDDetLayers") << "groupedCompatibleDets are currently given always in inside-out order";
  // this should be fixed either in RecoMTD/MeasurementDet/MTDDetLayerMeasurements or
  // RecoMTD/DetLayers/MTDRingForwardDoubleLayer

  detWithStates1 = theFrontLayer.compatibleDets(startingState, prop, est);
  detWithStates2 = theBackLayer.compatibleDets(startingState, prop, est);

  vector<DetGroup> result;
  if (!detWithStates1.empty())
    result.push_back(DetGroup(detWithStates1));
  if (!detWithStates2.empty())
    result.push_back(DetGroup(detWithStates2));
  LogTrace("MTDDetLayers") << "DoubleLayer Compatible dets: " << result.size();
  return result;
}

bool MTDRingForwardDoubleLayer::isCrack

(

const GlobalPoint &

gp

)

const

Definition at line 152 of file MTDRingForwardDoubleLayer.cc.

References cms::cuda::assert(), LogTrace, PV3DBase< T, PVType, FrameType >::perp(), alignCSCRings::r, mps_fire::result, MTDRingForwardLayer::rings(), ForwardDetRing::specificSurface(), and theBackLayer.

                                                                   {
  // approximate
  bool result = false;
  double r = gp.perp();
  const std::vector<const ForwardDetRing*>& backRings = theBackLayer.rings();
  if (backRings.size() > 1) {
    const MTDDetRing* innerRing = dynamic_cast<const MTDDetRing*>(backRings[0]);
    const MTDDetRing* outerRing = dynamic_cast<const MTDDetRing*>(backRings[1]);
    assert(innerRing && outerRing);
    float crackInner = innerRing->specificSurface().outerRadius();
    float crackOuter = outerRing->specificSurface().innerRadius();
    LogTrace("MTDDetLayers") << "In a crack:" << crackInner << " " << r << " " << crackOuter;
    if (r > crackInner && r < crackOuter)
      return true;
  }
  // non-overlapping rings
  return result;
}

bool MTDRingForwardDoubleLayer::isInsideOut

(

const TrajectoryStateOnSurface &

tsos

)

const

Definition at line 71 of file MTDRingForwardDoubleLayer.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), Basic3DVector< T >::dot(), TrajectoryStateOnSurface::globalMomentum(), and TrajectoryStateOnSurface::globalPosition().

Referenced by compatible().

                                                                                      {
  return tsos.globalPosition().basicVector().dot(tsos.globalMomentum().basicVector()) > 0;
}

virtual const std::vector<const ForwardDetRing*>& MTDRingForwardDoubleLayer::rings ( ) const

inlinevirtual

Return the vector of rings.

Definition at line 54 of file MTDRingForwardDoubleLayer.h.

References theRings.

{ return theRings; }

void MTDRingForwardDoubleLayer::selfTest

(

)

const

Definition at line 171 of file MTDRingForwardDoubleLayer.cc.

References funct::abs(), cms::cuda::assert(), MTDRingForwardLayer::basicComponents(), ETLDetId::kETLv1maxRing, theBackLayer, and theFrontLayer.

Referenced by MTDRingForwardDoubleLayer().

                                               {
  const std::vector<const GeomDet*>& frontDets = theFrontLayer.basicComponents();
  const std::vector<const GeomDet*>& backDets = theBackLayer.basicComponents();

  for (int iring = 0; iring <= ETLDetId::kETLv1maxRing; iring++) {
    float frontz(0.);
    float backz(1e3f);
    for (const auto& thisFront : frontDets) {
      if (static_cast<ETLDetId>(thisFront->geographicalId().rawId()).mtdRR() == iring) {
        float tmpz(std::abs(thisFront->surface().position().z()));
        if (tmpz > frontz) {
          frontz = tmpz;
        }
      }
    }
    for (const auto& thisBack : backDets) {
      if (static_cast<ETLDetId>(thisBack->geographicalId().rawId()).mtdRR() == iring) {
        float tmpz(std::abs(thisBack->surface().position().z()));
        if (tmpz < backz) {
          backz = tmpz;
        }
      }
    }
    assert(frontz < backz);
  }
}

SubDetector MTDRingForwardDoubleLayer::subDetector ( ) const

inlineoverride

Definition at line 49 of file MTDRingForwardDoubleLayer.h.

References MTDRingForwardLayer::subDetector(), and theBackLayer.

{ return theBackLayer.subDetector(); }

Member Data Documentation

MTDRingForwardLayer MTDRingForwardDoubleLayer::theBackLayer

private

Definition at line 68 of file MTDRingForwardDoubleLayer.h.

Referenced by backLayer(), compatible(), computeSurface(), groupedCompatibleDets(), isCrack(), selfTest(), and subDetector().

std::vector<const GeomDet*> MTDRingForwardDoubleLayer::theBasicComponents

private

Definition at line 71 of file MTDRingForwardDoubleLayer.h.

Referenced by basicComponents(), and MTDRingForwardDoubleLayer().

std::vector<const GeometricSearchDet*> MTDRingForwardDoubleLayer::theComponents

private

Definition at line 70 of file MTDRingForwardDoubleLayer.h.

Referenced by components(), and MTDRingForwardDoubleLayer().

MTDRingForwardLayer MTDRingForwardDoubleLayer::theFrontLayer

private

Definition at line 67 of file MTDRingForwardDoubleLayer.h.

Referenced by compatible(), computeSurface(), frontLayer(), groupedCompatibleDets(), and selfTest().

std::vector<const ForwardDetRing*> MTDRingForwardDoubleLayer::theRings

private

Definition at line 69 of file MTDRingForwardDoubleLayer.h.

Referenced by MTDRingForwardDoubleLayer(), and rings().