MuonTrackDT13ChamberResidual Class Reference

#include <MuonTrackDT13ChamberResidual.h>

Inheritance diagram for MuonTrackDT13ChamberResidual:

List of all members.

Public Member Functions
virtual void	addResidual (const TrajectoryStateOnSurface *tsos, const TransientTrackingRecHit *hit)
	MuonTrackDT13ChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
virtual void	setSegmentResidual (const reco::MuonChamberMatch *, const reco::MuonSegmentMatch *)

---

Detailed Description

Implementation of tracker muon chamber residuals for axial DT layers

Id:

MuonTrackDT13ChamberResidual.h,v 1.1 2011/10/12 23:32:08 khotilov Exp

Definition at line 13 of file MuonTrackDT13ChamberResidual.h.

---

Constructor & Destructor Documentation

MuonTrackDT13ChamberResidual::MuonTrackDT13ChamberResidual	(	edm::ESHandle< GlobalTrackingGeometry >	globalGeometry,
		AlignableNavigator *	navigator,
		DetId	chamberId,
		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
	)		[inline, virtual]

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 reco::MuonSegmentMatch::dtSegmentRef, reco::MuonSegmentMatch::dXdZ, reco::MuonChamberMatch::dXdZ, reco::MuonChamberMatch::dXdZErr, reco::MuonSegmentMatch::dXdZErr, reco::MuonSegmentMatch::dYdZ, reco::MuonChamberMatch::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::MuonSegmentMatch::x, reco::MuonChamberMatch::x, reco::MuonSegmentMatch::xErr, reco::MuonChamberMatch::xErr, reco::MuonSegmentMatch::y, and reco::MuonChamberMatch::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;
}