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

#include <MuonCSCChamberResidual.h>

Inheritance diagram for MuonCSCChamberResidual:
MuonHitsChamberResidual MuonChamberResidual

Public Member Functions

virtual void addResidual (const TrajectoryStateOnSurface *tsos, const TransientTrackingRecHit *hit)
 
 MuonCSCChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const 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 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, const 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_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_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 CSC

Id:

Definition at line 13 of file MuonCSCChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 9 of file MuonCSCChamberResidual.cc.

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

11  : MuonHitsChamberResidual(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::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:86
MuonHitsChamberResidual::MuonHitsChamberResidual
MuonHitsChamberResidual(edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const AlignableDetOrUnitPtr &chamberAlignable)
Definition: MuonHitsChamberResidual.cc:8

Member Function Documentation

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

Implements MuonChamberResidual.

Definition at line 19 of file MuonCSCChamberResidual.cc.

References angle(), assert(), funct::cos(), TrackingRecHit::geographicalId(), CSCLayer::geometry(), CSCGeometry::layer(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPosition(), TrackingRecHit::localPositionError(), MuonChamberResidual::m_chamberAlignable, MuonChamberResidual::m_globalGeometry, MuonHitsChamberResidual::m_hitx_1, MuonHitsChamberResidual::m_hitx_x, MuonHitsChamberResidual::m_hitx_xx, MuonHitsChamberResidual::m_hitx_xy, MuonHitsChamberResidual::m_hitx_y, MuonHitsChamberResidual::m_hity_1, MuonHitsChamberResidual::m_hity_x, MuonHitsChamberResidual::m_hity_xx, MuonHitsChamberResidual::m_hity_xy, MuonHitsChamberResidual::m_hity_y, MuonChamberResidual::m_individual_weight, MuonChamberResidual::m_individual_x, MuonChamberResidual::m_individual_y, MuonChamberResidual::m_localIDs, MuonChamberResidual::m_localResids, MuonChamberResidual::m_numHits, M_PI, MuonHitsChamberResidual::m_residual_1, MuonHitsChamberResidual::m_residual_x, MuonHitsChamberResidual::m_residual_xx, MuonHitsChamberResidual::m_residual_xy, MuonHitsChamberResidual::m_residual_y, MuonHitsChamberResidual::m_trackx_1, MuonHitsChamberResidual::m_trackx_x, MuonHitsChamberResidual::m_trackx_xx, MuonHitsChamberResidual::m_trackx_xy, MuonHitsChamberResidual::m_trackx_y, MuonHitsChamberResidual::m_tracky_1, MuonHitsChamberResidual::m_tracky_x, MuonHitsChamberResidual::m_tracky_xx, MuonHitsChamberResidual::m_tracky_xy, MuonHitsChamberResidual::m_tracky_y, CSCLayerGeometry::nearestStrip(), MuonChamberResidual::residual(), MuonHitsChamberResidual::segment_fit(), funct::sin(), CSCLayerGeometry::stripAngle(), Alignable::surface(), AlignableSurface::toLocal(), histoStyle::weight, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), create_public_lumi_plots::xy, PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

20 {
21  DetId id = hit->geographicalId();
22  const CSCGeometry *cscGeometry = dynamic_cast<const CSCGeometry*>(m_globalGeometry->slaveGeometry(id));
23  assert(cscGeometry);
24 
25  align::LocalPoint hitChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(hit->localPosition()));
26  align::LocalPoint tsosChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(tsos->localPosition()));
27 
28  int strip = cscGeometry->layer(id)->geometry()->nearestStrip(hit->localPosition());
29  double angle = cscGeometry->layer(id)->geometry()->stripAngle(strip) - M_PI/2.;
30  double sinAngle = sin(angle);
31  double cosAngle = cos(angle);
32 
33  double residual = cosAngle * (tsosChamberPos.x() - hitChamberPos.x()) + sinAngle * (tsosChamberPos.y() - hitChamberPos.y()); // yes, that's +sin()
34 
35  double xx = hit->localPositionError().xx();
36  double xy = hit->localPositionError().xy();
37  double yy = hit->localPositionError().yy();
38  double weight = 1. / (xx*cosAngle*cosAngle + 2.*xy*sinAngle*cosAngle + yy*sinAngle*sinAngle);
39 
40  double layerPosition = tsosChamberPos.z(); // the layer's position in the chamber's coordinate system
41  double layerHitPos = hitChamberPos.z();
42 
43  m_numHits++;
44 
45  // "x" is the layerPosition, "y" is the residual (this is a linear fit to residual versus layerPosition)
46  m_residual_1 += weight;
47  m_residual_x += weight * layerPosition;
48  m_residual_y += weight * residual;
49  m_residual_xx += weight * layerPosition * layerPosition;
50  m_residual_xy += weight * layerPosition * residual;
51 
52  // "x" is the layerPosition, "y" is chamberx (this is a linear fit to chamberx versus layerPosition)
53  m_trackx_1 += weight;
54  m_trackx_x += weight * layerPosition;
55  m_trackx_y += weight * tsosChamberPos.x();
56  m_trackx_xx += weight * layerPosition * layerPosition;
57  m_trackx_xy += weight * layerPosition * tsosChamberPos.x();
58 
59  // "x" is the layerPosition, "y" is chambery (this is a linear fit to chambery versus layerPosition)
60  m_tracky_1 += weight;
61  m_tracky_x += weight * layerPosition;
62  m_tracky_y += weight * tsosChamberPos.y();
63  m_tracky_xx += weight * layerPosition * layerPosition;
64  m_tracky_xy += weight * layerPosition * tsosChamberPos.y();
65 
66  m_hitx_1 += weight;
67  m_hitx_x += weight * layerHitPos;
68  m_hitx_y += weight * hitChamberPos.x();
69  m_hitx_xx += weight * layerHitPos * layerHitPos;
70  m_hitx_xy += weight * layerHitPos * hitChamberPos.x();
71 
72  m_hity_1 += weight;
73  m_hity_x += weight * layerHitPos;
74  m_hity_y += weight * hitChamberPos.y();
75  m_hity_xx += weight * layerHitPos * layerHitPos;
76  m_hity_xy += weight * layerHitPos * hitChamberPos.y();
77 
78  m_localIDs.push_back(id);
79  m_localResids.push_back(residual);
80  m_individual_x.push_back(layerPosition);
81  m_individual_y.push_back(residual);
82  m_individual_weight.push_back(weight);
83 
84  if (m_numHits>1) segment_fit();
85 }
LocalError::xx
float xx() const
Definition: LocalError.h:24
MuonChamberResidual::m_individual_y
std::vector< double > m_individual_y
Definition: MuonChamberResidual.h:95
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
assert
assert(m_qm.get())
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
create_public_lumi_plots.xy
tuple xy
Definition: create_public_lumi_plots.py:1810
LocalError::xy
float xy() const
Definition: LocalError.h:25
MuonChamberResidual::m_individual_x
std::vector< double > m_individual_x
Definition: MuonChamberResidual.h:94
MuonChamberResidual::m_localIDs
std::vector< DetId > m_localIDs
Definition: MuonChamberResidual.h:92
AlignableSurface::toLocal
align::RotationType toLocal(const align::RotationType &) const
Return in local frame a rotation given in global frame.
LocalError::yy
float yy() const
Definition: LocalError.h:26
MuonChamberResidual::m_localResids
std::vector< double > m_localResids
Definition: MuonChamberResidual.h:93
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
MuonChamberResidual::m_chamberAlignable
AlignableDetOrUnitPtr m_chamberAlignable
Definition: MuonChamberResidual.h:89
M_PI
#define M_PI
Definition: readTestVector.cc:25
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
CSCLayerGeometry::nearestStrip
int nearestStrip(const LocalPoint &lp) const
Definition: CSCLayerGeometry.h:94
DetId
Definition: DetId.h:18
TrackingRecHit::localPositionError
virtual LocalError localPositionError() const =0
CSCGeometry::layer
const CSCLayer * layer(CSCDetId id) const
Return the layer corresponding to given DetId.
Definition: CSCGeometry.cc:124
MuonChamberResidual::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:86
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:116
histoStyle.weight
int weight
Definition: histoStyle.py:50
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
TrackingRecHit::localPosition
virtual LocalPoint localPosition() const =0
MuonChamberResidual::m_individual_weight
std::vector< double > m_individual_weight
Definition: MuonChamberResidual.h:96
CSCLayerGeometry::stripAngle
float stripAngle(int strip) const
Definition: CSCLayerGeometry.cc:140
MuonChamberResidual::residual
double residual() const
Definition: MuonChamberResidual.h:55
CSCLayer::geometry
const CSCLayerGeometry * geometry() const
Definition: CSCLayer.h:47
angle
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11
virtual void MuonCSCChamberResidual::setSegmentResidual ( const reco::MuonChamberMatch *  ,
const reco::MuonSegmentMatch *   
)
inlinevirtual

Implements MuonChamberResidual.

Definition at line 24 of file MuonCSCChamberResidual.h.

24 {}
Generated for CMSSW Reference Manual by doxygen 1.8.5