DirectMuonNavigation Class Reference

#include <DirectMuonNavigation.h>

Public Member Functions
DirectMuonNavigation *	clone () const
std::vector< const DetLayer * >	compatibleEndcapLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const
	DirectMuonNavigation (const edm::ESHandle< MuonDetLayerGeometry > &)
	DirectMuonNavigation (const edm::ESHandle< MuonDetLayerGeometry > &, const edm::ParameterSet &)
	~DirectMuonNavigation ()

Private Member Functions
bool	checkCompatible (const FreeTrajectoryState &fts, const BarrelDetLayer *) const
bool	checkCompatible (const FreeTrajectoryState &fts, const ForwardDetLayer *) const
void	inOutBackward (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const
void	inOutBarrel (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const
void	inOutForward (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const
void	outInBackward (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const
void	outInBarrel (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const
void	outInForward (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const
bool	outward (const FreeTrajectoryState &fts) const

Private Attributes
float	epsilon_
bool	theBarrelFlag
bool	theEndcapFlag
edm::ESHandle< MuonDetLayerGeometry >	theMuonDetLayerGeometry

Detailed Description

Definition at line 20 of file DirectMuonNavigation.h.

Constructor & Destructor Documentation

DirectMuonNavigation::DirectMuonNavigation

(

const edm::ESHandle< MuonDetLayerGeometry > &

muonLayout

)

Definition at line 19 of file DirectMuonNavigation.cc.

Referenced by clone().

    : theMuonDetLayerGeometry(muonLayout), epsilon_(100.), theEndcapFlag(true), theBarrelFlag(true) {}

DirectMuonNavigation::DirectMuonNavigation	(	const edm::ESHandle< MuonDetLayerGeometry > &	muonLayout,
		const edm::ParameterSet &	par
	)

Definition at line 22 of file DirectMuonNavigation.cc.

    : theMuonDetLayerGeometry(muonLayout),
      epsilon_(100.),
      theEndcapFlag(par.getParameter<bool>("Endcap")),
      theBarrelFlag(par.getParameter<bool>("Barrel")) {}

DirectMuonNavigation::~DirectMuonNavigation ( )

inline

Definition at line 30 of file DirectMuonNavigation.h.

References checkCompatible(), compatibleEndcapLayers(), compatibleLayers(), inOutBackward(), inOutBarrel(), inOutForward(), outInBackward(), outInBarrel(), outInForward(), and outward().

{}

Member Function Documentation

bool DirectMuonNavigation::checkCompatible ( const FreeTrajectoryState &  fts, const BarrelDetLayer *  dl  ) const

private

Definition at line 208 of file DirectMuonNavigation.cc.

References epsilon_, FreeTrajectoryState::momentum(), outward(), PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), CosmicsPD_Skims::radius, submit::rm, slope, BarrelDetLayer::specificSurface(), PV3DBase< T, PVType, FrameType >::z(), and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by inOutBackward(), inOutBarrel(), inOutForward(), outInBackward(), outInBarrel(), outInForward(), and ~DirectMuonNavigation().

                                                                                                         {
  float z0 = fts.position().z();
  float r0 = fts.position().perp();
  float zm = fts.momentum().z();
  float rm = fts.momentum().perp();
  float slope = zm / rm;
  if (!outward(fts))
    slope = -slope;
  const BoundCylinder& bc = dl->specificSurface();
  float radius = bc.radius();
  float length = bc.bounds().length() / 2.;

  float z1 = slope * (radius - r0) + z0;
  return (fabs(z1) <= fabs(length) + epsilon_);
}

bool DirectMuonNavigation::checkCompatible ( const FreeTrajectoryState &  fts, const ForwardDetLayer *  dl  ) const

private

Definition at line 224 of file DirectMuonNavigation.cc.

References epsilon_, FreeTrajectoryState::momentum(), outward(), PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), diffTwoXMLs::r1, submit::rm, slope, ForwardDetLayer::specificSurface(), z, PV3DBase< T, PVType, FrameType >::z(), and HLTMuonOfflineAnalyzer_cfi::z0.

                                                                                                          {
  float z0 = fts.position().z();
  float r0 = fts.position().perp();
  float zm = fts.momentum().z();
  float rm = fts.momentum().perp();
  float slope = rm / zm;

  if (!outward(fts))
    slope = -slope;

  const BoundDisk& bd = dl->specificSurface();

  float outRadius = bd.outerRadius();
  float inRadius = bd.innerRadius();
  float z = bd.position().z();

  float r1 = slope * (z - z0) + r0;
  return (r1 >= inRadius - epsilon_ && r1 <= outRadius + epsilon_);
}

DirectMuonNavigation* DirectMuonNavigation::clone ( void  ) const

inline

Definition at line 27 of file DirectMuonNavigation.h.

References DirectMuonNavigation().

{ return new DirectMuonNavigation(*this); }

vector< const DetLayer * > DirectMuonNavigation::compatibleEndcapLayers	(	const FreeTrajectoryState &	fts,
		PropagationDirection	timeDirection
	)		const

Definition at line 93 of file DirectMuonNavigation.cc.

References alongMomentum, inOutForward(), FreeTrajectoryState::momentum(), oppositeToMomentum, outInBackward(), convertSQLitetoXML_cfg::output, groupFilesInBlocks::reverse, and PV3DBase< T, PVType, FrameType >::z().

Referenced by CosmicMuonTrajectoryBuilder::build(), CosmicMuonTrajectoryBuilder::trajectories(), and ~DirectMuonNavigation().

                                                                                                     {
  float zm = fts.momentum().z();

  vector<const DetLayer*> output;

  // collect all endcap layers on 2 sides
  outInBackward(fts, output);
  inOutForward(fts, output);

  // check direction FTS to get a correct order of DetLayers
  if ((zm > 0 && dir == oppositeToMomentum) || (zm < 0 && dir == alongMomentum))
    std::reverse(output.begin(), output.end());

  return output;
}

vector< const DetLayer * > DirectMuonNavigation::compatibleLayers	(	const FreeTrajectoryState &	fts,
		PropagationDirection	timeDirection
	)		const

Definition at line 30 of file DirectMuonNavigation.cc.

References inOutBackward(), inOutBarrel(), inOutForward(), FreeTrajectoryState::momentum(), oppositeToMomentum, outInBackward(), outInBarrel(), outInForward(), convertSQLitetoXML_cfg::output, outward(), FreeTrajectoryState::position(), groupFilesInBlocks::reverse, theBarrelFlag, theEndcapFlag, PV3DBase< T, PVType, FrameType >::z(), and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by CosmicMuonTrajectoryBuilder::build(), DTChamberEfficiency::compatibleLayers(), StandAloneMuonFilter::compatibleLayers(), StandAloneMuonTrajectoryBuilder::propagateTheSeedTSOS(), CosmicMuonTrajectoryBuilder::trajectories(), and ~DirectMuonNavigation().

                                                                                               {
  float z0 = fts.position().z();
  float zm = fts.momentum().z();

  bool inOut = outward(fts);

  vector<const DetLayer*> output;

  // check direction and position of FTS to get a correct order of DetLayers

  if (inOut) {
    if ((zm * z0) >= 0) {
      if (theBarrelFlag)
        inOutBarrel(fts, output);
      if (theEndcapFlag) {
        if (z0 >= 0)
          inOutForward(fts, output);
        else
          inOutBackward(fts, output);
      }
    } else {
      if (theEndcapFlag) {
        if (z0 >= 0)
          outInForward(fts, output);
        else
          outInBackward(fts, output);
      }
      if (theBarrelFlag)
        inOutBarrel(fts, output);
    }
  } else {
    if ((zm * z0) >= 0) {
      if (theBarrelFlag)
        outInBarrel(fts, output);
      if (theEndcapFlag) {
        if (z0 >= 0)
          inOutForward(fts, output);
        else
          inOutBackward(fts, output);
      }
    } else {
      if (theEndcapFlag) {
        if (z0 >= 0)
          outInForward(fts, output);
        else
          outInBackward(fts, output);
      }
      if (theBarrelFlag)
        outInBarrel(fts, output);
    }
  }

  if (dir == oppositeToMomentum)
    std::reverse(output.begin(), output.end());

  return output;
}

void DirectMuonNavigation::inOutBackward ( const FreeTrajectoryState &  fts, std::vector< const DetLayer * > &  output  ) const

private

Definition at line 176 of file DirectMuonNavigation.cc.

References MuonDetLayerGeometry::allBackwardLayers(), checkCompatible(), generateEDF::cont, and theMuonDetLayerGeometry.

Referenced by compatibleLayers(), and ~DirectMuonNavigation().

                                                                                                              {
  bool cont = false;
  const vector<const DetLayer*>& backward = theMuonDetLayerGeometry->allBackwardLayers();

  for (vector<const DetLayer*>::const_iterator iter_E = backward.begin(); iter_E != backward.end(); iter_E++) {
    if (cont)
      output.push_back((*iter_E));
    else if (checkCompatible(fts, dynamic_cast<const ForwardDetLayer*>(*iter_E))) {
      output.push_back((*iter_E));
      cont = true;
    }
  }
}

void DirectMuonNavigation::inOutBarrel ( const FreeTrajectoryState &  fts, std::vector< const DetLayer * > &  output  ) const

private

Definition at line 110 of file DirectMuonNavigation.cc.

References MuonDetLayerGeometry::allBarrelLayers(), Reference_intrackfit_cff::barrel, checkCompatible(), generateEDF::cont, and theMuonDetLayerGeometry.

Referenced by compatibleLayers(), and ~DirectMuonNavigation().

                                                                                                            {
  bool cont = false;
  const vector<const DetLayer*>& barrel = theMuonDetLayerGeometry->allBarrelLayers();

  for (vector<const DetLayer*>::const_iterator iter_B = barrel.begin(); iter_B != barrel.end(); iter_B++) {
    if (cont)
      output.push_back((*iter_B));
    else if (checkCompatible(fts, dynamic_cast<const BarrelDetLayer*>(*iter_B))) {
      output.push_back((*iter_B));
      cont = true;
    }
  }
}

void DirectMuonNavigation::inOutForward ( const FreeTrajectoryState &  fts, std::vector< const DetLayer * > &  output  ) const

private

Definition at line 144 of file DirectMuonNavigation.cc.

References MuonDetLayerGeometry::allForwardLayers(), checkCompatible(), generateEDF::cont, and theMuonDetLayerGeometry.

Referenced by compatibleEndcapLayers(), compatibleLayers(), and ~DirectMuonNavigation().

                                                                                                             {
  const vector<const DetLayer*>& forward = theMuonDetLayerGeometry->allForwardLayers();
  bool cont = false;
  for (vector<const DetLayer*>::const_iterator iter_E = forward.begin(); iter_E != forward.end(); iter_E++) {
    if (cont)
      output.push_back((*iter_E));
    else if (checkCompatible(fts, dynamic_cast<const ForwardDetLayer*>(*iter_E))) {
      output.push_back((*iter_E));
      cont = true;
    }
  }
}

void DirectMuonNavigation::outInBackward ( const FreeTrajectoryState &  fts, std::vector< const DetLayer * > &  output  ) const

private

Definition at line 190 of file DirectMuonNavigation.cc.

References MuonDetLayerGeometry::allBackwardLayers(), checkCompatible(), generateEDF::cont, and theMuonDetLayerGeometry.

Referenced by compatibleEndcapLayers(), compatibleLayers(), and ~DirectMuonNavigation().

                                                                                                              {
  bool cont = false;
  const vector<const DetLayer*>& backward = theMuonDetLayerGeometry->allBackwardLayers();

  vector<const DetLayer*>::const_iterator rbegin = backward.end();
  rbegin--;
  vector<const DetLayer*>::const_iterator rend = backward.begin();
  rend--;
  for (vector<const DetLayer*>::const_iterator iter_E = rbegin; iter_E != rend; iter_E--) {
    if (cont)
      output.push_back((*iter_E));
    else if (checkCompatible(fts, dynamic_cast<const ForwardDetLayer*>(*iter_E))) {
      output.push_back((*iter_E));
      cont = true;
    }
  }
}

void DirectMuonNavigation::outInBarrel ( const FreeTrajectoryState &  fts, std::vector< const DetLayer * > &  output  ) const

private

Definition at line 124 of file DirectMuonNavigation.cc.

References MuonDetLayerGeometry::allBarrelLayers(), Reference_intrackfit_cff::barrel, checkCompatible(), generateEDF::cont, and theMuonDetLayerGeometry.

Referenced by compatibleLayers(), and ~DirectMuonNavigation().

                                                                                                            {
  // default barrel layers are in out, reverse order
  const vector<const DetLayer*>& barrel = theMuonDetLayerGeometry->allBarrelLayers();

  bool cont = false;
  vector<const DetLayer*>::const_iterator rbegin = barrel.end();
  rbegin--;
  vector<const DetLayer*>::const_iterator rend = barrel.begin();
  rend--;

  for (vector<const DetLayer*>::const_iterator iter_B = rbegin; iter_B != rend; iter_B--) {
    if (cont)
      output.push_back((*iter_B));
    else if (checkCompatible(fts, dynamic_cast<const BarrelDetLayer*>(*iter_B))) {
      output.push_back((*iter_B));
      cont = true;
    }
  }
}

void DirectMuonNavigation::outInForward ( const FreeTrajectoryState &  fts, std::vector< const DetLayer * > &  output  ) const

private

Definition at line 157 of file DirectMuonNavigation.cc.

References MuonDetLayerGeometry::allForwardLayers(), checkCompatible(), generateEDF::cont, and theMuonDetLayerGeometry.

Referenced by compatibleLayers(), and ~DirectMuonNavigation().

                                                                                                             {
  // default forward layers are in out, reverse order

  bool cont = false;
  const vector<const DetLayer*>& forward = theMuonDetLayerGeometry->allForwardLayers();
  vector<const DetLayer*>::const_iterator rbegin = forward.end();
  rbegin--;
  vector<const DetLayer*>::const_iterator rend = forward.begin();
  rend--;
  for (vector<const DetLayer*>::const_iterator iter_E = rbegin; iter_E != rend; iter_E--) {
    if (cont)
      output.push_back((*iter_E));
    else if (checkCompatible(fts, dynamic_cast<const ForwardDetLayer*>(*iter_E))) {
      output.push_back((*iter_E));
      cont = true;
    }
  }
}

bool DirectMuonNavigation::outward ( const FreeTrajectoryState &  fts ) const

private

Definition at line 244 of file DirectMuonNavigation.cc.

References FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by checkCompatible(), compatibleLayers(), and ~DirectMuonNavigation().

                                                                       {
  //  return (fts.position().basicVector().dot(fts.momentum().basicVector())>0);

  float x0 = fts.position().x();
  float y0 = fts.position().y();

  float xm = fts.momentum().x();
  float ym = fts.momentum().y();

  return ((x0 * xm + y0 * ym) > 0);
}

Member Data Documentation

float DirectMuonNavigation::epsilon_

private

Definition at line 53 of file DirectMuonNavigation.h.

Referenced by checkCompatible().

bool DirectMuonNavigation::theBarrelFlag

private

Definition at line 55 of file DirectMuonNavigation.h.

Referenced by compatibleLayers().

bool DirectMuonNavigation::theEndcapFlag

private

Definition at line 54 of file DirectMuonNavigation.h.

Referenced by compatibleLayers().

edm::ESHandle<MuonDetLayerGeometry> DirectMuonNavigation::theMuonDetLayerGeometry

private

Definition at line 52 of file DirectMuonNavigation.h.

Referenced by inOutBackward(), inOutBarrel(), inOutForward(), outInBackward(), outInBarrel(), and outInForward().