#include <MuonHitsChamberResidual.h>

Inheritance diagram for MuonHitsChamberResidual:

Public Member Functions
	MuonHitsChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
void	segment_fit ()
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

Detailed Description

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

Id:: MuonHitsChamberResidual.h,v 1.1 2011/10/12 23:32:07 khotilov Exp

Author:: V. Khotilovich - Texas A&M University <khotilov@cern.ch>

Id:: MuonHitsChamberResidual.cc,v 1.1 2011/10/12 23:32:08 khotilov Exp

Author:: V. Khotilovich - Texas A&M University <khotilov@cern.ch>

Definition at line 22 of file MuonHitsChamberResidual.h.

Constructor & Destructor Documentation

MuonHitsChamberResidual::MuonHitsChamberResidual	(	edm::ESHandle< GlobalTrackingGeometry >	globalGeometry,
		AlignableNavigator *	navigator,
		DetId	chamberId,
		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(), MuonDT13ChamberResidual::addResidual(), and MuonDT2ChamberResidual::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;
}