Regressor function object for JShower3EZ fit using Abstract minimiser, that just computes the chi2 without a fit. More...

#include <JShower3EZRegressor.hh>

Inheritance diagram for JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >:

Public Types
typedef JTOOLS::JSplineFunction1S_t	JFunction1D_t

typedef JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > >	JPDFMaplist_t

typedef JPHYSICS::JPDFTable< JFunction1D_t, JPDFMaplist_t >	JPDF_t

typedef JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > >	JNPEMaplist_t

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

typedef JAbstractMinimiser< JShower3EZ >	minimiser_type

typedef JRegressor< JShower3EZ, JAbstractMinimiser >	regressor_type

typedef minimiser_type::result_type	result_type

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

double	operator() (const JShower3EZ &shower, const JPMTW0 &pmt) const
	Fit function. More...

JNPE_t::result_type	getH0 (const double R_Hz) const
	Get background hypothesis value for time integrated PDF. More...

JNPE_t::result_type	getH1 (const double D, const double cd, const double theta, const double phi, const double E) const
	Get signal hypothesis value for time integrated PDF. More...

JNPE_t::result_type	getH1 (const JShower3EZ &shower, const JPMTW0 &pmt) const
	Get signal hypothesis value for time integrated PDF. More...

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

template<class JFunction_t , class T >
result_type	operator() (const JFunction_t &fit, T __begin, T __end)
	Get chi2. More...

Public Attributes
JNPE_t	npe [NUMBER_OF_PDFS]
	PDF. More...

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

JModel_t	value
	model value More...

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]

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

Detailed Description

Regressor function object for JShower3EZ fit using Abstract minimiser, that just computes the chi2 without a fit.

Definition at line 478 of file JShower3EZRegressor.hh.

Member Typedef Documentation

◆ JFunction1D_t

typedef JTOOLS::JSplineFunction1S_t JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JFunction1D_t

Definition at line 483 of file JShower3EZRegressor.hh.

◆ JPDFMaplist_t

typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JPDFMaplist_t

Definition at line 487 of file JShower3EZRegressor.hh.

◆ JPDF_t

typedef JPHYSICS::JPDFTable<JFunction1D_t, JPDFMaplist_t> JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JPDF_t

Definition at line 488 of file JShower3EZRegressor.hh.

◆ JNPEMaplist_t

typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JNPEMaplist_t

Definition at line 493 of file JShower3EZRegressor.hh.

◆ JNPE_t

typedef JPHYSICS::JNPETable<double, double, JNPEMaplist_t> JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JNPE_t

Definition at line 494 of file JShower3EZRegressor.hh.

◆ minimiser_type

typedef JAbstractMinimiser <JShower3EZ > JFIT::JAbstractRegressor< JShower3EZ , JAbstractMinimiser >::minimiser_type

inherited

Definition at line 80 of file JRegressor.hh.

◆ regressor_type

typedef JRegressor<JShower3EZ , JAbstractMinimiser > JFIT::JAbstractRegressor< JShower3EZ , JAbstractMinimiser >::regressor_type

inherited

Definition at line 81 of file JRegressor.hh.

◆ result_type

typedef minimiser_type::result_type JFIT::JAbstractRegressor< JShower3EZ , JAbstractMinimiser >::result_type

inherited

Definition at line 82 of file JRegressor.hh.

Constructor & Destructor Documentation

◆ JRegressor()

JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::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<JShower3Z, JGandalf>::pdf_t.

Parameters

fileDescriptor PDF file descriptor

Definition at line 505 of file JShower3EZRegressor.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);
  
           npe[ i ] = JNPE_t(pdf);
         }
         catch(const JException& error) {
           FATAL(error.what() << endl);
         }
       }
  
       // Add PDFs
       for (int i = 1; i < NUMBER_OF_PDFS; i += 2) {
   
         npe[ i ].add(npe[i-1]);
  
         JNPE_t buffer;
         
         npe[i-1].swap(buffer);
       }    
     }

Member Function Documentation

◆ operator()() [1/3]

double JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::operator()	(	const JShower3EZ &	shower,
		const JPMTW0 &	pmt
	)		const

inline

Fit function.

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

Parameters

shower	shower
pmt	pmt

Returns: chi2

Definition at line 555 of file JShower3EZRegressor.hh.

     {
       using namespace JPP;
       
       JPosition3D  D(pmt.getPosition());
       JDirection3D U(pmt.getDirection());
  
       D.sub(shower.getPosition());
  
       const double z  = D.getDot(shower.getDirection());
       const double x  = D.getX();
       const double y  = D.getY();
       const double cd = z/D.getLength();                   // cosine angle between shower direction and PMT position
  
       U.rotate(JRotation3Z(-atan2(y,x)));                  // rotate PMT axis to x-z plane
  
       const double theta = getPMTAngle(U.getTheta());
       const double phi   = getPMTAngle(U.getPhi());
  
       JNPE_t::result_type H0 = getH0(pmt.getR());                                    // background hypothesis value for time integrated PDF. 
       JNPE_t::result_type H1 = getH1(D.getLength(), cd, theta, phi, shower.getE());  // signal hypothesis value for time integrated PDF. 
  
       if (get_value(H1) >= Vmax_npe) {
         H1 *= Vmax_npe / get_value(H1);
       }
  
       H1 += H0; // now H1 is signal + background
  
       const bool hit = pmt.getN() != 0;
       const double u = getChi2(get_value(H1), hit) - getChi2(get_value(H0), hit); // -log likelihood ratio
  
       return estimator->getRho(u);
     }

◆ getH0()

JNPE_t::result_type JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::getH0 ( const double R_Hz ) const

inline

