CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions
MuonTrackCSCChamberResidual Class Reference

#include <MuonTrackCSCChamberResidual.h>

Inheritance diagram for MuonTrackCSCChamberResidual:
MuonChamberResidual

Public Member Functions

virtual void addResidual (const TrajectoryStateOnSurface *tsos, const TransientTrackingRecHit *hit)
 
 MuonTrackCSCChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
 
virtual void setSegmentResidual (const reco::MuonChamberMatch *, const reco::MuonSegmentMatch *)
 
- 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, 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_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:
MuonTrackCSCChamberResidual.h,v 1.1 2011/10/12 23:32:07 khotilov Exp

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 9 of file MuonTrackCSCChamberResidual.cc.

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

11  : MuonChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
12 {
13  m_type = MuonChamberResidual::kCSC;
14  align::GlobalVector zDirection(0., 0., 1.);
15  m_sign = m_globalGeometry->idToDet(m_chamberId)->toLocal(zDirection).z() > 0. ? 1. : -1.;
16 }
MuonChamberResidual::MuonChamberResidual
MuonChamberResidual(edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable)
Definition: MuonChamberResidual.cc:14
MuonChamberResidual::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:86

Member Function Documentation

virtual void MuonTrackCSCChamberResidual::addResidual ( const TrajectoryStateOnSurface *  tsos,
const TransientTrackingRecHit *  hit 
)
inlinevirtual

Implements MuonChamberResidual.

Definition at line 20 of file MuonTrackCSCChamberResidual.h.

20 {}
void MuonTrackCSCChamberResidual::setSegmentResidual ( const reco::MuonChamberMatch *  trk,
const reco::MuonSegmentMatch *  seg 
)
virtual

Implements MuonChamberResidual.

Definition at line 19 of file MuonTrackCSCChamberResidual.cc.

References reco::MuonSegmentMatch::cscSegmentRef, reco::MuonChamberMatch::dXdZ, reco::MuonSegmentMatch::dXdZ, reco::MuonChamberMatch::dXdZErr, reco::MuonSegmentMatch::dXdZErr, reco::MuonChamberMatch::dYdZ, reco::MuonSegmentMatch::dYdZ, 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, 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.

20 {
21  CSCDetId id(trk->id.rawId());
22 
23  CSCSegmentRef segmentCSC = seg->cscSegmentRef;
24  if (segmentCSC.get() != 0)
25  {
26  const CSCSegment* segment = segmentCSC.get();
27  m_numHits = segment->nRecHits();
28  m_ndof = segment->degreesOfFreedom();
29  m_chi2 = segment->chi2();
30  //std::cout<<"csc seg position = "<<segment->localPosition()<<" numhits="<<m_numHits<<" id: "<<id<<std::endl;
31  }
32 
33  align::LocalPoint l_seg(seg->x, seg->y, 0.);
34  align::LocalPoint l_trk(trk->x, trk->y, 0.);
35  //align::GlobalPoint g_seg = m_globalGeometry->idToDet(chamber)->toGlobal(l_seg);
36  //align::GlobalPoint g_trk = m_globalGeometry->idToDet(chamber)->toGlobal(l_trk);
37 
38  /*
39  double dphi = g_trk.phi() - g_seg.phi();
40  while (dphi > M_PI) dphi -= 2.*M_PI;
41  while (dphi < -M_PI) dphi += 2.*M_PI;
42  m_residual = - m_sign * g_trk.perp() * dphi; // coming from global, need to adjust the sign
43  std::cout<<"cscres="<<m_residual<<" dx="<<trk->x-seg->x<<" diff="<<trk->x-seg->x - m_residual<<std::endl;
44  */
45  m_residual = trk->x-seg->x;
46  m_residual_error = sqrt( pow(trk->xErr, 2) + pow(seg->xErr, 2) );
47  m_resslope = trk->dXdZ - seg->dXdZ;
48  m_resslope_error = sqrt( pow(trk->dXdZErr, 2) + pow(seg->dXdZErr, 2) );
49 
50  m_trackx = trk->x;
51  m_tracky = trk->y;
52  m_trackdxdz = trk->dXdZ;
53  m_trackdydz = trk->dYdZ;
54 
55  m_segx = seg->x;
56  m_segy = seg->y;
57  m_segdxdz = seg->dXdZ;
58  m_segdydz = seg->dYdZ;
59 
60  //double yresidual_error = sqrt( pow(trk->yErr, 2) + pow(seg->yErr, 2) );
61  //double yresslope_error = sqrt( pow(trk->dYdZErr, 2) + pow(seg->dYdZErr, 2) );
62  //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;
63  //std::cout<<"csc trk "<<m_trackx<<" "<<m_tracky<<" "<<m_trackdxdz<<" "<<m_trackdydz<<std::endl;
64  //std::cout<<"csc seg "<<m_segx<<" "<<m_segy<<" "<<m_segdxdz<<" "<<m_segdydz<<std::endl;
65 }
reco::MuonSegmentMatch::cscSegmentRef
CSCSegmentRef cscSegmentRef
Definition: MuonSegmentMatch.h:48
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
Generated for CMSSW Reference Manual by doxygen 1.8.5