TFoam is the main class of the multi-dimensional general purpose Monte Carlo event generator (integrator) FOAM.
FOAM is the simplified version of the multi-dimensional general purpose Monte Carlo event generator (integrator) FOAM. It creates hyper-rectangular "foam of cells", which is more dense around its peaks. See the following 2-dim. example of the map of 1000 cells for doubly peaked distribution:
FOAM is now fully integrated with the ROOT package. The important bonus of the ROOT use is persistency of the FOAM objects!
For more sophisticated problems full version of FOAM may be more appropriate: See full version of FOAM
Name | default | Description |
---|---|---|
kDim | 0 | Dimension of the integration space. Must be redefined! |
nCells | 1000 | No of allocated number of cells, |
nSampl | 200 | No. of MC events in the cell MC exploration |
nBin | 8 | No. of bins in edge-histogram in cell exploration |
OptRej | 1 | OptRej = 0, weighted; OptRej=1, wt=1 MC events |
OptDrive | 2 | Maximum weight reduction, =1 for variance reduction |
EvPerBin | 25 | Maximum number of the effective wt=1 events/bin, |
EvPerBin=0 deactivates this option | ||
Chat | 1 | =0,1,2 is the `‘chat level’' in the standard output |
MaxWtRej | 1.1 | Maximum weight used to get w=1 MC events |
The above can be redefined before calling Initialize() method, for instance FoamObject->SetkDim(15)
sets dimension of the distribution to 15. Only kDim
HAS TO BE redefined, the other parameters may be left at their defaults. nCell
may be increased up to about million cells for wildly peaked distributions. Increasing nSampl
sometimes helps, but it may cost CPU time. MaxWtRej
may need to be increased for wild a distribution, while using OptRej=0
.
Past versions of FOAM: August 2003, v.1.00; September 2003 v.1.01 Adopted starting from FOAM-2.06 by P. Sawicki
Users of FOAM are kindly requested to cite the following work: S. Jadach, Computer Physics Communications 152 (2003) 55.
Public Member Functions | |
TFoam () | |
Default constructor for streamer, user should not use it. | |
TFoam (const Char_t *) | |
User constructor, to be employed by the user. | |
TFoam (const TFoam &) | |
Copy Constructor NOT IMPLEMENTED (NEVER USED) | |
~TFoam () override | |
Default destructor. | |
virtual void | Carver (Int_t &, Double_t &, Double_t &) |
Internal method used by Initialize. | |
virtual Int_t | CellFill (Int_t, TFoamCell *) |
Internal method used by Initialize. | |
virtual void | CheckAll (Int_t) |
User utility, miscellaneous and debug. | |
virtual Int_t | Divide (TFoamCell *) |
Internal method used by Initialize. | |
virtual Double_t | Eval (Double_t *) |
Internal method. | |
virtual void | Explore (TFoamCell *Cell) |
Internal method used by Initialize. | |
virtual void | Finalize (Double_t &, Double_t &) |
May be called optionally by the user after the MC run. | |
virtual void | GenerCel2 (TFoamCell *&) |
Internal method. | |
virtual void | GetIntegMC (Double_t &, Double_t &) |
User method. | |
virtual void | GetIntNorm (Double_t &, Double_t &) |
User method. | |
virtual void | GetMCvect (Double_t *) |
User may get generated MC point/vector with help of this method. | |
virtual Double_t | GetMCwt () |
User may get weight MC weight using this method. | |
virtual void | GetMCwt (Double_t &) |
User may get weight MC weight using this method. | |
virtual Long_t | GetnCalls () const |
virtual Long_t | GetnEffev () const |
virtual Double_t | GetPrimary () const |
virtual void | GetPrimary (Double_t &prime) |
virtual TRandom * | GetPseRan () const |
virtual TFoamIntegrand * | GetRho () |
virtual Int_t | GetTotDim () const |
virtual const char * | GetVersion () const |
virtual void | GetWtParams (Double_t, Double_t &, Double_t &, Double_t &) |
May be called optionally after the MC run. | |
virtual void | Grow () |
Internal method used by Initialize. | |
virtual void | InitCells () |
Internal method used by Initialize. | |
virtual void | Initialize () |
Basic initialization of FOAM invoked by the user. | |
virtual void | Initialize (TRandom *, TFoamIntegrand *) |
TClass * | IsA () const override |
virtual void | LinkCells (void) |
virtual void | MakeActiveList () |
Internal method used by Initialize. | |
virtual void | MakeAlpha () |
Internal method used by Initialize. | |
virtual void | MakeEvent () |
User method. | |
virtual Double_t | MCgenerate (Double_t *MCvect) |
User method which generates MC event and returns MC weight. | |
virtual Long_t | PeekMax () |
Internal method used by Initialize. | |
virtual void | PrintCells () |
Prints geometry of ALL cells of the FOAM. | |
virtual void | ResetPseRan (TRandom *PseRan) |
User may optionally reset random number generator using this method. | |
virtual void | ResetRho (TFoamIntegrand *Rho) |
User may optionally reset the distribution using this method. | |
virtual void | RootPlot2dim (Char_t *) |
Debugging tool which plots 2-dimensional cells as rectangles in C++ format readable for root. | |
virtual void | SetChat (Int_t Chat) |
virtual void | SetEvPerBin (Int_t EvPerBin) |
virtual void | SetInhiDiv (Int_t, Int_t) |
This can be called before Initialize, after setting kDim It defines which variables are excluded in the process of the cell division. | |
virtual void | SetkDim (Int_t kDim) |
virtual void | SetMaxWtRej (Double_t MaxWtRej) |
virtual void | SetnBin (Int_t nBin) |
virtual void | SetnCells (Long_t nCells) |
virtual void | SetnSampl (Long_t nSampl) |
virtual void | SetOptDrive (Int_t OptDrive) |
virtual void | SetOptRej (Int_t OptRej) |
virtual void | SetPseRan (TRandom *PseRan) |
virtual void | SetRho (TFoamIntegrand *Rho) |
User may use this method to set the distribution object. | |
virtual void | SetRhoInt (Double_t(*fun)(Int_t, Double_t *)) |
User may use this method to set the distribution object as a global function pointer (and not as an interpreted function). | |
virtual void | SetXdivPRD (Int_t, Int_t, Double_t[]) |
This should be called before Initialize, after setting kDim It predefines values of the cell division for certain variable iDim. | |
void | Streamer (TBuffer &) override |
Stream an object of class TObject. | |
void | StreamerNVirtual (TBuffer &ClassDef_StreamerNVirtual_b) |
virtual void | Varedu (Double_t[5], Int_t &, Double_t &, Double_t &) |
Internal method used by Initialize. | |
![]() | |
TObject () | |
TObject constructor. | |
TObject (const TObject &object) | |
TObject copy ctor. | |
virtual | ~TObject () |
TObject destructor. | |
void | AbstractMethod (const char *method) const |
Use this method to implement an "abstract" method that you don't want to leave purely abstract. | |
virtual void | AppendPad (Option_t *option="") |
Append graphics object to current pad. | |
virtual void | Browse (TBrowser *b) |
Browse object. May be overridden for another default action. | |
ULong_t | CheckedHash () |
Check and record whether this class has a consistent Hash/RecursiveRemove setup (*) and then return the regular Hash value for this object. | |
virtual const char * | ClassName () const |
Returns name of class to which the object belongs. | |
virtual void | Clear (Option_t *="") |
virtual TObject * | Clone (const char *newname="") const |
Make a clone of an object using the Streamer facility. | |
virtual Int_t | Compare (const TObject *obj) const |
Compare abstract method. | |
virtual void | Copy (TObject &object) const |
Copy this to obj. | |
virtual void | Delete (Option_t *option="") |
Delete this object. | |
virtual Int_t | DistancetoPrimitive (Int_t px, Int_t py) |
Computes distance from point (px,py) to the object. | |
virtual void | Draw (Option_t *option="") |
Default Draw method for all objects. | |
virtual void | DrawClass () const |
Draw class inheritance tree of the class to which this object belongs. | |
virtual TObject * | DrawClone (Option_t *option="") const |
Draw a clone of this object in the current selected pad with: gROOT->SetSelectedPad(c1) . | |
virtual void | Dump () const |
Dump contents of object on stdout. | |
virtual void | Error (const char *method, const char *msgfmt,...) const |
Issue error message. | |
virtual void | Execute (const char *method, const char *params, Int_t *error=nullptr) |
Execute method on this object with the given parameter string, e.g. | |
virtual void | Execute (TMethod *method, TObjArray *params, Int_t *error=nullptr) |
Execute method on this object with parameters stored in the TObjArray. | |
virtual void | ExecuteEvent (Int_t event, Int_t px, Int_t py) |
Execute action corresponding to an event at (px,py). | |
virtual void | Fatal (const char *method, const char *msgfmt,...) const |
Issue fatal error message. | |
virtual TObject * | FindObject (const char *name) const |
Must be redefined in derived classes. | |
virtual TObject * | FindObject (const TObject *obj) const |
Must be redefined in derived classes. | |
virtual Option_t * | GetDrawOption () const |
Get option used by the graphics system to draw this object. | |
virtual const char * | GetIconName () const |
Returns mime type name of object. | |
virtual const char * | GetName () const |
Returns name of object. | |
virtual char * | GetObjectInfo (Int_t px, Int_t py) const |
Returns string containing info about the object at position (px,py). | |
virtual Option_t * | GetOption () const |
virtual const char * | GetTitle () const |
Returns title of object. | |
virtual UInt_t | GetUniqueID () const |
Return the unique object id. | |
virtual Bool_t | HandleTimer (TTimer *timer) |
Execute action in response of a timer timing out. | |
virtual ULong_t | Hash () const |
Return hash value for this object. | |
Bool_t | HasInconsistentHash () const |
Return true is the type of this object is known to have an inconsistent setup for Hash and RecursiveRemove (i.e. | |
virtual void | Info (const char *method, const char *msgfmt,...) const |
Issue info message. | |
virtual Bool_t | InheritsFrom (const char *classname) const |
Returns kTRUE if object inherits from class "classname". | |
virtual Bool_t | InheritsFrom (const TClass *cl) const |
Returns kTRUE if object inherits from TClass cl. | |
virtual void | Inspect () const |
Dump contents of this object in a graphics canvas. | |
void | InvertBit (UInt_t f) |
Bool_t | IsDestructed () const |
IsDestructed. | |
virtual Bool_t | IsEqual (const TObject *obj) const |
Default equal comparison (objects are equal if they have the same address in memory). | |
virtual Bool_t | IsFolder () const |
Returns kTRUE in case object contains browsable objects (like containers or lists of other objects). | |
R__ALWAYS_INLINE Bool_t | IsOnHeap () const |
virtual Bool_t | IsSortable () const |
R__ALWAYS_INLINE Bool_t | IsZombie () const |
virtual void | ls (Option_t *option="") const |
The ls function lists the contents of a class on stdout. | |
void | MayNotUse (const char *method) const |
Use this method to signal that a method (defined in a base class) may not be called in a derived class (in principle against good design since a child class should not provide less functionality than its parent, however, sometimes it is necessary). | |
virtual Bool_t | Notify () |
This method must be overridden to handle object notification (the base implementation is no-op). | |
void | Obsolete (const char *method, const char *asOfVers, const char *removedFromVers) const |
Use this method to declare a method obsolete. | |
void | operator delete (void *ptr) |
Operator delete. | |
void | operator delete (void *ptr, void *vp) |
Only called by placement new when throwing an exception. | |
void | operator delete[] (void *ptr) |
Operator delete []. | |
void | operator delete[] (void *ptr, void *vp) |
Only called by placement new[] when throwing an exception. | |
void * | operator new (size_t sz) |
void * | operator new (size_t sz, void *vp) |
void * | operator new[] (size_t sz) |
void * | operator new[] (size_t sz, void *vp) |
TObject & | operator= (const TObject &rhs) |
TObject assignment operator. | |
virtual void | Paint (Option_t *option="") |
This method must be overridden if a class wants to paint itself. | |
virtual void | Pop () |
Pop on object drawn in a pad to the top of the display list. | |
virtual void | Print (Option_t *option="") const |
This method must be overridden when a class wants to print itself. | |
virtual Int_t | Read (const char *name) |
Read contents of object with specified name from the current directory. | |
virtual void | RecursiveRemove (TObject *obj) |
Recursively remove this object from a list. | |
void | ResetBit (UInt_t f) |
virtual void | SaveAs (const char *filename="", Option_t *option="") const |
Save this object in the file specified by filename. | |
virtual void | SavePrimitive (std::ostream &out, Option_t *option="") |
Save a primitive as a C++ statement(s) on output stream "out". | |
void | SetBit (UInt_t f) |
void | SetBit (UInt_t f, Bool_t set) |
Set or unset the user status bits as specified in f. | |
virtual void | SetDrawOption (Option_t *option="") |
Set drawing option for object. | |
virtual void | SetUniqueID (UInt_t uid) |
Set the unique object id. | |
void | StreamerNVirtual (TBuffer &ClassDef_StreamerNVirtual_b) |
virtual void | SysError (const char *method, const char *msgfmt,...) const |
Issue system error message. | |
R__ALWAYS_INLINE Bool_t | TestBit (UInt_t f) const |
Int_t | TestBits (UInt_t f) const |
virtual void | UseCurrentStyle () |
Set current style settings in this object This function is called when either TCanvas::UseCurrentStyle or TROOT::ForceStyle have been invoked. | |
virtual void | Warning (const char *method, const char *msgfmt,...) const |
Issue warning message. | |
virtual Int_t | Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0) |
Write this object to the current directory. | |
virtual Int_t | Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0) const |
Write this object to the current directory. | |
Static Public Member Functions | |
static TClass * | Class () |
static const char * | Class_Name () |
static constexpr Version_t | Class_Version () |
static const char * | DeclFileName () |
![]() | |
static TClass * | Class () |
static const char * | Class_Name () |
static constexpr Version_t | Class_Version () |
static const char * | DeclFileName () |
static Longptr_t | GetDtorOnly () |
Return destructor only flag. | |
static Bool_t | GetObjectStat () |
Get status of object stat flag. | |
static void | SetDtorOnly (void *obj) |
Set destructor only flag. | |
static void | SetObjectStat (Bool_t stat) |
Turn on/off tracking of objects in the TObjectTable. | |
Protected Attributes | |
Double_t * | fAlpha |
[fDim] Internal parameters of the hyper-rectangle | |
TFoamCell ** | fCells = nullptr |
[fNCells] Array of ALL cells | |
std::vector< Long_t > | fCellsAct |
Index of active cells, constructed at the end of foam build-up. | |
Int_t | fChat |
Chat=0,1,2 chat level in output, Chat=1 normal level. | |
TString | fDate |
Release date of FOAM. | |
Int_t | fDim |
Dimension of the integration/simulation space. | |
Int_t | fEvPerBin |
Maximum number of effective (wt=1) events per bin. | |
TObjArray * | fHistDbg |
Histograms of wt, for debug. | |
TObjArray * | fHistEdg |
Histograms of wt, one for each cell edge. | |
TH1D * | fHistWt |
Histogram of the MC wt. | |
Int_t * | fInhiDiv |
! [fDim] Flags for inhibiting cell division | |
Int_t | fLastCe |
Index of the last cell. | |
Int_t * | fMaskDiv |
! [fDim] Dynamic Mask for cell division | |
Double_t | fMaxWtRej |
Maximum weight in rejection for getting wt=1 events. | |
Double_t | fMCerror |
and its error | |
TFoamMaxwt * | fMCMonit |
Monitor of the MC weight for measuring MC efficiency. | |
Double_t | fMCresult |
True Integral R from MC series. | |
Double_t * | fMCvect |
[fDim] Generated MC vector for the outside user | |
Double_t | fMCwt |
MC weight. | |
TMethodCall * | fMethodCall |
! ROOT's pointer to user-defined global distribution function | |
TString | fName |
Name of a given instance of the FOAM class. | |
Int_t | fNBin |
No. of bins in the edge histogram for cell MC exploration. | |
Long_t | fNCalls |
Total number of the function calls. | |
Int_t | fNCells |
Maximum number of cells. | |
Long_t | fNEffev |
Total number of effective events (wt=1) in the foam buildup. | |
Double_t | fNevGen |
Total number of the generated MC events. | |
Int_t | fNoAct |
Number of active cells. | |
Int_t | fNSampl |
No. of MC events, when dividing (exploring) cell. | |
Int_t | fOptDrive |
Optimization switch =1,2 for variance or maximum weight optimization. | |
Int_t | fOptPRD |
Option switch for predefined division, for quick check. | |
Int_t | fOptRej |
Switch =0 for weighted events; =1 for unweighted events in MC. | |
Double_t * | fPrimAcu |
[fNoAct] Array of cumulative probability of all active cells | |
Double_t | fPrime |
Primary integral R' (R=R'<wt> ) | |
TRandom * | fPseRan |
Pointer to user-defined generator of pseudorandom numbers. | |
TFoamIntegrand * | fRho |
! Pointer to the user-defined integrand function/distribution | |
Int_t | fRNmax |
Maximum No. of the rand. numb. requested at once. | |
Double_t * | fRvec |
[fRNmax] random number vector from r.n. generator fDim+1 maximum elements | |
Double_t | fSumOve |
Total Sum of overweighted events. | |
Double_t | fSumWt |
Double_t | fSumWt2 |
Total sum of wt and wt^2. | |
TString | fVersion |
Actual version of the FOAM like (1.01m) | |
Double_t | fWtMax |
Double_t | fWtMin |
Maximum/Minimum MC weight. | |
TFoamVect ** | fXdivPRD |
! Lists of division values encoded in one vector per direction | |
Private Member Functions | |
TFoamCell * | getCell (std::size_t i) const |
Double_t | Sqr (Double_t x) const |
Additional Inherited Members | |
![]() | |
enum | { kIsOnHeap = 0x01000000 , kNotDeleted = 0x02000000 , kZombie = 0x04000000 , kInconsistent = 0x08000000 , kBitMask = 0x00ffffff } |
enum | { kSingleKey = (1ULL << ( 0 )) , kOverwrite = (1ULL << ( 1 )) , kWriteDelete = (1ULL << ( 2 )) } |
enum | EDeprecatedStatusBits { kObjInCanvas = (1ULL << ( 3 )) } |
enum | EStatusBits { kCanDelete = (1ULL << ( 0 )) , kMustCleanup = (1ULL << ( 3 )) , kIsReferenced = (1ULL << ( 4 )) , kHasUUID = (1ULL << ( 5 )) , kCannotPick = (1ULL << ( 6 )) , kNoContextMenu = (1ULL << ( 8 )) , kInvalidObject = (1ULL << ( 13 )) } |
![]() | |
enum | { kOnlyPrepStep = (1ULL << ( 3 )) } |
![]() | |
virtual void | DoError (int level, const char *location, const char *fmt, va_list va) const |
Interface to ErrorHandler (protected). | |
void | MakeZombie () |
#include <TFoam.h>
TFoam::TFoam | ( | ) |
Internal method used by Initialize.
Determines the best edge-candidate and the position of the division plane for the future cell division, in the case of the optimization of the maximum weight. It exploits results of the cell MC exploration run stored in fHistEdg.
|
virtual |
Internal method used by Initialize.
It divides cell iCell into two daughter cells. The iCell is retained and tagged as inactive, daughter cells are appended at the end of the buffer. New vertex is added to list of vertices. List of active cells is updated, iCell removed, two daughters added and their properties set with help of MC sampling (TFoam_Explore) Returns Code RC=-1 of buffer limit is reached, fLastCe=fnBuf.
|
virtual |
Internal method used by Initialize.
It explores newly defined cell with help of special short MC sampling. As a result, estimates of true and drive volume is defined/determined Average and dispersion of the weight distribution will is found along each edge and the best edge (minimum dispersion, best maximum weight) is memorized for future use.
The optimal division point for eventual future cell division is determined/recorded. Recorded are also minimum and maximum weight etc. The volume estimate in all (inactive) parent cells is updated. Note that links to parents and initial volume = 1/2 parent has to be already defined prior to calling this routine.
|
virtual |
|
private |
User method.
It returns NORMALIZATION integral to be combined with the average weights and content of the histograms in order to get proper absolute normalization of the integrand and distributions. It can be called after initialization, before or during the MC run.
IntNorm | variable where the integral will be stored |
ErrInt | variable where the absolute error of the integral will be stored if internal rejection is active, otherwise set to 0. |
|
virtual |
|
virtual |
|
virtual |
|
inlinevirtual |
|
virtual |
|
virtual |
|
virtual |
|
virtual |
This method starts the process of the cell build-up. User must invoke Initialize with two arguments or Initialize without arguments. This is done BEFORE generating first MC event and AFTER allocating FOAM object and resetting (optionally) its internal parameters/switches. The overall operational scheme of the FOAM is the following:
InitCells initializes memory storage for cells and begins exploration process from the root cell. The empty cells are allocated/filled using CellFill. The procedure Grow which loops over cells, picks up the cell with the biggest `‘driver integral’', see Computer Physics Communications 152 152 (2003) 55 for explanations, with the help of PeekMax procedure. The chosen cell is split using Divide. Subsequently, the procedure Explore called by the Divide (and by InitCells for the root cell) does the most important job in the FOAM object build-up: it performs a small MC run for each newly allocated daughter cell. Explore calculates how profitable the future split of the cell will be and defines the optimal cell division geometry with the help of Carver or Varedu procedures, for maximum weight or variance optimization respectively. All essential results of the exploration are written into the explored cell object. At the very end of the foam build-up, Finally, MakeActiveList is invoked to create a list of pointers to all active cells, for the purpose of the quick access during the MC generation. The procedure Explore employs MakeAlpha to generate random coordinates inside a given cell with the uniform distribution. The above sequence of the procedure calls is depicted in the following figure:
|
inlineoverridevirtual |
|
virtual |
|
virtual |
|
virtual |
User method.
It generates randomly point/vector according to user-defined distribution. Prior initialization with help of Initialize() is mandatory. Generated MC point/vector is available using GetMCvect and the MC weight with GetMCwt. MC point is generated with wt=1 or with variable weight, see OptRej switch.
|
virtual |
|
virtual |
|
virtual |
User may optionally reset random number generator using this method.
Usually it is done when FOAM object is restored from the disk. IMPORTANT: this method deletes existing random number generator registered in the FOAM object. In particular such an object is created by the streamer during the disk-read operation.
|
virtual |
User may optionally reset the distribution using this method.
Usually it is done when FOAM object is restored from the disk. IMPORTANT: this method deletes existing distribution object registered in the FOAM object. In particular such an object is created by the streamer diring the disk-read operation. This method is used only in very special cases, because the distribution in most cases should be "owned" by the FOAM object and should not be replaced by another one after initialization.
|
virtual |
This can be called before Initialize, after setting kDim It defines which variables are excluded in the process of the cell division.
For example 'FoamX->SetInhiDiv(1, 1);' inhibits division of y-variable. The resulting map of cells in 2-dim. case will look as follows:
|
virtual |
This should be called before Initialize, after setting kDim It predefines values of the cell division for certain variable iDim.
For example setting 3 predefined division lines using:
results in the following 2-dim. pattern of the cells:
|
overridevirtual |
|
inline |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |
|
protected |