CMS 3D CMS Logo
List of all members | Public Member Functions
MuonCSCChamberResidual Class Reference

#include <MuonCSCChamberResidual.h>

Inheritance diagram for MuonCSCChamberResidual:
MuonHitsChamberResidual MuonChamberResidual

Public Member Functions

void addResidual (edm::ESHandle< Propagator > prop, const TrajectoryStateOnSurface *tsos, const TrackingRecHit *hit, double, double) override
 
 MuonCSCChamberResidual (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 CSC

Id

Definition at line 14 of file MuonCSCChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 11 of file MuonCSCChamberResidual.cc.

References GlobalTrackingGeometry::idToDet(), MuonChamberResidual::kCSC, MuonChamberResidual::m_chamberId, MuonChamberResidual::m_globalGeometry, MuonChamberResidual::m_sign, MuonChamberResidual::m_type, GeomDet::toLocal(), and PV3DBase< T, PVType, FrameType >::z().

13 : MuonHitsChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
14 {
15  m_type = MuonChamberResidual::kCSC;
16  align::GlobalVector zDirection(0., 0., 1.);
17  m_sign = m_globalGeometry->idToDet(m_chamberId)->toLocal(zDirection).z() > 0. ? 1. : -1.;
18 }
GeomDet::toLocal
LocalPoint toLocal(const GlobalPoint &gp) const
Conversion to the R.F. of the GeomDet.
Definition: GeomDet.h:69
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
MuonChamberResidual::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:92
GlobalTrackingGeometry::idToDet
const GeomDet * idToDet(DetId) const override
Definition: GlobalTrackingGeometry.cc:44
MuonHitsChamberResidual::MuonHitsChamberResidual
MuonHitsChamberResidual(edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const AlignableDetOrUnitPtr &chamberAlignable)
Definition: MuonHitsChamberResidual.cc:8

Member Function Documentation

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

Implements MuonChamberResidual.

Definition at line 21 of file MuonCSCChamberResidual.cc.

References angle(), funct::cos(), gather_cfg::cout, TrackingRecHit::geographicalId(), CSCLayer::geometry(), GlobalTrackingGeometry::idToDet(), 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(), GlobalTrackingGeometry::slaveGeometry(), digitizers_cfi::strip, CSCLayerGeometry::stripAngle(), Alignable::surface(), GeomDet::toGlobal(), AlignableSurface::toLocal(), mps_merge::weight, PV3DBase< T, PVType, FrameType >::x(), geometryCSVtoXML::xx, LocalError::xx(), geometryCSVtoXML::xy, LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), geometryCSVtoXML::yy, LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

