MuonDT13ChamberResidual Class Reference

#include <MuonDT13ChamberResidual.h>

Inheritance diagram for MuonDT13ChamberResidual:

Public Member Functions
void	addResidual (const TrajectoryStateOnSurface *tsos, const TransientTrackingRecHit *hit)
	MuonDT13ChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
double	signConvention (const unsigned int rawId=0) const
int	type () const
Public Member Functions inherited from MuonChamberResidual
AlignableDetOrUnitPtr	chamberAlignable () const
DetId	chamberId () const
double	chi2 () const
double	global_hitresid (int i) const
double	global_residual () const
double	global_resslope () const
GlobalPoint	global_trackpos ()
int	hitlayer (int i) const
double	hitposition (int i) const
double	hitresid (int i) const
DetId	localid (int i) const
	MuonChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
int	ndof () const
int	numHits () const
double	residual () const
double	residual_error () const
double	resslope () const
double	resslope_error () const
double	trackdxdz () const
double	trackdydz () const
double	trackx () const
double	tracky () const
virtual	~MuonChamberResidual ()

Additional Inherited Members
Public Types inherited from MuonChamberResidual
enum	{ kDT13, kDT2, kCSC }
Protected Attributes inherited from MuonChamberResidual
AlignableDetOrUnitPtr	m_chamberAlignable
DetId	m_chamberId
edm::ESHandle < GlobalTrackingGeometry >	m_globalGeometry
std::vector< double >	m_individual_weight
std::vector< double >	m_individual_x
std::vector< double >	m_individual_y
std::vector< DetId >	m_localIDs
std::vector< double >	m_localResids
AlignableNavigator *	m_navigator
int	m_numHits
double	m_residual_1
double	m_residual_x
double	m_residual_xx
double	m_residual_xy
double	m_residual_y
double	m_trackx_1
double	m_trackx_x
double	m_trackx_xx
double	m_trackx_xy
double	m_trackx_y
double	m_tracky_1
double	m_tracky_x
double	m_tracky_xx
double	m_tracky_xy
double	m_tracky_y

Detailed Description

Date:

2009/02/27 18:58:29

Revision:

1.1

Author

J. Pivarski - Texas A&M University pivar.nosp@m.ski@.nosp@m.physi.nosp@m.cs.t.nosp@m.amu.e.nosp@m.du

Definition at line 15 of file MuonDT13ChamberResidual.h.

Constructor & Destructor Documentation

MuonDT13ChamberResidual::MuonDT13ChamberResidual ( edm::ESHandle< GlobalTrackingGeometry >  globalGeometry, AlignableNavigator *  navigator, DetId  chamberId, AlignableDetOrUnitPtr  chamberAlignable  )

inline

Definition at line 17 of file MuonDT13ChamberResidual.h.

    : MuonChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
  {};

Member Function Documentation

void MuonDT13ChamberResidual::addResidual ( const TrajectoryStateOnSurface *  tsos, const TransientTrackingRecHit *  hit  )

virtual

Implements MuonChamberResidual.

Definition at line 3 of file MuonDT13ChamberResidual.cc.

References TrackingRecHit::geographicalId(), TrackingRecHit::localPosition(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPositionError(), MuonChamberResidual::m_chamberAlignable, MuonChamberResidual::m_globalGeometry, MuonChamberResidual::m_individual_weight, MuonChamberResidual::m_individual_x, MuonChamberResidual::m_individual_y, MuonChamberResidual::m_localIDs, MuonChamberResidual::m_localResids, MuonChamberResidual::m_numHits, MuonChamberResidual::m_residual_1, MuonChamberResidual::m_residual_x, MuonChamberResidual::m_residual_xx, MuonChamberResidual::m_residual_xy, MuonChamberResidual::m_residual_y, MuonChamberResidual::m_trackx_1, MuonChamberResidual::m_trackx_x, MuonChamberResidual::m_trackx_xx, MuonChamberResidual::m_trackx_xy, MuonChamberResidual::m_trackx_y, MuonChamberResidual::m_tracky_1, MuonChamberResidual::m_tracky_x, MuonChamberResidual::m_tracky_xx, MuonChamberResidual::m_tracky_xy, MuonChamberResidual::m_tracky_y, MuonChamberResidual::residual(), Alignable::surface(), AlignableSurface::toLocal(), CommonMethods::weight(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                  {
  DetId id = hit->geographicalId();

  LocalPoint hitChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(hit->localPosition()));
  LocalPoint tsosChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(tsos->localPosition()));

  double residual = tsosChamberPos.x() - hitChamberPos.x();  // residual is track minus hit
  double weight = 1. / hit->localPositionError().xx();  // weight linear fit by hit-only local error
  double layerPosition = tsosChamberPos.z();  // the layer's position in the chamber's coordinate system

  m_numHits++;

  // "x" is the layerPosition, "y" is the residual (this is a linear fit to residual versus layerPosition)
  m_residual_1 += weight;
  m_residual_x += weight * layerPosition;
  m_residual_y += weight * residual;
  m_residual_xx += weight * layerPosition * layerPosition;
  m_residual_xy += weight * layerPosition * residual;

  // "x" is the layerPosition, "y" is chamberx (this is a linear fit to chamberx versus layerPosition)
  m_trackx_1 += weight;
  m_trackx_x += weight * layerPosition;
  m_trackx_y += weight * tsosChamberPos.x();
  m_trackx_xx += weight * layerPosition * layerPosition;
  m_trackx_xy += weight * layerPosition * tsosChamberPos.x();

  // "x" is the layerPosition, "y" is chambery (this is a linear fit to chambery versus layerPosition)
  m_tracky_1 += weight;
  m_tracky_x += weight * layerPosition;
  m_tracky_y += weight * tsosChamberPos.y();
  m_tracky_xx += weight * layerPosition * layerPosition;
  m_tracky_xy += weight * layerPosition * tsosChamberPos.y();

  m_localIDs.push_back(id);
  m_localResids.push_back(tsos->localPosition().x() - hit->localPosition().x());
  m_individual_x.push_back(layerPosition);
  m_individual_y.push_back(residual);
  m_individual_weight.push_back(weight);
}

double MuonDT13ChamberResidual::signConvention ( const unsigned int  rawId = 0 ) const

inlinevirtual

Implements MuonChamberResidual.

Definition at line 28 of file MuonDT13ChamberResidual.h.

References funct::cos(), MuonChamberResidual::m_chamberId, MuonChamberResidual::m_globalGeometry, M_PI, and funct::sin().

                                                          {
    DetId id = m_chamberId;
    if (rawId != 0) id = DetId(rawId);
    double rphiAngle = atan2(m_globalGeometry->idToDet(id)->position().y(), m_globalGeometry->idToDet(id)->position().x()) + M_PI/2.;
    GlobalVector rphiDirection(cos(rphiAngle), sin(rphiAngle), 0.);
    return (m_globalGeometry->idToDet(id)->toLocal(rphiDirection).x() > 0. ? 1. : -1.);
  };

int MuonDT13ChamberResidual::type ( ) const

inlinevirtual

Implements MuonChamberResidual.

Definition at line 21 of file MuonDT13ChamberResidual.h.

References MuonChamberResidual::kDT13.

Referenced by python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), and python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties().

{ return MuonChamberResidual::kDT13; };