MultiGaussianState.h

Go to the documentation of this file.

#ifndef MultiGaussianState_H
#define MultiGaussianState_H

#include "TrackingTools/GsfTools/interface/SingleGaussianState.h"
#include <memory>

#include <vector>

// #include <iostream>
// #include <signal.h>

template <unsigned int N>
class MultiGaussianState {
public:
  typedef typename SingleGaussianState<N>::Vector Vector;
  typedef typename SingleGaussianState<N>::Matrix Matrix;
  typedef SingleGaussianState<N> SingleState;
  typedef std::shared_ptr<SingleState> SingleStatePtr;
//   typedef std::vector< std::shared_ptr<const SingleState> > SingleStateContainer;
  typedef std::vector< SingleStatePtr > SingleStateContainer;

public:

  MultiGaussianState() : theCombinedStateUp2Date(false) {
//     ++instances_;++maxInstances_;
//     std::cout << "MultiGaussianState() " << N << " " << instances_ << std::endl;
  }

  MultiGaussianState(const SingleStateContainer& stateV)
    : theComponents(stateV), theCombinedStateUp2Date(false) {
//     theComponents[0]->rescaleWeight(1.);
//     ++instances_;++maxInstances_;
//     std::cout << "MultiGaussianState(const SingleStateContainer&) " << N << " " 
//        << instances_ << std::endl;
  }

//   MultiGaussianState(const MultiGaussianState<N>& rhs) :
//     theComponents(rhs.theComponents), theCombinedState(rhs.theCombinedState),
//     theCombinedStateUp2Date(rhs.theCombinedStateUp2Date) {
//     ++instances_;++maxInstances_;
//     std::cout << "MultiGaussianState(const MultiGaussianState<N>&) " << N << " " 
//        << instances_ << std::endl;
//   }

  ~MultiGaussianState() {
//     --instances_;
//     std::cout << "~MultiGaussianState " << N << " " << instances_ << std::endl;
  }

//   /**
//    * Creates a new multi-state with the given components.
//    * For this base class, no information is passed from the initial
//    * instance.
//    */
//   virtual MultiGaussianState createState(
//  const std::vector<SingleGaussianState> & stateV) const {
//     return MultiGaussianState(stateV);
//   }

//   /**
//    * Creates a new single-state with the given information.
//    * For this base class, no information is passed from the initial
//    * instance.
//    */
//   virtual SingleGaussianState createSingleState (
//  const AlgebraicVector & aMean, const AlgebraicSymMatrix & aCovariance,
//  double aWeight = 1.) const {
//     return SingleGaussianState(aMean, aCovariance, aWeight);
//   }

  double weight() const;
  const Vector & mean() const;
  const Matrix & covariance() const;
  const Matrix & weightMatrix() const;
  inline const SingleStateContainer& components() const {return theComponents;}
  int dimension () const {
    return N;
  }
  void setWeight (double newWeight);
  void rescaleWeight (double scale);

// protected:
private:
  void checkCombinedState() const;

//   std::vector<SingleState> theComponents;
// should become a vector of pointers to const SingleState ...
  const SingleStateContainer theComponents;
  mutable SingleStatePtr theCombinedState;
  mutable bool theCombinedStateUp2Date;

// public:
//   static int instances_;
//   static int maxInstances_;
//   static int constructsCombinedState_;
};

//NOTE: Circular dependency between MultiGaussianState and
// MultiGaussianStateCombiner requires they be in the same
// headerfile

template <unsigned int N>
class MultiGaussianStateCombiner {

private:
  typedef SingleGaussianState<N> SingleState;
  typedef MultiGaussianState<N> MultiState;
  typedef typename MultiGaussianState<N>::SingleStatePtr SingleStatePtr;
  typedef typename MultiGaussianState<N>::SingleStateContainer VSC;
  
public:

//   typedef std::vector<SingleState> VSC;

  SingleStatePtr combine(const MultiState & theState) const;
  SingleStatePtr combine(const VSC& theComponents) const;

};


#include "TrackingTools/GsfTools/interface/MultiGaussianStateCombiner.icc"
#include "TrackingTools/GsfTools/interface/MultiGaussianState.icc"

//   template <unsigned int N> int MultiGaussianState<N>::instances_ = 0;
//   template <unsigned int N> int MultiGaussianState<N>::maxInstances_ = 0;
//   template <unsigned int N> int MultiGaussianState<N>::constructsCombinedState_ = 0;

#endif