1
/* This file is part of MAUS: http://micewww.pp.rl.ac.uk:8080/projects/maus
3
* MAUS is free software: you can redistribute it and/or modify
4
* it under the terms of the GNU General Public License as published by
5
* the Free Software Foundation, either version 3 of the License, or
6
* (at your option) any later version.
8
* MAUS is distributed in the hope that it will be useful,
9
* but WITHOUT ANY WARRANTY; without even the implied warranty of
10
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
11
* GNU General Public License for more details.
13
* You should have received a copy of the GNU General Public License
14
* along with MAUS. If not, see <http://www.gnu.org/licenses/>.
19
#ifndef STRAIGHTTRACK_HH
20
#define STRAIGHTTRACK_HH
33
#include "src/common_cpp/Recon/Kalman/KalmanPropagator.hh"
37
/** @class KalmanStraightPropagator
39
* @brief Models the propagation of particles in the absence of fields.
41
* Inherits from KalmanTrack, where most of the functionality is defined.
42
* The methods overwritten here are the ones related with particle propagation.
45
class KalmanStraightPropagator : public KalmanPropagator {
47
/** @brief Straight Track constructor.
49
KalmanStraightPropagator();
51
/** @brief Destructor.
53
virtual ~KalmanStraightPropagator();
55
/** @brief Calculates the propagator (F), using Taylor expansion at current site.
57
void UpdatePropagator(const KalmanState *old_site, const KalmanState *new_site);
59
/** @brief Calculates the projected state.
61
void CalculatePredictedState(const KalmanState *old_site, KalmanState *new_site);
63
/** @brief Returns the momentum of the track. In the straight track case, p is just a guess.
65
double GetTrackMomentum(const KalmanState *a_site = 0);
67
TMatrixD BuildQ(TMatrixD a, double L0, double w);
69
TMatrixD GetIntermediateState(const KalmanState *old_site, double delta_z, TMatrixD &F);