#include <MuonDT13ChamberResidual.h>

Inheritance diagram for MuonDT13ChamberResidual:

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

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

void	setSegmentResidual (const reco::MuonChamberMatch , const reco::MuonSegmentMatch ) override

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 axial DT layers

Id

Definition at line 15 of file MuonDT13ChamberResidual.h.

Constructor & Destructor Documentation

◆ MuonDT13ChamberResidual()

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

Definition at line 7 of file MuonDT13ChamberResidual.cc.

     : MuonHitsChamberResidual(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.;
 }

References funct::cos(), GlobalTrackingGeometry::idToDet(), MuonChamberResidual::kDT13, MuonChamberResidual::m_chamberId, MuonChamberResidual::m_globalGeometry, M_PI, MuonChamberResidual::m_sign, MuonChamberResidual::m_type, GeomDet::position(), funct::sin(), GeomDet::toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Member Function Documentation

◆ addResidual()

void MuonDT13ChamberResidual::addResidual	(	edm::ESHandle< Propagator >	prop,
		const TrajectoryStateOnSurface *	tsos,
		const TrackingRecHit *	hit,
		double	chamber_width,
		double	chamber_length
	)

overridevirtual

Implements MuonChamberResidual.

Definition at line 22 of file MuonDT13ChamberResidual.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 << " MuonDT13ChamberResidual 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.x() - hitChamberPos.x();  // residual is hit minus hit
   double weight = 1. / hit->localPositionError().xx();       // weight linear fit by hit-only local error
   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 * layerPosition;
   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();
 }

◆ setSegmentResidual()

void MuonDT13ChamberResidual::setSegmentResidual	(	const reco::MuonChamberMatch *	,
		const reco::MuonSegmentMatch *
	)

inlineoverridevirtual

Implements MuonChamberResidual.

Definition at line 31 of file MuonDT13ChamberResidual.h.

31 {}