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