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 (const MultiState & state);
  
  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:
  const MultiState theInitialState;
//   bool sortStates;
  double minFractionalWeight;

  bool combinationDone;

  double theValidWeightSum;
  SingleStateContainer theStates;
//   static TimingReport::Item * theTimerAdd;
//   static TimingReport::Item * theTimerComb;
  
};

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

#endif