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.

	~JHistGenerator ()

JPosition3D	getPosition ()
	Return a randomly generated position.

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

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 ;
    }