22 {
23  bool m_debug = false;
24 
25  if (m_debug) std::cout << "MuonCSCChamberResidual::addResidual 1" << std::endl;
26  DetId id = hit->geographicalId();
27  if (m_debug) std::cout << "MuonCSCChamberResidual::addResidual 2" << std::endl;
28  const CSCGeometry *cscGeometry = dynamic_cast<const CSCGeometry*>(m_globalGeometry->slaveGeometry(id));
29  if (m_debug) std::cout << "MuonCSCChamberResidual::addResidual 3" << std::endl;
30  assert(cscGeometry);
31 
32  if (m_debug) {
33  std::cout << " MuonCSCChamberResidual hit->localPosition() x: " << hit->localPosition().x() << " tsos->localPosition() x: " << tsos->localPosition().x() << std::endl;
34  std::cout << " hit->localPosition() y: " << hit->localPosition().y() << " tsos->localPosition() y: " << tsos->localPosition().y() << std::endl;
35  std::cout << " hit->localPosition() z: " << hit->localPosition().z() << " tsos->localPosition() z: " << tsos->localPosition().z() << std::endl;
36  }
37 
38  // hit->localPosition() is coordinate in local system of LAYER. Transfer it to coordiante in local system of chamber
39  align::LocalPoint hitChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(hit->localPosition()));
40  // TSOS->localPosition() is given in local system of CHAMBER (for segment-based reconstruction)
41  // align::LocalPoint tsosChamberPos = tsos->localPosition();
42  align::LocalPoint tsosChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(tsos->localPosition()));
43 
44  int strip = cscGeometry->layer(id)->geometry()->nearestStrip(hit->localPosition());
45  double angle = cscGeometry->layer(id)->geometry()->stripAngle(strip) - M_PI/2.;
46  double sinAngle = sin(angle);
47  double cosAngle = cos(angle);
48 
49  double residual = cosAngle * (tsosChamberPos.x() - hitChamberPos.x()) + sinAngle * (tsosChamberPos.y() - hitChamberPos.y()); // yes, that's +sin()
50 
51  if (m_debug) std::cout << " MuonCSCChamberResidual residual: " << residual << std::endl;
52 
53  double xx = hit->localPositionError().xx();
54  double xy = hit->localPositionError().xy();
55  double yy = hit->localPositionError().yy();
56  double weight = 1. / (xx*cosAngle*cosAngle + 2.*xy*sinAngle*cosAngle + yy*sinAngle*sinAngle);
57 
58  double layerPosition = tsosChamberPos.z(); // the layer's position in the chamber's coordinate system
59  double layerHitPos = hitChamberPos.z();
60 
61  m_numHits++;
62 
63  // "x" is the layerPosition, "y" is the residual (this is a linear fit to residual versus layerPosition)
64  m_residual_1 += weight;
65  m_residual_x += weight * layerPosition;
66  m_residual_y += weight * residual;
67  m_residual_xx += weight * layerPosition * layerPosition;
68  m_residual_xy += weight * layerPosition * residual;
69 
70  // "x" is the layerPosition, "y" is chamberx (this is a linear fit to chamberx versus layerPosition)
71  m_trackx_1 += weight;
72  m_trackx_x += weight * layerPosition;
73  m_trackx_y += weight * tsosChamberPos.x();
74  m_trackx_xx += weight * layerPosition * layerPosition;
75  m_trackx_xy += weight * layerPosition * tsosChamberPos.x();
76 
77  // "x" is the layerPosition, "y" is chambery (this is a linear fit to chambery versus layerPosition)
78  m_tracky_1 += weight;
79  m_tracky_x += weight * layerPosition;
80  m_tracky_y += weight * tsosChamberPos.y();
81  m_tracky_xx += weight * layerPosition * layerPosition;
82  m_tracky_xy += weight * layerPosition * tsosChamberPos.y();
83 
84  m_hitx_1 += weight;
85  m_hitx_x += weight * layerHitPos;
86  m_hitx_y += weight * hitChamberPos.x();
87  m_hitx_xx += weight * layerHitPos * layerHitPos;
88  m_hitx_xy += weight * layerHitPos * hitChamberPos.x();
89 
90  m_hity_1 += weight;
91  m_hity_x += weight * layerHitPos;
92  m_hity_y += weight * hitChamberPos.y();
93  m_hity_xx += weight * layerHitPos * layerHitPos;
94  m_hity_xy += weight * layerHitPos * hitChamberPos.y();
95 
96  m_localIDs.push_back(id);
97  m_localResids.push_back(residual); // FIXME Check if this method is needed
98  m_individual_x.push_back(layerPosition);
99  m_individual_y.push_back(residual);
100  m_individual_weight.push_back(weight);
101 
102  if (m_numHits>1) segment_fit();
103 }
LocalError::xx
float xx() const
Definition: LocalError.h:24
MuonChamberResidual::m_individual_y
std::vector< double > m_individual_y
Definition: MuonChamberResidual.h:101
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:54
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
weight
Definition: weight.py:1
LocalError::xy
float xy() const
Definition: LocalError.h:25
MuonChamberResidual::m_individual_x
std::vector< double > m_individual_x
Definition: MuonChamberResidual.h:100
MuonChamberResidual::m_localIDs
std::vector< DetId > m_localIDs
Definition: MuonChamberResidual.h:98
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:99
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
GlobalTrackingGeometry::slaveGeometry
const TrackingGeometry * slaveGeometry(DetId id) const
Return the pointer to the actual geometry for a given DetId.
Definition: GlobalTrackingGeometry.cc:55
MuonChamberResidual::m_chamberAlignable
AlignableDetOrUnitPtr m_chamberAlignable
Definition: MuonChamberResidual.h:95
TrackingRecHit::localPosition
virtual LocalPoint localPosition() const =0
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:135
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
CSCLayerGeometry::nearestStrip
int nearestStrip(const LocalPoint &lp) const
Definition: CSCLayerGeometry.h:94
DetId
Definition: DetId.h:18
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
CSCGeometry::layer
const CSCLayer * layer(CSCDetId id) const
Return the layer corresponding to given DetId.
Definition: CSCGeometry.cc:139
MuonChamberResidual::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:92
GlobalTrackingGeometry::idToDet
const GeomDet * idToDet(DetId) const override
Definition: GlobalTrackingGeometry.cc:44
gather_cfg.cout
cout
Definition: gather_cfg.py:144
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:116
TrackingRecHit::localPositionError
virtual LocalError localPositionError() const =0
mps_merge.weight
weight
Definition: mps_merge.py:87
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
MuonChamberResidual::m_individual_weight
std::vector< double > m_individual_weight
Definition: MuonChamberResidual.h:102
CSCLayerGeometry::stripAngle
float stripAngle(int strip) const
Definition: CSCLayerGeometry.cc:140
MuonChamberResidual::residual
double residual() const
Definition: MuonChamberResidual.h:57
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
void MuonCSCChamberResidual::setSegmentResidual ( const reco::MuonChamberMatch *  ,
const reco::MuonSegmentMatch *   
)
inlineoverridevirtual

Implements MuonChamberResidual.

Definition at line 25 of file MuonCSCChamberResidual.h.

25 {}
Generated for CMSSW Reference Manual by doxygen 1.8.11