MultiGaussianStateAssembler.h

Go to the documentation of this file.

#ifndef MultiGaussianStateAssembler_h_
#define MultiGaussianStateAssembler_h_

#include "TrackingTools/GsfTools/interface/SingleGaussianState.h"
#include "TrackingTools/GsfTools/interface/MultiGaussianState.h"

#include <vector>

template <unsigned int N>
class MultiGaussianStateAssembler {
private:
  typedef SingleGaussianState<N> SingleState;
  typedef MultiGaussianState<N> MultiState;
  typedef typename MultiGaussianState<N>::SingleStatePtr SingleStatePtr;
  typedef typename MultiGaussianState<N>::SingleStateContainer SingleStateContainer;

public:
  //
  // constructors
  //
  MultiGaussianStateAssembler() = default;

  void addState(const MultiState& state);
  void addState(const SingleStatePtr& state);

  MultiState combinedState();
  MultiState combinedState(const float weight);

private:
  void addStateVector(const SingleStateContainer&);

  bool prepareCombinedState();

  MultiState reweightedCombinedState(const double) const;

  void removeSmallWeights();

private:
  double minFractionalWeight = 1.e-16;

  bool combinationDone = false;

  double theValidWeightSum = 0;
  SingleStateContainer theStates;
};

#include "TrackingTools/GsfTools/interface/MultiGaussianStateAssembler.icc"

#endif