Implementation of the JGenerator interface. More...

#include <JMarkovGenerator.hh>

Inheritance diagram for JMARKOV::JHistGenerator:

Public Member Functions
	JHistGenerator (TH1 _hx, TH1 _hy, TH1 *_hz)
	Constructor. More...

	~JHistGenerator ()

JPosition3D	getPosition ()
	Return a randomly generated position. More...

double	getWeight (JPosition3D pos)
	return the weight (=probability density dP/dV) for the given position. More...

Public Attributes
TH1 *	hx

TH1 *	hy

TH1 *	hz

Detailed Description

Implementation of the JGenerator interface.

Generates positions by sampling from 3 histograms that represent the projection of the distribution along the X, Y and Z axis.

Definition at line 521 of file JMarkovGenerator.hh.

Constructor & Destructor Documentation

◆ JHistGenerator()

JMARKOV::JHistGenerator::JHistGenerator	(	TH1 *	_hx,
		TH1 *	_hy,
		TH1 *	_hz
	)

inline

Constructor.

The histograms are cloned and automatically normalized.

Note: empty bins should always be avoided!

Definition at line 530 of file JMarkovGenerator.hh.

                                                    {
       // copy histograms
       hx = (TH1*) _hx->Clone( _hx->GetName() ) ;
       hy = (TH1*) _hy->Clone( _hy->GetName() ) ;
       hz = (TH1*) _hz->Clone( _hz->GetName() ) ;
       // normalize histograms
       hx->Scale( 1.0/hx->Integral("width") ) ;
       hy->Scale( 1.0/hy->Integral("width") ) ;
       hz->Scale( 1.0/hz->Integral("width") ) ;
     }

◆ ~JHistGenerator()

JMARKOV::JHistGenerator::~JHistGenerator ( )

inline

Definition at line 541 of file JMarkovGenerator.hh.

                       {
       delete hx ;
       delete hy ;
       delete hz ;
     }

Member Function Documentation

◆ getPosition()

JPosition3D JMARKOV::JHistGenerator::getPosition ( )

inlinevirtual

Return a randomly generated position.

Implements JMARKOV::JGenerator.

Definition at line 547 of file JMarkovGenerator.hh.

                               {
       double x = hx->GetRandom() ;
       double y = hy->GetRandom() ;
       double z = hz->GetRandom() ;
       return JPosition3D(x,y,z) ;
     }

◆ getWeight()

double JMARKOV::JHistGenerator::getWeight ( JPosition3D pos )

inlinevirtual

return the weight (=probability density dP/dV) for the given position.

When the generator is properly normalized, the integral of this quantity over the whole space is 1.

Implements JMARKOV::JGenerator.

Definition at line 554 of file JMarkovGenerator.hh.

                                         {
       Int_t bin ;
       double w = 1 ;
       
       bin = hx->GetXaxis()->FindBin( pos.getX() ) ;
       if( bin==0 || bin==hx->GetNbinsX()+1 ) return 0.0 ; // out of range
       w *= hx->GetBinContent(bin) ;
  
       bin = hy->GetXaxis()->FindBin( pos.getY() ) ;
       if( bin==0 || bin==hy->GetNbinsX()+1 ) return 0.0 ; // out of range
       w *= hy->GetBinContent(bin) ;
  
       bin = hz->GetXaxis()->FindBin( pos.getZ() ) ;
       if( bin==0 || bin==hz->GetNbinsX()+1 ) return 0.0 ; // out of range
       w *= hz->GetBinContent(bin) ;
  
       return w ;
     }