d9/d2c/CylinderKalman_8h_source.html

#ifndef KALMAN_CYLINDERKALMAN_H

#define KALMAN_CYLINDERKALMAN_H


#include "HelixKalman.h"


#include <Eigen/Core>


#include <vector>


class HelixKalmanState;

class SimpleHit3D;


class CylinderKalman : public HelixKalman

{

  public:

    CylinderKalman(std::vector<float>& detector_radii, std::vector<float>& detector_material, float B);

    virtual ~CylinderKalman();


  protected:

    void calculateProjections(SimpleHit3D& hit, HelixKalmanState& state, Eigen::Matrix<float,2,5>& H, Eigen::Matrix<float,2,1>& ha);

    void calculateMeasurements(SimpleHit3D& hit, Eigen::Matrix<float,2,1>& m, Eigen::Matrix<float,2,2>& G);

    bool calculateScatteringVariance(HelixKalmanState& state, float& var);

    void updateIntersection(HelixKalmanState& state, int layer);

    void subtractProjections(Eigen::Matrix<float,2,1>& m, Eigen::Matrix<float,2,1>& ha, Eigen::Matrix<float,2,1>& diff);


  private:

    void calculate_dxda(SimpleHit3D& hit, HelixKalmanState& state, Eigen::Matrix<float,3,5>& dxda, float& x, float& y, float& z);

    std::vector<float> det_rad;

    std::vector<float> det_scatter_variance;

    unsigned int nlayers;

    float signk_store;

};


#endif