#include <MuonDT2ChamberResidual.h>

Inheritance diagram for MuonDT2ChamberResidual:

Public Member Functions
virtual void	addResidual (edm::ESHandle< Propagator > prop, const TrajectoryStateOnSurface tsos, const TrackingRecHit hit, double, double)

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

double	ChambW () 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, 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_chamber_length

double	m_chamber_width

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_Chambl

double	m_ChambW

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,
		AlignableDetOrUnitPtr	chamberAlignable
	)

Definition at line 7 of file MuonDT2ChamberResidual.cc.

References GlobalTrackingGeometry::idToDet(), MuonChamberResidual::kDT2, MuonChamberResidual::m_chamberId, MuonChamberResidual::m_globalGeometry, MuonChamberResidual::m_sign, MuonChamberResidual::m_type, GeomDet::toLocal(), and PV3DBase< T, PVType, FrameType >::y().

 : 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	(	edm::ESHandle< Propagator >	prop,
		const TrajectoryStateOnSurface *	tsos,
		const TrackingRecHit *	hit,
		double	chamber_width,
		double	chamber_length
	)

virtual

Implements MuonChamberResidual.

Definition at line 19 of file MuonDT2ChamberResidual.cc.

 {
     bool m_debug = false;
 
     m_chamber_width = chamber_width;
     m_chamber_length = chamber_length;
 
     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()));
 
     if (m_debug) {
         std::cout << " MuonDT2ChamberResidual hitChamberPos x: " <<  hitChamberPos.x() << " tsosChamberPos x: " << tsosChamberPos.x() << std::endl;
         std::cout << "                        hitChamberPos y: " <<  hitChamberPos.y() << " tsosChamberPos y: " << tsosChamberPos.y() << std::endl;
         std::cout << "                        hitChamberPos z: " <<  hitChamberPos.z() << " tsosChamberPos z: " << tsosChamberPos.z() << std::endl;
     }
 
     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()); //FIXME looks like this line is not used anywhere, moreover it is wrong for segment-based reconstruction, I changed it to the follwoing line
     m_localResids.push_back(residual);
     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