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