MuonTrackCSCChamberResidual Class Reference

#include <MuonTrackCSCChamberResidual.h>

Inheritance diagram for MuonTrackCSCChamberResidual:

Public Member Functions
void	addResidual (edm::ESHandle< Propagator > prop, const TrajectoryStateOnSurface *tsos, const TrackingRecHit *hit, double, double) override=0
	MuonTrackCSCChamberResidual (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 CSC

Id:

Definition at line 13 of file MuonTrackCSCChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 8 of file MuonTrackCSCChamberResidual.cc.

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

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

Member Function Documentation

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

overridepure virtual

Implements MuonChamberResidual.

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

overridevirtual

Implements MuonChamberResidual.

Definition at line 18 of file MuonTrackCSCChamberResidual.cc.

References CSCSegment::chi2(), reco::MuonSegmentMatch::cscSegmentRef, CSCSegment::degreesOfFreedom(), reco::MuonChamberMatch::dXdZ, reco::MuonSegmentMatch::dXdZ, reco::MuonChamberMatch::dXdZErr, reco::MuonSegmentMatch::dXdZErr, reco::MuonChamberMatch::dYdZ, reco::MuonSegmentMatch::dYdZ, gpuClustering::id, reco::MuonChamberMatch::id, 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, CSCSegment::nRecHits(), funct::pow(), DetId::rawId(), mathSSE::sqrt(), reco::MuonChamberMatch::x, reco::MuonSegmentMatch::x, reco::MuonChamberMatch::xErr, reco::MuonSegmentMatch::xErr, reco::MuonChamberMatch::y, and reco::MuonSegmentMatch::y.

                                                                                      {
  CSCDetId id(trk->id.rawId());

  CSCSegmentRef segmentCSC = seg->cscSegmentRef;
  if (segmentCSC.get() != nullptr) {
    const CSCSegment *segment = segmentCSC.get();
    m_numHits = segment->nRecHits();
    m_ndof = segment->degreesOfFreedom();
    m_chi2 = segment->chi2();
    //std::cout<<"csc seg position = "<<segment->localPosition()<<"  numhits="<<m_numHits<<"  id: "<<id<<std::endl;
  }

  align::LocalPoint l_seg(seg->x, seg->y, 0.);
  align::LocalPoint l_trk(trk->x, trk->y, 0.);
  //align::GlobalPoint g_seg = m_globalGeometry->idToDet(chamber)->toGlobal(l_seg);
  //align::GlobalPoint g_trk = m_globalGeometry->idToDet(chamber)->toGlobal(l_trk);

  /*
  double dphi = g_trk.phi() - g_seg.phi();
  while (dphi >  M_PI) dphi -= 2.*M_PI;
  while (dphi < -M_PI) dphi += 2.*M_PI;
  m_residual = - m_sign * g_trk.perp() * dphi; // coming from global, need to adjust the sign
  std::cout<<"cscres="<<m_residual<<"  dx="<<trk->x-seg->x<<"  diff="<<trk->x-seg->x - m_residual<<std::endl;
  */
  m_residual = trk->x - seg->x;
  m_residual_error = sqrt(pow(trk->xErr, 2) + pow(seg->xErr, 2));
  m_resslope = trk->dXdZ - seg->dXdZ;
  m_resslope_error = sqrt(pow(trk->dXdZErr, 2) + pow(seg->dXdZErr, 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;

  //double yresidual_error = sqrt( pow(trk->yErr, 2) + pow(seg->yErr, 2) );
  //double yresslope_error = sqrt( pow(trk->dYdZErr, 2) + pow(seg->dYdZErr, 2) );
  //std::cout<<"csc res "<<m_residual<<"+-"<<m_residual_error<<"  "<<m_resslope<<"+-"<<m_resslope_error<<"  "<<trk->y-seg->y<<"+-"<<yresidual_error<<"  "<<trk->dYdZ - seg->dYdZ<<"+-"<<yresslope_error<<std::endl;
  //std::cout<<"csc trk "<<m_trackx<<" "<<m_tracky<<" "<<m_trackdxdz<<" "<<m_trackdydz<<std::endl;
  //std::cout<<"csc seg "<<m_segx<<" "<<m_segy<<" "<<m_segdxdz<<" "<<m_segdydz<<std::endl;
}