EventShapeVariables.h

Go to the documentation of this file.

#ifndef EventShapeVariables_h
#define EventShapeVariables_h

#include "DataFormats/Math/interface/Vector3D.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"

#include "TMatrixDSym.h"
#include "TVectorD.h"

#include <vector>

class EventShapeVariables {

 public:
  explicit EventShapeVariables(const edm::View<reco::Candidate>& inputVectors);  
  explicit EventShapeVariables(const std::vector<math::XYZVector>& inputVectors);  
  explicit EventShapeVariables(const std::vector<math::RhoEtaPhiVector>& inputVectors);  
  explicit EventShapeVariables(const std::vector<math::RThetaPhiVector>& inputVectors);  
  ~EventShapeVariables(){};

  double isotropy(const unsigned int& numberOfSteps = 1000) const;
  
  double circularity(const unsigned int& numberOfSteps = 1000) const;

  double sphericity(double = 2.)  const;
  double aplanarity(double = 2.)  const;
  double C(double = 2.) const;
  double D(double = 2.) const;
  
 private:
  TMatrixDSym compMomentumTensor(double = 2.) const;
  TVectorD compEigenValues(double = 2.) const;

  std::vector<math::XYZVector> inputVectors_;
};

#endif