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

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

                                         :
       elog (Elog10),
       Emin (1.0),
       Emax (1.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;
         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 (buffer.is_valid(&JHead::cut_in) &&
                  buffer.is_valid(&JHead::livetime)) {
 
         Emin  = buffer.cut_in.Emin;
         Emax  = buffer.cut_in.Emax;
         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 67 of file JVolume.hh.

     {
       return alpha;
     }

double JAANET::JVolume::getWall ( ) const

inline

Get generation dependent weight.

Returns: weight

Definition at line 78 of file JVolume.hh.

     {
       return Wall;
     }

Double_t JAANET::JVolume::getXmin ( ) const

inline

Get minimal abscissa value.

Returns: abscissa value

Definition at line 89 of file JVolume.hh.

     {
       return getX(Emin);
     }

Double_t JAANET::JVolume::getXmax ( ) const

inline

Get maximal abscissa value.

Returns: abscissa value

Definition at line 100 of file JVolume.hh.

     {
       return getX(Emax);
     }

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 113 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 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 136 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 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 159 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 179 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 194 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 208 of file JVolume.hh.

     {
       Wall /= factor;
 
       return *this;
     }

Member Data Documentation

bool JAANET::JVolume::elog

protected

histogram option

Definition at line 217 of file JVolume.hh.

double JAANET::JVolume::Emin

protected

minimal energy

Definition at line 218 of file JVolume.hh.

double JAANET::JVolume::Emax

protected

maximal energy

Definition at line 219 of file JVolume.hh.

double JAANET::JVolume::alpha

protected

spectral index

Definition at line 220 of file JVolume.hh.

double JAANET::JVolume::Wall

protected

generation volume

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