Auxiliary class for histogramming of effective volume. More...

#include <JVolume.hh>

Public Member Functions
	JVolume (const Head &head, const bool Elog10=false)
	Constructor. More...

double	getAlpha () const
	Get spectral index of energy distribution. More...

double	getWall () const
	Get generation dependent weight. More...

Double_t	getXmin () const
	Get minimal abscissa value. More...

Double_t	getXmax () const
	Get maximal abscissa value. More...

Double_t	getX (const Double_t E, double constrain=false) const
	Get abscissa value. More...

Double_t	getE (const Double_t x, double constrain=false) const
	Get energy. More...

Double_t	getW (TAxis *axis, const Double_t E) const
	Get bin width corrected energy spectrum dependent weight. More...

double	getW (const double E) const
	Get generation dependent integral value of given energy. More...

JVolume &	mul (const double factor)
	Multiply weight. More...

JVolume &	div (const double factor)
	Divide weight. More...

Protected Attributes
bool	elog
	histogram option More...

JTOOLS::JRange< double >	E
	Energy range [GeV]. More...

double	alpha
	spectral index More...

double	Wall
	generation volume More...

Detailed Description

Auxiliary class for histogramming of effective volume.

In this, it is assumed that events are generated according $\frac{dN}{dE} \propto E^{\alpha}$ . The weight is expressed in $\mathrm{km}^{3}$ .

Definition at line 29 of file JVolume.hh.

Constructor & Destructor Documentation

JAANET::JVolume::JVolume	(	const Head &	head,
		const bool	Elog10 = `false`
	)

inline

Constructor.

Parameters

head	Monte Carlo run header
Elog10	application of log10 to energy

Definition at line 36 of file JVolume.hh.

                                         :
       elog  (Elog10),
       E(0.0, 0.0),
       alpha (0.0),
       Wall  (1.0)
     {
       const JHead buffer(head);
 
       if        (buffer.is_valid(&JHead::spectrum) && 
                  buffer.is_valid(&JHead::cut_nu)   && 
                  buffer.is_valid(&JHead::genvol)) {
 
         alpha  = buffer.spectrum.alpha;
         E      = buffer.cut_nu.E;
         Wall   = (1.0e-9 * (getW(E.getUpperLimit()) -
                             getW(E.getLowerLimit())) *
                   buffer.genvol.volume / buffer.genvol.numberOfEvents);    // km^3
 
       } else if (buffer.is_valid(&JHead::cut_in) &&
                  buffer.is_valid(&JHead::livetime)) {
 
         E      = buffer.cut_in.E;
         Wall   = 1.0 / buffer.livetime.numberOfSeconds;                     // Hz
       }
     }

Member Function Documentation

double JAANET::JVolume::getAlpha ( ) const

inline

Get spectral index of energy distribution.

Returns: spectral index

Definition at line 69 of file JVolume.hh.

     {
       return alpha;
     }

double JAANET::JVolume::getWall ( ) const

inline

Get generation dependent weight.

Returns: weight

Definition at line 80 of file JVolume.hh.

     {
       return Wall;
     }

Double_t JAANET::JVolume::getXmin ( ) const

inline

Get minimal abscissa value.

Returns: abscissa value

Definition at line 91 of file JVolume.hh.

     {
       return getX(E.getLowerLimit());
     }

Double_t JAANET::JVolume::getXmax ( ) const

inline

Get maximal abscissa value.

Returns: abscissa value

Definition at line 102 of file JVolume.hh.

     {
       return getX(E.getUpperLimit());
     }

Double_t JAANET::JVolume::getX	(	const Double_t	E,
		double	constrain = `false`
	)		const

inline

Get abscissa value.

Parameters

E	energy
constrain	constrain

Returns: abscissa value

Definition at line 115 of file JVolume.hh.

     {
       double x = (elog ? log10(E) : E);
 
       if (constrain) {
         if        (x < getXmin()) {
           return getXmin();
         } else if (x > getXmax()) {
           return getXmax();
         }
       }
 
       return x;
     }

Double_t JAANET::JVolume::getE	(	const Double_t	x,
		double	constrain = `false`
	)		const

inline

Get energy.

Parameters

x	abscissa value
constrain	constrain

Returns: energy

Definition at line 138 of file JVolume.hh.

     {
       const double Ex = (elog ? pow(10.0, x) : x);
 
       if (constrain)
         return E.constrain(Ex);
       else
         return Ex;
     }

Double_t JAANET::JVolume::getW	(	TAxis *	axis,
		const Double_t	E
	)		const

inline

Get bin width corrected energy spectrum dependent weight.

Parameters

axis	axis
E	energy

Returns: weight [km^3]

Definition at line 156 of file JVolume.hh.

     {
       const Int_t    index = axis->FindBin(getX(E));
 
       const Double_t xmin  = axis->GetBinLowEdge(index);
       const Double_t xmax  = axis->GetBinUpEdge (index);
       
       const Double_t Wmin  = getW(getE(xmin));
       const Double_t Wmax  = getW(getE(xmax));
 
       return Wall / (Wmax - Wmin);
     }

double JAANET::JVolume::getW ( const double E ) const

inline

Get generation dependent integral value of given energy.

Parameters

E energy

Returns: integral value

Definition at line 176 of file JVolume.hh.

     {
       if (alpha != -1.0)
         return pow(E, 1.0 + alpha) / (1.0 + alpha);
       else
         return log(E);
     }

JVolume& JAANET::JVolume::mul ( const double factor )

inline

Multiply weight.

Parameters

factor factor

Returns: this volume

Definition at line 191 of file JVolume.hh.

     {
       Wall *= factor;
 
       return *this;
     }

JVolume& JAANET::JVolume::div ( const double factor )

inline

Divide weight.

Parameters

factor factor

Returns: this volume

Definition at line 205 of file JVolume.hh.

     {
       Wall /= factor;
 
       return *this;
     }

Member Data Documentation

bool JAANET::JVolume::elog

protected

histogram option

Definition at line 214 of file JVolume.hh.

JTOOLS::JRange<double> JAANET::JVolume::E

protected

Energy range [GeV].

Definition at line 215 of file JVolume.hh.

double JAANET::JVolume::alpha

protected

spectral index

Definition at line 216 of file JVolume.hh.

double JAANET::JVolume::Wall

protected

generation volume

Definition at line 217 of file JVolume.hh.

The documentation for this struct was generated from the following file:

JVolume.hh

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation