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

#include <MuonDT13ChamberResidual.h>

Inheritance diagram for MuonDT13ChamberResidual:
MuonHitsChamberResidual MuonChamberResidual

Public Member Functions

virtual void addResidual (edm::ESHandle< Propagator > prop, const TrajectoryStateOnSurface *tsos, const TrackingRecHit *hit, double, double)
 
 MuonDT13ChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, 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 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 axial DT layers

Id:

Definition at line 15 of file MuonDT13ChamberResidual.h.

Constructor & Destructor Documentation

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

Definition at line 8 of file MuonDT13ChamberResidual.cc.

References funct::cos(), MuonChamberResidual::kDT13, MuonChamberResidual::m_chamberId, MuonChamberResidual::m_globalGeometry, M_PI, MuonChamberResidual::m_sign, MuonChamberResidual::m_type, and funct::sin().

11 : MuonHitsChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable)
12 {
13  m_type = MuonChamberResidual::kDT13;
14  double rphiAngle = atan2(m_globalGeometry->idToDet(m_chamberId)->position().y(), m_globalGeometry->idToDet(m_chamberId)->position().x()) + M_PI/2.;
15  align::GlobalVector rphiDirection(cos(rphiAngle), sin(rphiAngle), 0.);
16  m_sign = m_globalGeometry->idToDet(m_chamberId)->toLocal(rphiDirection).x() > 0. ? 1. : -1.;
17 }
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
MuonChamberResidual::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:92
MuonHitsChamberResidual::MuonHitsChamberResidual
MuonHitsChamberResidual(edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, const AlignableDetOrUnitPtr &chamberAlignable)
Definition: MuonHitsChamberResidual.cc:8

Member Function Documentation

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

Implements MuonChamberResidual.

Definition at line 21 of file MuonDT13ChamberResidual.cc.

References gather_cfg::cout, TrackingRecHit::geographicalId(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPosition(), TrackingRecHit::localPositionError(), MuonHitsChamberResidual::m_chamber_length, MuonHitsChamberResidual::m_chamber_width, 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, 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, MuonChamberResidual::residual(), MuonHitsChamberResidual::segment_fit(), Alignable::surface(), AlignableSurface::toLocal(), puppiForMET_cff::weight, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

22 {
23  bool m_debug = false;
24 
25  m_chamber_width = chamber_width;
26  m_chamber_length = chamber_length;
27 
28 
29  DetId id = hit->geographicalId();
30 
31  align::LocalPoint hitChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(hit->localPosition()));
32  align::LocalPoint tsosChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(tsos->localPosition()));
33 
34  if (m_debug) {
35  std::cout << " MuonDT13ChamberResidual hitChamberPos x: " << hitChamberPos.x() << " tsosChamberPos x: " << tsosChamberPos.x() << std::endl;
36  std::cout << " hitChamberPos y: " << hitChamberPos.y() << " tsosChamberPos y: " << tsosChamberPos.y() << std::endl;
37  std::cout << " hitChamberPos z: " << hitChamberPos.z() << " tsosChamberPos z: " << tsosChamberPos.z() << std::endl;
38  }
39 
40  double residual = tsosChamberPos.x() - hitChamberPos.x(); // residual is hit minus hit
41  double weight = 1. / hit->localPositionError().xx(); // weight linear fit by hit-only local error
42  double layerPosition = tsosChamberPos.z(); // the layer's position in the chamber's coordinate system
43  double layerHitPos = hitChamberPos.z();
44 
45  m_numHits++;
46 
47  // "x" is the layerPosition, "y" is the residual (this is a linear fit to residual versus layerPosition)
48  m_residual_1 += weight;
49  m_residual_x += weight * layerPosition;
50  m_residual_y += weight * residual;
51  m_residual_xx += weight * layerPosition * layerPosition;
52  m_residual_xy += weight * layerPosition * residual;
53 
54  // "x" is the layerPosition, "y" is chamberx (this is a linear fit to chamberx versus layerPosition)
55  m_trackx_1 += weight;
56  m_trackx_x += weight * layerPosition;
57  m_trackx_y += weight * tsosChamberPos.x();
58  m_trackx_xx += weight * layerPosition * layerPosition;
59  m_trackx_xy += weight * layerPosition * tsosChamberPos.x();
60 
61  // "x" is the layerPosition, "y" is chambery (this is a linear fit to chambery versus layerPosition)
62  m_tracky_1 += weight;
63  m_tracky_x += weight * layerPosition;
64  m_tracky_y += weight * tsosChamberPos.y();
65  m_tracky_xx += weight * layerPosition * layerPosition;
66  m_tracky_xy += weight * layerPosition * tsosChamberPos.y();
67 
68  m_hitx_1 += weight;
69  m_hitx_x += weight * layerHitPos;
70  m_hitx_y += weight * hitChamberPos.x();
71  m_hitx_xx += weight * layerHitPos * layerHitPos;
72  m_hitx_xy += weight * layerHitPos * hitChamberPos.x();
73 
74  m_hity_1 += weight;
75  m_hity_x += weight * layerHitPos;
76  m_hity_y += weight * hitChamberPos.y();
77  m_hity_xx += weight * layerHitPos * layerPosition;
78  m_hity_xy += weight * layerHitPos * hitChamberPos.y();
79 
80  m_localIDs.push_back(id);
81  // m_localResids.push_back(tsos->localPosition().x() - hit->localPosition().x()); //FIXME looks like this line is not used anywhere, moreover it is wrong for segment-based reconstruction, I changed it to the follwoing line
82  m_localResids.push_back(residual);
83  m_individual_x.push_back(layerPosition);
84  m_individual_y.push_back(residual);
85  m_individual_weight.push_back(weight);
86 
87  if (m_numHits>1) segment_fit();
88 }
LocalError::xx
float xx() const
Definition: LocalError.h:24
MuonChamberResidual::m_individual_y
std::vector< double > m_individual_y
Definition: MuonChamberResidual.h:101
puppiForMET_cff.weight
tuple weight
Definition: puppiForMET_cff.py:35
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
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.
MuonChamberResidual::m_localResids
std::vector< double > m_localResids
Definition: MuonChamberResidual.h:99
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
MuonChamberResidual::m_chamberAlignable
AlignableDetOrUnitPtr m_chamberAlignable
Definition: MuonChamberResidual.h:95
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:131
DetId
Definition: DetId.h:18
TrackingRecHit::localPositionError
virtual LocalError localPositionError() const =0
MuonChamberResidual::m_globalGeometry
edm::ESHandle< GlobalTrackingGeometry > m_globalGeometry
Definition: MuonChamberResidual.h:92
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:116
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:102
MuonChamberResidual::residual
double residual() const
Definition: MuonChamberResidual.h:57
virtual void MuonDT13ChamberResidual::setSegmentResidual ( const reco::MuonChamberMatch *  ,
const reco::MuonSegmentMatch *   
)
inlinevirtual

Implements MuonChamberResidual.

Definition at line 26 of file MuonDT13ChamberResidual.h.

26 {}
Generated for CMSSW Reference Manual by doxygen 1.8.5