Regressor function object for JShower3EZ fit using JSimplex minimiser. More...

#include <JShowerEnergyRegressor.hh>

Inheritance diagram for JFIT::JRegressor< JEnergy, JSimplex >:

Public Types
typedef JTOOLS::JSplineFunction1S_t	JFunction1D_t

typedef JTOOLS::JMAPLIST< JTOOLS::JPolint1FunctionalMap, JTOOLS::JPolint1FunctionalGridMap >::maplist	JMapList_t

typedef JPHYSICS::JPDFTable< JFunction1D_t, JMapList_t >	JPDF_t

typedef JPHYSICS::JNPETable< double, double, JMapList_t >	JNPE_t

typedef JSimplex< JEnergy >	minimiser_type

typedef JRegressor< JEnergy, JSimplex >	regressor_type

typedef minimiser_type::result_type	result_type

Public Member Functions
	JRegressor (const std::string &fileDescriptor)
	Parameterized constructor. More...

double	operator() (const JEnergy &x, const JShowerNPEHit &npe) const
	Fit function. More...

JShowerNPE	getNPE (const JAxis3D &axis, const double R_Hz) const
	Create data structure for handling light yields for PMT. More...

result_type	operator() (const JEnergy &value, T __begin, T __end)
	Global fit. More...

template<class JFunction_t , class T >
double	operator() (const JFunction_t &fit, T __begin, T __end)
	Multi-dimensional fit. More...

template<class JFunction_t , class T >
double	operator() (const JFunction_t &fit, T __begin, T __end, const JModel_t &step)
	1D fit. More...

Static Public Member Functions
static double	getNPE (const std::vector< JNPE_t > &NPE, const double D, const double cd)
	Get number of photo-electrons. More...

Public Attributes
std::vector< JNPE_t >	Y
	light from EM showers More...

JLANG::JSharedPointer< JMEstimator >	estimator
	M-Estimator function. More...

JModel_t	value

std::vector< JModel_t >	step

int	numberOfIterations

Static Public Attributes
static JTimeRange	T_ns
	Time window with respect to Cherenkov hypothesis [ns]. More...

static double	Vmax_npe = std::numeric_limits<double>::max()
	Maximal integral of PDF [npe]. More...

static const int	NUMBER_OF_PDFS = 2

static const JPDFType_t	pdf_t [NUMBER_OF_PDFS]
	PDF types. More...

static int	MAXIMUM_ITERATIONS = 1000
	maximal number of iterations More...

static double	EPSILON = 1.0e-4
	maximal distance to minimum More...

static int	debug = 0
	debug level (default is off). More...

Private Member Functions
template<class JFunction_t , class T >
double	evaluate (const JFunction_t &fit, T __begin, T __end) const
	Evaluate chi2 for given data set. More...

Private Attributes
JModel_t	p0

JModel_t	p1

JModel_t	wall

Detailed Description

Regressor function object for JShower3EZ fit using JSimplex minimiser.

Definition at line 50 of file JShowerEnergyRegressor.hh.

Member Typedef Documentation

◆ JFunction1D_t

typedef JTOOLS::JSplineFunction1S_t JFIT::JRegressor< JEnergy, JSimplex >::JFunction1D_t

Definition at line 55 of file JShowerEnergyRegressor.hh.

◆ JMapList_t

typedef JTOOLS::JMAPLIST<JTOOLS::JPolint1FunctionalMap, JTOOLS::JPolint1FunctionalGridMap>::maplist JFIT::JRegressor< JEnergy, JSimplex >::JMapList_t

Definition at line 57 of file JShowerEnergyRegressor.hh.

◆ JPDF_t

typedef JPHYSICS::JPDFTable<JFunction1D_t, JMapList_t> JFIT::JRegressor< JEnergy, JSimplex >::JPDF_t

Definition at line 58 of file JShowerEnergyRegressor.hh.

◆ JNPE_t

typedef JPHYSICS::JNPETable<double, double, JMapList_t> JFIT::JRegressor< JEnergy, JSimplex >::JNPE_t

Definition at line 60 of file JShowerEnergyRegressor.hh.

◆ minimiser_type

typedef JSimplex <JEnergy > JFIT::JAbstractRegressor< JEnergy , JSimplex >::minimiser_type

inherited

Definition at line 80 of file JRegressor.hh.

◆ regressor_type

typedef JRegressor<JEnergy , JSimplex > JFIT::JAbstractRegressor< JEnergy , JSimplex >::regressor_type

inherited

Definition at line 81 of file JRegressor.hh.

◆ result_type

typedef minimiser_type::result_type JFIT::JAbstractRegressor< JEnergy , JSimplex >::result_type

inherited

Definition at line 82 of file JRegressor.hh.

Constructor & Destructor Documentation

◆ JRegressor()

JFIT::JRegressor< JEnergy, JSimplex >::JRegressor ( const std::string & fileDescriptor )

inline

Parameterized constructor.

The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD which will be replaced by the corresponding PDF types listed in JRegressor<JEnergy, JSimplex>::pdf_t.

Parameters

fileDescriptor PDF file descriptor

Definition at line 71 of file JShowerEnergyRegressor.hh.

                                                :
       estimator(new JMEstimatorNull())
     {
       using namespace std;
       using namespace JPP;
       
       
       const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(JMATH::zero));
       
       for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
         
         try {
  
           JPDF_t pdf;
  
           const string file_name = getFilename(fileDescriptor, pdf_t[i]);
           
           NOTICE("loading PDF from file " << file_name << "... " << flush);
  
           pdf.load(file_name.c_str());
  
           NOTICE("OK" << endl);
  
           pdf.setExceptionHandler(supervisor);
                   
           Y.push_back(JNPE_t(pdf));
         }
         catch(const JException& error) {
           FATAL(error.what() << endl);
         }
       }
  
       // Add PDFs
  
       for (int i = 1; i < NUMBER_OF_PDFS; i += 2) {
         Y[i].add(Y[i-1]);
       }    
  
       Y.erase(Y.begin());
     }

Member Function Documentation

◆ operator()() [1/4]

double JFIT::JRegressor< JEnergy, JSimplex >::operator()	(	const JEnergy &	x,
		const JShowerNPEHit &	npe
	)		const

inline

Fit function.

This method is used to determine the chi2 of given PMT with respect to shower hypothesis.

Parameters

x	energy
npe	number of photoelectrons

Returns: chi2

Definition at line 120 of file JShowerEnergyRegressor.hh.

     {      
       using namespace JPP;
  
       const double E = x.getE();
       const double u = npe.getChi2(E);
  
       return estimator->getRho(u);
     }

◆ getNPE() [1/2]

JShowerNPE JFIT::JRegressor< JEnergy, JSimplex >::getNPE	(	const JAxis3D &	axis,
		const double	R_Hz
	)		const

inline

Create data structure for handling light yields for PMT.

Parameters

axis	PMT axis
R_Hz	singles rate [Hz]

Returns: light yields

Definition at line 137 of file JShowerEnergyRegressor.hh.

     {
       using namespace JPP;
  
       const JPosition3D  D(axis.getPosition());
       const JDirection3D U(axis.getDirection());      
  
       const double U_length = std::sqrt(U.getDX()*U.getDX() + 
                                         U.getDY()*U.getDY() + 
                                         U.getDZ()*U.getDZ());
  
       const double ct = U.getDot(D) / (D.getLength()*U_length);
       
       const double y = getNPE(Y, D.getLength(), ct);
  
       return JShowerNPE(getN(T_ns, R_Hz * 1.0e-9), y);      
     }

◆ getNPE() [2/2]

static double JFIT::JRegressor< JEnergy, JSimplex >::getNPE	(	const std::vector< JNPE_t > &	NPE,
		const double	D,
		const double	cd
	)

inlinestatic

Get number of photo-electrons.

Parameters

NPE	NPE tables
D	PMT distance from shower [m]
cd	cosine of the PMT angle wrt the photon direction

Returns: number of photo-electrons

Definition at line 165 of file JShowerEnergyRegressor.hh.

      {
        double npe = 0.0;
        
        for (std::vector<JNPE_t>::const_iterator i = NPE.begin(); i != NPE.end(); ++i) {
          
          if (D <= i->getXmax()) {
            
            try {
              
              const double y = get_value((*i)(std::max(D, i->getXmin()), cd));
              
              if (y > 0.0) {
                npe += y;
              }
            }
            catch(const JLANG::JException& error) {
              ERROR(error << std::endl);
            }
          }
        }
        
        return npe;
      }    

◆ operator()() [2/4]

result_type JFIT::JAbstractRegressor< JEnergy , JSimplex >::operator()	(	const JEnergy &	value,
		T	__begin,
		T	__end
	)

inlineinherited

Global fit.

Parameters

value	start value
__begin	begin of data set
__end	end of data set

Returns: chi2

Definition at line 94 of file JRegressor.hh.

     {
       static_cast<minimiser_type&>(*this).value = value;
  
       return static_cast<minimiser_type&>(*this)(static_cast<regressor_type&>(*this), __begin, __end);
     }

◆ operator()() [3/4]

template<class JModel_t >

template<class JFunction_t , class T >

double JFIT::JSimplex< JModel_t >::operator()	(	const JFunction_t &	fit,
		T	__begin,
		T	__end
	)

inlineinherited

Multi-dimensional fit.

The given fit function should return the equivalent of chi2 for the current value of the given model and a given data point.

