MuonTrackDT2ChamberResidual Class Reference

#include <MuonTrackDT2ChamberResidual.h>

Inheritance diagram for MuonTrackDT2ChamberResidual:

Public Member Functions
void	addResidual (edm::ESHandle< Propagator > prop, const TrajectoryStateOnSurface *tsos, const TrackingRecHit *hit, double, double) override=0
	MuonTrackDT2ChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
void	setSegmentResidual (const reco::MuonChamberMatch *, const reco::MuonSegmentMatch *) override
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 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 tracker muon chamber residuals for transverse DT layers

Id

Definition at line 13 of file MuonTrackDT2ChamberResidual.h.

Constructor & Destructor Documentation

MuonTrackDT2ChamberResidual()

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

Definition at line 7 of file MuonTrackDT2ChamberResidual.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().

    : MuonChamberResidual(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

addResidual()

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

overridepure virtual

Implements MuonChamberResidual.

setSegmentResidual()

void MuonTrackDT2ChamberResidual::setSegmentResidual ( const reco::MuonChamberMatch *  trk, const reco::MuonSegmentMatch *  seg  )

overridevirtual

Implements MuonChamberResidual.

Definition at line 17 of file MuonTrackDT2ChamberResidual.cc.

References cms::cuda::assert(), reco::MuonSegmentMatch::dtSegmentRef, reco::MuonChamberMatch::dXdZ, reco::MuonSegmentMatch::dXdZ, reco::MuonChamberMatch::dYdZ, reco::MuonSegmentMatch::dYdZ, reco::MuonChamberMatch::dYdZErr, reco::MuonSegmentMatch::dYdZErr, edm::Ref< C, T, F >::get(), DTRecSegment4D::hasZed(), MuonChamberResidual::m_chi2, MuonChamberResidual::m_ndof, MuonChamberResidual::m_numHits, MuonChamberResidual::m_residual, MuonChamberResidual::m_residual_error, MuonChamberResidual::m_resslope, MuonChamberResidual::m_resslope_error, MuonChamberResidual::m_segdxdz, MuonChamberResidual::m_segdydz, MuonChamberResidual::m_segx, MuonChamberResidual::m_segy, MuonChamberResidual::m_trackdxdz, MuonChamberResidual::m_trackdydz, MuonChamberResidual::m_trackx, MuonChamberResidual::m_tracky, funct::pow(), mathSSE::sqrt(), reco::MuonChamberMatch::x, reco::MuonSegmentMatch::x, reco::MuonChamberMatch::y, reco::MuonSegmentMatch::y, reco::MuonChamberMatch::yErr, and reco::MuonSegmentMatch::yErr.

                                                                                      {
  DTRecSegment4DRef segmentDT = seg->dtSegmentRef;
  if (segmentDT.get() != nullptr) {
    const DTRecSegment4D *segment = segmentDT.get();
    assert(segment->hasZed());
    const DTSLRecSegment2D *zSeg = (*segment).zSegment();
    m_numHits = zSeg->recHits().size();
    m_ndof = zSeg->degreesOfFreedom();
    m_chi2 = zSeg->chi2();
    //std::cout<<"z seg position = "<<zSeg->localPosition()<<"  numhits="<<m_numHits<<std::endl;
  }

  m_residual = trk->y - seg->y;
  m_residual_error = sqrt(pow(trk->yErr, 2) + pow(seg->yErr, 2));
  m_resslope = trk->dYdZ - seg->dYdZ;
  m_resslope_error = sqrt(pow(trk->dYdZErr, 2) + pow(seg->dYdZErr, 2));

  m_trackx = trk->x;
  m_tracky = trk->y;
  m_trackdxdz = trk->dXdZ;
  m_trackdydz = trk->dYdZ;

  m_segx = seg->x;
  m_segy = seg->y;
  m_segdxdz = seg->dXdZ;
  m_segdydz = seg->dYdZ;

  //std::cout<<"d2 res "<<m_residual<<"+-"<<m_residual_error<<"  "<<m_resslope<<"+-"<<m_resslope_error<<std::endl;
  //std::cout<<"d2 trk "<<m_trackx<<" "<<m_tracky<<" "<<m_trackdxdz<<" "<<m_trackdydz<<std::endl;
  //std::cout<<"d2 seg "<<m_segx<<" "<<m_segy<<" "<<m_segdxdz<<" "<<m_segdydz<<std::endl;
}