Get background hypothesis value for time integrated PDF.

Parameters

R_Hz rate [Hz]

Returns: hypothesis value

Definition at line 595 of file JShower3EZRegressor.hh.

     {
       return JNPE_t::result_type(R_Hz * 1e-9 * T_ns.getLength());
     }

◆ getH1() [1/2]

JNPE_t::result_type JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::getH1	(	const double	D,
		const double	cd,
		const double	theta,
		const double	phi,
		const double	E
	)		const

inline

Get signal hypothesis value for time integrated PDF.

Parameters

D	PMT distance from shower [m]
cd	cosine angle between shower direction and PMT position
theta	PMT zenith angle [deg]
phi	PMT azimuth angle [deg]
E	shower energy [GeV]

Returns: hypothesis value

Definition at line 610 of file JShower3EZRegressor.hh.

     {
       JNPE_t::result_type h1 = JMATH::zero;
         
       for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
  
         if (!npe[i].empty() && D <= npe[i].getXmax()) {
  
           try {
  
             JNPE_t::result_type y1 = E * npe[i](std::max(D, npe[i].getXmin()), cd, theta, phi);
  
             // safety measure
  
             if(y1 < 0){
               y1 = 0;  
             }
             
             h1 += y1;
           
           }
           catch(JLANG::JException& error) {
             ERROR(error << std::endl);
           }
         }
       }
  
       return h1;
     }

◆ getH1() [2/2]

JNPE_t::result_type JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::getH1	(	const JShower3EZ &	shower,
		const JPMTW0 &	pmt
	)		const

inline

Get signal hypothesis value for time integrated PDF.

Parameters

shower	shower
pmt	pmt

Returns: hypothesis value

Definition at line 651 of file JShower3EZRegressor.hh.

     {
       using namespace JPP; 
     
       JPosition3D  D(pmt.getPosition());
       JDirection3D U(pmt.getDirection());
  
  
       const double z  = D.getDot(shower.getDirection());
       const double x  = D.getX();
       const double y  = D.getY();
       const double cd = z/D.getLength();                   // cosine angle between shower direction and PMT position
  
       U.rotate(JRotation3Z(-atan2(y,x)));                  // rotate PMT axis to x-z plane
  
       const double theta = getPMTAngle(U.getTheta());
       const double phi   = getPMTAngle(U.getPhi());
  
       JNPE_t::result_type H1 = getH1(D.getLength(), cd, theta, phi, 1.0);  // signal hypothesis value for time integrated PDF. E=1 because it is linear with E.
  
       if (get_value(H1) >= Vmax_npe) {
         H1 *= Vmax_npe / get_value(H1);
       }
  
     return H1;
     }

◆ operator()() [2/3]

result_type JFIT::JAbstractRegressor< JShower3EZ , JAbstractMinimiser >::operator()	(	const JShower3EZ &	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/3]

template<class JModel_t >

template<class JFunction_t , class T >

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

inlineinherited

Get chi2.

Parameters

fit	fit function
__begin	begin of data
__end	end of data

Returns: chi2

Definition at line 48 of file JRegressor.hh.

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

Member Data Documentation

◆ T_ns

JTimeRange JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::T_ns

static

Time window with respect to Cherenkov hypothesis [ns].

Definition at line 677 of file JShower3EZRegressor.hh.

◆ Vmax_npe

double JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::Vmax_npe = std::numeric_limits<double>::max()

static

Maximal integral of PDF [npe].

Definition at line 678 of file JShower3EZRegressor.hh.

◆ NUMBER_OF_PDFS

const int JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::NUMBER_OF_PDFS = 2

static

Definition at line 680 of file JShower3EZRegressor.hh.

◆ pdf_t

const JPDFType_t JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::pdf_t

static

Initial value:

= { DIRECT_LIGHT_FROM_EMSHOWER,

SCATTERED_LIGHT_FROM_EMSHOWER }

JPHYSICS::SCATTERED_LIGHT_FROM_EMSHOWER

@ SCATTERED_LIGHT_FROM_EMSHOWER

scattered light from EM shower

Definition: JPDFTypes.hh:38

JPHYSICS::DIRECT_LIGHT_FROM_EMSHOWER

@ DIRECT_LIGHT_FROM_EMSHOWER

direct light from EM shower

Definition: JPDFTypes.hh:37

Definition at line 682 of file JShower3EZRegressor.hh.

◆ npe

JNPE_t JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::npe[NUMBER_OF_PDFS]

PDF.

Definition at line 684 of file JShower3EZRegressor.hh.

◆ estimator

JLANG::JSharedPointer<JMEstimator> JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::estimator

M-Estimator function.

Definition at line 686 of file JShower3EZRegressor.hh.

◆ value

template<class JModel_t >

JModel_t JFIT::JAbstractMinimiser< JModel_t >::value

inherited

model value

Definition at line 59 of file JRegressor.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:

JShower3EZRegressor.hh

Public Types

Public Member Functions

Public Attributes

Static Public Attributes

Detailed Description

Member Typedef Documentation

◆ JFunction1D_t

◆ JPDFMaplist_t

◆ JPDF_t

◆ JNPEMaplist_t

◆ JNPE_t

◆ minimiser_type

◆ regressor_type

◆ result_type

Constructor & Destructor Documentation

◆ JRegressor()

Member Function Documentation

◆ operator()() [1/3]

◆ getH0()

◆ getH1() [1/2]

◆ getH1() [2/2]

◆ operator()() [2/3]

◆ operator()() [3/3]

Member Data Documentation

◆ T_ns

◆ Vmax_npe

◆ NUMBER_OF_PDFS

◆ pdf_t

◆ npe

◆ estimator

◆ value

◆ debug