#include <MuonDT2ChamberResidual.h>

Inheritance diagram for MuonDT2ChamberResidual:

Public Member Functions
virtual void	addResidual (const TrajectoryStateOnSurface tsos, const TransientTrackingRecHit hit)

	MuonDT2ChamberResidual (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 MuonHitsChamberResidual
	MuonHitsChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const AlignableDetOrUnitPtr &chamberAlignable)

void	segment_fit ()

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 MuonHitsChamberResidual
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

Detailed Description

Implementation of muon chamber residuals for transverse DT layers

Id:

Definition at line 13 of file MuonDT2ChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file MuonDT2ChamberResidual.cc.

References MuonChamberResidual::kDT2, MuonChamberResidual::m_chamberId, MuonChamberResidual::m_globalGeometry, MuonChamberResidual::m_sign, and MuonChamberResidual::m_type.

   : MuonHitsChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
 {
   m_type = MuonChamberResidual::kDT2; 
   align::GlobalVector zDirection(0., 0., 1.);
   m_sign = m_globalGeometry->idToDet(m_chamberId)->toLocal(zDirection).y() > 0. ? 1. : -1.; 
 }

Member Function Documentation

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

virtual

Implements MuonChamberResidual.

Definition at line 16 of file MuonDT2ChamberResidual.cc.

 {
   DetId id = hit->geographicalId();
 
   align::LocalPoint hitChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(hit->localPosition()));
   align::LocalPoint tsosChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(tsos->localPosition()));
 
   double residual = tsosChamberPos.y() - hitChamberPos.y();  // residual is track minus hit
   double weight = 1. / hit->localPositionError().xx();  // weight linear fit by hit-only local error (yes, xx: layer x is chamber y)
   double layerPosition = tsosChamberPos.z();  // the layer's position in the chamber's coordinate system
   double layerHitPos = hitChamberPos.z();
 
   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_hitx_1 += weight;
   m_hitx_x += weight * layerHitPos;
   m_hitx_y += weight * hitChamberPos.x();
   m_hitx_xx += weight * layerHitPos * layerHitPos;
   m_hitx_xy += weight * layerHitPos * hitChamberPos.x();
 
   m_hity_1 += weight;
   m_hity_x += weight * layerHitPos;
   m_hity_y += weight * hitChamberPos.y();
   m_hity_xx += weight * layerHitPos * layerHitPos;
   m_hity_xy += weight * layerHitPos * hitChamberPos.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);
   
   if (m_numHits>1) segment_fit();
 }

virtual void MuonDT2ChamberResidual::setSegmentResidual	(	const reco::MuonChamberMatch *	,
		const reco::MuonSegmentMatch *
	)

inlinevirtual

Implements MuonChamberResidual.

Definition at line 24 of file MuonDT2ChamberResidual.h.

24 {}

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation