Commit de79bb99 authored by Paul Schuetze's avatar Paul Schuetze
Browse files

Give multiplet object a member that states the position at the scatterer

parent cc0c16df
......@@ -33,6 +33,9 @@ void Multiplet::fit() {
// FIXME: Currently asking for direction of "". Should be the last detector plane -> Would enable using more generic
// tracks
m_positionAtScatterer = ((m_downstream->intercept(m_scattererPosition) -
(ROOT::Math::XYZPoint(0, 0, 0) - m_upstream->intercept(m_scattererPosition))) /
2.);
m_offsetAtScatterer = m_downstream->intercept(m_scattererPosition) - m_upstream->intercept(m_scattererPosition);
// Calculate the angle
......@@ -43,6 +46,7 @@ void Multiplet::fit() {
m_kinkAtScatterer = ROOT::Math::XYVector(slopeXdown - slopeXup, slopeYdown - slopeYup);
this->calculateChi2();
m_isFitted = true;
}
ROOT::Math::XYZPoint Multiplet::intercept(double z) const {
......
......@@ -82,6 +82,7 @@ namespace corryvreckan {
void calculateChi2();
double m_scattererPosition;
ROOT::Math::XYZVector m_positionAtScatterer;
ROOT::Math::XYVector m_offsetAtScatterer;
ROOT::Math::XYVector m_kinkAtScatterer;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment