#include <MuonHitsChamberResidual.h>

Inheritance diagram for MuonHitsChamberResidual:

Public Member Functions
	MuonHitsChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const AlignableDetOrUnitPtr &chamberAlignable)

void	segment_fit ()

Public Member Functions inherited from MuonChamberResidual
virtual void	addResidual (const TrajectoryStateOnSurface , const TransientTrackingRecHit )=0

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 void	setSegmentResidual (const reco::MuonChamberMatch , const reco::MuonSegmentMatch )=0

virtual double	signConvention () const

double	trackdxdz () const

double	trackdydz () const

double	trackx () const

double	tracky () const

int	type () const

virtual	~MuonChamberResidual ()

Protected Attributes
double	m_hitx_1

double	m_hitx_x

double	m_hitx_xx

double	m_hitx_xy

double	m_hitx_y

double	m_hity_1

double	m_hity_x

double	m_hity_xx

double	m_hity_xy

double	m_hity_y

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

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

Additional Inherited Members
Public Types inherited from MuonChamberResidual
enum	{ kDT13, kDT2, kCSC }

Detailed Description

Second level abstraction class for muon chamber residulas: for alignment using individual rechits it implements linear segment fit of hits.

Id:

Author: V. Khotilovich - Texas A&M University khoti.nosp@m.lov@.nosp@m.cern..nosp@m.ch

Id:

Author: V. Khotilovich - Texas A&M University khoti.nosp@m.lov@.nosp@m.cern..nosp@m.ch

Definition at line 22 of file MuonHitsChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 8 of file MuonHitsChamberResidual.cc.

   : MuonChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
   , m_residual_1(0.)
   , m_residual_x(0.)  , m_residual_y(0.)
   , m_residual_xx(0.) , m_residual_xy(0.)
   , m_trackx_1(0.)
   , m_trackx_x(0.)  , m_trackx_y(0.)
   , m_trackx_xx(0.) , m_trackx_xy(0.)
   , m_tracky_1(0.)
   , m_tracky_x(0.)  , m_tracky_y(0.)
   , m_tracky_xx(0.) , m_tracky_xy(0.)
   , m_hitx_1(0.)
   , m_hitx_x(0.)  , m_hitx_y(0.)
   , m_hitx_xx(0.) , m_hitx_xy(0.)
   , m_hity_1(0.)
   , m_hity_x(0.)  , m_hity_y(0.)
   , m_hity_xx(0.) , m_hity_xy(0.)
 {}

Member Function Documentation

void MuonHitsChamberResidual::segment_fit ( )

Definition at line 31 of file MuonHitsChamberResidual.cc.

Referenced by MuonCSCChamberResidual::addResidual(), MuonDT2ChamberResidual::addResidual(), and MuonDT13ChamberResidual::addResidual().

 {
   assert(m_numHits > 1);
   
   double delta = m_residual_1 * m_residual_xx - m_residual_x * m_residual_x;
   m_residual = (m_residual_xx * m_residual_y - m_residual_x * m_residual_xy) / delta;
   
   delta = m_residual_1 * m_residual_xx - m_residual_x * m_residual_x;
   m_residual_error = sqrt(m_residual_xx / delta);
   
   delta = m_residual_1 * m_residual_xx - m_residual_x * m_residual_x;
   m_resslope = (m_residual_1 * m_residual_xy - m_residual_x * m_residual_y) / delta;
   
   delta = m_residual_1 * m_residual_xx - m_residual_x * m_residual_x;
   m_resslope_error = sqrt(m_residual_1 / delta);
   
   m_ndof = m_individual_x.size() - 2;
   
   m_chi2 = 0.;
   double a = m_residual;
   double b = m_resslope;
   std::vector< double >::const_iterator x = m_individual_x.begin();
   std::vector< double >::const_iterator y = m_individual_y.begin();
   std::vector< double >::const_iterator w = m_individual_weight.begin();
   for (; x != m_individual_x.end(); ++x, ++y, ++w)   m_chi2 += pow((*y) - a - b * (*x), 2) * (*w);
 
   delta = m_trackx_1 * m_trackx_xx - m_trackx_x * m_trackx_x;
   m_trackdxdz = (m_trackx_1 * m_trackx_xy - m_trackx_x * m_trackx_y) / delta;
   
   delta = m_tracky_1 * m_tracky_xx - m_tracky_x * m_tracky_x;
   m_trackdydz = (m_tracky_1 * m_tracky_xy - m_tracky_x * m_tracky_y) / delta;
   
   delta = m_trackx_1 * m_trackx_xx - m_trackx_x * m_trackx_x;
   m_trackx = (m_trackx_xx * m_trackx_y - m_trackx_x * m_trackx_xy) / delta;
   
   delta = m_tracky_1 * m_tracky_xx - m_tracky_x * m_tracky_x;
   m_tracky = (m_tracky_xx * m_tracky_y - m_tracky_x * m_tracky_xy) / delta;
 
   delta = m_hitx_1 * m_hitx_xx - m_hitx_x * m_hitx_x;
   m_segdxdz = (m_hitx_1 * m_hitx_xy - m_hitx_x * m_hitx_y) / delta;
 
   delta = m_hity_1 * m_hity_xx - m_hity_x * m_hity_x;
   m_segdydz = (m_hity_1 * m_hity_xy - m_hity_x * m_hity_y) / delta;
 
   delta = m_hitx_1 * m_hitx_xx - m_hitx_x * m_hitx_x;
   m_segx = (m_hitx_xx * m_hitx_y - m_hitx_x * m_hitx_xy) / delta;
 
   delta = m_hity_1 * m_hity_xx - m_hity_x * m_hity_x;
   m_segy = (m_hity_xx * m_hity_y - m_hity_x * m_hity_xy) / delta;
 }