/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/RecoVertex/GaussianSumVertexFit/src/BasicMultiVertexState.cc

Go to the documentation of this file.

#include "RecoVertex/GaussianSumVertexFit/interface/BasicMultiVertexState.h"
#include "RecoVertex/VertexPrimitives/interface/VertexException.h"

using namespace std;

BasicMultiVertexState::
BasicMultiVertexState(const vector<VertexState>& vsComp) :
  valid(true), theComponents(vsComp), theCombinedStateUp2Date( false) {}


GlobalPoint BasicMultiVertexState::position() const
{
  checkCombinedState();
  return theCombinedState.position();
}


GlobalError BasicMultiVertexState::error() const
{
  checkCombinedState();
  return theCombinedState.error();
}


GlobalWeight BasicMultiVertexState::weight() const
{
  checkCombinedState();
  return theCombinedState.weight();
}


AlgebraicVector3 BasicMultiVertexState::weightTimesPosition() const
{
  checkCombinedState();
  return theCombinedState.weightTimesPosition();
}


// RefCountedVertexSeed BasicMultiVertexState::seedWithoutTracks() const
// {
//   checkCombinedState();
//   return theCombinedState.seedWithoutTracks();
// }

double BasicMultiVertexState::weightInMixture() const {
  if (!valid) throw VertexException("BasicSingleVertexState::invalid");
  if (theComponents.empty()) {
    cout << "Asking for weight of empty MultiVertexState, returning zero!" << endl;
    throw VertexException("Asking for weight of empty MultiVertexState, returning zero!");
    return 0.;
  }

  double weight = 0.;
  for (vector<VertexState>::const_iterator it = theComponents.begin(); 
        it != theComponents.end(); it++) {
    weight += it->weightInMixture();
  }
  return weight;
}

void BasicMultiVertexState::checkCombinedState() const
{
  if (!valid) throw VertexException("BasicSingleVertexState::invalid");
  if (theCombinedStateUp2Date) return;

  theCombinedState = theCombiner.combine(theComponents);
  theCombinedStateUp2Date = true;
}