ACDCGenCell.h

// -*- C++ -*-
//
// ACDCGenCell.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ACDCGenCell_H
#define ACDCGenCell_H

#include "ACDCGenConfig.h"
#include "ACDCTraits.h"

namespace ACDCGenerator {

struct ACDCGenCellInfo;

class ACDCGenCell {

public:

  inline ACDCGenCell(double newG);

  inline ACDCGenCell(double newG, double newV);

  inline ~ACDCGenCell();

  template <typename RndType>
  inline ACDCGenCell * generate(DVector & lo, DVector & up, RndType * rnd);

  inline ACDCGenCell * generate(DVector & lo, DVector & up, DVector & rndv);

  inline ACDCGenCell * getCell(DVector & lo, const DVector & x, DVector &  up);

  inline void smooth(double frac);

  inline bool isSplit() const;

  inline double doMaxInt(double rescale = 1.0);

  inline double maxInt() const;

  inline void splitme(double lo, double newDiv, double up, DimType newDim);

  inline void g(double newG);

  inline double g() const;

  inline double v() const;

  inline DimType dim() const;

  inline double div() const;

  inline ACDCGenCell * upper() const;

  inline ACDCGenCell * lower() const;

  inline int nBins() const;

  inline int depth() const;

  inline void extract(DVector & lo, DVector & up,
                      vector<ACDCGenCellInfo> & v) const;

  inline long getIndex(const ACDCGenCell * c) const;

  inline long getIndex(const ACDCGenCell * c, long & indx) const;
   
  inline ACDCGenCell * getCell(long i);

  inline ACDCGenCell * getCell(long i, long & indx);

public:

  double theG;

  double theV;

  ACDCGenCell * theUpper;

  ACDCGenCell * theLower;

  double theDivision;

  DimType theSplitDimension;

private:

  ACDCGenCell();

  ACDCGenCell(const ACDCGenCell &);

  ACDCGenCell & operator=(const ACDCGenCell &) = delete;

};


struct ACDCGenCellInfo {

  typedef vector<ACDCGenCellInfo>::size_type Index;

  double g;

  double v;

  DVector up;
  DVector lo;

  Index iup;

  Index ilo;

};

}

#include "ACDCGenCell.icc"

#endif