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...

double	Emin
	minimal energy More...

double	Emax
	maximal energy 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 26 of file JVolume.hh.

Constructor & Destructor Documentation

JGIZMO::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 33 of file JVolume.hh.

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

Member Function Documentation

double JGIZMO::JVolume::getAlpha ( ) const

inline

Get spectral index of energy distribution.

Returns: spectral index

Definition at line 69 of file JVolume.hh.

     {
       return alpha;
     }

double JGIZMO::JVolume::getWall ( ) const

inline

Get generation dependent weight.

Returns: weight

Definition at line 80 of file JVolume.hh.

     {
       return Wall;
     }

Double_t JGIZMO::JVolume::getXmin ( ) const

inline

Get minimal abscissa value.

Returns: abscissa value

Definition at line 91 of file JVolume.hh.

     {
       return getX(Emin);
     }

Double_t JGIZMO::JVolume::getXmax ( ) const

inline

Get maximal abscissa value.

Returns: abscissa value

Definition at line 102 of file JVolume.hh.

     {
       return getX(Emax);
     }

Double_t JGIZMO::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.

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

Double_t JGIZMO::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 E = (elog ? pow(10.0, x) : x);
 
       if (constrain) {
         if        (E < Emin) { 
           return Emin;
         } else if (E > Emax) {
           return Emax;
         }
       }
 
       return E;
     }

Double_t JGIZMO::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 161 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 JGIZMO::JVolume::getW ( const double E ) const

inline

Get generation dependent integral value of given energy.

Parameters

E energy

Returns: integral value

Definition at line 181 of file JVolume.hh.

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

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

inline

Multiply weight.

Parameters

factor factor

Returns: this volume

Definition at line 196 of file JVolume.hh.

     {
       Wall *= factor;
 
       return *this;
     }

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

inline

Divide weight.

Parameters

factor factor

Returns: this volume

Definition at line 210 of file JVolume.hh.

     {
       Wall /= factor;
 
       return *this;
     }

Member Data Documentation

bool JGIZMO::JVolume::elog

protected

histogram option

Definition at line 219 of file JVolume.hh.

double JGIZMO::JVolume::Emin

protected

minimal energy

Definition at line 220 of file JVolume.hh.

double JGIZMO::JVolume::Emax

protected

maximal energy

Definition at line 221 of file JVolume.hh.

double JGIZMO::JVolume::alpha

protected

spectral index

Definition at line 222 of file JVolume.hh.

double JGIZMO::JVolume::Wall

protected

generation volume

Definition at line 223 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