MuonTrackDT13ChamberResidual Class Reference

#include <MuonTrackDT13ChamberResidual.h>

Inheritance diagram for MuonTrackDT13ChamberResidual:

Public Member Functions
virtual void	addResidual (const TrajectoryStateOnSurface *tsos, const TransientTrackingRecHit *hit)
	MuonTrackDT13ChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const AlignableDetOrUnitPtr &chamberAlignable)
virtual void	setSegmentResidual (const reco::MuonChamberMatch *, const reco::MuonSegmentMatch *)
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
align::GlobalPoint	global_stubpos ()
align::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, const AlignableDetOrUnitPtr &chamberAlignable)
int	ndof () const
int	numHits () const
double	residual () const
double	residual_error () const
double	resslope () const
double	resslope_error () const
double	segdxdz () const
double	segdydz () const
double	segx () const
double	segy () const
virtual double	signConvention () const
double	trackdxdz () const
double	trackdydz () const
double	trackx () const
double	tracky () const
int	type () 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
double	m_chi2
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_ndof
int	m_numHits
double	m_residual
double	m_residual_error
double	m_resslope
double	m_resslope_error
double	m_segdxdz
double	m_segdydz
double	m_segx
double	m_segy
double	m_sign
double	m_trackdxdz
double	m_trackdydz
double	m_trackx
double	m_tracky
int	m_type

Detailed Description

Implementation of tracker muon chamber residuals for axial DT layers

Id:

Definition at line 13 of file MuonTrackDT13ChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 8 of file MuonTrackDT13ChamberResidual.cc.

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

  : MuonChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
{
  m_type = MuonChamberResidual::kDT13;
  double rphiAngle = atan2(m_globalGeometry->idToDet(m_chamberId)->position().y(), m_globalGeometry->idToDet(m_chamberId)->position().x()) + M_PI/2.;
  align::GlobalVector rphiDirection(cos(rphiAngle), sin(rphiAngle), 0.);
  m_sign = m_globalGeometry->idToDet(m_chamberId)->toLocal(rphiDirection).x() > 0. ? 1. : -1.;
}

Member Function Documentation

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

inlinevirtual

Implements MuonChamberResidual.

Definition at line 20 of file MuonTrackDT13ChamberResidual.h.

{}

void MuonTrackDT13ChamberResidual::setSegmentResidual ( const reco::MuonChamberMatch *  trk, const reco::MuonSegmentMatch *  seg  )

virtual

Implements MuonChamberResidual.

Definition at line 19 of file MuonTrackDT13ChamberResidual.cc.

References assert(), reco::MuonSegmentMatch::dtSegmentRef, reco::MuonChamberMatch::dXdZ, reco::MuonSegmentMatch::dXdZ, reco::MuonChamberMatch::dXdZErr, reco::MuonSegmentMatch::dXdZErr, reco::MuonChamberMatch::dYdZ, reco::MuonSegmentMatch::dYdZ, edm::Ref< C, T, F >::get(), DTRecSegment4D::hasPhi(), MuonChamberResidual::m_chi2, MuonChamberResidual::m_ndof, MuonChamberResidual::m_numHits, MuonChamberResidual::m_residual, MuonChamberResidual::m_residual_error, MuonChamberResidual::m_resslope, MuonChamberResidual::m_resslope_error, MuonChamberResidual::m_segdxdz, MuonChamberResidual::m_segdydz, MuonChamberResidual::m_segx, MuonChamberResidual::m_segy, MuonChamberResidual::m_trackdxdz, MuonChamberResidual::m_trackdydz, MuonChamberResidual::m_trackx, MuonChamberResidual::m_tracky, DTRecSegment4D::phiSegment(), funct::pow(), mathSSE::sqrt(), reco::MuonChamberMatch::x, reco::MuonSegmentMatch::x, reco::MuonChamberMatch::xErr, reco::MuonSegmentMatch::xErr, reco::MuonChamberMatch::y, and reco::MuonSegmentMatch::y.

{
  DTRecSegment4DRef segmentDT = seg->dtSegmentRef;
  if (segmentDT.get() != 0)
  {
    const DTRecSegment4D* segment = segmentDT.get();
    assert(segment->hasPhi());
    const DTChamberRecSegment2D* phiSeg = segment->phiSegment();
    m_numHits = phiSeg->recHits().size();
    m_ndof = phiSeg->degreesOfFreedom();
    m_chi2 = phiSeg->chi2();
    //std::cout<<"phi seg position = "<<phiSeg->localPosition()<<"  numhits="<<m_numHits<<std::endl;
  }

  m_residual = trk->x - seg->x;
  m_residual_error = sqrt( pow(trk->xErr, 2) + pow(seg->xErr, 2) );
  m_resslope = trk->dXdZ - seg->dXdZ;
  m_resslope_error = sqrt( pow(trk->dXdZErr, 2) + pow(seg->dXdZErr, 2) );
  
  m_trackx = trk->x;
  m_tracky = trk->y;
  m_trackdxdz = trk->dXdZ;
  m_trackdydz = trk->dYdZ;

  m_segx = seg->x;
  m_segy = seg->y;
  m_segdxdz = seg->dXdZ;
  m_segdydz = seg->dYdZ;
  
  //std::cout<<"d13 res "<<m_residual<<"+-"<<m_residual_error<<"  "<<m_resslope<<"+-"<<m_resslope_error<<std::endl;
  //std::cout<<"d13 trk "<<m_trackx<<" "<<m_tracky<<" "<<m_trackdxdz<<" "<<m_trackdydz<<std::endl;
  //std::cout<<"d13 seg "<<m_segx<<" "<<m_segy<<" "<<m_segdxdz<<" "<<m_segdydz<<std::endl;
}