|
|
|
#include <MuonCSCChamberResidual.h>
Public Member Functions | |
| void | addResidual (const TrajectoryStateOnSurface *tsos, const TransientTrackingRecHit *hit) |
| MuonCSCChamberResidual (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, AlignableNavigator *navigator, DetId chamberId, AlignableDetOrUnitPtr chamberAlignable) | |
| double | signConvention (const unsigned int rawId=0) const |
| int | type () const |
| MuonCSCChamberResidual::MuonCSCChamberResidual | ( | edm::ESHandle< GlobalTrackingGeometry > | globalGeometry, |
| AlignableNavigator * | navigator, | ||
| DetId | chamberId, | ||
| AlignableDetOrUnitPtr | chamberAlignable | ||
| ) | [inline] |
Definition at line 14 of file MuonCSCChamberResidual.h.
: MuonChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable) {};
| void MuonCSCChamberResidual::addResidual | ( | const TrajectoryStateOnSurface * | tsos, |
| const TransientTrackingRecHit * | hit | ||
| ) | [virtual] |
Implements MuonChamberResidual.
Definition at line 4 of file MuonCSCChamberResidual.cc.
References angle(), funct::cos(), TrackingRecHit::geographicalId(), CSCLayer::geometry(), CSCGeometry::layer(), TrackingRecHit::localPosition(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPositionError(), MuonChamberResidual::m_chamberAlignable, MuonChamberResidual::m_globalGeometry, MuonChamberResidual::m_individual_weight, MuonChamberResidual::m_individual_x, MuonChamberResidual::m_individual_y, MuonChamberResidual::m_localIDs, MuonChamberResidual::m_localResids, MuonChamberResidual::m_numHits, M_PI, MuonChamberResidual::m_residual_1, MuonChamberResidual::m_residual_x, MuonChamberResidual::m_residual_xx, MuonChamberResidual::m_residual_xy, MuonChamberResidual::m_residual_y, MuonChamberResidual::m_trackx_1, MuonChamberResidual::m_trackx_x, MuonChamberResidual::m_trackx_xx, MuonChamberResidual::m_trackx_xy, MuonChamberResidual::m_trackx_y, MuonChamberResidual::m_tracky_1, MuonChamberResidual::m_tracky_x, MuonChamberResidual::m_tracky_xx, MuonChamberResidual::m_tracky_xy, MuonChamberResidual::m_tracky_y, CSCLayerGeometry::nearestStrip(), MuonChamberResidual::residual(), funct::sin(), strip(), CSCLayerGeometry::stripAngle(), CommonMethods::weight(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().
{
DetId id = hit->geographicalId();
const CSCGeometry *cscGeometry = dynamic_cast<const CSCGeometry*>(m_globalGeometry->slaveGeometry(id));
assert(cscGeometry);
LocalPoint hitChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(hit->localPosition()));
LocalPoint tsosChamberPos = m_chamberAlignable->surface().toLocal(m_globalGeometry->idToDet(id)->toGlobal(tsos->localPosition()));
int strip = cscGeometry->layer(id)->geometry()->nearestStrip(hit->localPosition());
double angle = cscGeometry->layer(id)->geometry()->stripAngle(strip) - M_PI/2.;
double sinAngle = sin(angle);
double cosAngle = cos(angle);
double residual = cosAngle * (tsosChamberPos.x() - hitChamberPos.x()) + sinAngle * (tsosChamberPos.y() - hitChamberPos.y()); // yes, that's +sin()
double xx = hit->localPositionError().xx();
double xy = hit->localPositionError().xy();
double yy = hit->localPositionError().yy();
double weight = 1. / (xx*cosAngle*cosAngle + 2.*xy*sinAngle*cosAngle + yy*sinAngle*sinAngle);
double layerPosition = tsosChamberPos.z(); // the layer's position in the chamber's coordinate system
m_numHits++;
// "x" is the layerPosition, "y" is the residual (this is a linear fit to residual versus layerPosition)
m_residual_1 += weight;
m_residual_x += weight * layerPosition;
m_residual_y += weight * residual;
m_residual_xx += weight * layerPosition * layerPosition;
m_residual_xy += weight * layerPosition * residual;
// "x" is the layerPosition, "y" is chamberx (this is a linear fit to chamberx versus layerPosition)
m_trackx_1 += weight;
m_trackx_x += weight * layerPosition;
m_trackx_y += weight * tsosChamberPos.x();
m_trackx_xx += weight * layerPosition * layerPosition;
m_trackx_xy += weight * layerPosition * tsosChamberPos.x();
// "x" is the layerPosition, "y" is chambery (this is a linear fit to chambery versus layerPosition)
m_tracky_1 += weight;
m_tracky_x += weight * layerPosition;
m_tracky_y += weight * tsosChamberPos.y();
m_tracky_xx += weight * layerPosition * layerPosition;
m_tracky_xy += weight * layerPosition * tsosChamberPos.y();
m_localIDs.push_back(id);
m_localResids.push_back(residual);
m_individual_x.push_back(layerPosition);
m_individual_y.push_back(residual);
m_individual_weight.push_back(weight);
}
| double MuonCSCChamberResidual::signConvention | ( | const unsigned int | rawId = 0 | ) | const [inline, virtual] |
Implements MuonChamberResidual.
Definition at line 24 of file MuonCSCChamberResidual.h.
References MuonChamberResidual::m_chamberId, and MuonChamberResidual::m_globalGeometry.
{
DetId id = m_chamberId;
if (rawId != 0) id = DetId(rawId);
GlobalVector zDirection(0., 0., 1.);
return (m_globalGeometry->idToDet(id)->toLocal(zDirection).z() > 0. ? 1. : -1.);
};
| int MuonCSCChamberResidual::type | ( | ) | const [inline, virtual] |
Implements MuonChamberResidual.
Definition at line 18 of file MuonCSCChamberResidual.h.
References MuonChamberResidual::kCSC.
{ return MuonChamberResidual::kCSC; };
1.7.3