Parameters

fit	fit function
__begin	begin of data
__end	end of data

Returns: chi2

Definition at line 71 of file JSimplex.hh.

     {
       using namespace std;
       using namespace JPP;
  
       double chi2_old = evaluate(fit, __begin, __end);
  
       const int N = step.size();
  
       if (N != 0) {
  
         p0   = value;
         p1   = value;
         wall = value;
  
         double chi2[N];
  
         numberOfIterations = 0;
  
         while (numberOfIterations != MAXIMUM_ITERATIONS) {
  
           DEBUG("old: " << FIXED(12,5) << chi2_old << endl);
  
           p0 = value;
  
           for (int i = 0; i != N; ++i) {
  
             DEBUG("step: " << i << ' ' << setw(5) << numberOfIterations << endl);
  
             chi2[i] = (*this)(fit, __begin, __end, step[i]);
           }
  
           // overall step direction of last iteration
  
           wall  = value;
           wall -= p0;
  
           const double chi2_new = (*this)(fit, __begin, __end, wall);
  
           DEBUG("new: " << FIXED(12,5) << chi2_new << endl);
  
  
           if (fabs(chi2_old - chi2_new) < EPSILON*fabs(chi2_old)) {
             return chi2_new;
           }
  
           // double overall step
  
           wall   = value;
           wall  -= p0;
           value += wall;
  
           const double fe = evaluate(fit, __begin, __end);
  
           value -= wall;
  
  
           for (int i = N-1; i != 0; --i) {
             chi2[i] = chi2[i-1] - chi2[i];
           }
  
           chi2[0] = chi2_old - chi2[0];
  
  
           double df = 0.0;
  
           for (int i = 0; i != N; ++i) {
             if (chi2[i] > df) {
               df = chi2[i];
             }
           }
  
           const double fn = chi2_new;
           const double f0 = chi2_old;
           const double ff = f0 - fn - df;
  
           // shift step directions
  
           if (fe < f0 && 2.0*(f0 - 2.0*fn + fe)*ff*ff < (f0-fe)*(f0-fe)*df) {
  
             for (int i = 0; i != N - 1; ++i) {
               step[i] = step[i+1];
             }
  
             step[N-1] = wall;
           }
  
           chi2_old = chi2_new;
         }
       }
  
       return chi2_old;
     }

◆ operator()() [4/4]

template<class JModel_t >

template<class JFunction_t , class T >

double JFIT::JSimplex< JModel_t >::operator()	(	const JFunction_t &	fit,
		T	__begin,
		T	__end,
		const JModel_t &	step
	)

inlineinherited

1D fit.

The given fit function should return the equivalent of chi2 for the current value of the given model and a given data point.

Parameters

fit	fit function
__begin	begin of data
__end	end of data
step	step direction

Definition at line 178 of file JSimplex.hh.

     {
       using namespace std;
       using namespace JPP;
  
       double lambda = 0.5;                                     // control parameter
       double factor = 1.0;                                     // multiplication factor
  
       double chi2_old = evaluate(fit, __begin, __end);
  
       for (int i = 0; numberOfIterations != MAXIMUM_ITERATIONS; ++numberOfIterations, ++i) {
  
         p1      = step;
         p1     *= lambda;
         value  += p1;
  
         const double chi2_new = evaluate(fit, __begin, __end);
  
         DEBUG("step: " << setw(3) << i << ' ' << FIXED(12,5) << chi2_old << ' ' << FIXED(12,5) << chi2_new << ' ' << FIXED(5,2) << lambda << endl);
  
         if (fabs(chi2_old - chi2_new) < EPSILON*fabs(chi2_old)) {
  
           if (chi2_new > chi2_old) {
  
             p1      = step;
             p1     *= lambda;
             value  -= p1;                                      // undo last step
  
             return chi2_old;
  
           } else {
  
             return chi2_new;
           }
         }
  
         if (chi2_new < chi2_old) {
  
           chi2_old = chi2_new;
  
         } else {
             
           p1      = step;
           p1     *= lambda;
           value  -=  p1;                                       // step back
           lambda  = -lambda;                                   // change direction
           
           if (i != 0) {
             factor = 0.5;                                      // reduce step size
           }
         }
         
         lambda = factor * lambda;
       }
  
       return chi2_old;
     }

◆ evaluate()

template<class JModel_t >

template<class JFunction_t , class T >

double JFIT::JSimplex< JModel_t >::evaluate	(	const JFunction_t &	fit,
		T	__begin,
		T	__end
	)		const

inlineprivateinherited

Evaluate chi2 for given data set.

Parameters

fit	fit function
__begin	begin of data
__end	end of data

