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