Definition at line 253 of file JSimplex.hh.

     {
       double chi2 = 0.0;
  
       for (T hit = __begin; hit != __end; ++hit) {
         chi2 += fit(value, *hit);
       }
  
       return chi2;
     }

Member Data Documentation

◆ T_ns

JTimeRange JFIT::JRegressor< JEnergy, JSimplex >::T_ns

static

Time window with respect to Cherenkov hypothesis [ns].

Default values.

Definition at line 193 of file JShowerEnergyRegressor.hh.

◆ Vmax_npe

double JFIT::JRegressor< JEnergy, JSimplex >::Vmax_npe = std::numeric_limits<double>::max()

static

Maximal integral of PDF [npe].

Definition at line 194 of file JShowerEnergyRegressor.hh.

◆ NUMBER_OF_PDFS

const int JFIT::JRegressor< JEnergy, JSimplex >::NUMBER_OF_PDFS = 2

static

Definition at line 196 of file JShowerEnergyRegressor.hh.

◆ pdf_t

const JPDFType_t JFIT::JRegressor< JEnergy, JSimplex >::pdf_t

static

Initial value:

= { DIRECT_LIGHT_FROM_BRIGHT_POINT,

SCATTERED_LIGHT_FROM_BRIGHT_POINT }

JPHYSICS::SCATTERED_LIGHT_FROM_BRIGHT_POINT

@ SCATTERED_LIGHT_FROM_BRIGHT_POINT

scattered light from bright point

Definition: JPDFTypes.hh:43

JPHYSICS::DIRECT_LIGHT_FROM_BRIGHT_POINT

@ DIRECT_LIGHT_FROM_BRIGHT_POINT

direct light from bright point

Definition: JPDFTypes.hh:42

PDF types.

Definition at line 198 of file JShowerEnergyRegressor.hh.

◆ Y

std::vector<JNPE_t> JFIT::JRegressor< JEnergy, JSimplex >::Y

light from EM showers

Definition at line 200 of file JShowerEnergyRegressor.hh.

◆ estimator

JLANG::JSharedPointer<JMEstimator> JFIT::JRegressor< JEnergy, JSimplex >::estimator

M-Estimator function.

Definition at line 202 of file JShowerEnergyRegressor.hh.

◆ MAXIMUM_ITERATIONS

template<class JModel_t >

int JFIT::JSimplex< JModel_t >::MAXIMUM_ITERATIONS = 1000

staticinherited

maximal number of iterations

maximal number of iterations.

Definition at line 237 of file JSimplex.hh.

◆ EPSILON

template<class JModel_t >

double JFIT::JSimplex< JModel_t >::EPSILON = 1.0e-4

staticinherited

maximal distance to minimum

maximal distance to minimum.

Definition at line 238 of file JSimplex.hh.

◆ value

template<class JModel_t >

JModel_t JFIT::JSimplex< JModel_t >::value

inherited

Definition at line 240 of file JSimplex.hh.

◆ step

template<class JModel_t >

std::vector<JModel_t> JFIT::JSimplex< JModel_t >::step

inherited

Definition at line 241 of file JSimplex.hh.

◆ numberOfIterations

template<class JModel_t >

int JFIT::JSimplex< JModel_t >::numberOfIterations

inherited

Definition at line 242 of file JSimplex.hh.

◆ p0

template<class JModel_t >

JModel_t JFIT::JSimplex< JModel_t >::p0

mutableprivateinherited

Definition at line 264 of file JSimplex.hh.

◆ p1

template<class JModel_t >

JModel_t JFIT::JSimplex< JModel_t >::p1

mutableprivateinherited

Definition at line 265 of file JSimplex.hh.

◆ wall

template<class JModel_t >

JModel_t JFIT::JSimplex< JModel_t >::wall

mutableprivateinherited

Definition at line 266 of file JSimplex.hh.

◆ debug

template<class T >

int JEEP::JMessage< T >::debug = 0

staticinherited

debug level (default is off).

Definition at line 45 of file JMessage.hh.

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

JShowerEnergyRegressor.hh

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ JFunction1D_t

◆ JMapList_t

◆ JPDF_t

◆ JNPE_t

◆ minimiser_type

◆ regressor_type

◆ result_type

Constructor & Destructor Documentation

◆ JRegressor()

Member Function Documentation

◆ operator()() [1/4]

◆ getNPE() [1/2]

◆ getNPE() [2/2]

◆ operator()() [2/4]

◆ operator()() [3/4]

◆ operator()() [4/4]

◆ evaluate()

Member Data Documentation

◆ T_ns

◆ Vmax_npe

◆ NUMBER_OF_PDFS

◆ pdf_t

◆ Y

◆ estimator

◆ MAXIMUM_ITERATIONS

◆ EPSILON

◆ value

◆ step

◆ numberOfIterations

◆ p0

◆ p1

◆ wall

◆ debug